This is the mail archive of the
glibc-cvs@sourceware.org
mailing list for the glibc project.
GNU C Library master sources branch master updated. glibc-2.22-67-g739babd
- From: jsm28 at sourceware dot org
- To: glibc-cvs at sourceware dot org
- Date: 14 Aug 2015 17:15:33 -0000
- Subject: GNU C Library master sources branch master updated. glibc-2.22-67-g739babd
This is an automated email from the git hooks/post-receive script. It was
generated because a ref change was pushed to the repository containing
the project "GNU C Library master sources".
The branch, master has been updated
via 739babd775d4b69525e3494cad4310742a9b360a (commit)
from db7f8c8fe0575fa4c441c2f04e666a01cb06f0cc (commit)
Those revisions listed above that are new to this repository have
not appeared on any other notification email; so we list those
revisions in full, below.
- Log -----------------------------------------------------------------
http://sourceware.org/git/gitweb.cgi?p=glibc.git;a=commitdiff;h=739babd775d4b69525e3494cad4310742a9b360a
commit 739babd775d4b69525e3494cad4310742a9b360a
Author: Joseph Myers <joseph@codesourcery.com>
Date: Fri Aug 14 17:15:06 2015 +0000
Fix fma spurious underflows (bug 18824).
Various fma implementations have logic that, when computing fma (x, y,
z) where z is large (so care needs taking to avoid internal overflow)
but x * y is small, scale x * y up instead of down to avoid internal
underflows resulting from scaling down. (In these cases, x * y is
small enough that only its sign actually matters rather than the exact
value.)
The threshold for scaling up instead of down was correct for "if the
unscaled values were multiplied, the low part of the multiplication
could underflow", and the scaling was sufficient to ensure that the
low part of the multiplication did not underflow (given that cases of
very small x * y - less than half the least subnormal - were
previously dealt with). However, the choice in the functions wasn't
between scaling up or no scaling, but between scaling up and scaling
down (scaling down actually being needed when x * y isn't so small
compared to z and so the exact value does matter). Thus a larger
threshold is needed to ensure that scaling down doesn't produce values
the multiplication of whose low parts underflows. This patch
increases the thresholds accordingly.
Tested for x86_64, x86 and mips64 (with the MIPS version of s_fmal.c
removed so that the ldbl-128 version gets tested instead of the
soft-fp one).
[BZ #18824]
* sysdeps/ieee754/dbl-64/s_fma.c (__fma): Increase threshold for
scaling x * y up instead of down.
* sysdeps/ieee754/ldbl-128/s_fmal.c (__fmal): Likewise.
* sysdeps/ieee754/ldbl-96/s_fmal.c (__fmal): Likewise.
* math/auto-libm-test-in: Add more tests of fma.
* math/auto-libm-test-out: Regenerated.
diff --git a/ChangeLog b/ChangeLog
index 9a74f25..f510bea 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,13 @@
+2015-08-14 Joseph Myers <joseph@codesourcery.com>
+
+ [BZ #18824]
+ * sysdeps/ieee754/dbl-64/s_fma.c (__fma): Increase threshold for
+ scaling x * y up instead of down.
+ * sysdeps/ieee754/ldbl-128/s_fmal.c (__fmal): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_fmal.c (__fmal): Likewise.
+ * math/auto-libm-test-in: Add more tests of fma.
+ * math/auto-libm-test-out: Regenerated.
+
2015-08-14 Paul Pluzhnikov <ppluzhnikov@google.com>
* sysdeps/x86_64/fpu/libm-test-ulps: Regenerated.
diff --git a/NEWS b/NEWS
index 5ca3de8..088969c 100644
--- a/NEWS
+++ b/NEWS
@@ -10,7 +10,7 @@ Version 2.23
* The following bugs are resolved with this release:
16517, 16519, 16520, 16734, 17905, 18086, 18265, 18480, 18525, 18618,
- 18647, 18661, 18674, 18778, 18781, 18787, 18789, 18790, 18820.
+ 18647, 18661, 18674, 18778, 18781, 18787, 18789, 18790, 18820, 18824.
Version 2.22
diff --git a/math/auto-libm-test-in b/math/auto-libm-test-in
index d3a4834..4480371 100644
--- a/math/auto-libm-test-in
+++ b/math/auto-libm-test-in
@@ -1794,6 +1794,24 @@ fma 0x1.0000000000000000000000000001p-16382 -0x1.0000000000000000000000000001p-6
fma 0x1.0000000000000000000000000001p-16382 0x1.0000000000000000000000000001p-66 -0x1p16319
fma 0x1.0000000000000000000000000001p-16382 -0x1.0000000000000000000000000001p-66 -0x1p16319
+# Bug 6801: errno setting may be missing.
+fma 0x1.fffffep-126 0x1.fffffep25 0x1.fffffep127 missing-errno
+fma 0x1.fffffep-126 -0x1.fffffep25 0x1.fffffep127
+fma 0x1.fffffep-126 0x1.fffffep25 -0x1.fffffep127
+fma 0x1.fffffep-126 -0x1.fffffep25 -0x1.fffffep127 missing-errno
+fma 0x1.fffffffffffffp-1022 0x1.fffffffffffffp54 0x1.fffffffffffffp1023 missing-errno
+fma 0x1.fffffffffffffp-1022 -0x1.fffffffffffffp54 0x1.fffffffffffffp1023
+fma 0x1.fffffffffffffp-1022 0x1.fffffffffffffp54 -0x1.fffffffffffffp1023
+fma 0x1.fffffffffffffp-1022 -0x1.fffffffffffffp54 -0x1.fffffffffffffp1023 missing-errno
+fma 0x1.fffffffffffffffep-16382 0x1.fffffffffffffffep65 0x1.fffffffffffffffep16383 missing-errno
+fma 0x1.fffffffffffffffep-16382 -0x1.fffffffffffffffep65 0x1.fffffffffffffffep16383
+fma 0x1.fffffffffffffffep-16382 0x1.fffffffffffffffep65 -0x1.fffffffffffffffep16383
+fma 0x1.fffffffffffffffep-16382 -0x1.fffffffffffffffep65 -0x1.fffffffffffffffep16383 missing-errno
+fma 0x1.ffffffffffffffffffffffffffffp-16382 0x1.ffffffffffffffffffffffffffffp114 0x1.ffffffffffffffffffffffffffffp16383 missing-errno
+fma 0x1.ffffffffffffffffffffffffffffp-16382 -0x1.ffffffffffffffffffffffffffffp114 0x1.ffffffffffffffffffffffffffffp16383
+fma 0x1.ffffffffffffffffffffffffffffp-16382 0x1.ffffffffffffffffffffffffffffp114 -0x1.ffffffffffffffffffffffffffffp16383
+fma 0x1.ffffffffffffffffffffffffffffp-16382 -0x1.ffffffffffffffffffffffffffffp114 -0x1.ffffffffffffffffffffffffffffp16383 missing-errno
+
hypot 0 0
hypot 0 -0
hypot -0 0
diff --git a/math/auto-libm-test-out b/math/auto-libm-test-out
index 1900e78..d8ede58 100644
--- a/math/auto-libm-test-out
+++ b/math/auto-libm-test-out
@@ -132191,6 +132191,262 @@ fma 0x1.0000000000000000000000000001p-16382 -0x1.0000000000000000000000000001p-6
= fma tonearest ldbl-128 0x4.0000000000000000000000000004p-16384L -0x4.0000000000000000000000000004p-68L -0x8p+16316L : -0x8p+16316L : inexact
= fma towardzero ldbl-128 0x4.0000000000000000000000000004p-16384L -0x4.0000000000000000000000000004p-68L -0x8p+16316L : -0x8p+16316L : inexact
= fma upward ldbl-128 0x4.0000000000000000000000000004p-16384L -0x4.0000000000000000000000000004p-68L -0x8p+16316L : -0x8p+16316L : inexact
+fma 0x1.fffffep-126 0x1.fffffep25 0x1.fffffep127 missing-errno
+= fma downward flt-32 0x7.fffff8p-128f 0x3.fffffcp+24f 0xf.fffffp+124f : 0xf.fffffp+124f : inexact
+= fma tonearest flt-32 0x7.fffff8p-128f 0x3.fffffcp+24f 0xf.fffffp+124f : 0xf.fffffp+124f : inexact
+= fma towardzero flt-32 0x7.fffff8p-128f 0x3.fffffcp+24f 0xf.fffffp+124f : 0xf.fffffp+124f : inexact
+= fma upward flt-32 0x7.fffff8p-128f 0x3.fffffcp+24f 0xf.fffffp+124f : plus_infty : inexact overflow errno-erange errno-erange-ok
+= fma downward dbl-64 0x7.fffff8p-128 0x3.fffffcp+24 0xf.fffffp+124 : 0xf.fffffp+124 : inexact
+= fma tonearest dbl-64 0x7.fffff8p-128 0x3.fffffcp+24 0xf.fffffp+124 : 0xf.fffffp+124 : inexact
+= fma towardzero dbl-64 0x7.fffff8p-128 0x3.fffffcp+24 0xf.fffffp+124 : 0xf.fffffp+124 : inexact
+= fma upward dbl-64 0x7.fffff8p-128 0x3.fffffcp+24 0xf.fffffp+124 : 0xf.fffff00000008p+124 : inexact
+= fma downward ldbl-96-intel 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma tonearest ldbl-96-intel 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma towardzero ldbl-96-intel 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma upward ldbl-96-intel 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffff0000000001p+124L : inexact
+= fma downward ldbl-96-m68k 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma tonearest ldbl-96-m68k 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma towardzero ldbl-96-m68k 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma upward ldbl-96-m68k 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffff0000000001p+124L : inexact
+= fma downward ldbl-128 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma tonearest ldbl-128 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma towardzero ldbl-128 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma upward ldbl-128 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffff00000000000000000000008p+124L : inexact
+= fma downward ldbl-128ibm 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma tonearest ldbl-128ibm 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma towardzero ldbl-128ibm 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma upward ldbl-128ibm 0x7.fffff8p-128L 0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffff000000000000000000004p+124L : inexact
+fma 0x1.fffffep-126 -0x1.fffffep25 0x1.fffffep127
+= fma downward flt-32 0x7.fffff8p-128f -0x3.fffffcp+24f 0xf.fffffp+124f : 0xf.ffffep+124f : inexact
+= fma tonearest flt-32 0x7.fffff8p-128f -0x3.fffffcp+24f 0xf.fffffp+124f : 0xf.fffffp+124f : inexact
+= fma towardzero flt-32 0x7.fffff8p-128f -0x3.fffffcp+24f 0xf.fffffp+124f : 0xf.ffffep+124f : inexact
+= fma upward flt-32 0x7.fffff8p-128f -0x3.fffffcp+24f 0xf.fffffp+124f : 0xf.fffffp+124f : inexact
+= fma downward dbl-64 0x7.fffff8p-128 -0x3.fffffcp+24 0xf.fffffp+124 : 0xf.ffffefffffff8p+124 : inexact
+= fma tonearest dbl-64 0x7.fffff8p-128 -0x3.fffffcp+24 0xf.fffffp+124 : 0xf.fffffp+124 : inexact
+= fma towardzero dbl-64 0x7.fffff8p-128 -0x3.fffffcp+24 0xf.fffffp+124 : 0xf.ffffefffffff8p+124 : inexact
+= fma upward dbl-64 0x7.fffff8p-128 -0x3.fffffcp+24 0xf.fffffp+124 : 0xf.fffffp+124 : inexact
+= fma downward ldbl-96-intel 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.ffffeffffffffffp+124L : inexact
+= fma tonearest ldbl-96-intel 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma towardzero ldbl-96-intel 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.ffffeffffffffffp+124L : inexact
+= fma upward ldbl-96-intel 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma downward ldbl-96-m68k 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.ffffeffffffffffp+124L : inexact
+= fma tonearest ldbl-96-m68k 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma towardzero ldbl-96-m68k 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.ffffeffffffffffp+124L : inexact
+= fma upward ldbl-96-m68k 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma downward ldbl-128 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.ffffeffffffffffffffffffffff8p+124L : inexact
+= fma tonearest ldbl-128 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma towardzero ldbl-128 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.ffffeffffffffffffffffffffff8p+124L : inexact
+= fma upward ldbl-128 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma downward ldbl-128ibm 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.ffffeffffffffffffffffffffcp+124L : inexact
+= fma tonearest ldbl-128ibm 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+= fma towardzero ldbl-128ibm 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.ffffeffffffffffffffffffffcp+124L : inexact
+= fma upward ldbl-128ibm 0x7.fffff8p-128L -0x3.fffffcp+24L 0xf.fffffp+124L : 0xf.fffffp+124L : inexact
+fma 0x1.fffffep-126 0x1.fffffep25 -0x1.fffffep127
+= fma downward flt-32 0x7.fffff8p-128f 0x3.fffffcp+24f -0xf.fffffp+124f : -0xf.fffffp+124f : inexact
+= fma tonearest flt-32 0x7.fffff8p-128f 0x3.fffffcp+24f -0xf.fffffp+124f : -0xf.fffffp+124f : inexact
+= fma towardzero flt-32 0x7.fffff8p-128f 0x3.fffffcp+24f -0xf.fffffp+124f : -0xf.ffffep+124f : inexact
+= fma upward flt-32 0x7.fffff8p-128f 0x3.fffffcp+24f -0xf.fffffp+124f : -0xf.ffffep+124f : inexact
+= fma downward dbl-64 0x7.fffff8p-128 0x3.fffffcp+24 -0xf.fffffp+124 : -0xf.fffffp+124 : inexact
+= fma tonearest dbl-64 0x7.fffff8p-128 0x3.fffffcp+24 -0xf.fffffp+124 : -0xf.fffffp+124 : inexact
+= fma towardzero dbl-64 0x7.fffff8p-128 0x3.fffffcp+24 -0xf.fffffp+124 : -0xf.ffffefffffff8p+124 : inexact
+= fma upward dbl-64 0x7.fffff8p-128 0x3.fffffcp+24 -0xf.fffffp+124 : -0xf.ffffefffffff8p+124 : inexact
+= fma downward ldbl-96-intel 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma tonearest ldbl-96-intel 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma towardzero ldbl-96-intel 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.ffffeffffffffffp+124L : inexact
+= fma upward ldbl-96-intel 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.ffffeffffffffffp+124L : inexact
+= fma downward ldbl-96-m68k 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma tonearest ldbl-96-m68k 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma towardzero ldbl-96-m68k 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.ffffeffffffffffp+124L : inexact
+= fma upward ldbl-96-m68k 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.ffffeffffffffffp+124L : inexact
+= fma downward ldbl-128 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma tonearest ldbl-128 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma towardzero ldbl-128 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.ffffeffffffffffffffffffffff8p+124L : inexact
+= fma upward ldbl-128 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.ffffeffffffffffffffffffffff8p+124L : inexact
+= fma downward ldbl-128ibm 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma tonearest ldbl-128ibm 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma towardzero ldbl-128ibm 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.ffffeffffffffffffffffffffcp+124L : inexact
+= fma upward ldbl-128ibm 0x7.fffff8p-128L 0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.ffffeffffffffffffffffffffcp+124L : inexact
+fma 0x1.fffffep-126 -0x1.fffffep25 -0x1.fffffep127 missing-errno
+= fma downward flt-32 0x7.fffff8p-128f -0x3.fffffcp+24f -0xf.fffffp+124f : minus_infty : inexact overflow errno-erange errno-erange-ok
+= fma tonearest flt-32 0x7.fffff8p-128f -0x3.fffffcp+24f -0xf.fffffp+124f : -0xf.fffffp+124f : inexact
+= fma towardzero flt-32 0x7.fffff8p-128f -0x3.fffffcp+24f -0xf.fffffp+124f : -0xf.fffffp+124f : inexact
+= fma upward flt-32 0x7.fffff8p-128f -0x3.fffffcp+24f -0xf.fffffp+124f : -0xf.fffffp+124f : inexact
+= fma downward dbl-64 0x7.fffff8p-128 -0x3.fffffcp+24 -0xf.fffffp+124 : -0xf.fffff00000008p+124 : inexact
+= fma tonearest dbl-64 0x7.fffff8p-128 -0x3.fffffcp+24 -0xf.fffffp+124 : -0xf.fffffp+124 : inexact
+= fma towardzero dbl-64 0x7.fffff8p-128 -0x3.fffffcp+24 -0xf.fffffp+124 : -0xf.fffffp+124 : inexact
+= fma upward dbl-64 0x7.fffff8p-128 -0x3.fffffcp+24 -0xf.fffffp+124 : -0xf.fffffp+124 : inexact
+= fma downward ldbl-96-intel 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffff0000000001p+124L : inexact
+= fma tonearest ldbl-96-intel 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma towardzero ldbl-96-intel 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma upward ldbl-96-intel 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma downward ldbl-96-m68k 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffff0000000001p+124L : inexact
+= fma tonearest ldbl-96-m68k 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma towardzero ldbl-96-m68k 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma upward ldbl-96-m68k 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma downward ldbl-128 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffff00000000000000000000008p+124L : inexact
+= fma tonearest ldbl-128 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma towardzero ldbl-128 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma upward ldbl-128 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma downward ldbl-128ibm 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffff000000000000000000004p+124L : inexact
+= fma tonearest ldbl-128ibm 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma towardzero ldbl-128ibm 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+= fma upward ldbl-128ibm 0x7.fffff8p-128L -0x3.fffffcp+24L -0xf.fffffp+124L : -0xf.fffffp+124L : inexact
+fma 0x1.fffffffffffffp-1022 0x1.fffffffffffffp54 0x1.fffffffffffffp1023 missing-errno
+= fma downward dbl-64 0x7.ffffffffffffcp-1024 0x7.ffffffffffffcp+52 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffff8p+1020 : inexact
+= fma tonearest dbl-64 0x7.ffffffffffffcp-1024 0x7.ffffffffffffcp+52 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffff8p+1020 : inexact
+= fma towardzero dbl-64 0x7.ffffffffffffcp-1024 0x7.ffffffffffffcp+52 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffff8p+1020 : inexact
+= fma upward dbl-64 0x7.ffffffffffffcp-1024 0x7.ffffffffffffcp+52 0xf.ffffffffffff8p+1020 : plus_infty : inexact overflow errno-erange errno-erange-ok
+= fma downward ldbl-96-intel 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma tonearest ldbl-96-intel 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-96-intel 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma upward ldbl-96-intel 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff801p+1020L : inexact
+= fma downward ldbl-96-m68k 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma tonearest ldbl-96-m68k 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-96-m68k 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma upward ldbl-96-m68k 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff801p+1020L : inexact
+= fma downward ldbl-128 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma tonearest ldbl-128 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-128 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma upward ldbl-128 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8000000000000008p+1020L : inexact
+= fma downward ldbl-128ibm 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma tonearest ldbl-128ibm 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-128ibm 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma upward ldbl-128ibm 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff80000000000004p+1020L : inexact
+fma 0x1.fffffffffffffp-1022 -0x1.fffffffffffffp54 0x1.fffffffffffffp1023
+= fma downward dbl-64 0x7.ffffffffffffcp-1024 -0x7.ffffffffffffcp+52 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffp+1020 : inexact
+= fma tonearest dbl-64 0x7.ffffffffffffcp-1024 -0x7.ffffffffffffcp+52 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffff8p+1020 : inexact
+= fma towardzero dbl-64 0x7.ffffffffffffcp-1024 -0x7.ffffffffffffcp+52 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffffp+1020 : inexact
+= fma upward dbl-64 0x7.ffffffffffffcp-1024 -0x7.ffffffffffffcp+52 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffff8p+1020 : inexact
+= fma downward ldbl-96-intel 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff7ffp+1020L : inexact
+= fma tonearest ldbl-96-intel 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-96-intel 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff7ffp+1020L : inexact
+= fma upward ldbl-96-intel 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma downward ldbl-96-m68k 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff7ffp+1020L : inexact
+= fma tonearest ldbl-96-m68k 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-96-m68k 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff7ffp+1020L : inexact
+= fma upward ldbl-96-m68k 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma downward ldbl-128 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff7ffffffffffffff8p+1020L : inexact
+= fma tonearest ldbl-128 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-128 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff7ffffffffffffff8p+1020L : inexact
+= fma upward ldbl-128 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma downward ldbl-128ibm 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff7ffffffffffffcp+1020L : inexact
+= fma tonearest ldbl-128ibm 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-128ibm 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff7ffffffffffffcp+1020L : inexact
+= fma upward ldbl-128ibm 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffff8p+1020L : inexact
+fma 0x1.fffffffffffffp-1022 0x1.fffffffffffffp54 -0x1.fffffffffffffp1023
+= fma downward dbl-64 0x7.ffffffffffffcp-1024 0x7.ffffffffffffcp+52 -0xf.ffffffffffff8p+1020 : -0xf.ffffffffffff8p+1020 : inexact
+= fma tonearest dbl-64 0x7.ffffffffffffcp-1024 0x7.ffffffffffffcp+52 -0xf.ffffffffffff8p+1020 : -0xf.ffffffffffff8p+1020 : inexact
+= fma towardzero dbl-64 0x7.ffffffffffffcp-1024 0x7.ffffffffffffcp+52 -0xf.ffffffffffff8p+1020 : -0xf.ffffffffffffp+1020 : inexact
+= fma upward dbl-64 0x7.ffffffffffffcp-1024 0x7.ffffffffffffcp+52 -0xf.ffffffffffff8p+1020 : -0xf.ffffffffffffp+1020 : inexact
+= fma downward ldbl-96-intel 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma tonearest ldbl-96-intel 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-96-intel 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff7ffp+1020L : inexact
+= fma upward ldbl-96-intel 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff7ffp+1020L : inexact
+= fma downward ldbl-96-m68k 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma tonearest ldbl-96-m68k 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-96-m68k 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff7ffp+1020L : inexact
+= fma upward ldbl-96-m68k 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff7ffp+1020L : inexact
+= fma downward ldbl-128 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma tonearest ldbl-128 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-128 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff7ffffffffffffff8p+1020L : inexact
+= fma upward ldbl-128 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff7ffffffffffffff8p+1020L : inexact
+= fma downward ldbl-128ibm 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma tonearest ldbl-128ibm 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-128ibm 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff7ffffffffffffcp+1020L : inexact
+= fma upward ldbl-128ibm 0x7.ffffffffffffcp-1024L 0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff7ffffffffffffcp+1020L : inexact
+fma 0x1.fffffffffffffp-1022 -0x1.fffffffffffffp54 -0x1.fffffffffffffp1023 missing-errno
+= fma downward dbl-64 0x7.ffffffffffffcp-1024 -0x7.ffffffffffffcp+52 -0xf.ffffffffffff8p+1020 : minus_infty : inexact overflow errno-erange errno-erange-ok
+= fma tonearest dbl-64 0x7.ffffffffffffcp-1024 -0x7.ffffffffffffcp+52 -0xf.ffffffffffff8p+1020 : -0xf.ffffffffffff8p+1020 : inexact
+= fma towardzero dbl-64 0x7.ffffffffffffcp-1024 -0x7.ffffffffffffcp+52 -0xf.ffffffffffff8p+1020 : -0xf.ffffffffffff8p+1020 : inexact
+= fma upward dbl-64 0x7.ffffffffffffcp-1024 -0x7.ffffffffffffcp+52 -0xf.ffffffffffff8p+1020 : -0xf.ffffffffffff8p+1020 : inexact
+= fma downward ldbl-96-intel 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff801p+1020L : inexact
+= fma tonearest ldbl-96-intel 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-96-intel 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma upward ldbl-96-intel 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma downward ldbl-96-m68k 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff801p+1020L : inexact
+= fma tonearest ldbl-96-m68k 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-96-m68k 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma upward ldbl-96-m68k 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma downward ldbl-128 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8000000000000008p+1020L : inexact
+= fma tonearest ldbl-128 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-128 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma upward ldbl-128 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma downward ldbl-128ibm 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff80000000000004p+1020L : inexact
+= fma tonearest ldbl-128ibm 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma towardzero ldbl-128ibm 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+= fma upward ldbl-128ibm 0x7.ffffffffffffcp-1024L -0x7.ffffffffffffcp+52L -0xf.ffffffffffff8p+1020L : -0xf.ffffffffffff8p+1020L : inexact
+fma 0x1.fffffffffffffffep-16382 0x1.fffffffffffffffep65 0x1.fffffffffffffffep16383 missing-errno
+= fma downward ldbl-96-intel 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma tonearest ldbl-96-intel 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma towardzero ldbl-96-intel 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma upward ldbl-96-intel 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : plus_infty : inexact overflow errno-erange errno-erange-ok
+= fma downward ldbl-96-m68k 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma tonearest ldbl-96-m68k 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma towardzero ldbl-96-m68k 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma upward ldbl-96-m68k 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : plus_infty : inexact overflow errno-erange errno-erange-ok
+= fma downward ldbl-128 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma tonearest ldbl-128 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma towardzero ldbl-128 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma upward ldbl-128 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffff0000000000008p+16380L : inexact
+fma 0x1.fffffffffffffffep-16382 -0x1.fffffffffffffffep65 0x1.fffffffffffffffep16383
+= fma downward ldbl-96-intel 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.ffffffffffffffep+16380L : inexact
+= fma tonearest ldbl-96-intel 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma towardzero ldbl-96-intel 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.ffffffffffffffep+16380L : inexact
+= fma upward ldbl-96-intel 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma downward ldbl-96-m68k 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.ffffffffffffffep+16380L : inexact
+= fma tonearest ldbl-96-m68k 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma towardzero ldbl-96-m68k 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.ffffffffffffffep+16380L : inexact
+= fma upward ldbl-96-m68k 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma downward ldbl-128 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.ffffffffffffffeffffffffffff8p+16380L : inexact
+= fma tonearest ldbl-128 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+= fma towardzero ldbl-128 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.ffffffffffffffeffffffffffff8p+16380L : inexact
+= fma upward ldbl-128 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact
+fma 0x1.fffffffffffffffep-16382 0x1.fffffffffffffffep65 -0x1.fffffffffffffffep16383
+= fma downward ldbl-96-intel 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma tonearest ldbl-96-intel 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma towardzero ldbl-96-intel 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.ffffffffffffffep+16380L : inexact
+= fma upward ldbl-96-intel 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.ffffffffffffffep+16380L : inexact
+= fma downward ldbl-96-m68k 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma tonearest ldbl-96-m68k 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma towardzero ldbl-96-m68k 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.ffffffffffffffep+16380L : inexact
+= fma upward ldbl-96-m68k 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.ffffffffffffffep+16380L : inexact
+= fma downward ldbl-128 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma tonearest ldbl-128 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma towardzero ldbl-128 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.ffffffffffffffeffffffffffff8p+16380L : inexact
+= fma upward ldbl-128 0x7.fffffffffffffff8p-16384L 0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.ffffffffffffffeffffffffffff8p+16380L : inexact
+fma 0x1.fffffffffffffffep-16382 -0x1.fffffffffffffffep65 -0x1.fffffffffffffffep16383 missing-errno
+= fma downward ldbl-96-intel 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : minus_infty : inexact overflow errno-erange errno-erange-ok
+= fma tonearest ldbl-96-intel 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma towardzero ldbl-96-intel 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma upward ldbl-96-intel 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma downward ldbl-96-m68k 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : minus_infty : inexact overflow errno-erange errno-erange-ok
+= fma tonearest ldbl-96-m68k 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma towardzero ldbl-96-m68k 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma upward ldbl-96-m68k 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma downward ldbl-128 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffff0000000000008p+16380L : inexact
+= fma tonearest ldbl-128 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma towardzero ldbl-128 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+= fma upward ldbl-128 0x7.fffffffffffffff8p-16384L -0x3.fffffffffffffffcp+64L -0xf.fffffffffffffffp+16380L : -0xf.fffffffffffffffp+16380L : inexact
+fma 0x1.ffffffffffffffffffffffffffffp-16382 0x1.ffffffffffffffffffffffffffffp114 0x1.ffffffffffffffffffffffffffffp16383 missing-errno
+= fma downward ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L 0x7.fffffffffffffffffffffffffffcp+112L 0xf.fffffffffffffffffffffffffff8p+16380L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact
+= fma tonearest ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L 0x7.fffffffffffffffffffffffffffcp+112L 0xf.fffffffffffffffffffffffffff8p+16380L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact
+= fma towardzero ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L 0x7.fffffffffffffffffffffffffffcp+112L 0xf.fffffffffffffffffffffffffff8p+16380L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact
+= fma upward ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L 0x7.fffffffffffffffffffffffffffcp+112L 0xf.fffffffffffffffffffffffffff8p+16380L : plus_infty : inexact overflow errno-erange errno-erange-ok
+fma 0x1.ffffffffffffffffffffffffffffp-16382 -0x1.ffffffffffffffffffffffffffffp114 0x1.ffffffffffffffffffffffffffffp16383
+= fma downward ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L -0x7.fffffffffffffffffffffffffffcp+112L 0xf.fffffffffffffffffffffffffff8p+16380L : 0xf.fffffffffffffffffffffffffffp+16380L : inexact
+= fma tonearest ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L -0x7.fffffffffffffffffffffffffffcp+112L 0xf.fffffffffffffffffffffffffff8p+16380L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact
+= fma towardzero ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L -0x7.fffffffffffffffffffffffffffcp+112L 0xf.fffffffffffffffffffffffffff8p+16380L : 0xf.fffffffffffffffffffffffffffp+16380L : inexact
+= fma upward ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L -0x7.fffffffffffffffffffffffffffcp+112L 0xf.fffffffffffffffffffffffffff8p+16380L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact
+fma 0x1.ffffffffffffffffffffffffffffp-16382 0x1.ffffffffffffffffffffffffffffp114 -0x1.ffffffffffffffffffffffffffffp16383
+= fma downward ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L 0x7.fffffffffffffffffffffffffffcp+112L -0xf.fffffffffffffffffffffffffff8p+16380L : -0xf.fffffffffffffffffffffffffff8p+16380L : inexact
+= fma tonearest ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L 0x7.fffffffffffffffffffffffffffcp+112L -0xf.fffffffffffffffffffffffffff8p+16380L : -0xf.fffffffffffffffffffffffffff8p+16380L : inexact
+= fma towardzero ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L 0x7.fffffffffffffffffffffffffffcp+112L -0xf.fffffffffffffffffffffffffff8p+16380L : -0xf.fffffffffffffffffffffffffffp+16380L : inexact
+= fma upward ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L 0x7.fffffffffffffffffffffffffffcp+112L -0xf.fffffffffffffffffffffffffff8p+16380L : -0xf.fffffffffffffffffffffffffffp+16380L : inexact
+fma 0x1.ffffffffffffffffffffffffffffp-16382 -0x1.ffffffffffffffffffffffffffffp114 -0x1.ffffffffffffffffffffffffffffp16383 missing-errno
+= fma downward ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L -0x7.fffffffffffffffffffffffffffcp+112L -0xf.fffffffffffffffffffffffffff8p+16380L : minus_infty : inexact overflow errno-erange errno-erange-ok
+= fma tonearest ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L -0x7.fffffffffffffffffffffffffffcp+112L -0xf.fffffffffffffffffffffffffff8p+16380L : -0xf.fffffffffffffffffffffffffff8p+16380L : inexact
+= fma towardzero ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L -0x7.fffffffffffffffffffffffffffcp+112L -0xf.fffffffffffffffffffffffffff8p+16380L : -0xf.fffffffffffffffffffffffffff8p+16380L : inexact
+= fma upward ldbl-128 0x7.fffffffffffffffffffffffffffcp-16384L -0x7.fffffffffffffffffffffffffffcp+112L -0xf.fffffffffffffffffffffffffff8p+16380L : -0xf.fffffffffffffffffffffffffff8p+16380L : inexact
hypot 0 0
= hypot downward flt-32 0x0p+0f 0x0p+0f : 0x0p+0f : inexact-ok
= hypot tonearest flt-32 0x0p+0f 0x0p+0f : 0x0p+0f : inexact-ok
diff --git a/sysdeps/ieee754/dbl-64/s_fma.c b/sysdeps/ieee754/dbl-64/s_fma.c
index 716b412..278b690 100644
--- a/sysdeps/ieee754/dbl-64/s_fma.c
+++ b/sysdeps/ieee754/dbl-64/s_fma.c
@@ -117,7 +117,7 @@ __fma (double x, double y, double z)
very small, adjust them up to avoid spurious underflows,
rather than down. */
if (u.ieee.exponent + v.ieee.exponent
- <= IEEE754_DOUBLE_BIAS + DBL_MANT_DIG)
+ <= IEEE754_DOUBLE_BIAS + 2 * DBL_MANT_DIG)
{
if (u.ieee.exponent > v.ieee.exponent)
u.ieee.exponent += 2 * DBL_MANT_DIG + 2;
diff --git a/sysdeps/ieee754/ldbl-128/s_fmal.c b/sysdeps/ieee754/ldbl-128/s_fmal.c
index b13178f..5abc910 100644
--- a/sysdeps/ieee754/ldbl-128/s_fmal.c
+++ b/sysdeps/ieee754/ldbl-128/s_fmal.c
@@ -121,7 +121,7 @@ __fmal (long double x, long double y, long double z)
very small, adjust them up to avoid spurious underflows,
rather than down. */
if (u.ieee.exponent + v.ieee.exponent
- <= IEEE854_LONG_DOUBLE_BIAS + LDBL_MANT_DIG)
+ <= IEEE854_LONG_DOUBLE_BIAS + 2 * LDBL_MANT_DIG)
{
if (u.ieee.exponent > v.ieee.exponent)
u.ieee.exponent += 2 * LDBL_MANT_DIG + 2;
diff --git a/sysdeps/ieee754/ldbl-96/s_fmal.c b/sysdeps/ieee754/ldbl-96/s_fmal.c
index eec5a02..1232c9e 100644
--- a/sysdeps/ieee754/ldbl-96/s_fmal.c
+++ b/sysdeps/ieee754/ldbl-96/s_fmal.c
@@ -119,7 +119,7 @@ __fmal (long double x, long double y, long double z)
very small, adjust them up to avoid spurious underflows,
rather than down. */
if (u.ieee.exponent + v.ieee.exponent
- <= IEEE854_LONG_DOUBLE_BIAS + LDBL_MANT_DIG)
+ <= IEEE854_LONG_DOUBLE_BIAS + 2 * LDBL_MANT_DIG)
{
if (u.ieee.exponent > v.ieee.exponent)
u.ieee.exponent += 2 * LDBL_MANT_DIG + 2;
-----------------------------------------------------------------------
Summary of changes:
ChangeLog | 10 ++
NEWS | 2 +-
math/auto-libm-test-in | 18 +++
math/auto-libm-test-out | 256 +++++++++++++++++++++++++++++++++++++
sysdeps/ieee754/dbl-64/s_fma.c | 2 +-
sysdeps/ieee754/ldbl-128/s_fmal.c | 2 +-
sysdeps/ieee754/ldbl-96/s_fmal.c | 2 +-
7 files changed, 288 insertions(+), 4 deletions(-)
hooks/post-receive
--
GNU C Library master sources