levelRieszElement.hh

#ifndef LevelRieszElement_hh
#define LevelRieszElement_hh
 
#include "../hp1D.hh"
#include "../hp2D.hh"
#include "../formula.hh"
#include "../function.hh"
 
using concepts::Real;
using concepts::Real2d;
 
namespace hp1D {
 
  //*************************************************** LevelRieszElement *****
  class LevelRieszElement : public concepts::LinearForm<concepts::Real> {
    public:
 
    LevelRieszElement(const hp2D::Element<Real>& elm,
                      const concepts::ElementFormulaContainer<Real>& levelFunction,
                      const concepts::ElementFormulaContainer<Real2d>& levelFunctionGrad,
                      const concepts::ElementFormulaContainer<Real>& observable);
 
    ~LevelRieszElement();
 
    void operator()(const concepts::Element<Real>& elm,
                    concepts::ElementMatrix<Real>& em) const;
 
    protected:
 
      virtual std::ostream& info(std::ostream& os) const;
 
    private:
      const hp2D::Element<Real>& elm2D_;
      const concepts::ElementFormulaContainer<Real> levelFunction_;
      const concepts::ElementFormulaContainer<Real2d> levelFunctionGrad_;
      const concepts::ElementFormulaContainer<Real> observable_;
 
 
  };//class LevelRiesz
 
}// namespace hp1D
 
#endif    //levelRiesz_hh