From 5e908464b9603e9cd95047d2edd00de5f78fedaf Mon Sep 17 00:00:00 2001 From: Jakub Jelinek Date: Wed, 13 Oct 2010 22:27:03 -0400 Subject: Implement accurate fma. --- sysdeps/ieee754/dbl-64/s_fmaf.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'sysdeps/ieee754/dbl-64/s_fmaf.c') diff --git a/sysdeps/ieee754/dbl-64/s_fmaf.c b/sysdeps/ieee754/dbl-64/s_fmaf.c index 9721876601..cd16cd1dce 100644 --- a/sysdeps/ieee754/dbl-64/s_fmaf.c +++ b/sysdeps/ieee754/dbl-64/s_fmaf.c @@ -39,7 +39,7 @@ __fmaf (float x, float y, float z) fesetround (FE_TOWARDZERO); /* Perform addition with round to odd. */ u.d = temp + (double) z; - if ((u.ieee.mantissa1 & 1) == 0 && u.ieee.exponent != 0xff) + if ((u.ieee.mantissa1 & 1) == 0 && u.ieee.exponent != 0x7ff) u.ieee.mantissa1 |= fetestexcept (FE_INEXACT) != 0; feupdateenv (&env); /* And finally truncation with round to nearest. */ -- cgit 1.4.1