interpolate.h

Go to the documentation of this file.

 
#pragma once
#include <vector>
 
namespace LIBRPA
{
namespace utils
{
 
class CubicSpline
{
private:
    const std::vector<double> d_xs;
    const std::vector<double> d_ys;
    std::vector<double> d_a;
    std::vector<double> d_b;
    std::vector<double> d_c;
    std::vector<double> d_d;
public:
    // default constructor
    CubicSpline(const std::vector<double>& xs, const std::vector<double>& ys);
    // destructor
    ~CubicSpline() {};
    // copy & move constructor
    CubicSpline(const CubicSpline& s) = delete;
    CubicSpline(CubicSpline&& s) = delete;
    // copy & move assignment operator
    CubicSpline& operator=(const CubicSpline& s) = delete;
    CubicSpline& operator=(CubicSpline&& s) = delete;
 
    double operator()(double x) const;
    std::vector<double> operator()(const std::vector<double>& xs) const;
}; /* end of class CubicSpline */
 
std::vector<double> interp_cubic_spline(const std::vector<double> &xs,
                                        const std::vector<double> &ys,
                                        const std::vector<double> &xs_new);
 
} /* end of namepsace utils */
 
} /* end of namespace LIBRPA */