levelRiesz.hh

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