spline_1d_interpolation.hpp

Go to the documentation of this file.

#ifndef SPLINE_1D_INTERPOLATION_HPP
#define SPLINE_1D_INTERPOLATION_HPP

#include "uniform_range.hpp"
#include "binary_search.hpp"
template<class Device>
struct SplinePrecomputations{
    View<double*, CLayout, Device> dx;
    View<double*, CLayout, Device> dx_inv;
    View<double*, CLayout, Device> l_inv;
    View<double*, CLayout, Device> mu;

    View<double*, CLayout, Device> b;
    View<double*, CLayout, Device> c;
    View<double*, CLayout, Device> d;

    View<int*, CLayout, Device> xmap;

    SplinePrecomputations(const View<double*, CLayout, Device>& x, const UniformRange& range, double period);

    SplinePrecomputations(){}

    static int get_np(int n, double period){
        bool periodic = (period>0.0);
        constexpr int N_WRAP = 3;
        int offset = (periodic ? N_WRAP : 0);
        return (offset + n + offset);
    }
};

template<class Device>
void spline_interpolation_to_uniform(const SplinePrecomputations<Device>& p, const View<double*, CLayout, Device>& x0, const View<double*, CLayout, Device>& y0, const UniformRange& range, const View<double*, CLayout, Device>& y, double period);


template<class Device>
View<double*, CLayout, Device> linear_interpolation_to_uniform(const View<double*, CLayout, Device>& x0, const View<double*, CLayout, Device>& y0, const UniformRange& range, const View<double*, CLayout, Device>& y, double period);
#endif