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