We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 3e3e1f5 commit d7733b1Copy full SHA for d7733b1
1 file changed
src/PoissonRandom.jl
@@ -95,7 +95,7 @@ function procf(λ, K::Int, s::Float64)
95
c0 = 1 - b1 + 3 * b2 - 15 * c3
96
97
if K < 10
98
- px = -λ
+ px = -float(λ)
99
log_py = K * log(λ) - loggamma(K + 1) # log(K!) via loggamma
100
py = exp(log_py)
101
else
0 commit comments