]> git.ipfire.org Git - thirdparty/glibc.git/commitdiff
Fix sign of remquo zero remainder in round-downward mode (bug 17987).
authorJoseph Myers <joseph@codesourcery.com>
Tue, 17 Feb 2015 00:41:50 +0000 (00:41 +0000)
committerJoseph Myers <joseph@codesourcery.com>
Tue, 17 Feb 2015 00:41:50 +0000 (00:41 +0000)
Various remquo implementations produce a zero remainder with the wrong
sign (a zero remainder should always have the sign of the first
argument, as specified in IEEE 754) in round-downward mode, resulting
from the sign of 0 - 0.  This patch checks for zero results and fixes
their sign accordingly.

Tested for x86_64, x86, mips64 and powerpc.

[BZ #17987]
* sysdeps/ieee754/dbl-64/s_remquo.c (__remquo): Ensure sign of
zero result does not depend on the sign resulting from
subtraction.
* sysdeps/ieee754/dbl-64/wordsize-64/s_remquo.c (__remquo):
Likewise.
* sysdeps/ieee754/flt-32/s_remquof.c (__remquof): Likewise.
* sysdeps/ieee754/ldbl-128/s_remquol.c (__remquol): Likewise.
* sysdeps/ieee754/ldbl-128ibm/s_remquol.c (__remquol): Likewise.
* sysdeps/ieee754/ldbl-96/s_remquol.c (__remquol): Likewise.
* math/libm-test.inc (remquo_test_data): Add more tests.

ChangeLog
NEWS
math/libm-test.inc
sysdeps/ieee754/dbl-64/s_remquo.c
sysdeps/ieee754/dbl-64/wordsize-64/s_remquo.c
sysdeps/ieee754/flt-32/s_remquof.c
sysdeps/ieee754/ldbl-128/s_remquol.c
sysdeps/ieee754/ldbl-128ibm/s_remquol.c
sysdeps/ieee754/ldbl-96/s_remquol.c

index 540e997931a4a1598049e7ada235456ed4108344..f74ef2378cb54bea802a0dfae6c47e7d637755f8 100644 (file)
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,17 @@
+2015-02-16  Joseph Myers  <joseph@codesourcery.com>
+
+       [BZ #17987]
+       * sysdeps/ieee754/dbl-64/s_remquo.c (__remquo): Ensure sign of
+       zero result does not depend on the sign resulting from
+       subtraction.
+       * sysdeps/ieee754/dbl-64/wordsize-64/s_remquo.c (__remquo):
+       Likewise.
+       * sysdeps/ieee754/flt-32/s_remquof.c (__remquof): Likewise.
+       * sysdeps/ieee754/ldbl-128/s_remquol.c (__remquol): Likewise.
+       * sysdeps/ieee754/ldbl-128ibm/s_remquol.c (__remquol): Likewise.
+       * sysdeps/ieee754/ldbl-96/s_remquol.c (__remquol): Likewise.
+       * math/libm-test.inc (remquo_test_data): Add more tests.
+
 2015-02-16  Paul Eggert  <eggert@cs.ucla.edu>
 
        * manual/time.texi (TZ Variable): glibc no longer comes with tzdata.
diff --git a/NEWS b/NEWS
index 4f4f740af0773804c94337aee0b0c29b5edc1d54..1c5e5909c4dd1d8ab6fb704296ca883bab4f3956 100644 (file)
--- a/NEWS
+++ b/NEWS
@@ -10,7 +10,7 @@ Version 2.22
 * The following bugs are resolved with this release:
 
   4719, 15467, 15790, 16560, 17569, 17792, 17912, 17932, 17944, 17949,
-  17964, 17965, 17967, 17969, 17978.
+  17964, 17965, 17967, 17969, 17978, 17987.
 \f
 Version 2.21
 
index 1cf80e3379e4af31b5430821d3c6df2ff6e77b16..85d8b676f796cdc06eeffb73eef395a322412b6b 100644 (file)
@@ -8818,6 +8818,16 @@ static const struct test_ffI_f1_data remquo_test_data[] =
     TEST_ffI_f1 (remquo, -1, -max_value / 4, -1, 0, NO_INEXACT_EXCEPTION),
     TEST_ffI_f1 (remquo, -1, max_value / 8, -1, 0, NO_INEXACT_EXCEPTION),
     TEST_ffI_f1 (remquo, -1, -max_value / 8, -1, 0, NO_INEXACT_EXCEPTION),
+
+    TEST_ffI_f1 (remquo, max_value, max_value / 2, plus_zero, 2, NO_INEXACT_EXCEPTION),
+    TEST_ffI_f1 (remquo, max_value, -max_value / 2, plus_zero, -2, NO_INEXACT_EXCEPTION),
+    TEST_ffI_f1 (remquo, -max_value, max_value / 2, minus_zero, -2, NO_INEXACT_EXCEPTION),
+    TEST_ffI_f1 (remquo, -max_value, -max_value / 2, minus_zero, 2, NO_INEXACT_EXCEPTION),
+
+    TEST_ffI_f1 (remquo, 2, 1, plus_zero, 2, NO_INEXACT_EXCEPTION),
+    TEST_ffI_f1 (remquo, 2, -1, plus_zero, -2, NO_INEXACT_EXCEPTION),
+    TEST_ffI_f1 (remquo, -2, 1, minus_zero, -2, NO_INEXACT_EXCEPTION),
+    TEST_ffI_f1 (remquo, -2, -1, minus_zero, 2, NO_INEXACT_EXCEPTION),
   };
 
 static void
index e07efa89f83bc5093c926478ec8069f58befa9f1..b081d26b0dd460c0bfb631728e814e90cd43529e 100644 (file)
@@ -101,6 +101,9 @@ __remquo (double x, double y, int *quo)
 
   *quo = qs ? -cquo : cquo;
 
+  /* Ensure correct sign of zero result in round-downward mode.  */
+  if (x == 0.0)
+    x = 0.0;
   if (sx)
     x = -x;
   return x;
index ab561780c88fe1d269fa351104ac1aa6fe348765..304e9d0fb4438e05a47d42ae3f3e371db1586aa6 100644 (file)
@@ -100,6 +100,9 @@ __remquo (double x, double y, int *quo)
 
   *quo = qs ? -cquo : cquo;
 
+  /* Ensure correct sign of zero result in round-downward mode.  */
+  if (x == 0.0)
+    x = 0.0;
   if (sx)
     x = -x;
   return x;
index 04944fdcd386eece7613898641749f700bab4d09..36cf359b9d4e90a8dd655873d3d873cdef28240a 100644 (file)
@@ -100,6 +100,9 @@ __remquof (float x, float y, int *quo)
 
   *quo = qs ? -cquo : cquo;
 
+  /* Ensure correct sign of zero result in round-downward mode.  */
+  if (x == 0.0f)
+    x = 0.0f;
   if (sx)
     x = -x;
   return x;
index da175a14f67d6b63a8fd803c8ae10155e3da269e..82cdf1a5894ee2205ab7744ff3cf89b95fc90956 100644 (file)
@@ -102,6 +102,9 @@ __remquol (long double x, long double y, int *quo)
 
   *quo = qs ? -cquo : cquo;
 
+  /* Ensure correct sign of zero result in round-downward mode.  */
+  if (x == 0.0L)
+    x = 0.0L;
   if (sx)
     x = -x;
   return x;
index ac228796b0d0013288c78b85691b3c4680760634..57c50590c90455b6f4bf7548de8df63b833a4a5e 100644 (file)
@@ -107,6 +107,9 @@ __remquol (long double x, long double y, int *quo)
 
   *quo = qs ? -cquo : cquo;
 
+  /* Ensure correct sign of zero result in round-downward mode.  */
+  if (x == 0.0L)
+    x = 0.0L;
   if (sx)
     x = -x;
   return x;
index 1462b54456900610921a9f6e2aea3c8e7cc28f9e..8d6b10ed619ec9b4b2bb387ce166e3f662c98c8d 100644 (file)
@@ -101,6 +101,9 @@ __remquol (long double x, long double p, int *quo)
 
   *quo = qs ? -cquo : cquo;
 
+  /* Ensure correct sign of zero result in round-downward mode.  */
+  if (x == 0.0L)
+    x = 0.0L;
   if (sx)
     x = -x;
   return x;