bilinearFormHelper.hh

Go to the documentation of this file.

 
#ifndef hp2dbilinearformhelper_hh
#define hp2dbilinearformhelper_hh
 
#include "linearFormHelper.hh"
#include "hp2D/arrayElementFormula.hh"
 
namespace concepts {
 
  // forward declarations
  class SharedJacobianDet;
 
  template<uint dim>
  class SharedJacobianAdj;
 
  template<class F>
  class Array;
 
  template<class F, class G>
  class ElementFormula;
 
  template<class F, uint dim>
  class Point;
}
 
namespace hp2D {
 
  // forward declarations
  template<class F>
  class BaseQuad;
 
  template<class F>
  class ArrayElementFormula;
 
  using concepts::Real;
 
  // *********************************************** BilinearFormHelper_0_0 **
 
  template<class F>
  class BilinearFormHelper_0_0 : public LinearFormHelper_0<F> {
  public:
    BilinearFormHelper_0_0
    (const concepts::ElementFormulaContainer<F> 
     frm = concepts::ElementFormulaContainer<F>())
      : LinearFormHelper_0<F>(frm) {}
    
    const concepts::ElementFormulaContainer<F> formula() const { 
      return this->frm_; }
  };
 
  // ******************************************** BilinearFormHelper_0_1_Part **
 
  template<class F>
  class BilinearFormHelper_0_1_Part {
  public:
    BilinearFormHelper_0_1_Part(const concepts::ElementFormulaContainer<F> frm);
 
    virtual ~BilinearFormHelper_0_1_Part();
 
    void data(const concepts::RCP<concepts::SharedJacobianAdj<2> > d);
    concepts::RCP<concepts::SharedJacobianAdj<2> > data() const;
  protected:
    void computeIntermediate_(const BaseQuad<concepts::Real>& elm, const int i) const;
    mutable ArrayElementFormula<concepts::Point<F,2> > intermediateVector_;
    concepts::ElementFormulaContainer<F> frm_;
 
    concepts::RCP<concepts::SharedJacobianAdj<2> > sharedData_;
  };
 
  // ************************************************ BilinearFormHelper_1_1 **
 
  template<class F, class G = typename concepts::Realtype<F>::type>
  class BilinearFormHelper_1_1 {
  public:
    BilinearFormHelper_1_1(const concepts::ElementFormulaContainer<F> frm);
    BilinearFormHelper_1_1
    (const concepts::ElementFormulaContainer<concepts::Mapping<G,2> > frmM);
 
    virtual ~BilinearFormHelper_1_1();
 
    void data(const concepts::RCP<concepts::SharedJacobianAdj<2> > d);
    concepts::RCP<concepts::SharedJacobianAdj<2> > data() const;
  protected:
    void computeIntermediate_(const BaseQuad<Real>& elm,
                              const int i = -1, const int j = -1) const;
    mutable concepts::Array<F> intermediateValue_;
    mutable concepts::Array<concepts::Mapping<G,2> > intermediateMatrix_;
    concepts::ElementFormulaContainer<F> frm_;
    concepts::ElementFormulaContainer<concepts::Mapping<G,2> > frmM_;
  private:
    void computeJacobianMatrix_(const BaseQuad<Real>& elm,
                                concepts::Array<concepts::Mapping<Real,2> >& J, 
                                concepts::Array<Real>& detJ_inv) const;
    void computeadjJ_adjJT_rank1_
    (concepts::Array<concepts::Mapping<Real,2> >& intermediateMatrix,
     const int i, const int j) const;
 
    concepts::RCP<concepts::SharedJacobianAdj<2> > sharedData_;
  };
 
  // ************************************************ BilinearFormHelper_2_1 **
 
  template<class F>
  class BilinearFormHelper_2_1 {
  public:
    BilinearFormHelper_2_1(const concepts::ElementFormulaContainer<F> frm1,
         const concepts::ElementFormulaContainer<F> frm2);
 
    BilinearFormHelper_2_1
    (const concepts::ElementFormulaContainer<concepts::Point<F, 2> > frm);
 
    virtual ~BilinearFormHelper_2_1();
  protected:
    void computeIntermediate_(const BaseQuad<Real>& elm) const;
    mutable ArrayElementFormula<concepts::Point<F,2> > intermediateVector_;
 
    concepts::ElementFormulaContainer<concepts::Point<F, 2> > frm_;
  };
 
  // ************************************************ BilinearFormHelper_2_2 **
 
  template<class F>
  class BilinearFormHelper_2_2 {
  public:
    BilinearFormHelper_2_2(const concepts::ElementFormulaContainer<F> 
                           frm = concepts::ElementFormulaContainer<F>());
    virtual ~BilinearFormHelper_2_2();
  protected:
    void computeIntermediate_(const BaseQuad<Real>& elm) const;
    mutable concepts::Array<F> intermediateValue_;
    concepts::ElementFormulaContainer<F> frm_;
  };
 
} // namespace hp2D
 
#endif // hp2dbilinearformhelper_hh