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);
 
    using Base::assemble;
    virtual void assemble(bool /* saveEliminationMatrix */);
 
    virtual void assemble() { assemble(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