frm_sum.hh

Go to the documentation of this file.

 
#ifndef sumformula_hh
#define sumformula_hh
 
#include "basics/typedefs.hh"
#include "basics/output.hh"
#include "basics/vectorsMatricesForward.hh"
#include "formulaContainer.hh"
#include "exceptions.hh"
 
namespace concepts {
 
 
 // ************************************************************ Frm_Sum **
 
  template<class F, class H = F, class J = F>
  class Frm_Sum : public Formula<F> {
  public:
    Frm_Sum(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_Sum<F,H,J>* clone() const {
      return new Frm_Sum<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_Sum<Real> operator+(const Formula<Real>& frm, const Real a);
 
   Frm_Sum<Cmplx,Real> operator+(const Formula<Real>& frm, const Cmplx a);
 
   Frm_Sum<Cmplx,Cmplx,Real> operator+(const Formula<Cmplx>& frm, const Real a);
 
   Frm_Sum<Cmplx> operator+(const Formula<Cmplx>& frm, const Cmplx a);
 
   //                             ************
 
   Frm_Sum<Real> operator+(const Formula<Real>& frm1, const Formula<Real>& frm2);
 
   Frm_Sum<Cmplx,Real> operator+(const Formula<Real>& frm1, const Formula<Cmplx>& frm2);
 
   Frm_Sum<Cmplx,Cmplx,Real> operator+(const Formula<Cmplx>& frm1, const Formula<Real>& frm2);
 
   Frm_Sum<Cmplx> operator+(const Formula<Cmplx>& frm1, const Formula<Cmplx>& frm2);
 
 
} // namespace concepts
 
 
#endif // sumformula_hh