about summary refs log tree commit diff
path: root/sysdeps/ieee754/ldbl-96/s_fma.c
diff options
context:
space:
mode:
authorRichard Henderson <rth@twiddle.net>2014-08-01 12:13:50 -1000
committerRichard Henderson <rth@twiddle.net>2014-08-01 12:13:50 -1000
commit4896f049205013528f1f0df7f7cbb52e8c5664d1 (patch)
tree7f397f3b733485f3a4859eed19251c64d2eca64d /sysdeps/ieee754/ldbl-96/s_fma.c
parent6f64e7d75fdab0be50edc814a41ddc43ef4ec8ef (diff)
downloadglibc-4896f049205013528f1f0df7f7cbb52e8c5664d1.tar.gz
glibc-4896f049205013528f1f0df7f7cbb52e8c5664d1.tar.xz
glibc-4896f049205013528f1f0df7f7cbb52e8c5664d1.zip
Force eval for fma implementations
Diffstat (limited to 'sysdeps/ieee754/ldbl-96/s_fma.c')
-rw-r--r--sysdeps/ieee754/ldbl-96/s_fma.c11
1 files changed, 6 insertions, 5 deletions
diff --git a/sysdeps/ieee754/ldbl-96/s_fma.c b/sysdeps/ieee754/ldbl-96/s_fma.c
index 354be4eae9..a27102eff1 100644
--- a/sysdeps/ieee754/ldbl-96/s_fma.c
+++ b/sysdeps/ieee754/ldbl-96/s_fma.c
@@ -64,16 +64,17 @@ __fma (double x, double y, double z)
   t1 = m1 - t1;
   t2 = z - t2;
   long 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)
     {
       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;
     }