gsThermoAssembler.h

Go to the documentation of this file.

#pragma once

#include <gsElasticity/gsElasticityAssembler.h>

namespace gismo
{

template <class T>
class gsThermoAssembler : public gsElasticityAssembler<T>
{
public:
    typedef gsElasticityAssembler<T> Base;

    gsThermoAssembler(const gsMultiPatch<T> & patches,
                      const gsMultiBasis<T> & bases,
                      const gsBoundaryConditions<T> & b_conditions,
                      const gsFunction<T> & body_force,
                      const gsFunctionSet<T> & temperature_field);

    virtual void assemble(bool saveEliminationMatrix = false);

    void assembleThermo();

protected:
    void findNonDirichletSides();

protected:
    const gsFunctionSet<T> & m_temperatureField;
    bool assembledElasticity;
    std::vector<std::pair<index_t,boxSide> > nonDirichletSides;
    gsMatrix<T> elastRhs;

    using Base::m_pde_ptr;
    using Base::m_options;

}; // class definition ends
} // namespace ends

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