qgrid.h

Go to the documentation of this file.

//
// APFEL++ 2017
//
// Author: Valerio Bertone: valerio.bertone@cern.ch
//
 
#pragma once
 
#include <vector>
#include <tuple>
#include <iostream>
#include <functional>
 
namespace apfel
{
  template<class T>
  class QGrid
  {
  public:
    QGrid() = delete;
 
    QGrid(int                                  const& nQ,
          double                               const& QMin,
          double                               const& QMax,
          int                                  const& InterDegree,
          std::vector<double>                  const& Thresholds,
          std::function<double(double const&)> const& TabFunc,
          std::function<double(double const&)> const& InvTabFunc);
 
    QGrid(int                 const& nQ,
          double              const& QMin,
          double              const& QMax,
          int                 const& InterDegree,
          std::vector<double> const& Thresholds,
          double              const& Lambda = 0.25);
 
    QGrid(std::vector<double> const& Qg,
          int                 const& InterDegree);
 
    T Evaluate(double const& Q) const;
 
    T Derive(double const& Q) const;
 
    T Integrate(double const& Qa, double const& Qb) const;
 
    bool operator == (QGrid const& sg) const;
    bool operator != (QGrid const& sg) const;
 
    int                                         nQ()                 const { return _nQ; }          
    int                                         InterDegree()        const { return _InterDegree; } 
    double                                      QMin()               const { return _QMin; }        
    double                                      QMax()               const { return _QMax; }        
    std::function<double(double const&)> const& TabFunc()            const { return _TabFunc; }     
    std::vector<double>                  const& GetThresholds()      const { return _Thresholds;}   
    std::vector<double>                  const& GetQGrid()           const { return _Qg; }          
    std::vector<double>                  const& GetFQGrid()          const { return _fQg; }         
    std::vector<int>                     const& GetThesholdIndices() const { return _nQg; }         
    std::vector<T>                       const& GetQGridValues()     const { return _GridValues; }  
 
    double Interpolant(int const& tQ, int const& tau, double const& fq) const;
 
    double DerInterpolant(int const& tQ, int const& tau, double const& Q) const;
 
    double IntInterpolant(int const& tQ, int const& tau, double const& Qa, double const& Qb) const;
 
    std::tuple<int, int, int> SumBounds(double const& Q) const;
 
    void Print() const { std::cout << *this << std::endl; }
 
  protected:
    int                                  _nQ;           
    double                               _QMin;         
    double                               _QMax;         
    int                                  _InterDegree;  
    std::vector<double>                  _Thresholds;   
    std::function<double(double const&)> _TabFunc;      
    std::vector<double>                  _Qg;           
    std::vector<double>                  _fQg;          
    std::vector<int>                     _nQg;          
    std::vector<T>                       _GridValues;   
 
    template<class U>
    friend std::ostream& operator << (std::ostream& os, QGrid<U> const& dt);
  };
 
  template<class T>
  inline std::ostream& operator << (std::ostream& os, QGrid<T> const& Qg)
  {
    os << "QGrid: " << &Qg << "\n";
    os << "nQ                = " << Qg._nQ << "\n";
    os << "QMin              = " << Qg._QMin << "\n";
    os << "QMax              = " << Qg._QMax << "\n";
    os << "InterDegree       = " << Qg._InterDegree << "\n";
    os << "Thresholds        = ";
    for (const auto &v: Qg._Thresholds) os << v << " ";
    os << "\n";
    os << "Threshold indices = ";
    for (const auto &v: Qg._nQg) os << v << " ";
    os << "\n";
    os << "Qg                = ";
    for (const auto &v: Qg._Qg) os << v << " ";
    os << "\n";
    return os;
  }
}