Failure: Test: y1 (0X1.27E204P+99) == -0X1.FFFD6EP-51 Result: is: 0.00000000000000000000e+00 0x0.00000000000000000000p+0 should be: -8.88161002585147315788e-16 -0x1.fffd6e00000000000000p-51 difference: 8.88161002585147315788e-16 0x1.fffd6e00000000000000p-51 ulp : 16776887.0000 max.ulp : 0.0000
libm bugs remain suspended until somebody comes up with a patch. This doesn't mean that the expected behavior here is correct.
See http://sourceware.org/ml/libc-alpha/2012-02/msg00512.html for analysis.
Fixed by: commit c36e1d2369f149a64c759f3ad6f4cc5a03b03b62 Author: Joseph Myers <joseph@codesourcery.com> Date: Fri Mar 16 20:08:02 2012 +0000 Disable Bessel function TLOSS errors in POSIX mode.