URL:
<
http://savannah.gnu.org/bugs/?24211>
Summary: gsl_interp_accel_find() boundary condition error
Project: GNU Scientific Library
Submitted by: bjg
Submitted on: Fri 05 Sep 2008 04:53:17 PM BST
Category: Runtime error
Severity: 4 - Important
Operating System:
Status: Confirmed
Assigned to: bjg
Open/Closed: Open
Release: 1.11
Discussion Lock: Any
_______________________________________________________
Details:
From: Thomas Weber <thomas.weber.mail@.........>
To: bug-gsl@gnu.org
Subject: [Bug-gsl] Bug in gsl_interp_accel_find()
Date: Wed, 03 Sep 2008 14:18:15 +0200
[1 <text/plain (7bit)>]
Hi,
sorry, just realized that there's a separate list for bug reports
(originally posted to
http://lists.gnu.org/archive/html/help-gsl/2008-09/msg00003.html)
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.
Thomas
_______________________________________________________
File Attachments:
-------------------------------------------------------
Date: Fri 05 Sep 2008 04:53:17 PM BST Name: break-accel.c Size: 513B By:
bjg
<
http://savannah.gnu.org/bugs/download.php?file_id=16433>
-------------------------------------------------------
Date: Fri 05 Sep 2008 04:53:17 PM BST Name: interp.diff Size: 503B By:
bjg
<
http://savannah.gnu.org/bugs/download.php?file_id=16434>
_______________________________________________________
Reply to this item at:
<
http://savannah.gnu.org/bugs/?24211>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
