[Bug-gsl] bug in periodic csplines?

[Stefan Jahn]

hello!

i think i found a small bug in periodic cubic splines calculation.  for
the 3-point calculation you need to apply this patch:

cspline.c: cspline_init_periodic patch

    - const double h2 = xa[3] - xa[2];

    - g[1] = 3.0 * ((ya[1] - ya[2]) / h2 - (ya[2] - ya[1]) / h1);
    + g[1] = 3.0 * ((ya[1] - ya[2]) / h0 - (ya[2] - ya[1]) / h1);

can anyone verify that?  please answer privately as well, because i am
not subscribed to the list.

cheers, stefan.