formula/frm_sum.hh

Go to the documentation of this file.

#ifndef sumformula_hh
#define sumformula_hh

#include "elementFormulaRCP.hh"

namespace concepts {

 
 // ************************************************************ Frm_Sum **

  template<class F, class H = F, class J = F>
  class Frm_Sum : public Formula<F> {
  public:
    Frm_Sum(const RCP<const Formula<H> > frm1,
            const RCP<const Formula<J> > frm2) 
      : frm1_(frm1), frm2_(frm2) 
    {
      conceptsAssert(frm1.get(), Assertion());
      conceptsAssert(frm2.get(), Assertion());
    }
    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 << "Frm_Sum(" << *frm1_ << " + " << *frm2_ << ")";
    }
  private:
    RCP<const Formula<H> > frm1_;
    RCP<const Formula<J> > frm2_;
  };


  // ************************************************************ operator+ **

//   Frm_Sum<Real> 
//   operator+(const RCP<const Formula<Real> >  frm1, const RCP<const Formula<Real> > frm2);

//   Frm_Sum<Cmplx,Cmplx,Real> 
//   operator+(const RCP<const Formula<Cmplx> > frm1, const RCP<const Formula<Real> > frm2);

//   Frm_Sum<Cmplx,Real,Cmplx> 
//   operator+(const RCP<const Formula<Real> >  frm1, const RCP<const Formula<Cmplx> > frm2);

//   Frm_Sum<Cmplx,Cmplx,Cmplx> 
//   operator+(const RCP<const Formula<Cmplx> > frm1, const RCP<const Formula<Cmplx> > frm2);


} // namespace concepts


#endif // sumformula_hh