frmE_inverse.hh

Go to the documentation of this file.

 
#ifndef inverseelementformula_hh
#define inverseelementformula_hh
 
#include "basics/debug.hh"
//#include "elementFormulaContainer.hh"
#include "exceptions.hh"
 
#define InvElemFormConstr_D 0
#define InvElemFormDestr_D  0
#define InvElemFormClone_D  0
 
namespace concepts {
 
  // forward declaration
  template<class F, class G>
  class ElementFormulaContainer;
 
 
  // *********************************************************** FrmE_Inverse **
 
  template<typename F, typename G = typename Realtype<F>::type>
  class FrmE_Inverse : public ElementFormula<F,G> {
  public:
    FrmE_Inverse(const ElementFormulaContainer<F,G> frm) 
      : frm_(frm)
    {
      if (frm.empty()) throw(PointerToEmptyElementFormula());
    }
    FrmE_Inverse(const FrmE_Inverse<F,G>& frm) 
      : frm_(frm.frm_) {}
 
    virtual ~FrmE_Inverse() {
      DEBUGL(InvElemFormDestr_D, "Delete " << *this)
    }
    
    virtual F operator() (const ElementWithCell<G>& elm, 
                          const Real p,    const Real t = 0.0) const {
      F v = frm_(elm, p,t);
 
      return 1.0 / v;
    }
    virtual F operator() (const ElementWithCell<G>& elm, 
                          const Real2d& p, const Real t = 0.0) const {
      F v = frm_(elm, p,t);
      
      return 1.0 / v;
    }
    virtual F operator() (const ElementWithCell<G>& elm, 
                          const Real3d& p, const Real t = 0.0) const {
      F v = frm_(elm, p,t);
      
      return 1.0 / v;
    }
    virtual FrmE_Inverse<F,G>* clone() const {
      return new FrmE_Inverse<F,G>(frm_);
    }
  protected:
    virtual std::ostream& info(std::ostream& os) const {
      return os << concepts::typeOf(*this) << "(1 / " << frm_ << ')';
    }
  private:
    ElementFormulaContainer<F,G> frm_;
  };
 
} // namespace concepts
 
 
#endif // inverseelementformula_hh