diff options
| author | Joseph Myers <joseph@codesourcery.com> | 2015-06-24 15:12:03 +0000 |
|---|---|---|
| committer | Joseph Myers <joseph@codesourcery.com> | 2015-06-24 15:12:03 +0000 |
| commit | 8475ab16848cf2adb6ce3446749c2f1d10233286 (patch) | |
| tree | ce32b378b0e67c8d52f6a6a63347ae4b3711d159 /sysdeps/ieee754 | |
| parent | 36870482d2a33e3b703bd74b52063594389e827a (diff) | |
| download | glibc-8475ab16848cf2adb6ce3446749c2f1d10233286.tar.gz glibc-8475ab16848cf2adb6ce3446749c2f1d10233286.tar.xz glibc-8475ab16848cf2adb6ce3446749c2f1d10233286.zip | |
Fix ldbl-128 expl missing underflows (bug 18586).
Similar to various other bugs in this area, the ldbl-128 expl implementation does not raise the underflow exception for all subnormal results, if the scaling down is exact although the actual result is inexact. This patch fixes this by forcing the exception in this case (the tests that failed before and pass after the test are already in the testsuite). Tested for mips64. [BZ #18586] * sysdeps/ieee754/ldbl-128/e_expl.c (__ieee754_expl): Force underflow exception for small results.
Diffstat (limited to 'sysdeps/ieee754')
| -rw-r--r-- | sysdeps/ieee754/ldbl-128/e_expl.c | 10 |
1 files changed, 9 insertions, 1 deletions
diff --git a/sysdeps/ieee754/ldbl-128/e_expl.c b/sysdeps/ieee754/ldbl-128/e_expl.c index b4b789658a..1cd095cb77 100644 --- a/sysdeps/ieee754/ldbl-128/e_expl.c +++ b/sysdeps/ieee754/ldbl-128/e_expl.c @@ -230,7 +230,15 @@ __ieee754_expl (long double x) if (!unsafe) return result; else - return result * scale_u.d; + { + result *= scale_u.d; + if (result < LDBL_MIN) + { + long double force_underflow = result * result; + math_force_eval (force_underflow); + } + return result; + } } /* Exceptional cases: */ else if (isless (x, himark)) |