interpolator.h

Go to the documentation of this file.

//
// APFEL++ 2017
//
// Author: Valerio Bertone: valerio.bertone@cern.ch
//
 
#pragma once
 
#include "apfel/grid.h"
 
#include <array>
 
namespace apfel
{
  class Interpolator
  {
  public:
    virtual ~Interpolator() = default;
 
    Interpolator(Grid const& gr);
 
    Interpolator(Grid                             const& gr,
                 std::vector<std::vector<double>> const& distsubgrid,
                 std::vector<double>              const& distjointgrid);
 
 
    double Evaluate(double const& x) const;
 
    double Evaluate(double const& x, int const& ig) const;
 
    double Derive(double const& x) const;
 
    double Integrate(double const& a, double const& b) const;
 
    virtual double InterpolantLog(int const& beta, double const& lnx, SubGrid const& sg) const = 0;
 
    virtual double Interpolant(int const& beta, double const& x, SubGrid const& sg) const = 0;
 
    virtual double DerInterpolant(int const&, double const&, SubGrid const&) const { return 0; };
 
    virtual double IntInterpolant(int const&, double const&, double const&, SubGrid const&) const { return 0; };
 
    virtual std::array<int, 2> SumBounds(double const& x, SubGrid const& sg) const = 0;
 
    Grid                             const& GetGrid()                  const { return _grid; }                    
    std::vector<std::vector<double>> const& GetDistributionSubGrid()   const { return _distributionSubGrid; }     
    std::vector<double>              const& GetDistributionJointGrid() const { return _distributionJointGrid; }   
    void                                    Print()                    const { std::cout << *this << std::endl; } 
 
  protected:
    Grid                             const& _grid;                  
    std::vector<std::vector<double>>        _distributionSubGrid;   
    std::vector<double>                     _distributionJointGrid; 
 
    friend std::ostream& operator << (std::ostream& os, Interpolator const& sg);
  };
 
  std::ostream& operator << (std::ostream& os, Interpolator const& in);
}