#include <interpolation.h>

template<typename T>
class smash::InterpolateDataLinear< T >

Represent a piecewise linear interpolation.

Parameters

T	Type of interpolated values.

Definition at line 62 of file interpolation.h.

Public Member Functions
	InterpolateDataLinear (const std::vector< T > &x, const std::vector< T > &y)
	Interpolate function f given discrete samples f(x_i) = y_i. More...

T	operator() (T x) const
	Calculate spline interpolation at x. More...

Private Attributes
std::vector< T >	x_
	x_i More...

std::vector< InterpolateLinear< T > >	f_
	Piecewise linear interpolation using f(x_i) More...

Constructor & Destructor Documentation

◆ InterpolateDataLinear()

template<typename T >

smash::InterpolateDataLinear< T >::InterpolateDataLinear	(	const std::vector< T > &	x,
		const std::vector< T > &	y
	)

Interpolate function f given discrete samples f(x_i) = y_i.

Parameters

x	x-values.
y	y-values.

Returns: The interpolation function.

Piecewise linear interpolation is used. Values outside the given samples will use the outmost linear interpolation.

Definition at line 195 of file interpolation.h.

                                                                        {
   assert(x.size() == y.size());
   const size_t n = x.size();
   const auto p = generate_sort_permutation(
       x, [&](T const& a, T const& b) { return a < b; });
   x_ = apply_permutation(x, p);
   check_duplicates(x_, "InterpolateDataLinear");
   std::vector<T> y_sorted = std::move(apply_permutation(y, p));
   f_.reserve(n - 1);
   for (size_t i = 0; i < n - 1; i++) {
     f_.emplace_back(
         InterpolateLinear<T>(x_[i], y_sorted[i], x_[i + 1], y_sorted[i + 1]));
   }
 }

Member Function Documentation

◆ operator()()

template<typename T >

T smash::InterpolateDataLinear< T >::operator() ( T x ) const

Calculate spline interpolation at x.

Parameters

x	Interpolation argument.

Returns: Interpolated value.

Definition at line 240 of file interpolation.h.

                                                  {
   // Find the piecewise linear interpolation corresponding to x0.
   size_t i = find_index(x_, x0);
   if (i >= f_.size()) {
     // We don't have a linear interpolation beyond the last point in x_.
     // Use the last linear interpolation instead.
     i = f_.size() - 1;
   }
   return f_[i](x0);
 }