FM / interpolators

/usr/include/boost/math/interpolators

Name	Size	Mode	Actions
detail/	-	0755	rm
barycentric_rational.hpp	3805	0644	edit dl rm
cardinal_cubic_b_spline.hpp	3715	0644	edit dl rm
cardinal_quadratic_b_spline.hpp	2347	0644	edit dl rm
cardinal_quintic_b_spline.hpp	2675	0644	edit dl rm
cardinal_trigonometric.hpp	1443	0644	edit dl rm
catmull_rom.hpp	8767	0644	edit dl rm
cubic_b_spline.hpp	3665	0644	edit dl rm
cubic_hermite.hpp	3470	0644	edit dl rm
makima.hpp	5978	0644	edit dl rm
pchip.hpp	3720	0644	edit dl rm
quintic_hermite.hpp	3639	0644	edit dl rm
septic_hermite.hpp	3697	0644	edit dl rm
vector_barycentric_rational.hpp	2866	0644	edit dl rm
whittaker_shannon.hpp	1291	0644	edit dl rm

Edit: /usr/include/boost/math/interpolators/cardinal_quintic_b_spline.hpp (2675B)

#ifndef BOOST_MATH_INTERPOLATORS_CARDINAL_QUINTIC_B_SPLINE_HPP
#define BOOST_MATH_INTERPOLATORS_CARDINAL_QUINTIC_B_SPLINE_HPP
#include <memory>
#include <limits>
#include <boost/math/interpolators/detail/cardinal_quintic_b_spline_detail.hpp>

namespace boost{ namespace math{ namespace interpolators {

template <class Real>
class cardinal_quintic_b_spline
{
public:
    // If you don't know the value of the derivative at the endpoints, leave them as nans and the routine will estimate them.
    // y[0] = y(a), y[n - 1] = y(b), step_size = (b - a)/(n -1).
    cardinal_quintic_b_spline(const Real* const y,
                                size_t n,
                                Real t0 /* initial time, left endpoint */,
                                Real h  /*spacing, stepsize*/,
                                std::pair<Real, Real> left_endpoint_derivatives = {std::numeric_limits<Real>::quiet_NaN(), std::numeric_limits<Real>::quiet_NaN()},
                                std::pair<Real, Real> right_endpoint_derivatives = {std::numeric_limits<Real>::quiet_NaN(), std::numeric_limits<Real>::quiet_NaN()})
     : impl_(std::make_shared<detail::cardinal_quintic_b_spline_detail<Real>>(y, n, t0, h, left_endpoint_derivatives, right_endpoint_derivatives))
    {}

// Oh the bizarre error messages if we template this on a RandomAccessContainer:
    cardinal_quintic_b_spline(std::vector<Real> const & y,
                                Real t0 /* initial time, left endpoint */,
                                Real h  /*spacing, stepsize*/,
                                std::pair<Real, Real> left_endpoint_derivatives = {std::numeric_limits<Real>::quiet_NaN(), std::numeric_limits<Real>::quiet_NaN()},
                                std::pair<Real, Real> right_endpoint_derivatives = {std::numeric_limits<Real>::quiet_NaN(), std::numeric_limits<Real>::quiet_NaN()})
     : impl_(std::make_shared<detail::cardinal_quintic_b_spline_detail<Real>>(y.data(), y.size(), t0, h, left_endpoint_derivatives, right_endpoint_derivatives))
    {}

Real operator()(Real t) const {
        return impl_->operator()(t);
    }

Real prime(Real t) const {
       return impl_->prime(t);
    }

Real double_prime(Real t) const {
        return impl_->double_prime(t);
    }

Real t_max() const {
        return impl_->t_max();
    }

private:
    std::shared_ptr<detail::cardinal_quintic_b_spline_detail<Real>> impl_;
};

}}}
#endif

// Copyright Nick Thompson, 2019 // Use, modification and distribution are subject to the // Boost Software License, Version 1.0. // (See accompanying file LICENSE_1_0.txt // or copy at http://www.boost.org/LICENSE_1_0.txt) #ifndef BOOST_MATH_INTERPOLATORS_CARDINAL_QUINTIC_B_SPLINE_HPP #define BOOST_MATH_INTERPOLATORS_CARDINAL_QUINTIC_B_SPLINE_HPP #include #include #include namespace boost{ namespace math{ namespace interpolators { template class cardinal_quintic_b_spline { public: // If you don't know the value of the derivative at the endpoints, leave them as nans and the routine will estimate them. // y[0] = y(a), y[n - 1] = y(b), step_size = (b - a)/(n -1). cardinal_quintic_b_spline(const Real* const y, size_t n, Real t0 /* initial time, left endpoint */, Real h /*spacing, stepsize*/, std::pair left_endpoint_derivatives = {std::numeric_limits::quiet_NaN(), std::numeric_limits::quiet_NaN()}, std::pair right_endpoint_derivatives = {std::numeric_limits::quiet_NaN(), std::numeric_limits::quiet_NaN()}) : impl_(std::make_shared>(y, n, t0, h, left_endpoint_derivatives, right_endpoint_derivatives)) {} // Oh the bizarre error messages if we template this on a RandomAccessContainer: cardinal_quintic_b_spline(std::vector const & y, Real t0 /* initial time, left endpoint */, Real h /*spacing, stepsize*/, std::pair left_endpoint_derivatives = {std::numeric_limits::quiet_NaN(), std::numeric_limits::quiet_NaN()}, std::pair right_endpoint_derivatives = {std::numeric_limits::quiet_NaN(), std::numeric_limits::quiet_NaN()}) : impl_(std::make_shared>(y.data(), y.size(), t0, h, left_endpoint_derivatives, right_endpoint_derivatives)) {} Real operator()(Real t) const { return impl_->operator()(t); } Real prime(Real t) const { return impl_->prime(t); } Real double_prime(Real t) const { return impl_->double_prime(t); } Real t_max() const { return impl_->t_max(); } private: std::shared_ptr> impl_; }; }}} #endif