This commit is contained in:
Mark Dickinson 2011-06-25 12:02:53 +02:00
commit 5b0c22ced3
1 changed files with 3 additions and 5 deletions

View File

@ -427,11 +427,9 @@ class Random(_random.Random):
# lambd: rate lambd = 1/mean
# ('lambda' is a Python reserved word)
random = self.random
u = random()
while u <= 1e-7:
u = random()
return -_log(u)/lambd
# we use 1-random() instead of random() to preclude the
# possibility of taking the log of zero.
return -_log(1.0 - self.random())/lambd
## -------------------- von Mises distribution --------------------