bugGNU Scientific Library - Bugs: bug #52359, Unexpected results in airy_Ai...


You are not allowed to post comments on this tracker with your current authentication level.

bug #52359: Unexpected results in airy_Ai function

Submitter:  Patrick Alken <psa>
Submitted:  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

Tue 07 Nov 2017 03:32:37 PM UTC, original submission:  

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.

Patrick Alken <psa>
Group administrator


(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attached Files
file #42343:  airy_divbyzero.c added by psa (678B - text/x-csrc)


Depends on the following items: None found

Items that depend on this one: None found


Carbon-Copy List
  • -email is unavailable- added by psa (Submitted the item)

    Follows 1 latest change.

    Date Changed by Updated Field Previous Value => Replaced by
    2017-11-07 psa Attached File- Added airy_divbyzero.c, #42343

    Back to the top

    Powered by Savane 3.13-b921.
    Corresponding source code