gsCDRAssembler.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsAssembler/gsPoissonAssembler.h>
#include <gsPde/gsConvDiffRePde.h>
 
namespace gismo
{
 
 
struct stabilizerCDR
{
    enum method
    {
        SUPG = 1, 
        none = 0  
    };
};
 
 
template<class T>
class gsCDRAssembler : public gsPoissonAssembler<T>
{
public:
    typedef gsPoissonAssembler<T> Base;
 
public:
 
    static gsOptionList defaultOptions();
 
    gsCDRAssembler(const gsConvDiffRePde<T> & pde,
                   const gsMultiBasis<T>    & bases,
                   stabilizerCDR::method      flagStabilization = stabilizerCDR::none)
    {
        // enrich options in constructor, refresh to apply options
 
        // 0: no stabilization
        // 1: SUPG
        m_options = defaultOptions();
        m_options.addInt("Stabilization", "Choice of stabilization method; 0 := no; 1 := SUPG;", flagStabilization);
        Base::initialize(pde, bases, m_options);
    }
 
    gsCDRAssembler(const gsConvDiffRePde<T> & pde,
                   const gsMultiBasis<T>    & bases,
                   dirichlet::strategy        dirStrategy,
                   iFace::strategy            intStrategy = iFace::glue,
                   stabilizerCDR::method      flagStabilization = stabilizerCDR::none)
    {
        m_options = defaultOptions();
        m_options.setInt("DirichletStrategy", dirStrategy);
        m_options.setInt("InterfaceStrategy", intStrategy);
 
        // 0: no stabilization
        // 1: SUPG
        m_options.addInt("Stabilization", "Choice of stabilization method; 0 := no; 1 := SUPG;", flagStabilization);
 
        Base::initialize(pde, bases, m_options);
    }
 
    gsCDRAssembler(gsMultiPatch<T> const         & patches,
                   gsMultiBasis<T> const         & bases,
                   gsBoundaryConditions<T> const & bconditions,
                   const gsFunction<T>           & rhs,
                   const gsFunction<T>           & coeff_A,
                   const gsFunction<T>           & coeff_b,
                   const gsFunction<T>           & coeff_c,
                   dirichlet::strategy             dirStrategy = dirichlet::elimination,
                   iFace::strategy                 intStrategy = iFace::glue,
                   stabilizerCDR::method           flagStabilization = stabilizerCDR::none)
    {
        m_options = defaultOptions();
        m_options.setInt("DirichletStrategy", dirStrategy);
        m_options.setInt("InterfaceStrategy", intStrategy);
 
        // 0: no stabilization
        // 1: SUPG
        m_options.addInt("Stabilization", "Choice of stabilization method; 0 := no; 1 := SUPG;", flagStabilization);
 
        typename gsPde<T>::Ptr pde(new gsConvDiffRePde<T>
                                       (patches, bconditions, &coeff_A, &coeff_b, &coeff_c, &rhs));
        Base::initialize(pde, bases, m_options);
    }
 
    using Base::assemble;
    void assemble();
 
protected:
 
    // Members from gsAssembler
    using Base::m_pde_ptr;
    using Base::m_bases;
    using Base::m_ddof;
    using Base::m_options;
    using Base::m_system;
};
 
} // namespace gismo
 
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsCDRAssembler.hpp)
#endif