From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: by sourceware.org (Postfix, from userid 1827) id CE4D9385840C; Mon, 24 Oct 2022 14:00:33 +0000 (GMT) DKIM-Filter: OpenDKIM Filter v2.11.0 sourceware.org CE4D9385840C DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=sourceware.org; s=default; t=1666620036; bh=hEgOm8h7UGLqHZKpfCuwRbI/HF08LfNkGnuHa2K4a5A=; h=From:To:Subject:Date:From; b=OscbSR/gXklkuqkJTdSoLGX+W0fHTYZFpMkcroUbFOLLUYZy6zw1IuZjc9IkZ835Y 7ABPCESwlsslkhwo3Bwc8Ua316mmK59fetRSCBzCnelKetprIQL8y4dBdX28w5hDAM 6amC0eumEh5b9aOvRQLDQh2+ANIsdG5GsvGOsbD4= Content-Type: text/plain; charset="us-ascii" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit From: Tulio Magno Quites Machado Filho To: glibc-cvs@sourceware.org Subject: [glibc] Fix BZ #29463 in the ibm128 implementation of y1l too X-Act-Checkin: glibc X-Git-Author: Michael Hudson-Doyle X-Git-Refname: refs/heads/master X-Git-Oldrev: b87a70e5e26e2656637da5f58ba4d26777d9e37b X-Git-Newrev: b6e37b7805b0182c3e25cdab39ebf5f001c04d05 Message-Id: <20221024140036.CE4D9385840C@sourceware.org> Date: Mon, 24 Oct 2022 14:00:33 +0000 (GMT) List-Id: https://sourceware.org/git/gitweb.cgi?p=glibc.git;h=b6e37b7805b0182c3e25cdab39ebf5f001c04d05 commit b6e37b7805b0182c3e25cdab39ebf5f001c04d05 Author: Michael Hudson-Doyle Date: Mon Aug 22 14:05:04 2022 +1200 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 Reviewed-by: Aurelien Jarno Reviewed-by: Tulio Magno Quites Machado Filho Diff: --- sysdeps/ieee754/ldbl-128ibm/e_j1l.c | 3 +++ 1 file changed, 3 insertions(+) diff --git a/sysdeps/ieee754/ldbl-128ibm/e_j1l.c b/sysdeps/ieee754/ldbl-128ibm/e_j1l.c index f85ba94466..0a5fe68342 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; }