[Bug-gsl] [bug #24211] gsl_interp_accel_find() boundary condition error

[Brian Gough]

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 <address@hidden>
To: address@hidden
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/