diff options
| author | Richard Henderson <rth@twiddle.net> | 2014-08-01 12:13:50 -1000 |
|---|---|---|
| committer | Richard Henderson <rth@twiddle.net> | 2014-08-01 12:13:50 -1000 |
| commit | 4896f049205013528f1f0df7f7cbb52e8c5664d1 (patch) | |
| tree | 7f397f3b733485f3a4859eed19251c64d2eca64d /sysdeps/ieee754/dbl-64 | |
| parent | 6f64e7d75fdab0be50edc814a41ddc43ef4ec8ef (diff) | |
| download | glibc-4896f049205013528f1f0df7f7cbb52e8c5664d1.tar.gz glibc-4896f049205013528f1f0df7f7cbb52e8c5664d1.tar.xz glibc-4896f049205013528f1f0df7f7cbb52e8c5664d1.zip | |
Force eval for fma implementations
Diffstat (limited to 'sysdeps/ieee754/dbl-64')
| -rw-r--r-- | sysdeps/ieee754/dbl-64/s_fma.c | 11 |
1 files changed, 6 insertions, 5 deletions
diff --git a/sysdeps/ieee754/dbl-64/s_fma.c b/sysdeps/ieee754/dbl-64/s_fma.c index 389acd4cde..77065aaba1 100644 --- a/sysdeps/ieee754/dbl-64/s_fma.c +++ b/sysdeps/ieee754/dbl-64/s_fma.c @@ -198,16 +198,17 @@ __fma (double x, double y, double z) t1 = m1 - t1; t2 = z - t2; double a2 = t1 + t2; + /* Ensure the arithmetic is not scheduled after feclearexcept call. */ + math_force_eval (m2); + math_force_eval (a2); feclearexcept (FE_INEXACT); - /* If the result is an exact zero, ensure it has the correct - sign. */ + /* If the result is an exact zero, ensure it has the correct sign. */ if (a1 == 0 && m2 == 0) { libc_feupdateenv (&env); - /* Ensure that round-to-nearest value of z + m1 is not - reused. */ - asm volatile ("" : "=m" (z) : "m" (z)); + /* Ensure that round-to-nearest value of z + m1 is not reused. */ + z = math_opt_barrier (z); return z + m1; } |