NESSi  v1.0.2
The NonEquilibrium Systems Simulation Library

◆ 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().

375  {
376 
377  assert(tstp>=InterpolationOrder);
378  assert(tinter<=tstp);
379  assert(tinter>=(tstp-InterpolationOrder));
380 
381  if (f.size1() == 1)
382  return interpolation_dispatch<T, 1>(tstp,tinter, f, integration::I<T>(InterpolationOrder));
383  else
384  return interpolation_dispatch<T, LARGESIZE>(tstp,tinter, f, integration::I<T>(InterpolationOrder));
385 }
+ Here is the call graph for this function: