From 8ea587db2b0bb5cd3e527982f7a1499748a987c3 Mon Sep 17 00:00:00 2001
From: Stefan Liebler <stli@linux.vnet.ibm.com>
Date: Wed, 16 Apr 2014 12:52:52 +0200
Subject: [BZ #16824] Fix failing y1 due to too large ulps in downward/upward
 rounding mode.

---
 sysdeps/ieee754/ldbl-128/e_j1l.c | 1 +
 1 file changed, 1 insertion(+)

(limited to 'sysdeps/ieee754')

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;
-- 
cgit 1.4.1