boundaries.hh

Go to the documentation of this file.

 
#ifndef hp2dEdgeBoundaries_hh
#define hp2dEdgeBoundaries_hh
 
#include <memory>
#include "basics/typedefs.hh"
#include "function/linearForm.hh"
#include "space.hh"
#include "space/dirichlet.hh"
#include "hp2Dedge/bilinearForm.hh"
#include "hp2Dedge/quad.hh"
 
 
namespace concepts {
  // forward declarations
  template<class F>
  class Formula;
 
  class BoundaryConditions;
 
  template<class F>
  class Element;
 
  template<class F>
  class ElementMatrix;
 
  template<class F>
  class Array;
}
 
namespace hp2Dedge {
  // forward declarations
  template<class F>
  class Quad;
 
  using concepts::Real;
 
  // ************************************************* DirichletCoefficients **
 
  class DirichletCoefficients : public concepts::Dirichlet<Real> {
  public:
    DirichletCoefficients(Space& spc);
    ~DirichletCoefficients();
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  };
 
  // ******************************************************* DirichletRotRot **
 
  class DirichletRotRot : public concepts::LinearForm<Real> {
  public:
    DirichletRotRot(DirichletCoefficients* dc = 0);
 
    virtual ~DirichletRotRot();
 
    void operator()(const concepts::Element<Real>& elm,
                    concepts::ElementMatrix<Real>& em) const;
    void operator()(const Quad<>& elm, concepts::ElementMatrix<Real>& em) const;
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    concepts::Array<Real> jacobianReciprocal_, function_;
    mutable concepts::Array<Real> tmp_;
    DirichletCoefficients* dc_;
  };
 
  // ***************************************************** DirichletIdentity **
 
  class DirichletIdentity : public concepts::LinearForm<Real> {
  public:
    DirichletIdentity(DirichletCoefficients* dc = 0);
 
    virtual ~DirichletIdentity();
 
    void operator()(const concepts::Element<Real>& elm,
        concepts::ElementMatrix<Real>& em) const;
    void operator()(const Quad<>& elm, concepts::ElementMatrix<Real>& em) const;
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    mutable concepts::Array<Real> jacobian_, function_;
    concepts::Array<concepts::MapReal2d> jacobianInv_;
    concepts::Array<Real> tmp_;
    DirichletCoefficients* dc_;
  };
 
} // namespace hp2Dedge
 
#endif // hp2dEdgeBoundaries_hh