function.hh

Go to the documentation of this file.

 
#ifndef hpFunction2dedge_hh
#define hpFunction2dedge_hh
 
#include "hp2D/formula.hh"
#include "hp2D/functionBase.hh"
#include "quad.hh"
 
namespace hp2Dedge {
 
  using concepts::Real;
  using concepts::Real2d;
 
  // ***************************************************************** Value **
 
  template<typename F = Real, typename G = Real,
           typename H = typename concepts::Realtype<F>::type>
  class Value : public hp2D::ElementFunction<2, F, Quad<H> > {
  public:
    Value();
    Value(const concepts::ElementFormula<G,H>& frm);
    virtual Value<F,G,H>* clone() const {
      if (f_.get()) return new Value(*f_.get()); 
      return new Value();       
    }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    std::unique_ptr<const concepts::ElementFormula<G,H> > f_;
    virtual void compute_(const Quad<H>& elm, const concepts::Array<F>& coeff,
        concepts::Array<F>& val, const uint i[2]) const;
    virtual void compute_(const Quad<H>& elm, const uint j[3],
        concepts::Array<F>& val, const uint i[2]) const;
    virtual void compute_(const Quad<H>& elm, const concepts::Array<F>& coeff,
        concepts::Array<F>& val,
        const Real2d& p, const Real t = 0.0) const;
  };
 
  // ******************************************************************* Rot **
 
  template<typename F = Real, typename G = Real>
  class Rot : public  hp2D::ElementFunction<1, F, Quad<G> > {
  public:
    virtual Rot<F,G>* clone() const { return new Rot(); }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    virtual void compute_(const Quad<>& elm, const concepts::Array<F>& coeff,
        F& val, const uint i[2]) const;
    virtual void compute_(const Quad<>& elm, const uint j[3],
        F& val, const uint i[2]) const;
    virtual void compute_(const Quad<>& elm, const concepts::Array<F>& coeff,
        F& val,
        const Real2d& p, const Real t = 0.0) const;
  };
 
} // namespace hp2Dedge
 
 
#endif // hpRealunction2dedge_hh