linearForm2D.hh

Go to the documentation of this file.

 
#ifndef linearForm2D_hh
#define linearForm2D_hh
 
#include "basics/typedefs.hh"
#include "basics/vectorsMatrices.hh"
#include "toolbox/array.hh"
#include "formula/formula.hh"
#include "function/linearForm.hh"
 
namespace linearFEM {
 
  using concepts::Real;
 
  class Quad;
  class Triangle;
 
  // *************************************************************** Riesz2d **
 
  class Riesz2d : public concepts::LinearForm<Real> {
  public:
    Riesz2d(const concepts::Formula<Real>& frm) :
      frm_(frm.clone()), shpfct_(0) {}
    virtual void operator()(const concepts::Element<Real>& elm,
          concepts::ElementMatrix<Real>& em) const;
    void operator()(const linearFEM::Quad& elm,
        concepts::ElementMatrix<Real>& em) const;
    void operator()(const linearFEM::Triangle& elm,
        concepts::ElementMatrix<Real>& em) const;
  private:
    std::unique_ptr<const concepts::Formula<Real> > frm_;
    mutable concepts::Array<concepts::Real> shpfct_;
    static const Real weights_[4];
    static const concepts::Real2d coord_[4];
  };
 
} // namespace linearFEM
 
#endif // linearForm2D_hh