diff options
| author | Stefan Liebler <stli@linux.vnet.ibm.com> | 2014-04-16 12:52:52 +0200 |
|---|---|---|
| committer | Andreas Krebbel <krebbel@linux.vnet.ibm.com> | 2014-04-16 13:03:46 +0200 |
| commit | 8ea587db2b0bb5cd3e527982f7a1499748a987c3 (patch) | |
| tree | 5510f3fafaf3a170906e7edcdd06c42f7001745a /sysdeps/ieee754 | |
| parent | 9860b0450275ad2b69cb9360fd01d5c122a65fc5 (diff) | |
[BZ #16824] Fix failing y1 due to too large ulps in downward/upward rounding mode.
Diffstat (limited to 'sysdeps/ieee754')
| -rw-r--r-- | sysdeps/ieee754/ldbl-128/e_j1l.c | 1 |
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 70a1c86fd2..1264c95d49 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; |