[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Help-gsl] Bug in gsl_interp_accel_find()
|
From: |
Thomas Weber |
|
Subject: |
Re: [Help-gsl] Bug in gsl_interp_accel_find() |
|
Date: |
Wed, 03 Sep 2008 12:36:12 +0200 |
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