linDG3D/linearForm.hh

Go to the documentation of this file.

/* Class for DG linear forms
 */

#ifndef dgLinearForms_hh
#define dgLinearForms_hh

#include "basics/typedefs.hh"
#include "formula/formula.hh"
#include "formula/boundary.hh"
#include "geometry/boundaryConditions.hh"
#include "timestepping/vectors.hh"
#include "element.hh"

#include <vector>

#define DEBUG_FILLDATA 0

namespace linDG3D {
  using concepts::Real;

  // ************************************************************** ConstRhs **

  class ConstRhs : public timestepping::TimeLinearForm {
  public:
    ConstRhs(Real c = 1.0) : c_(c) {}
    virtual void operator()(const concepts::Element<Real>& elm,
                            concepts::ElementMatrix<Real>& em);
    void operator()(const FvdgElement& elm,
                    concepts::ElementMatrix<Real>& em);
  private:
    Real c_;
  };

  // *********************************************************** FvdgTimeRhs **
  class FvdgTimeRhs : public timestepping::TimeLinearForm {
  public:
    FvdgTimeRhs(concepts::Formula<Real>& frm) :
      timestepping::TimeLinearForm(),
      frm_(frm) {}
    void operator() (const concepts::Element<Real> &elmBase,
                     concepts::ElementMatrix<Real> &em);
  private:
    const concepts::Formula<Real>& frm_;
  };

  //******************************************************** FvdgTimeRhsDiff **
  class FvdgTimeRhsDiff : public timestepping::TimeLinearForm {
  public:
    FvdgTimeRhsDiff(concepts::Formula<Real>& frm, Real eps,
                    concepts::BoundaryConditions* bc, bool symm = false) :
      timestepping::TimeLinearForm(), frm_(frm), eps_(eps), bc_(bc) {
        symm_ = symm ? -1 : 1;
      }
    void operator() (const concepts::Element<Real> &elmBase,
                     concepts::ElementMatrix<Real> &em);
  private:
    const concepts::Formula<Real>& frm_;
    const Real eps_;
    const concepts::BoundaryConditions* bc_;
    int symm_;
    };

} // namespace linDG3D

#endif //dgLinearForms_hh