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
    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