sourceFunction.hh

#ifndef _sourceFunction_hh
#define _sourceFunction_hh
 
#include "formula.hh"
#include "function/elementFormula.hh"
#include "hp2D.hh"
#include <math.h>
//#include "auxiliaryFunctions.hh"
 
#define PI 3.14159265
 
using concepts::Real;
using concepts::Real2d;
using concepts::Real3d;
using concepts::Cmplx;
using concepts::besselJnu;
using concepts::besselYnu;
 
namespace concepts{
    
 
 
    
class SourceFunctionF0_x : public ElementFormula<Cmplx>{
  public:
    SourceFunctionF0_x( Real lambda, Real kn, Real R1, Real coef=1.0 ):
            lam_(lambda), kn_(kn), R1_(R1), C_(coef){};
    virtual Cmplx operator() (const ElementWithCell<Real>& elm, const Real p,
                             const Real t = 0.0) const;
    virtual Cmplx operator() (const ElementWithCell<Real>& elm, const Real2d& p,
                             const Real t = 0.0) const;
    virtual Cmplx operator() (const ElementWithCell<Real>& elm, const Real3d& p,
                             const Real t = 0.0) const;
    virtual SourceFunctionF0_x* clone() const;
 
 
    virtual ~SourceFunctionF0_x() {} ;
 
  private:
    Real lam_;
    Real kn_;
    Real R1_;
    Real C_;
};
 
 
class SourceFunctionF0_y : public ElementFormula<Cmplx>{
  public:
    SourceFunctionF0_y( Real lambda, Real kn, Real R1, Real coef=1.0 ):
            lam_(lambda), kn_(kn), R1_(R1), C_(coef) {};
    virtual Cmplx operator() (const ElementWithCell<Real>& elm, const Real p,
                             const Real t = 0.0) const;
    virtual Cmplx operator() (const ElementWithCell<Real>& elm, const Real2d& p,
                             const Real t = 0.0) const;
    virtual Cmplx operator() (const ElementWithCell<Real>& elm, const Real3d& p,
                             const Real t = 0.0) const;
    virtual SourceFunctionF0_y* clone() const;
 
 
    virtual ~SourceFunctionF0_y() {} ;
 
  private:
    Real lam_;
    Real kn_;
    Real R1_;
    Real C_;
};
 
 
}
 
 
#endif //_sourceFunction_hh