linearForm.hh

Go to the documentation of this file.

 
#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