From efb734887e10cb5adb0e722012fd65df8d1f8ff5 Mon Sep 17 00:00:00 2001 From: Aurelien Jarno Date: Fri, 1 Jun 2012 19:01:17 +0000 Subject: [PATCH] Add a barrier in the double fma implementation. --- ChangeLog | 5 +++++ sysdeps/ieee754/dbl-64/s_fma.c | 4 +++- 2 files changed, 8 insertions(+), 1 deletion(-) diff --git a/ChangeLog b/ChangeLog index 8cc25ef..f794e9f 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,3 +1,8 @@ +2012-06-01 Aurelien Jarno + + * sysdeps/ieee754/dbl-64/s_fma.c (__fma): Ensure a1 + u.d + computation is not scheduled after fetestexcept. + 2012-06-01 H.J. Lu [BZ #14117] diff --git a/sysdeps/ieee754/dbl-64/s_fma.c b/sysdeps/ieee754/dbl-64/s_fma.c index ab20a80..ce3bd36 100644 --- a/sysdeps/ieee754/dbl-64/s_fma.c +++ b/sysdeps/ieee754/dbl-64/s_fma.c @@ -1,5 +1,5 @@ /* Compute x * y + z as ternary operation. - Copyright (C) 2010, 2011 Free Software Foundation, Inc. + Copyright (C) 2010-2012 Free Software Foundation, Inc. This file is part of the GNU C Library. Contributed by Jakub Jelinek , 2010. @@ -158,6 +158,8 @@ __fma (double x, double y, double z) if ((u.ieee.mantissa1 & 1) == 0) u.ieee.mantissa1 |= libc_fetestexcept (FE_INEXACT) != 0; v.d = a1 + u.d; + /* Ensure the addition is not scheduled after fetestexcept call. */ + math_force_eval (v.d); } /* Reset rounding mode and test for inexact simultaneously. */ -- 2.7.4