From: Mark Dickinson <dickinsm@gmail.com>
Date: Sun, 10 Feb 2013 14:13:40 +0000 (+0000)
Subject: Issue #17149: Fix random.vonmisesvariate to always return results in [0, 2*math.pi].
X-Git-Tag: v2.7.4rc1~128
X-Git-Url: https://granicus.if.org/sourcecode?a=commitdiff_plain;h=9aaeb5e0c8b5946b305590eb85312c282a457098;p=python

Issue #17149: Fix random.vonmisesvariate to always return results in [0, 2*math.pi].
---

diff --git a/Lib/random.py b/Lib/random.py
index 36b956540b..bfa710be86 100644
--- a/Lib/random.py
+++ b/Lib/random.py
@@ -475,9 +475,9 @@ class Random(_random.Random):
 
         u3 = random()
         if u3 > 0.5:
-            theta = (mu % TWOPI) + _acos(f)
+            theta = (mu + _acos(f)) % TWOPI
         else:
-            theta = (mu % TWOPI) - _acos(f)
+            theta = (mu - _acos(f)) % TWOPI
 
         return theta
 
diff --git a/Lib/test/test_random.py b/Lib/test/test_random.py
index caa4f289ad..45391174e5 100644
--- a/Lib/test/test_random.py
+++ b/Lib/test/test_random.py
@@ -533,6 +533,20 @@ class TestDistributions(unittest.TestCase):
             self.assertAlmostEqual(s1/N, mu, 2)
             self.assertAlmostEqual(s2/(N-1), sigmasqrd, 2)
 
+    def test_von_mises_range(self):
+        # Issue 17149: von mises variates were not consistently in the
+        # range [0, 2*PI].
+        g = random.Random()
+        N = 100
+        for mu in 0.0, 0.1, 3.1, 6.2:
+            for kappa in 0.0, 2.3, 500.0:
+                for _ in range(N):
+                    sample = g.vonmisesvariate(mu, kappa)
+                    self.assertTrue(
+                        0 <= sample <= random.TWOPI,
+                        msg=("vonmisesvariate({}, {}) produced a result {} out"
+                             " of range [0, 2*pi]").format(mu, kappa, sample))
+
 class TestModule(unittest.TestCase):
     def testMagicConstants(self):
         self.assertAlmostEqual(random.NV_MAGICCONST, 1.71552776992141)
diff --git a/Misc/NEWS b/Misc/NEWS
index 39298b5278..8fa5cab674 100644
--- a/Misc/NEWS
+++ b/Misc/NEWS
@@ -202,6 +202,9 @@ Core and Builtins
 Library
 -------
 
+- Issue #17149: Fix random.vonmisesvariate to always return results in
+  the range [0, 2*math.pi].
+
 - Issue #1470548: XMLGenerator now works with UTF-16 and UTF-32 encodings.
 
 - Issue #6975: os.path.realpath() now correctly resolves multiple nested