formulaLiCo.hh

Go to the documentation of this file.

 
#ifndef formulaLiCo_hh
#define formulaLiCo_hh
 
#include "formula.hh"
 
namespace concepts {
 
  class FormulaLiCo : public Formula<Real> {
  public:
    FormulaLiCo(const Formula<Real>& frmA, const Formula<Real>& frmB, Real cA, Real cB)
      : frmA_(frmA), frmB_(frmB), cA_(cA), cB_(cB) { }
    virtual Formula<Real> * clone () const {
      return new FormulaLiCo(frmA_, frmB_, cA_, cB_);
    }
 
    virtual Real operator() (const Real p, const Real t=0.0) const {
      return cA_*frmA_(p, t) + cB_*frmB_(p, t);
    }
    virtual Real operator() (const Real2d &p, const Real t=0.0) const {
      return cA_*frmA_(p, t) + cB_*frmB_(p, t);
    }
    virtual Real operator() (const Real3d &p, const Real t=0.0) const {
      return cA_*frmA_(p, t) + cB_*frmB_(p, t);
    }
  private:
    const Formula<Real>& frmA_;
    const Formula<Real>& frmB_;
    const Real cA_;
    const Real cB_;
  };
 
} // namespace concepts
 
#endif //formulaLiCo_hh