[help-gsl] Re: Bug in gsl_interp_accel_find() (2008-09-03)
|
From: |
Thomas Weber <thomas.weber.mail@.........> |
|
Subject: |
Re: Bug in gsl_interp_accel_find() |
|
Date: |
Wed, 03 Sep 2008 12:36:12 +0200 |
|
To: |
help-gsl@gnu.org |
Am Mittwoch, den 03.09.2008, 10:57 +0200 schrieb Thomas Weber:
> Hi,
>
> consider the attached example code; the last assertion fails. Note that
> idx2 == 0, but idx3 == 1.
> The reason for that is that gsl_interp_accel_find() checks only for
> strictly smaller/bigger values.
>
> (Untested) patch attached.
Related: the function
int
gsl_interp_accel_reset (gsl_interp_accel * a)
(in interpolation/accel.c)
is not documented at all.
Thomas
