bilinearForm.hh

Go to the documentation of this file.

 
#ifndef linDG2DbilinearForm_hh
#define linDG2DbilinearForm_hh
 
#include "basics/typedefs.hh"
#include "operator/bilinearForm.hh"
#include "geometry/boundaryConditions.hh"
 
using concepts::Real;
 
namespace linDG2D {
 
  // ************************************************************ IdentityBf **
 
  class IdentityBf : public concepts::BilinearForm<Real> {
  public:
    IdentityBf() {}
    virtual void operator()(const concepts::Element<Real>& elmV,
                            const concepts::Element<Real>& elmU,
                            concepts::ElementMatrix<Real>& em) const;
    virtual IdentityBf* clone() const { return new IdentityBf(); }
  };
 
  // ********************************************************** LaplaceVolBf **
 
  class LaplaceVolBf : public concepts::BilinearForm<Real> {
  public:
    LaplaceVolBf() {}
    virtual void operator()(const concepts::Element<Real>& elmV,
                            const concepts::Element<Real>& elmU,
                            concepts::ElementMatrix<Real>& em) const;
    virtual LaplaceVolBf* clone() const { return new LaplaceVolBf(); }
  };
 
  // ********************************************************* LaplaceFluxBf **
 
  class LaplaceFluxBf : public concepts::BilinearForm<Real> {
  public:
    LaplaceFluxBf(const concepts::BoundaryConditions* bc) 
      : bc_(bc) {}
    virtual void operator()(const concepts::Element<Real>& elmV,
                            const concepts::Element<Real>& elmU,
                            concepts::ElementMatrix<Real>& em) const{
      conceptsAssert(false, concepts::Assertion());
    }
    virtual void operator()(const concepts::Element<Real>& elmV,
                            const concepts::Element<Real>& elmU,
                            concepts::ElementMatrix<Real>& em,
                            const concepts::ElementPair<Real>& ep) const;
    virtual LaplaceFluxBf* clone() const { return new LaplaceFluxBf(bc_); }
  private:
    const concepts::BoundaryConditions* bc_;
  };
 
  // ********************************************************* LaplaceStabBf **
 
  class LaplaceStabBf : public concepts::BilinearForm<Real> {
  public:
    LaplaceStabBf(Real omega, const concepts::BoundaryConditions* bc)
      : omega_(omega), bc_(bc) {}
    virtual void operator()(const concepts::Element<Real>& elmV,
                            const concepts::Element<Real>& elmU,
                            concepts::ElementMatrix<Real>& em) const {
      conceptsAssert(false, concepts::Assertion());
    }
    virtual void operator()(const concepts::Element<Real>& elmV,
                            const concepts::Element<Real>& elmU,
                            concepts::ElementMatrix<Real>& em,
                            const concepts::ElementPair<Real>& ep) const;
    virtual LaplaceStabBf* clone() const {
      return new LaplaceStabBf(omega_, bc_); }
  private:
    const Real omega_;
    const concepts::BoundaryConditions* bc_;
  };
 
} // namespace linDG2D
 
#endif // linDG2DbilinearForm_hh