tracesEval.hh

Go to the documentation of this file.

 
#ifndef TRACESEVAL_HH
#define TRACESEVAL_HH
 
#include "function/elementFormula.hh"
#include "basics/vectorsMatricesForward.hh"
 
namespace concepts {
  // forward declarations
  template<class F>
  class ElementWithCell;
 
  template<class F, class G>
  class ElementFunction;
}
 
namespace hp3D {
 
using concepts::Real;
 
// forward declarations
class Space;
class TraceSpace;
template<class F>
class Element;
 
template<uint dim, class F = Real, class G = F>
class ElementFormulaVectorOnTrace :
        public concepts::ElementFormula<concepts::Point<G,dim> > {
 
public:
  ElementFormulaVectorOnTrace(const TraceSpace &tspc,
                              const Space &spc,
                              const concepts::Vector<F> &vec,
                              const concepts::ElementFunction<G, Real> &func);
 
  virtual concepts::Point<G,dim> operator()(const concepts::ElementWithCell<Real> &elm,
                                            const concepts::Real p,
                                            const concepts::Real t = 0.0) const;
  virtual concepts::Point<G,dim> operator()(const concepts::ElementWithCell<Real> &elm,
                                            const concepts::Real2d &p,
                                            const concepts::Real t = 0.0) const;
  virtual concepts::Point<G,dim> operator()(const concepts::ElementWithCell<Real> &elm,
                                            const concepts::Real3d &p,
                                            const concepts::Real t = 0.0) const;
 
  virtual ElementFormulaVectorOnTrace<dim,F,G>* clone() const
  {
      return new ElementFormulaVectorOnTrace(tspc_, spc_, vec_, *func_);
  }
protected:
  virtual std::ostream& info(std::ostream &os) const;
 
private:
 
  const Element<Real>& getUnderlyingElementAndIndex_(const concepts::ElementWithCell<Real> &elm,
                                                     uint &face) const;
 
  const TraceSpace &tspc_;
  const Space &spc_;
  const concepts::Vector<F> &vec_;
  std::unique_ptr<const concepts::ElementFunction<G, Real> > func_;
  concepts::ElementFormulaVector<dim, F, G, Real> frmE_;
 
};
 
} // namespace hp3D
 
#endif // TRACESEVAL_HH