diff options
| author | John David Anglin <danglin@gcc.gnu.org> | 2025-02-25 15:57:53 -0500 |
|---|---|---|
| committer | John David Anglin <danglin@gcc.gnu.org> | 2025-02-25 16:00:49 -0500 |
| commit | 523f85558152a1b9cced6d669f758c27677775ba (patch) | |
| tree | ee9bafc289f8630f19311fe1b1389aa3fdd8ca0c | |
| parent | d6c156c326999f144cb5b73d29982108d549ad8a (diff) | |
| download | glibc-523f85558152a1b9cced6d669f758c27677775ba.tar.xz glibc-523f85558152a1b9cced6d669f758c27677775ba.zip | |
math: Add optimization barrier to ensure a1 + u.d is not reused [BZ #30664]
A number of fma tests started to fail on hppa when gcc was changed to
use Ranger rather than EVRP. Eventually I found that the value of
a1 + u.d in this is block of code was being computed in FE_TOWARDZERO
mode and not the original rounding mode:
if (TININESS_AFTER_ROUNDING)
{
w.d = a1 + u.d;
if (w.ieee.exponent == 109)
return w.d * 0x1p-108;
}
This caused the exponent value to be wrong and the wrong return path
to be used.
Here we add an optimization barrier after the rounding mode is reset
to ensure that the previous value of a1 + u.d is not reused.
Signed-off-by: John David Anglin <dave.anglin@bell.net>
| -rw-r--r-- | sysdeps/ieee754/dbl-64/s_fma.c | 3 |
1 files changed, 3 insertions, 0 deletions
diff --git a/sysdeps/ieee754/dbl-64/s_fma.c b/sysdeps/ieee754/dbl-64/s_fma.c index c5f5abdc68..79a3cd721d 100644 --- a/sysdeps/ieee754/dbl-64/s_fma.c +++ b/sysdeps/ieee754/dbl-64/s_fma.c @@ -244,6 +244,9 @@ __fma (double x, double y, double z) /* Reset rounding mode and test for inexact simultaneously. */ int j = libc_feupdateenv_test (&env, FE_INEXACT) != 0; + /* Ensure value of a1 + u.d is not reused. */ + a1 = math_opt_barrier (a1); + if (__glibc_likely (adjust == 0)) { if ((u.ieee.mantissa1 & 1) == 0 && u.ieee.exponent != 0x7ff) |