gsBasePde.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsPde/gsPde.h>
#include <gsCore/gsPiecewiseFunction.h>
 
namespace gismo
{
 
template<class T>
class gsBasePde : public gsPde<T>
{
 
public:
 
    gsBasePde( ) { }
 
    gsBasePde(const gsMultiPatch<T>         &domain,
              const gsBoundaryConditions<T> &bc,
              const gsPiecewiseFunction<T>  &rhs)
    : gsPde<T>(domain,bc), m_rhs(rhs)
    {
        m_unknownDim.setOnes(m_rhs.targetDim());
    }
 
    virtual int numRhs() const { return 1; }
    virtual std::ostream &print(std::ostream &os) const { return os; }
    const gsFunction<T> *    rhs()      const { return &m_rhs.piece(0); }
    virtual int numUnknowns() const     {return m_rhs.targetDim();}
 
protected:
    using gsPde<T>::m_unknownDim;
    gsPiecewiseFunction<T> m_rhs;
};
 
} // namespace gismo