linearForm1D.hh

Go to the documentation of this file.

 
#pragma once
 
#include "basics/typedefs.hh"
#include "basics/vectorsMatrices.hh"
#include "toolbox/array.hh"
#include "formula/formula.hh"
#include "function/linearForm.hh"
#include "integration/quadrature.hh"
 
namespace linearFEM {
 
  using concepts::Real;
 
  // *************************************************************** Riesz1d **
 
  class Riesz1d : public concepts::LinearForm<Real> {
  public:
    Riesz1d(const concepts::Formula<Real>& frm, const uint gauss_p = 1) :
      frm_(frm.clone()), quad_(gauss_p) {}
    virtual void operator()(const concepts::Element<Real>& elm,
          concepts::ElementMatrix<Real>& em) const;
  private:
    std::unique_ptr<const concepts::Formula<Real> > frm_;
    const concepts::Quadrature<4> quad_;
  };
 
} // namespace linearFEM