[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 |
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.
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Bug-gsl] bug in periodic csplines?,
Stefan Jahn <=