interpolation()

template<typename T >

cdmatrix cntr::interpolation	(	int	tstp,
		double	tinter,
		function< T > &	f,
		int	InterpolationOrder
	)

k-th order polynomial interpolation of the contour function

Purpose

k-th order polynomial interpolation of the contour function to arbitrary point 'tinter' on time interval [tstp-k,tstp]. One uses information at times t=(n-j) [j=0...k]

Parameters

tstp	[int] Time step.
tinter	[double] interpolation time point
f	[function] the contour function to be interpolated
InterpolationOrder	[int] Order of interpolation

Definition at line 374 of file cntr_utilities_impl.hpp.

References cntr::function< T >::size1().

                                                  {

    assert(tstp>=InterpolationOrder);
    assert(tinter<=tstp);
    assert(tinter>=(tstp-InterpolationOrder));

    if (f.size1() == 1)
        return interpolation_dispatch<T, 1>(tstp,tinter, f, integration::I<T>(InterpolationOrder));
    else
        return interpolation_dispatch<T, LARGESIZE>(tstp,tinter, f, integration::I<T>(InterpolationOrder));
}

Here is the call graph for this function: