hp2Dedge/linearForm.hh

Go to the documentation of this file.

/* Linear forms for hp 2D edgeFEM
 */

#ifndef hp2dEdgeLinearform_hh
#define hp2dEdgeLinearform_hh

#include <memory>
#include "basics/typedefs.hh"
#include "formula/boundary.hh"
#include "formula/elementFormulaContainer.hh"
#include "function/linearForm.hh"
#include "hp2D/formula.hh"
#include "hp2Dedge/quad.hh"

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

  class BoundaryConditions;

  template<class F>
  class Element;

  template<class F>
  class ElementMatrix;

  template<class F>
  class Array;
}

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

  using concepts::Real;

  // *************************************************************** Neumann **

  class Neumann : public concepts::Neumann<Real> {
  public:
    Neumann(const concepts::BoundaryConditions* bc);
    virtual void operator()(const concepts::Element<Real>& elm,
                            concepts::ElementMatrix<Real>& em);

    virtual void operator()(const Quad<Real>& elm,
                            concepts::ElementMatrix<Real>& em);
  private:
    virtual void part(const Edge<Real>& edge, const concepts::Boundary bd,
                      concepts::ElementMatrix<Real>& em);
  }; 

  // ***************************************************************** Riesz **

  class Riesz : public concepts::LinearForm<Real> {
  public:
    Riesz(const concepts::ElementFormulaContainer<Real> frm1,
          const concepts::ElementFormulaContainer<Real> frm2,
          const concepts::BoundaryConditions* bc = 0);
    Riesz(const concepts::ElementFormulaContainer<concepts::Real2d> frm,
          const concepts::BoundaryConditions* bc = 0);
    virtual ~Riesz();

    void operator()(const concepts::Element<Real>& elm,
                    concepts::ElementMatrix<Real>& em);
    void operator()(const Quad<Real>& elm, concepts::ElementMatrix<Real>& em);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    concepts::ElementFormulaContainer<concepts::Real2d> frm_;
    concepts::Array<Real> tmp_;
    Neumann neumann_;
  };

  // ************************************************************* EdgeRiesz **

  class EdgeRiesz : public concepts::LinearForm<Real> {
  public:
    EdgeRiesz(const concepts::Boundary& bd); 
    virtual ~EdgeRiesz();

    void operator()(const concepts::Element<Real>& elm,
                    concepts::ElementMatrix<Real>& em);
    void operator()(const Edge<Real>& elm, concepts::ElementMatrix<Real>& em);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    const concepts::Boundary bd_;
  };

} // namespace hp2Dedge

#endif // hp2dEdgeLinearform_hh