bilinearForm.hh

Go to the documentation of this file.

 
#ifndef linDG1DbilinearForm_hh
#define linDG1DbilinearForm_hh
 
#include "basics/typedefs.hh"
#include "operator/bilinearForm.hh"
#include "geometry/boundaryConditions.hh"
 
namespace linDG1D {
 
  using concepts::Real;
 
  // *********************************************************** BoundaryInt **
 
  class BoundaryInt : public concepts::BilinearForm<Real> {
  public:
    BoundaryInt(concepts::BoundaryConditions* bc = 0) : bc_(bc) {}
    virtual void operator()(const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY,
          concepts::ElementMatrix<Real>& em) const;
    virtual void operator()(const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY,
          concepts::ElementMatrix<Real>& em,
          const concepts::ElementPair<Real>& ep) const;
    virtual BoundaryInt* clone() const { return new BoundaryInt(bc_); }
  private:
    concepts::BoundaryConditions* bc_;
  };
 
  // ******************************************************* BoundaryIntStab **
 
  class BoundaryIntStab : public concepts::BilinearForm<Real> {
  public:
    virtual void operator()(const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY,
          concepts::ElementMatrix<Real>& em) const{
      conceptsAssert(false, concepts::Assertion());
    }
    virtual void operator()(const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY,
          concepts::ElementMatrix<Real>& em,
          const concepts::ElementPair<Real>& ep) const;
    virtual BoundaryIntStab* clone() const { return new BoundaryIntStab(); }
  };
 
} // namespace linDG1D
 
#endif // linDG1DbilinearForm_hh