hp2D/bilinearFormHelper.hh

Go to the documentation of this file.

/* @file bilinearFormHelper.hh Helper classes for linear forms in hp2D
 */

#ifndef hp2dbilinearformhelper_hh
#define hp2dbilinearformhelper_hh

#include "linearFormHelper.hh"

namespace concepts {

  // forward declarations
  class SharedJacobianDet;

  template<class F>
  class Array;

  template<class F, class G>
  class ElementFormula;

  template<class F, int 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 **

  template<class F>
  class BilinearFormHelper_0_1 : public LinearFormHelper_1<F> {
  public:
    BilinearFormHelper_0_1(const concepts::ElementFormulaContainer<F> frm1,
                           const concepts::ElementFormulaContainer<F> frm2)
      : LinearFormHelper_1<F>(frm1, frm2) {}

    BilinearFormHelper_0_1
    (const concepts::ElementFormulaContainer<concepts::Point<F,2> >& frm) 
      : LinearFormHelper_1<F>(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);
    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);
    concepts::Array<F> intermediateValue_;
    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);
    void computeadjJ_adjJT_rank1_
    (concepts::Array<concepts::Mapping<Real,2> >& intermediateMatrix,
     const int i, const int j);

    concepts::RCP<concepts::SharedJacobianAdj<2> > sharedData_;
  };

  // ************************************************ 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);
    concepts::Array<F> intermediateValue_;
    concepts::ElementFormulaContainer<F> frm_;
  };

} // namespace hp2D

#endif // hp2dbilinearformhelper_hh