frmE_conjugate.hh

Go to the documentation of this file.

 
#ifndef FORMULA_FRME_CONJUGATE_HH_
#define FORMULA_FRME_CONJUGATE_HH_
 
#include "basics/debug.hh"
#include "basics/output.hh"
#include "toolbox/sharedPointer.hh"
#include "elementFormulaContainer.hh"
#include "exceptions.hh"
 
#define ConjElemFormDestr_D 0
#define ConjElemFormClone_D 0
 
namespace concepts {
 
 
  // *********************************************************** FrmE_Product **
 
  template<typename F,
           typename G = typename Realtype<F>::type>
  class FrmE_Conjugate : public ElementFormula<F,G> {
  public:
    FrmE_Conjugate(const ElementFormulaContainer<F,G> frm) : frm_(frm) {
      if (frm.empty()) throw(PointerToEmptyElementFormula());
    }
 
    virtual ~FrmE_Conjugate() {
      DEBUGL(ConjElemFormDestr_D, "Delete " << *this)
    }
    virtual F operator() (const ElementWithCell<G>& elm,
                          const Real p,    const Real t = 0.0) const
    {
      return std::conj(frm_(elm, p, t));
    }
    virtual F operator() (const ElementWithCell<G>& elm,
                          const Real2d& p, const Real t = 0.0) const
    {
      return std::conj(frm_(elm, p, t));
    }
    virtual F operator() (const ElementWithCell<G>& elm,
                          const Real3d& p, const Real t = 0.0) const
    {
      return std::conj(frm_(elm, p, t));
    }
    virtual FrmE_Conjugate<F,G>* clone() const {
      DEBUGL(ConjElemFormClone_D, *this);
      return new FrmE_Conjugate<F,G>(frm_);
    }
  protected:
    virtual std::ostream& info(std::ostream& os) const {
      return os << concepts::typeOf(*this) << "(" << frm_ << ")";
    }
  private:
    ElementFormulaContainer<F,G> frm_;
  };
 
} // namespace concepts
 
namespace std {
 
  // ************************************************************ std::conj **
 
  template<class F>
  concepts::ElementFormulaContainer<F>
  conj(const concepts::ElementFormulaContainer<F> frm)
  {
    return FrmE_Conjugate(frm);
  }
 
} // namespace std
 
#endif /* FORMULA_FRME_CONJUGATE_HH_ */