Re: [Bug-gsl] Interpolation bug?

[Brian Gough]

At Thu, 25 Sep 2008 10:51:43 -0400,
Andrew W. Steiner wrote:
> The function accel_find(0.9) gives 0, and accel_find(1.0) also gives 0,
> then accel_find(1.1) gives 1, and then accel_find(1.0) gives 1, instead of
> zero, because the cache value is different. I think the second call to
> accel_find() is in error, and should actually give 1 rather than 0. Of
> course, this is floating point arithmetic, but possibly there is a 
> modification
> to accel_find() which will make more sense most of the time.

Thanks for the bug report.  I think this may be the same as reported
in https://savannah.gnu.org/bugs/index.php?24211

Could you check the latest version in the GIT repository and see if
it fixes the problem? 

-- 
Brian Gough

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