diff options
| author | Joseph Myers <joseph@codesourcery.com> | 2015-10-21 16:40:20 +0000 |
|---|---|---|
| committer | Joseph Myers <joseph@codesourcery.com> | 2015-10-21 16:40:20 +0000 |
| commit | 0a90a8f2b9af6dd226558df5d3825e0e9ccd4ce3 (patch) | |
| tree | 3089be81a979063289fe998a6f14681fa00ecd5a /sysdeps/ieee754/ldbl-128/e_j0l.c | |
| parent | 1245e51596059a9844cff0938b10009c3b156115 (diff) | |
| download | glibc-0a90a8f2b9af6dd226558df5d3825e0e9ccd4ce3.tar.gz glibc-0a90a8f2b9af6dd226558df5d3825e0e9ccd4ce3.tar.xz glibc-0a90a8f2b9af6dd226558df5d3825e0e9ccd4ce3.zip | |
Fix ldbl-128 j0l spurious underflows (bug 19156).
My recent addition of more tests for j0 showed up that the ldbl-128 implementation of j0l produces spurious underflow exceptions for arguments close to 0 (when the result is very close to 1). This patch fixes this by just returning the argument in that case. Tested for mips64 (where it fixes the recently-added tests that were previously failing). [BZ #19156] * sysdeps/ieee754/ldbl-128/e_j0l.c (__ieee754_j0l): Return 1 for arguments very close to 0.
Diffstat (limited to 'sysdeps/ieee754/ldbl-128/e_j0l.c')
| -rw-r--r-- | sysdeps/ieee754/ldbl-128/e_j0l.c | 2 |
1 files changed, 2 insertions, 0 deletions
diff --git a/sysdeps/ieee754/ldbl-128/e_j0l.c b/sysdeps/ieee754/ldbl-128/e_j0l.c index 1320de7044..c208916a79 100644 --- a/sysdeps/ieee754/ldbl-128/e_j0l.c +++ b/sysdeps/ieee754/ldbl-128/e_j0l.c @@ -693,6 +693,8 @@ __ieee754_j0l (long double x) xx = fabsl (x); if (xx <= 2.0L) { + if (xx < 0x1p-57L) + return 1.0L; /* 0 <= x <= 2 */ z = xx * xx; p = z * z * neval (z, J0_2N, NJ0_2N) / deval (z, J0_2D, NJ0_2D); |