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 "function/linearForm.hh"
#include "integration.hh"
#include "timestepping/vectors.hh"
 
#include "triangle.hh"
#include "meshInfo.hh"
 
#define DEBUG_FILLDATA 0
 
using concepts::Real;
 
namespace linDG2D {
 
  class LaplaceTimeLf : public timestepping::TimeLinearForm {
  public:
    LaplaceTimeLf(concepts::Formula<Real>& frm, 
                      concepts::BoundaryConditions* bc,
                      MeshInfo& mi,
                      Real omega, Real tau,
                      uint nQuadPts) :
      timestepping::TimeLinearForm(), frm_(frm), bc_(bc) , mi_(mi), 
      omega_(omega), tau_(tau), nQuadPts_(nQuadPts) { }
    void operator() (const concepts::Element<Real> &elmBase,
                     concepts::ElementMatrix<Real> &em);
  private:
    const concepts::Formula<Real>& frm_;
    const concepts::BoundaryConditions* bc_;
    const MeshInfo& mi_;
    const Real omega_;
    const Real tau_;
    const uint nQuadPts_;
  };
 
  class ScalarProductLf : public concepts::LinearForm<Real> {
  public:
    ScalarProductLf(concepts::Formula<Real>& frm, uint nQuadPts) 
      : frm_(frm), nQuadPts_(nQuadPts) { };
    void operator() (const concepts::Element<Real> &elmBase,
                     concepts::ElementMatrix<Real> &em);
  private:
    const concepts::Formula<Real>& frm_;
    const uint nQuadPts_;
  };
 
} // namespace linDG2D
 
#endif //dgLinearForms_hh