aboutsummaryrefslogtreecommitdiff
path: root/sysdeps
diff options
context:
space:
mode:
authorStefan Liebler <stli@linux.vnet.ibm.com>2014-04-16 12:52:52 +0200
committerAndreas Krebbel <krebbel@linux.vnet.ibm.com>2014-04-16 13:03:46 +0200
commit8ea587db2b0bb5cd3e527982f7a1499748a987c3 (patch)
tree5510f3fafaf3a170906e7edcdd06c42f7001745a /sysdeps
parent9860b0450275ad2b69cb9360fd01d5c122a65fc5 (diff)
downloadglibc-8ea587db2b0bb5cd3e527982f7a1499748a987c3.zip
glibc-8ea587db2b0bb5cd3e527982f7a1499748a987c3.tar.gz
glibc-8ea587db2b0bb5cd3e527982f7a1499748a987c3.tar.bz2
[BZ #16824] Fix failing y1 due to too large ulps in downward/upward rounding mode.
Diffstat (limited to 'sysdeps')
-rw-r--r--sysdeps/ieee754/ldbl-128/e_j1l.c1
1 files changed, 1 insertions, 0 deletions
diff --git a/sysdeps/ieee754/ldbl-128/e_j1l.c b/sysdeps/ieee754/ldbl-128/e_j1l.c
index 70a1c86..1264c95 100644
--- a/sysdeps/ieee754/ldbl-128/e_j1l.c
+++ b/sysdeps/ieee754/ldbl-128/e_j1l.c
@@ -856,6 +856,7 @@ __ieee754_y1l (long double x)
if (xx <= 2.0L)
{
/* 0 <= x <= 2 */
+ SET_RESTORE_ROUNDL (FE_TONEAREST);
z = xx * xx;
p = xx * neval (z, Y0_2N, NY0_2N) / deval (z, Y0_2D, NY0_2D);
p = -TWOOPI / xx + p;