linearFormHelper.hh

Go to the documentation of this file.

 
#ifndef hp2dlinearformhelper_hh
#define hp2dlinearformhelper_hh
 
#include "formula/elementFormulaContainer.hh"
 
namespace concepts {
 
  // forward declarations
  class SharedJacobianDet;
 
  template<class F>
  class Array;
 
  template<class F, uint dim>
  class Point;
}
 
namespace hp2D {
 
  // forward declarations
  template<class F>
  class BaseQuad;
 
  template<class F>
  class ArrayElementFormula;
 
  using concepts::Real;
 
  // *************************************************** LinearFormHelper_0 **
 
  template<class F>
  class LinearFormHelper_0 {
  public:
    LinearFormHelper_0
    (const concepts::ElementFormulaContainer<F> 
     frm = concepts::ElementFormulaContainer<F>());
 
    virtual ~LinearFormHelper_0();
 
    void data(const concepts::RCP<concepts::SharedJacobianDet> d);
    concepts::RCP<concepts::SharedJacobianDet> data() const;
  protected:
    void computeIntermediate_(const BaseQuad<concepts::Real>& elm) const;
 
    mutable concepts::Array<F> intermediateValue_;
 
    concepts::ElementFormulaContainer<F> frm_;
    concepts::RCP<concepts::SharedJacobianDet> sharedData_;
  };
 
  // **************************************************** LinearFormHelper_1 **
 
  template<class F>
  class LinearFormHelper_1 {
  public:
    LinearFormHelper_1(const concepts::ElementFormulaContainer<F> frm1,
                       const concepts::ElementFormulaContainer<F> frm2);
 
    LinearFormHelper_1
    (const concepts::ElementFormulaContainer<concepts::Point<F, 2> > frm);
    virtual ~LinearFormHelper_1();
  protected:
    void computeIntermediate_(const BaseQuad<concepts::Real>& elm) const;
 
    mutable ArrayElementFormula<concepts::Point<F,2> > intermediateVector_;
 
    concepts::ElementFormulaContainer<concepts::Point<F, 2> > frm_;
  };
 
  // **************************************************** LinearFormHelper_2 **
 
  template<class F>
  class LinearFormHelper_2 {
  public:
    LinearFormHelper_2(const concepts::ElementFormulaContainer<F> frm);
 
    virtual ~LinearFormHelper_2() {}
  protected:
    void computeIntermediate_(const BaseQuad<concepts::Real>& elm) const;
 
    mutable ArrayElementFormula<F> intermediateValue_;
 
    concepts::ElementFormulaContainer<F> frm_;
  };
 
 
} // namespace hp2D
 
#endif // hp2dlinearformhelper_hh