hp2D/bf_identity.hh

Go to the documentation of this file.

/* @bf_identity.hh Bilinear forms for Identity operator for hp 2D FEM
 */

#ifndef hp2dbf_identity_hh
#define hp2dbf_identity_hh

#include "basics/typedefs.hh"
#include "operator/bilinearForm.hh"
#include "geometry/sharedJacobian.hh"
#include "hp2D/bilinearFormHelper.hh"


namespace concepts {
  // forward declarations
  template<class F>
  class Array;

  template<class F, int dim>
  class Point;

  template<class F>
  class Element;

  template<class F>
  class ElementMatrix;

  template<typename F, typename G>
  class ElementFormula;
}

namespace vectorial {
  // forward declarations
  template<class F, class G>
  class BilinearForm;
}

namespace hp2D {
  // forward declarations
  template<class F>
  class Quad;

  class InfiniteLaguerreQuad;

  using concepts::Real;


  // ************************************************************** Identity **

  template<class F = Real>
  class Identity : public concepts::BilinearForm<F,Real>,
                   public BilinearFormHelper_0_0<F> {
  public:
    Identity(const concepts::ElementFormulaContainer<F> 
             frm = concepts::ElementFormulaContainer<F>(), 
             bool all = false);
    virtual ~Identity();
    virtual Identity<F>* clone() const;

    void operator()(const concepts::Element<Real>& elmX,
                    const concepts::Element<Real>& elmY, 
                    concepts::ElementMatrix<F>& em); 
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    concepts::ElementMatrix<Real> mass1D_;
    bool all_;

    bool operator()(const Quad<Real>* elmX,
                    const Quad<Real>* elmY, 
                    concepts::ElementMatrix<F>& em);
    bool operator()(const InfiniteLaguerreQuad* elmX,
                    const InfiniteLaguerreQuad* elmY, 
                    concepts::ElementMatrix<F>& em);
  };


  // ********************************************************** setupIdentity **

  template<class F>
  void setupIdentity
  (vectorial::BilinearForm<F,typename concepts::Realtype<F>::type>& bf,
   const concepts::ElementFormulaContainer<F> 
   frm = concepts::ElementFormulaContainer<F>());

} // namespace hp2D

#endif // hp2dbf_identity_hh