formula/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 << "FrmE_Inverse(1 / " << frm_ << ')';
    }
  private:
    ElementFormulaContainer<F,G> frm_;
  };

} // namespace concepts


#endif // inverseelementformula_hh