gsPoissonAssembler.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsAssembler/gsAssembler.h>
#include <gsPde/gsPoissonPde.h>
 
 
namespace gismo
{
 
template <class T>
class gsPoissonAssembler : public gsAssembler<T>
{
public:
    typedef gsAssembler<T> Base;
 
public:
 
    static gsOptionList defaultOptions();
 
    gsPoissonAssembler()
    { }
 
    gsPoissonAssembler( const gsPoissonPde<T>          & pde,
                        const gsMultiBasis<T>          & bases)
    {
        m_options = defaultOptions();
        Base::initialize(pde, bases, m_options);
    }
 
 
    gsPoissonAssembler( const gsPoissonPde<T>          & pde,
                        const gsMultiBasis<T>          & bases,
                        dirichlet::strategy           dirStrategy,
                        iFace::strategy               intStrategy = iFace::glue)
    {
        m_options = defaultOptions();
        m_options.setInt("DirichletStrategy", dirStrategy);
        m_options.setInt("InterfaceStrategy", intStrategy);
 
        Base::initialize(pde, bases, m_options);
    }
 
    gsPoissonAssembler( gsMultiPatch<T> const         & patches,
                        gsMultiBasis<T> const         & basis,
                        gsBoundaryConditions<T> const & bconditions,
                        const gsFunction<T>           & rhs,
                        dirichlet::strategy           dirStrategy = dirichlet::elimination,
                        iFace::strategy               intStrategy = iFace::glue)
    {
        m_options = defaultOptions();
        m_options.setInt("DirichletStrategy", dirStrategy);
        m_options.setInt("InterfaceStrategy", intStrategy);
 
        typename gsPde<T>::Ptr pde( new gsPoissonPde<T>(patches,bconditions,rhs) );
        Base::initialize(pde, basis, m_options);
    }
 
    virtual gsAssembler<T>* clone() const
    {
        return new gsPoissonAssembler<T>(*this);
    }
 
    virtual gsAssembler<T>* create() const
    {
        return new gsPoissonAssembler<T>();
    }
 
    // Refresh routine
    virtual void refresh();
 
    // Main assembly routine
    using Base::assemble;
    virtual void assemble();
 
    gsEigen::SparseSelfAdjointView< typename gsSparseMatrix<T>::Base, Lower> fullMatrix()
    {
        return m_system.matrix().template selfadjointView<Lower>();
    }
 
 
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(gsPoissonAssembler.hpp)
#endif