PreviousIndexLast

[bug-gsl] Re: Bug in gsl_interp_accel_find() (2008-09-05)

From: Brian Gough <bjg@gnu.org>
Subject: Re: Bug in gsl_interp_accel_find()
Date: Fri, 05 Sep 2008 16:48:43 +0100
To: Thomas Weber <thomas.weber.mail@.........>
Cc: bug-gsl@gnu.org

At Wed, 03 Sep 2008 14:18:15 +0200,
Thomas Weber wrote:
> 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.

Thanks for the excellent bug report. I have committed your patch for
the next release and added some extra tests for the accelerator.

--
Brian Gough

GNU Scientific Library -
http://www.gnu.org/software/gsl/