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 d9e3c0e commit 7a82d20Copy full SHA for 7a82d20
src/extras/random.jl
@@ -82,16 +82,15 @@ end
82
# normal variates
83
84
function randn(rng::AbstractRNG, ::Type{DoubleFloat{T}}) where {T<:IEEEFloat}
85
- urand1, urand2 = rand(rng, DoubleFloat{T}, 2)
86
- urand1 = urand1 + urand1 - 1
87
- urand2 = urand2 + urand2 - 1
88
- s = urand1*urand1 + urand2*urand2
+ urand1, urand2, s = ntuple(i -> zero(DoubleFloat{T}), Val{3}())
89
90
- while s >= 1 || s === 0
+ while true
91
urand1, urand2 = rand(rng, DoubleFloat{T}, 2)
92
urand1 = urand1 + urand1 - 1
93
urand2 = urand2 + urand2 - 1
94
s = urand1*urand1 + urand2*urand2
+
+ (s >= 1 || s === 0) || break
95
end
96
97
s = sqrt( -log(s) / s )
0 commit comments