Re: [Bug-gsl] bug in akima interpolation

[Brian Gough]

ABX writes:
 > GSL:
 > 
 >     m[-2] = 3.0 * m[0] - 2.0 * m[1]; m[-1] = 2.0 * m[0] -
 >     m[1]; m[size - 1] = 2.0 * m[size - 2] - m[size - 3];
 >     m[size] = 3.0 * m[size - 2] - 2.0 * m[size - 3];
 >   

Thanks for the bug report.  The code in GSL uses a variant of
akima splines called non-rounded akima splines, but this is
not documented in the manual.  I will make a note to document
this properly and see if we can also add the normal rounded
akima splines for the next release.

best regards,

-- 
Brian Gough

----------------------------------------------------------------------
Network Theory Ltd             Phone: +44 (0)117 3179309
15 Royal Park                    Fax: +44 (0)117 9048108
Bristol BS8 3AL                  WWW: http://www.network-theory.co.uk/
United Kingdom                 Email: address@hidden     
----------------------------------------------------------------------