Получения RNG семян в Джулии

В Джулии, семя для глобальной RNG может быть установлено сПолучения RNG семян в Джулии

srand(SEED)

Как можно получить семена или текущее состояние глобальной RNG, например снова в более позднем пункте?

Цель состоит в том, чтобы получить состояние RNG в любой момент времени и воссоздать его на другом сеансе, не зная начального семени или всех вызовов RNG, которые произошли в то же время.

Например, R позволяет получить доступ к текущему семени через

.Random.seed

Я надеялся, что эквивалентный способ будет существовать в Джулию.

источник

2014-10-17 Julian

[Read The Fabulous Руководство] (http://docs.julialang.org/en/release-0.3/stdlib/base/#random- номера). – jurgemaister

@jurgemaister Эта часть относится к установке семени.Я не мог найти упоминания о том, как вернуть его. – Julian

@jurgemaister: вы, кажется, не понимаете разницу между установкой семени и потом выяснением того, на что было установлено семя. – DSM

Очевидным решением является сохранение начального значения до того, как вы позвоните по телефону srand(seed).

В качестве альтернативы, если вы знаете используемый RNG и не криптографически безопасны, вы можете рассчитать значения из псевдослучайных чисел, которые он производит. Например, см Cracking a linear congruential generator

источник

2014-10-17 12:14:23 rossum

Генерация случайных чисел в Джулии использует библиотеку Мерсерн Твистер. – jurgemaister

Затем см. [Реверсивная функция темпера MERENNE Twister RNG] (http://b10l.com/reversing-the-mersenne-twister-rng-temper-function/) и тому подобное. – rossum

Base.Random.RANDOM_SEED является вашим другом для получения семян:

julia> srand(10) 

julia> Base.Random.RANDOM_SEED 
1-element Array{Uint32,1}: 
0x0000000a 

julia> srand(1) 

julia> Base.Random.RANDOM_SEED 
1-element Array{Uint32,1}: 
0x00000001 

julia> srand(0xF) 

julia> Base.Random.RANDOM_SEED 
1-element Array{Uint32,1}: 
0x0000000f

Это не документировано, но the source достаточно легко читать. Я не уверен, как получить текущее состояние RNG, но оно может быть в модуле dSFMT.

источник

2014-10-17 15:02:43 IainDunning

Теперь это доступно через 'Base.Random.GLOBAL_RNG.seed', а не' Base.Random.RANDOM_SEED'. – Chris

Использование специализированного MersenneTwister с явной переменной (вместо скрытой глобальной, предоставляемой функциями случайных значений по умолчанию) может быть обеспечена функциональность вам потребуется:

newmt = Base.Random.MersenneTwister(123) 
randvec = [rand(newmt) for i=1:100] 
# save state now 
savestate = deepcopy(newmt.state) 
randvec2 = [rand(newmt) for i=1:20] 
# rewind state to old state 
newmt.state = savestate 
randvec3 = [rand(newmt) for i=1:20] 
if randvec2==randvec3 
    println("works!") 
end

deepcopy бросил меня на секунду. Кроме того, было бы легче получить доступ к глобальному случайному состоянию генератора, но это может потребовать ccall ИНГОВ библиотеки libdSFMT (см random.jl и dSFMT.jl в Base.

источник

2014-10-19 22:29:11

Для лучшего контроля при использовании генератора случайных функций и, в общем,

RND = srand(0) 
function coolfun() 
    println(RND.idx) 
    output = srand(RND, 100) 
    ... 
end

источник

2017-12-19 12:29:19

Вы должны получить семена, как этот

reinterpret(Int32, Base.Random.GLOBAL_RNG.seed)

Тест:

julia> srand(123456789); 

julia> reinterpret(Int32, Base.Random.GLOBAL_RNG.seed) 
1-element Array{Int32,1}: 
123456789

Для сохранения восстановления полного состояния rng вы можете сделать простую вещь и просто сохранить весь объект Base.Random.GLOBAL_RNG. Простым способом будет использование пакета JLD.

В моем личном пакете я вручную сохраняю/читаю состояние rng на HDF5, см. here.

EDIT: Это, конечно, более явный вариант @ IainDunning Ответим

источник

2017-12-20 07:00:47 crstnbr

Получения RNG семян в Джулии

ответ

Смежные вопросы