gsALE.h

Go to the documentation of this file.

 
#pragma once
#include <gsIO/gsOptionList.h>
#include <gsElasticity/gsBaseUtils.h>
#include <gsCore/gsMultiPatch.h>
#include <gsElasticity/gsIterative.h>
#include <gsElasticity/gsBaseAssembler.h>
 
namespace gismo
{
 
template <class T>
class gsBaseAssembler;
template <class T>
class gsIterative;
 
template <class T>
class gsALE
{
public:
    gsALE(gsMultiPatch<T> & geometry, const gsMultiPatch<T> & displacement,
          const gsBoundaryInterface & interfaceStr2Mesh, ale_method::method method);
 
    static gsOptionList defaultOptions();
 
    gsOptionList & options() { return m_options; }
 
    index_t numDofs() const {return assembler->numDofs();}
 
    void constructSolution(gsMultiPatch<T> & solution) const;
 
    index_t updateMesh();
 
    void saveState();
 
    void recoverState();
 
    const gsBoundaryInterface & interface() { return m_interface;}
 
protected:
    void initialize();
 
    index_t linearMethod();
 
    index_t linearIncrementalMethod();
 
    index_t nonlinearMethod();
 
protected:
    const gsMultiPatch<T> & disp;
    const gsBoundaryInterface & m_interface;
    ale_method::method methodALE;
    gsOptionList m_options;
    typename gsBaseAssembler<T>::uPtr assembler;
    typename gsIterative<T>::uPtr solverNL;
    gsMultiPatch<T> ALEdisp;
    bool initialized;
 
 
    bool hasSavedState;
    gsMultiPatch<T> ALEdispSaved;
};
 
} // namespace ends
 
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsALE.hpp)
#endif