@@ -214,12 +214,12 @@ rng_native_52(::MersenneTwister) = Float64
214
214
# ### helper functions
215
215
216
216
# precondition: !mt_empty(r)
217
- rand_inbounds (r:: MersenneTwister , :: Close1Open2_64 ) = mt_pop! (r)
218
- rand_inbounds (r:: MersenneTwister , :: CloseOpen_64 = CloseOpen ()) =
219
- rand_inbounds (r, Close1Open2 ()) - 1.0
217
+ rand_inbounds (r:: MersenneTwister , :: CloseOpen12_64 ) = mt_pop! (r)
218
+ rand_inbounds (r:: MersenneTwister , :: CloseOpen01_64 = CloseOpen01 ()) =
219
+ rand_inbounds (r, CloseOpen12 ()) - 1.0
220
220
221
221
rand_inbounds (r:: MersenneTwister , :: UInt52Raw{T} ) where {T<: BitInteger } =
222
- reinterpret (UInt64, rand_inbounds (r, Close1Open2 ())) % T
222
+ reinterpret (UInt64, rand_inbounds (r, CloseOpen12 ())) % T
223
223
224
224
function rand (r:: MersenneTwister , x:: SamplerTrivial{UInt52Raw{UInt64}} )
225
225
reserve_1 (r)
238
238
239
239
# ### floats
240
240
241
- rand (r:: MersenneTwister , sp:: SamplerTrivial{Close1Open2_64 } ) =
241
+ rand (r:: MersenneTwister , sp:: SamplerTrivial{CloseOpen12_64 } ) =
242
242
(reserve_1 (r); rand_inbounds (r, sp[]))
243
243
244
244
# ### integers
@@ -323,7 +323,7 @@ function _rand_max383!(r::MersenneTwister, A::UnsafeView{Float64}, I::FloatInter
323
323
@gc_preserve r unsafe_copyto! (A. ptr+ m* sizeof (Float64), pointer (r. vals), n- m)
324
324
r. idx = n- m
325
325
end
326
- if I isa CloseOpen
326
+ if I isa CloseOpen01
327
327
for i= 1 : n
328
328
A[i] -= 1.0
329
329
end
@@ -332,10 +332,10 @@ function _rand_max383!(r::MersenneTwister, A::UnsafeView{Float64}, I::FloatInter
332
332
end
333
333
334
334
335
- fill_array! (s:: DSFMT_state , A:: Ptr{Float64} , n:: Int , :: CloseOpen_64 ) =
335
+ fill_array! (s:: DSFMT_state , A:: Ptr{Float64} , n:: Int , :: CloseOpen01_64 ) =
336
336
dsfmt_fill_array_close_open! (s, A, n)
337
337
338
- fill_array! (s:: DSFMT_state , A:: Ptr{Float64} , n:: Int , :: Close1Open2_64 ) =
338
+ fill_array! (s:: DSFMT_state , A:: Ptr{Float64} , n:: Int , :: CloseOpen12_64 ) =
339
339
dsfmt_fill_array_close1_open2! (s, A, n)
340
340
341
341
@@ -386,10 +386,10 @@ mask128(u::UInt128, ::Type{Float32}) =
386
386
(u & 0x007fffff007fffff007fffff007fffff ) | 0x3f8000003f8000003f8000003f800000
387
387
388
388
for T in (Float16, Float32)
389
- @eval function rand! (r:: MersenneTwister , A:: Array{$T} , :: SamplerTrivial{Close1Open2 {$T}} )
389
+ @eval function rand! (r:: MersenneTwister , A:: Array{$T} , :: SamplerTrivial{CloseOpen12 {$T}} )
390
390
n = length (A)
391
391
n128 = n * sizeof ($ T) ÷ 16
392
- _rand! (r, A, 2 * n128, Close1Open2 ())
392
+ _rand! (r, A, 2 * n128, CloseOpen12 ())
393
393
@gc_preserve A begin
394
394
A128 = UnsafeView {UInt128} (pointer (A), n128)
395
395
for i in 1 : n128
@@ -414,8 +414,8 @@ for T in (Float16, Float32)
414
414
A
415
415
end
416
416
417
- @eval function rand! (r:: MersenneTwister , A:: Array{$T} , :: SamplerTrivial{CloseOpen {$T}} )
418
- rand! (r, A, Close1Open2 ($ T))
417
+ @eval function rand! (r:: MersenneTwister , A:: Array{$T} , :: SamplerTrivial{CloseOpen01 {$T}} )
418
+ rand! (r, A, CloseOpen12 ($ T))
419
419
I32 = one (Float32)
420
420
for i in eachindex (A)
421
421
@inbounds A[i] = Float32 (A[i])- I32 # faster than "A[i] -= one(T)" for T==Float16
@@ -430,7 +430,7 @@ function rand!(r::MersenneTwister, A::UnsafeView{UInt128}, ::SamplerType{UInt128
430
430
n:: Int = length (A)
431
431
i = n
432
432
while true
433
- rand! (r, UnsafeView {Float64} (A. ptr, 2 i), Close1Open2 ())
433
+ rand! (r, UnsafeView {Float64} (A. ptr, 2 i), CloseOpen12 ())
434
434
n < 5 && break
435
435
i = 0
436
436
while n- i >= 5
0 commit comments