Assuming I have an array of doubles, what's a good algorithm to sample this series using Akima interpolation? I'm too stupid to translate that mathematical description into code.
// values is an array of doubles
// idx is the index of the left-hand value for the current interpolation
// t is the normalized parameter between values[idx] and values[idx+1]
// Don't worry about array bounds, I'll handle that separately.
public double InterpolateAkima(double[] values, int idx, double t)
{
...?
}
Repost and expansion on my answer to another SO question that was closed as a duplicate of this question - as suggested by a comment on that question.
Akima's original paper: ``A new method of interpolation and smooth curve fitting based on local procedures'', Journal of ACM 17, 4 (1970), 589-602
http://www.leg.ufpr.br/lib/exe/fetch.php/wiki:internas:biblioteca:akima.pdf
C implementation
https://github.com/ampl/gsl/blob/master/interpolation/akima.c
C# implementation
https://gist.github.com/dreikanter/3526685
Delphi implementation (see procedure BuildAkimaSpline in delphi/src/spline3.pas)
http://www.alglib.net/translator/re/alglib-2.6.0.delphi.zip
Akima's Fortran 66 implementation
http://cran.r-project.org/web/packages/akima/
Fortran 90 implementation
http://miyoshi.googlecode.com/svn-history/r72/trunk/common/common.f90
Java implementation
Lisp implementation "for AutoCAD 2d-Polylines"
http://autocad.xarch.at/code/candido/akima.lsp
Matlab implementation
http://www.mathworks.se/matlabcentral/fileexchange/1814-akima-interpolation
Pascal implementation (program description)
http://jean-pierre.moreau.pagesperso-orange.fr/Pascal/akima_pas.txt
Python implementation
http://www.lfd.uci.edu/~gohlke/code/akima.py.html
VB6 implementation (see subroutine BuildAkimaSpline in vb6/src/spline3.bas)
http://www.alglib.net/translator/re/alglib-2.6.0.vb6.zip
http://www.koders.com/cpp/fid1393B9D668316C1700966643DE0609660B9CB13A.aspx?s=%22Brian+Smith%22