[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Bug-gsl] bug in periodic csplines?

From: Stefan Jahn
Subject: [Bug-gsl] bug in periodic csplines?
Date: Thu, 27 Oct 2005 09:10:33 +0200 (CEST)
User-agent: SquirrelMail/1.4.5


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.

reply via email to

[Prev in Thread] Current Thread [Next in Thread]