frm_product.hh

Go to the documentation of this file.

 
#ifndef productformula_hh
#define productformula_hh
 
#include "basics/typedefs.hh"
#include "basics/output.hh"
#include "basics/vectorsMatricesForward.hh"
#include "formulaContainer.hh"
#include "exceptions.hh"
 
namespace concepts {
 
 
 // ************************************************************ Frm_Product **
 
  template<class F, class H = F, class J = F>
  class Frm_Product : public Formula<F> {
  public:
    Frm_Product(const FormulaContainer<H> frm1,
                const FormulaContainer<J> frm2)
      : frm1_(frm1), frm2_(frm2) 
    {
      if (frm1.empty()) throw(PointerToEmptyElementFormula());
      if (frm2.empty()) throw(PointerToEmptyElementFormula());
    }
    virtual F operator() (const Real p,    const Real t = 0.0) const {
      return frm1_(p,t) * frm2_(p,t);
    }
    virtual F operator() (const Real2d& p, const Real t = 0.0) const {
      return frm1_(p,t) * frm2_(p,t);
    }
    virtual F operator() (const Real3d& p, const Real t = 0.0) const {
      return frm1_(p,t) * frm2_(p,t);
    }
    virtual Frm_Product<F,H,J>* clone() const {
      return new Frm_Product<F,H,J>(frm1_, frm2_);
    }
  protected:
    virtual std::ostream& info(std::ostream& os) const {
      return os << concepts::typeOf(*this)<< "(" << frm1_ << " * " << frm2_ << ")";
    }
  private:
    const FormulaContainer<H> frm1_;
    const FormulaContainer<J> frm2_;
  };
 
  // ************************************************************ operator* **
 
  Frm_Product<Real> operator*(const Formula<Real>& frm, const Real a);
 
  Frm_Product<Cmplx,Real> operator*(const Formula<Real>& frm, const Cmplx a);
 
  Frm_Product<Cmplx,Cmplx,Real> operator*(const Formula<Cmplx>& frm, const Real a);
 
  Frm_Product<Cmplx> operator*(const Formula<Cmplx>& frm, const Cmplx a);
 
  //                             ************
 
  Frm_Product<Real> operator*(const Formula<Real>& frm1, const Formula<Real>& frm2);
 
  Frm_Product<Cmplx,Real> operator*(const Formula<Real>& frm1, const Formula<Cmplx>& frm2);
 
  Frm_Product<Cmplx,Cmplx,Real> operator*(const Formula<Cmplx>& frm1, const Formula<Real>& frm2);
 
  Frm_Product<Cmplx> operator*(const Formula<Cmplx>& frm1, const Formula<Cmplx>& frm2);
 
} // namespace concepts
 
 
#endif // productformula_hh