bilinearForm1D.hh

Go to the documentation of this file.

 
#ifndef bilinearLinFEM1d_hh
#define bilinearLinFEM1d_hh
 
#include <memory>
#include "basics/typedefs.hh"
#include "basics/vectorsMatrices.hh"
#include "formula/formula.hh"
#include "operator/bilinearForm.hh"
#include "integration/quadrature.hh"
 
 
namespace linearFEM {
 
  using concepts::Real;
 
  class Line;
 
  // ************************************************************* Laplace1d **
 
  class Laplace1d : public concepts::BilinearForm<Real> {
  public:
    virtual void operator()(const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY,
          concepts::ElementMatrix<Real>& em) const;
    virtual Laplace1d* clone() const { return new Laplace1d(); }
  };
 
  // ************************************************************ Identity1d **
 
  class Identity1d : public concepts::BilinearForm<Real> {
  public:
    virtual void operator()(const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY,
          concepts::ElementMatrix<Real>& em) const;
    virtual Identity1d* clone() const { return new Identity1d(); }
  };
 
  // *********************************************************** CIdentity1d **
 
  class CIdentity1d : public concepts::BilinearForm<Real> {
  public:
    CIdentity1d(const concepts::Formula<Real>& frm,
                const uint gauss_p = 1) :
      frm_(frm.clone()), quad_(gauss_p) {}
    virtual void operator()(const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY,
          concepts::ElementMatrix<Real>& em) const;
    virtual CIdentity1d* clone() const {
      return new CIdentity1d(*frm_.get(), quad_.n()); }
  private:
    std::unique_ptr<const concepts::Formula<Real> > frm_;
    const concepts::Quadrature<4> quad_;
  };
 
} // namespace linearFEM
 
#endif // bilinearLinFEM1d_hh