@@ -7,7 +7,7 @@ export srand,
7
7
randn, randn!,
8
8
randbool,
9
9
AbstractRNG, RNG, MersenneTwister,
10
- randmtzig_exprnd
10
+ randexp
11
11
12
12
abstract AbstractRNG
13
13
@@ -984,23 +984,23 @@ randn(dims::Int...) = randn!(Array(Float64, dims...))
984
984
randn (rng:: MersenneTwister , dims:: Dims ) = randn! (rng, Array (Float64, dims))
985
985
randn (rng:: MersenneTwister , dims:: Int... ) = randn! (rng, Array (Float64, dims... ))
986
986
987
- @inline function randmtzig_exprnd (rng:: MersenneTwister = GLOBAL_RNG)
987
+ @inline function randexp (rng:: MersenneTwister = GLOBAL_RNG)
988
988
@inbounds begin
989
989
ri = rand_ui52 (rng)
990
990
idx = ri & 0xFF
991
991
x = ri* we[idx+ 1 ]
992
992
ri < ke[idx+ 1 ] && return x # 98.9% of the time we return here 1st try
993
- return randmtzig_exprnd_unlikely (rng, idx, x)
993
+ return randexp_unlikely (rng, idx, x)
994
994
end
995
995
end
996
996
997
- function randmtzig_exprnd_unlikely (rng, idx, x)
997
+ function randexp_unlikely (rng, idx, x)
998
998
@inbounds if idx == 0
999
999
return ziggurat_exp_r - log (rand (rng))
1000
1000
elseif (fe[idx] - fe[idx+ 1 ])* rand (rng) + fe[idx+ 1 ] < exp (- x)
1001
1001
return x # return from the triangular area
1002
1002
else
1003
- return randmtzig_exprnd (rng)
1003
+ return randexp (rng)
1004
1004
end
1005
1005
end
1006
1006
0 commit comments