[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Bug in gsl_interp_accel_find()
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Bug in gsl_interp_accel_find() |
Date: |
Fri, 05 Sep 2008 16:48:43 +0100 |
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/22.2 Mule/5.0 (SAKAKI) |
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/