formula/piecewiseFormula.hh

Go to the documentation of this file.

/* piecewise formula class
*/

#ifndef piecewiseformula_hh
#define piecewiseformula_hh

#include "elementFormula.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 << "PiecewiseFormulaBase(default value = " << dflt_value()  << ")";
    return os;
  }

} // namespace concepts

#endif // piecewiseformula_hh