[Numpy-svn] r4554 - trunk/numpy/random/mtrand

[numpy-svn at scipy.org]

Author: rkern
Date: 2007-12-05 23:18:54 -0600 (Wed, 05 Dec 2007)
New Revision: 4554

Modified:
   trunk/numpy/random/mtrand/distributions.c
Log:
Use a correct upper bound for the inversion search in binomial distributions.

Modified: trunk/numpy/random/mtrand/distributions.c
===================================================================
--- trunk/numpy/random/mtrand/distributions.c	2007-12-05 19:59:29 UTC (rev 4553)
+++ trunk/numpy/random/mtrand/distributions.c	2007-12-06 05:18:54 UTC (rev 4554)
@@ -389,7 +389,7 @@
         state->q = q = 1.0 - p;
         state->r = qn = exp(n * log(q));
         state->c = np = n*p;
-        state->m = bound = min(n, np + 10.0*sqrt(np));
+        state->m = bound = min(n, np + 10.0*sqrt(np*q + 1));
     } else
     {
         q = state->q;