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