From 11e77ffcbc4a46001ba6c509f069692b4c5d0735 Mon Sep 17 00:00:00 2001 From: Mark Dickinson Date: Sat, 4 Dec 2010 12:38:19 +0000 Subject: [PATCH] Merged revisions 87032 via svnmerge from svn+ssh://pythondev@svn.python.org/python/branches/py3k ........ r87032 | mark.dickinson | 2010-12-04 12:25:30 +0000 (Sat, 04 Dec 2010) | 3 lines Issue #10596: Fix float.__mod__ to have the same behaviour as float.__divmod__ with respect to signed zeros. ........ --- Lib/test/test_float.py | 20 ++++++++++++++++++++ Misc/NEWS | 4 ++++ Objects/floatobject.c | 18 ++++++++++++++---- 3 files changed, 38 insertions(+), 4 deletions(-) diff --git a/Lib/test/test_float.py b/Lib/test/test_float.py index fd8436a8fa..6d90637298 100644 --- a/Lib/test/test_float.py +++ b/Lib/test/test_float.py @@ -231,6 +231,26 @@ class GeneralFloatCases(unittest.TestCase): self.assertTrue(d == d, "{%r : None} not equal to itself" % f) + @requires_IEEE_754 + def test_float_mod(self): + # Check behaviour of % operator for IEEE 754 special cases. + # In particular, check signs of zeros. + mod = operator.mod + + self.assertEqualAndEqualSign(mod(-1.0, 1.0), 0.0) + self.assertEqualAndEqualSign(mod(-1e-100, 1.0), 1.0) + self.assertEqualAndEqualSign(mod(-0.0, 1.0), 0.0) + self.assertEqualAndEqualSign(mod(0.0, 1.0), 0.0) + self.assertEqualAndEqualSign(mod(1e-100, 1.0), 1e-100) + self.assertEqualAndEqualSign(mod(1.0, 1.0), 0.0) + + self.assertEqualAndEqualSign(mod(-1.0, -1.0), -0.0) + self.assertEqualAndEqualSign(mod(-1e-100, -1.0), -1e-100) + self.assertEqualAndEqualSign(mod(-0.0, -1.0), -0.0) + self.assertEqualAndEqualSign(mod(0.0, -1.0), -0.0) + self.assertEqualAndEqualSign(mod(1e-100, -1.0), -1.0) + self.assertEqualAndEqualSign(mod(1.0, -1.0), -0.0) + class FormatFunctionsTestCase(unittest.TestCase): diff --git a/Misc/NEWS b/Misc/NEWS index 1c7abf7efc..76ccf7ab0d 100644 --- a/Misc/NEWS +++ b/Misc/NEWS @@ -10,6 +10,10 @@ What's New in Python 3.1.4? Core and Builtins ----------------- +- Issue #10596: Fix float.__mod__ to have the same behaviour as + float.__divmod__ with respect to signed zeros. -4.0 % 4.0 should be + 0.0, not -0.0. + Library ------- diff --git a/Objects/floatobject.c b/Objects/floatobject.c index 865c960b88..6a2af743af 100644 --- a/Objects/floatobject.c +++ b/Objects/floatobject.c @@ -600,10 +600,20 @@ float_rem(PyObject *v, PyObject *w) #endif PyFPE_START_PROTECT("modulo", return 0) mod = fmod(vx, wx); - /* note: checking mod*wx < 0 is incorrect -- underflows to - 0 if wx < sqrt(smallest nonzero double) */ - if (mod && ((wx < 0) != (mod < 0))) { - mod += wx; + if (mod) { + /* ensure the remainder has the same sign as the denominator */ + if ((wx < 0) != (mod < 0)) { + mod += wx; + } + } + else { + /* the remainder is zero, and in the presence of signed zeroes + fmod returns different results across platforms; ensure + it has the same sign as the denominator; we'd like to do + "mod = wx * 0.0", but that may get optimized away */ + mod *= mod; /* hide "mod = +0" from optimizer */ + if (wx < 0.0) + mod = -mod; } PyFPE_END_PROTECT(mod) return PyFloat_FromDouble(mod); -- 2.40.0