[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Bug-gsl] [bug #52359] Unexpected results in airy_Ai function

From: Patrick Alken
Subject: [Bug-gsl] [bug #52359] Unexpected results in airy_Ai function
Date: Tue, 7 Nov 2017 10:32:38 -0500 (EST)
User-agent: Mozilla/5.0 (X11; Linux x86_64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/61.0.3163.100 Safari/537.36


                 Summary: Unexpected results in airy_Ai function
                 Project: GNU Scientific Library
            Submitted by: psa
            Submitted on: Tue 07 Nov 2017 03:32:37 PM UTC
                Category: Runtime error
                Severity: 3 - Normal
        Operating System: 
                  Status: None
             Assigned to: None
             Open/Closed: Open
         Discussion Lock: Any



The GSL’s airy_Ai function gsl_sf_airy_Ai_e  triggers a floating-point
division-by-zero exception with the input -1.842761151977744. The exception
can be trapped by using “feenableexcept”, but it will disappear if you
slightly disturb the input, say, to -1.84276115198. By digging into the issue,
I have observed that in the airy_mod_phase function that airy_Ai invokes, the
variable result_m is divided whereas it has vanished following a non-trivial
computation (with loops, etc., in function cheb_eval_mode_e). Please check the
code airy_divbyzero.c below.


File Attachments:

Date: Tue 07 Nov 2017 03:32:37 PM UTC  Name: airy_divbyzero.c  Size: 678B  
By: psa



Reply to this item at:


  Message sent via/by Savannah

reply via email to

[Prev in Thread] Current Thread [Next in Thread]