bilinearForm2D.hh

Go to the documentation of this file.

 
#ifndef bilinearLinFEM_hh
#define bilinearLinFEM_hh
 
#include "basics/typedefs.hh"
#include "basics/vectorsMatrices.hh"
#include "operator/bilinearForm.hh"
#include "toolbox/array.hh"
 
namespace linearFEM {
 
  using concepts::Real;
 
  class Quad;
  class Triangle;
 
  // ************************************************************* Laplace2d **
 
  class Laplace2d : public concepts::BilinearForm<Real> {
  public:
    virtual void operator()(const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY,
          concepts::ElementMatrix<Real>& em) const;
    void operator()(const linearFEM::Quad& elmX,
        const linearFEM::Quad& elmY,
        concepts::ElementMatrix<Real>& em) const;
    void operator()(const linearFEM::Triangle& elmX,
        const linearFEM::Triangle& elmY,
        concepts::ElementMatrix<Real>& em) const;
    virtual Laplace2d* clone() const { return new Laplace2d(); }
  private:
    mutable concepts::Array<concepts::Real2d> shpfct_;
  };
 
  // ************************************************************ Identity2d **
 
  class Identity2d : public concepts::BilinearForm<Real> {
  public:
    virtual void operator()(const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY,
          concepts::ElementMatrix<Real>& em) const;
    void operator()(const linearFEM::Quad& elmX,
        const linearFEM::Quad& elmY,
        concepts::ElementMatrix<Real>& em) const;
    void operator()(const linearFEM::Triangle& elmX,
        const linearFEM::Triangle& elmY,
        concepts::ElementMatrix<Real>& em) const;
    virtual Identity2d* clone() const { return new Identity2d(); }
  private:
    mutable concepts::Array<Real> shpfct_;
  };
 
} // namespace linearFEM
 
#endif // bilinearLinFEM_hh