@@ -284,12 +284,12 @@ rng_native_52(::MersenneTwister) = Float64
284
284
# ### helper functions
285
285
286
286
# precondition: !mt_empty(r)
287
- rand_inbounds (r:: MersenneTwister , :: Close1Open2_64 ) = mt_pop! (r)
288
- rand_inbounds (r:: MersenneTwister , :: CloseOpen_64 = CloseOpen ()) =
289
- rand_inbounds (r, Close1Open2 ()) - 1.0
287
+ rand_inbounds (r:: MersenneTwister , :: CloseOpen12_64 ) = mt_pop! (r)
288
+ rand_inbounds (r:: MersenneTwister , :: CloseOpen01_64 = CloseOpen01 ()) =
289
+ rand_inbounds (r, CloseOpen12 ()) - 1.0
290
290
291
291
rand_inbounds (r:: MersenneTwister , :: UInt52Raw{T} ) where {T<: BitInteger } =
292
- reinterpret (UInt64, rand_inbounds (r, Close1Open2 ())) % T
292
+ reinterpret (UInt64, rand_inbounds (r, CloseOpen12 ())) % T
293
293
294
294
function rand (r:: MersenneTwister , x:: SamplerTrivial{UInt52Raw{UInt64}} )
295
295
reserve_1 (r)
308
308
309
309
# ### floats
310
310
311
- rand (r:: MersenneTwister , sp:: SamplerTrivial{Close1Open2_64 } ) =
311
+ rand (r:: MersenneTwister , sp:: SamplerTrivial{CloseOpen12_64 } ) =
312
312
(reserve_1 (r); rand_inbounds (r, sp[]))
313
313
314
314
# ### integers
@@ -377,7 +377,7 @@ function _rand_max383!(r::MersenneTwister, A::UnsafeView{Float64}, I::FloatInter
377
377
@gc_preserve r unsafe_copyto! (A. ptr+ m* sizeof (Float64), pointer (r. vals), n- m)
378
378
r. idxF = n- m
379
379
end
380
- if I isa CloseOpen
380
+ if I isa CloseOpen01
381
381
for i= 1 : n
382
382
A[i] -= 1.0
383
383
end
@@ -386,10 +386,10 @@ function _rand_max383!(r::MersenneTwister, A::UnsafeView{Float64}, I::FloatInter
386
386
end
387
387
388
388
389
- fill_array! (s:: DSFMT_state , A:: Ptr{Float64} , n:: Int , :: CloseOpen_64 ) =
389
+ fill_array! (s:: DSFMT_state , A:: Ptr{Float64} , n:: Int , :: CloseOpen01_64 ) =
390
390
dsfmt_fill_array_close_open! (s, A, n)
391
391
392
- fill_array! (s:: DSFMT_state , A:: Ptr{Float64} , n:: Int , :: Close1Open2_64 ) =
392
+ fill_array! (s:: DSFMT_state , A:: Ptr{Float64} , n:: Int , :: CloseOpen12_64 ) =
393
393
dsfmt_fill_array_close1_open2! (s, A, n)
394
394
395
395
@@ -440,10 +440,10 @@ mask128(u::UInt128, ::Type{Float32}) =
440
440
(u & 0x007fffff007fffff007fffff007fffff ) | 0x3f8000003f8000003f8000003f800000
441
441
442
442
for T in (Float16, Float32)
443
- @eval function rand! (r:: MersenneTwister , A:: Array{$T} , :: SamplerTrivial{Close1Open2 {$T}} )
443
+ @eval function rand! (r:: MersenneTwister , A:: Array{$T} , :: SamplerTrivial{CloseOpen12 {$T}} )
444
444
n = length (A)
445
445
n128 = n * sizeof ($ T) ÷ 16
446
- _rand! (r, A, 2 * n128, Close1Open2 ())
446
+ _rand! (r, A, 2 * n128, CloseOpen12 ())
447
447
@gc_preserve A begin
448
448
A128 = UnsafeView {UInt128} (pointer (A), n128)
449
449
for i in 1 : n128
@@ -468,8 +468,8 @@ for T in (Float16, Float32)
468
468
A
469
469
end
470
470
471
- @eval function rand! (r:: MersenneTwister , A:: Array{$T} , :: SamplerTrivial{CloseOpen {$T}} )
472
- rand! (r, A, Close1Open2 ($ T))
471
+ @eval function rand! (r:: MersenneTwister , A:: Array{$T} , :: SamplerTrivial{CloseOpen01 {$T}} )
472
+ rand! (r, A, CloseOpen12 ($ T))
473
473
I32 = one (Float32)
474
474
for i in eachindex (A)
475
475
@inbounds A[i] = Float32 (A[i])- I32 # faster than "A[i] -= one(T)" for T==Float16
@@ -484,7 +484,7 @@ function rand!(r::MersenneTwister, A::UnsafeView{UInt128}, ::SamplerType{UInt128
484
484
n:: Int = length (A)
485
485
i = n
486
486
while true
487
- rand! (r, UnsafeView {Float64} (A. ptr, 2 i), Close1Open2 ())
487
+ rand! (r, UnsafeView {Float64} (A. ptr, 2 i), CloseOpen12 ())
488
488
n < 5 && break
489
489
i = 0
490
490
while n- i >= 5
0 commit comments