bf_advection.hh

Go to the documentation of this file.

 
#ifndef HP2DEDGE_BF_ADVECTION_HH_
#define HP2DEDGE_BF_ADVECTION_HH_
 
#include <memory>
#include "operator/bilinearForm.hh"
#include "basics/typedefs.hh"
#include "basics/vectorsMatrices.hh"
#include "basics/vectorsMatricesForward.hh"
#include "formula/elementFormula.hh"
 
//#include "formula/boundary.hh"
#include "hp2D/bilinearFormHelper.hh"
/*
#include "geometry/sharedJacobian.hh"
#include "hp2D/formula.hh"
#include "hp2D/arrayElementFormula.hh"
#include "toolbox/sharedPointer.hh"
*/
 
namespace concepts
{
// forward declarations
  template<class F>
  class Element;
 
  template<class F>
  class ElementMatrix;
}
 
namespace hp2Dedge
{
  // forward declarations
  template<class F>
  class Quad;
 
  using concepts::Real;
  using concepts::Cmplx;
 
  // ************************************************************* Advection **
 
  template<class F = Real>
  class Advection : public concepts::BilinearForm<F,Real>,
                    public hp2D::BilinearFormHelper_2_1<F>
  {
  public:
    Advection(const concepts::ElementFormulaContainer<F> frm1,
              const concepts::ElementFormulaContainer<F> frm2)
      : hp2D::BilinearFormHelper_2_1<F>(frm1, frm2)
    {}
 
    Advection
    (const concepts::ElementFormulaContainer<concepts::Point<F, 2> > frm)
      : hp2D::BilinearFormHelper_2_1<F>(frm)
    {}
 
    virtual ~Advection() {}
 
    virtual Advection<F>* clone() const {
      return new Advection<F>(this->frm_);
    }
 
    virtual void operator()(const concepts::Element<Real>& elmX,
                            const concepts::Element<Real>& elmY,
                            concepts::ElementMatrix<F>& em) const;
 
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    void operator()(const Quad<Real>& elmX,
                    const Quad<Real>& elmY,
                    concepts::ElementMatrix<F>& em) const;
  };
 
} // namespace concepts
 
#endif /* HP2DEDGE_BF_ADVECTION_HH_ */