GNU 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 |
Release: |
Attached Files
Depends on the following items: None found
Items that depend on this one: None found
Carbon-Copy List
Follows 1 latest change.
Date | Changed by | Updated Field | Previous Value | => | Replaced by |
---|---|---|---|---|---|
2017-11-07 | psa | Attached File | - | ![]() |
Added airy_divbyzero.c, #42343 |
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.