interpolation.hh

Go to the documentation of this file.

#pragma once
 
// standard library
#include <memory>
#include <vector>
 
// external libraries
#include <boost/math/interpolators/barycentric_rational.hpp>
#include <gsl/gsl_spline.h>
 
namespace DiFfRG
{
  namespace Interpolation
  {
    class Barycentric
    {
    public:
      Barycentric(std::vector<double> &x, std::vector<double> &y, uint order = 1, bool copy = true);
 
      Barycentric(const std::vector<double> &x, const std::vector<double> &y, uint order = 1);
 
      Barycentric();
 
      bool check_consistency(double tolerance) const;
 
      double operator()(double x) const;
 
      double value(double x) const;
 
      double derivative(double x) const;
 
    private:
      std::unique_ptr<boost::math::interpolators::barycentric_rational<double>> interpolator;
    };
 
    class CubicSpline
    {
    public:
      CubicSpline(const std::vector<double> &x, const std::vector<double> &y);
 
      ~CubicSpline();
 
      double operator()(double x) const;
 
      double value(double x) const;
 
      double derivative(double x) const;
 
    private:
      gsl_interp_accel *acc;
      gsl_spline *spline;
    };
  } // namespace Interpolation
} // namespace DiFfRG