ogataquadrature.h

Go to the documentation of this file.

//
// APFEL++ 2017
//
// Author: Valerio Bertone: valerio.bertone@cern.ch
//
 
#pragma once
 
#include <vector>
#include <functional>
 
namespace apfel
{
  class OgataQuadrature
  {
  public:
    OgataQuadrature(int    const& nu = 0,
                    double const& CutOff = 1e-5,
                    double const& h = 0.001,
                    int    const& nZeroMax = 1000);
 
    void InitialiseWeights(int const& nZeroMax);
 
    template<typename T>
    T transform(std::function<T(double const&)> const& func, double const& qT, bool const& Dynh = true, int const& nmax = 1000, int const& period = 10) const;
 
    double GetBesselOrder() const { return _nu; }
 
    double GetCutOff() const { return _CutOff; }
 
    double GetStepParameter() const { return _h; }
 
    std::vector<double> GetCoordinates() const { return _xf; }
 
    std::vector<double> GetWeights() const { return _weights; }
 
    void SetStepParameter(double const& h) { _h = h; }
 
    void JnuZerosGenerator(int const& nu) const;
 
  private:
    int                 const _nu;       
    double              const _CutOff;   
    double                    _h;        
    int                 const _nZeroMax; 
    std::vector<double>       _xf;       
    std::vector<double>       _weights;  
  };
}