gsNsAssembler.h

Go to the documentation of this file.

#pragma once

#include <gsElasticity/gsBaseAssembler.h>
#include <gsElasticity/gsBaseUtils.h>

namespace gismo
{

template <class T>
class gsNsAssembler : public gsBaseAssembler<T>
{
public:
    typedef gsBaseAssembler<T> Base;

    gsNsAssembler(const gsMultiPatch<T> & patches,
                  const gsMultiBasis<T> & basisVel,
                  const gsMultiBasis<T> & basisPres,
                  const gsBoundaryConditions<T> & bconditions,
                  const gsFunction<T> & body_force);

    static gsOptionList defaultOptions();

    virtual void refresh();

    //--------------------- SYSTEM ASSEMBLY ----------------------------------//

    virtual void assemble(bool saveEliminationMatrix = false);

    virtual bool assemble(const gsMatrix<T> & solutionVector,
                          const std::vector<gsMatrix<T> > & fixedDoFs);

    virtual void assemble(const gsMultiPatch<T> & velocity, const gsMultiPatch<T> & pressure);

    //--------------------- SOLUTION CONSTRUCTION ----------------------------------//

    virtual void constructSolution(const gsMatrix<T> & solVector,
                                   const std::vector<gsMatrix<T> > & fixedDoFs,
                                   gsMultiPatch<T> & velocity) const;

    virtual void constructSolution(const gsMatrix<T> & solVector,
                                   const std::vector<gsMatrix<T> > & fixedDoFs,
                                   gsMultiPatch<T> & velocity, gsMultiPatch<T> & pressure) const;

    virtual void constructPressure(const gsMatrix<T> & solVector,
                                   const std::vector<gsMatrix<T> > & fixedDoFs,
                                   gsMultiPatch<T> & pressure) const;

    //--------------------- SPECIALS ----------------------------------//

    virtual gsMatrix<T> computeForce(const gsMultiPatch<T> & velocity, const gsMultiPatch<T> & pressure,
                                     const std::vector<std::pair<index_t,boxSide> > & bdrySides, bool split = false) const;

protected:
    virtual void reserve();

protected:

    short_t m_dim;

    using Base::m_pde_ptr;
    using Base::m_bases;
    using Base::m_ddof;
    using Base::m_options;
    using Base::m_system;
};

} // namespace gismo ends

#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsNsAssembler.hpp)
#endif