piecewiseFormula.hh

Go to the documentation of this file.

 
#ifndef piecewiseformula_hh
#define piecewiseformula_hh
 
#include "space/element.hh"
#include "elementFormula.hh"
#include "basics/output.hh"
 
 
namespace concepts {
 
  // ************************************************** PiecewiseFormulaBase **
 
  template<typename F>
  class PiecewiseFormulaBase : 
    public ElementFormula<F, typename concepts::Realtype<F>::type> {
  public:
    typedef typename concepts::Realtype<F>::type G; 
 
    PiecewiseFormulaBase(const F dflt_value = (F)0);
    virtual ~PiecewiseFormulaBase();
    virtual PiecewiseFormulaBase<F>* clone() const = 0;
 
    virtual F operator()(const ElementWithCell<G>& elm, const Real p,
                         const Real t = 0.0) const 
    {
      return operator()(elm.cell().connector(), elm.elemMap(p), t);
    }
    virtual F operator()(const ElementWithCell<G>& elm, const Real2d& p,
                         const Real t = 0.0) const
    {
      return operator()(elm.cell().connector(), elm.elemMap(p), t);
    }
 
    virtual F operator()(const ElementWithCell<G>& elm, const Real3d& p,
                         const Real t = 0.0) const
    {
      return operator()(elm.cell().connector(), elm.elemMap(p), t);
    }
 
    virtual F operator()(const Connector& cntr, const Real p,
                         const Real t = 0.0) const = 0;
    virtual F operator()(const Connector& cntr, const Real2d& p,
                         const Real t = 0.0) const = 0;
    virtual F operator()(const Connector& cntr, const Real3d& p,
                         const Real t = 0.0) const = 0;
    inline virtual const F& dflt_value() const { return dflt_value_;}
    inline virtual F& dflt_value() { return dflt_value_;}
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    F dflt_value_;
  };
 
  template<typename F>
  PiecewiseFormulaBase<F>::PiecewiseFormulaBase(const F dflt_value) : 
    dflt_value_(dflt_value) {
  }
  
  template<typename F>
  PiecewiseFormulaBase<F>::~PiecewiseFormulaBase() {}
 
  template<typename F>
  std::ostream& PiecewiseFormulaBase<F>::info(std::ostream& os) const {
    os << concepts::typeOf(*this) << "(default value = " << dflt_value()  << ")";
    return os;
  }
 
} // namespace concepts
 
#endif // piecewiseformula_hh