gsBiharmonicPde.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsPde/gsPde.h>
 
namespace gismo
{
 
template<class T>
class gsBiharmonicPde : public gsPde<T>
{
protected:
    gsBiharmonicPde( ) { }
    using gsPde<T>::m_domain;
public:
    gsBiharmonicPde(
        const gsMultiPatch<T>         &domain,
        const gsBoundaryConditions<T> &bc1,
        const gsBoundaryConditions<T> &bc2,
        const gsPiecewiseFunction<T>  &rhs
         )
    : gsPde<T>(domain,bc1), m_rhs(rhs), m_boundary_conditions_second(bc2)
    {
        m_unknownDim.setOnes(1);
    }
 
 
    ~gsBiharmonicPde( )
    { }
 
    virtual int numRhs() const
    {
        return m_rhs.piece(0).targetDim();
    }
 
    const gsFunction<T> *    rhs()      const { return &m_rhs.piece(0); }
 
    const gsBoundaryConditions<T> & bcFirstKind()  const {return this->bc();}
 
    const gsBoundaryConditions<T> & bcSecondKind() const {return m_boundary_conditions_second;}
 
    gsBoundaryConditions<T> &boundaryConditionsSecond() {return m_boundary_conditions_second;}
 
    const gsBoundaryConditions<T> &boundaryConditionsSecond() const {return m_boundary_conditions_second;}
 
    std::ostream &print(std::ostream &os) const
    {
        os<<"Biharmonic's equation  -\u0394^2 u = f ,  with:\n";
            os<<"Source function f= "<< m_rhs.piece(0) <<".\n";
            return os; 
        }
protected:
    using gsPde<T>::m_unknownDim;
 
    gsPiecewiseFunction<T> m_rhs;
    gsBoundaryConditions<T> m_boundary_conditions_second;
}; // class gsBiharmonicPde
 
} // namespace gismo