linearForm.hh

Go to the documentation of this file.

 
#pragma once
 
#include <memory>
#include "basics/typedefs.hh"
#include "function/linearForm.hh"
#include "formula.hh"
 
namespace concepts {
  // forward declarations
  template<class F>
  class Formula;
 
  class BoundaryConditions;
}
 
namespace hp1D {
 
  using concepts::Real;
 
  // *************************************************************** Neumann **
 
  class Neumann : public concepts::Neumann<Real> {
  public:
    Neumann(const concepts::BoundaryConditions* bc); 
    virtual ~Neumann();
 
    void operator()(const concepts::Element<Real>& elm,
                    concepts::ElementMatrix<Real>& em) const;
  }; 
 
  // ***************************************************************** Riesz **
 
  template<class F = Real>
  class Riesz : public concepts::LinearForm<F>, public LinearFormHelper<0,F>, public concepts::LinearFormChoice {
  public:
    Riesz(const concepts::ElementFormulaContainer<F> frm,
          const concepts::BoundaryConditions* bc = 0); 
    virtual ~Riesz();
 
    void operator()(const concepts::Element<Real>& elm,
                    concepts::ElementMatrix<F>& em) const;
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    Neumann neumann_;
 
  };
 
 
  // ********************************************** LinearFormGradInterp_Grad **
 
  template<typename F = Real>
  class LinearFormGradInterp_Grad : public concepts::LinearForm<F>,
                                    public LinearFormHelper<1,F> {
  public:
    LinearFormGradInterp_Grad(const concepts::ElementFormulaContainer<F> frm); 
    virtual ~LinearFormGradInterp_Grad();
 
    void operator()(const concepts::Element<Real>& elm,
                    concepts::ElementMatrix<F>& em) const;
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  };
 
  // ********************************************** GradLinearForm **
 
  template<typename F = Real>
  class GradLinearForm : public concepts::LinearForm<F>,
                         public LinearFormHelper<1,F> {
  public:
    GradLinearForm (const concepts::ElementFormulaContainer<F> frm); 
    virtual ~GradLinearForm();
 
    void operator()(const concepts::Element<Real>& elm,
                    concepts::ElementMatrix<F>& em) const;
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  };
 
  // ********************************************** PointLinearForm **
 
  template<typename F = Real>
  class PointEvaluation : public concepts::LinearForm<F> {
  public:
    PointEvaluation(const concepts::Set<concepts::Attribute> attrib,
                    const concepts::ElementFormulaContainer<F> frm = concepts::ConstFormula<F>(1.0))
      : attrib_(attrib), frm_(frm) {};
    void operator()(const concepts::Element<Real>& elm,
                    concepts::ElementMatrix<F>& em) const;
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    const concepts::Set<concepts::Attribute> attrib_;
    const concepts::ElementFormulaContainer<F> frm_;
 
  };
  
} // namespace hp1D