aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/ieee754/ldbl-128ibm/e_j1l.c
diff options
context:
space:
mode:
authorMichael Hudson-Doyle <michael.hudson@canonical.com>2022-08-22 14:05:04 +1200
committerTulio Magno Quites Machado Filho <tuliom@linux.ibm.com>2022-10-24 10:59:20 -0300
commitb6e37b7805b0182c3e25cdab39ebf5f001c04d05 (patch)
treed57e26f3abb4295763b9d371b3952574e1ebdfac /sysdeps/ieee754/ldbl-128ibm/e_j1l.c
parentb87a70e5e26e2656637da5f58ba4d26777d9e37b (diff)
downloadglibc-b6e37b7805b0182c3e25cdab39ebf5f001c04d05.zip
glibc-b6e37b7805b0182c3e25cdab39ebf5f001c04d05.tar.gz
glibc-b6e37b7805b0182c3e25cdab39ebf5f001c04d05.tar.bz2
Fix BZ #29463 in the ibm128 implementation of y1l too
Avoid moving code across SET_RESTORE_ROUNDL in order to fix [BZ #29463]. Tested-by: Aurelien Jarno <aurelien@aurel32.net> Reviewed-by: Aurelien Jarno <aurelien@aurel32.net> Reviewed-by: Tulio Magno Quites Machado Filho <tuliom@linux.ibm.com>
Diffstat (limited to 'sysdeps/ieee754/ldbl-128ibm/e_j1l.c')
-rw-r--r--sysdeps/ieee754/ldbl-128ibm/e_j1l.c3
1 files changed, 3 insertions, 0 deletions
diff --git a/sysdeps/ieee754/ldbl-128ibm/e_j1l.c b/sysdeps/ieee754/ldbl-128ibm/e_j1l.c
index f85ba94..0a5fe68 100644
--- a/sysdeps/ieee754/ldbl-128ibm/e_j1l.c
+++ b/sysdeps/ieee754/ldbl-128ibm/e_j1l.c
@@ -792,10 +792,13 @@ __ieee754_y1l (long double x)
{
/* 0 <= x <= 2 */
SET_RESTORE_ROUNDL (FE_TONEAREST);
+ xx = math_opt_barrier (xx);
+ x = math_opt_barrier (x);
z = xx * xx;
p = xx * neval (z, Y0_2N, NY0_2N) / deval (z, Y0_2D, NY0_2D);
p = -TWOOPI / xx + p;
p = TWOOPI * __ieee754_logl (x) * __ieee754_j1l (x) + p;
+ math_force_eval (p);
return p;
}