linearForm.hh

Go to the documentation of this file.

 
#ifndef hp3dlinearform_hh
#define hp3dlinearform_hh
 
#include <memory>
#include "formula/elementFormulaContainer.hh"
#include "function/linearForm.hh"
#include "basics/typedefs.hh"
 
namespace concepts {
  // forward declarations
  template<class F>
  class Formula;
 
  class BoundaryConditions;
 
  template<class F>
  class Element;
 
  template<class F>
  class ElementMatrix;
 
  template<class F>
  class Array;
}
 
namespace hp3D {
  // forward declarations
  class Hexahedron;
 
  using concepts::Real;
 
  // ***************************************************************** Riesz **
 
  class Riesz : public concepts::LinearForm<Real> {
  public:
    Riesz(const concepts::ElementFormulaContainer<Real, Real>& frm,
    concepts::BoundaryConditions* bc = 0); 
    virtual ~Riesz();
 
    void operator()(const concepts::Element<Real>& elm,
        concepts::ElementMatrix<Real>& em) const;
    void operator()(const Hexahedron& elm, concepts::ElementMatrix<Real>& em) const;
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    concepts::ElementFormulaContainer<Real, Real> frm_;
 
    mutable concepts::Array<Real> jacobian_;
 
    concepts::BoundaryConditions* bc_;
  }; 
 
} // namespace hp3D
 
#endif // hp3dlinearform_hh