|Summary:||expm1() raises bogus invalid exception|
|Product:||glibc||Reporter:||Michael Kerrisk <mtk.manpages>|
|Component:||math||Assignee:||Andreas Jaeger <aj>|
Description Michael Kerrisk 2008-08-04 04:41:39 UTC
For certain large values of x, expm1() returns a bogus FE_INVALID exception and returns a NaN instead of positive invinity.
Comment 2 Michael Kerrisk 2008-08-04 04:49:16 UTC
$ /tmp/mt_expm1 1e5 errno == 0 fetestexcept() says: FE_INVALID FE_OVERFLOW FE_INEXACT expm1(1.00000000000000000e+05)=nan 0 FE_OVERFLOW nan $ /tmp/mt_expm1 1.00199970127e5 errno == 0 fetestexcept() says: FE_INVALID FE_OVERFLOW FE_INEXACT expm1(1.00199970126999993e+05)=nan 0 FE_OVERFLOW nan $ /tmp/mt_expm1 1.00199970128e5 errno == 0 fetestexcept() says: FE_OVERFLOW FE_INEXACT expm1(1.00199970128000001e+05)=inf 0 FE_OVERFLOW +inf
Comment 3 Michael Kerrisk 2008-08-04 04:50:14 UTC
Comment #2 should have started with thw words: "The first two runs below demonstrate the problem:"
Comment 4 Joseph Myers 2012-02-29 22:33:42 UTC
This appears to have been fixed, probably by various changes in 2009 (there is a test in the testsuite of expm1 overflow, which will verify there is no such inappropriate exception, and that test was added in 2009). However, the OVERFLOW exception for expm1 seems to have disappeared on x86, so I've filed bug 13787 for that problem.