Как я могу получить текущее семя генератора случайных чисел NumPy?

Question

Следующие импортирует NumPy и устанавливает семя.

import numpy as np 
np.random.seed(42) 

Однако, я не заинтересован в установке семени, но больше в его чтении. random.get_state(), похоже, не содержит семена. documentation не дает явного ответа.

Как извлечь текущее семя, используемое numpy.random, если я не установил его вручную?

Я хочу использовать текущее семя для переноса на следующую итерацию процесса.

Answer 1

Короткий ответ: вы просто не можете (по крайней мере, не вообще).

Mersenne Twister ГСЧ, используемый NumPy имеет 2 -1 возможные внутренние состояния, в то время как один 64-битное целое имеет только 2 возможные значения. Поэтому невозможно сопоставить каждое состояние RNG с уникальным целым семенем.

Вы можете получить и установить внутреннее состояние ГСЧ непосредственно с помощью np.random.get_state и np.random.set_state. Вывод get_state является кортежем, вторым элементом которого является массив из (624,) из 32-битных целых чисел. Этот массив имеет более чем достаточно бит для представления каждого возможного внутреннего состояния RNG (2 ^{624 * 32}> 2 -1).

Кортеж, возвращенный get_state, может использоваться как семена, чтобы создать воспроизводимые последовательности случайных чисел. Например:

import numpy as np 

# randomly initialize the RNG from some platform-dependent source of entropy 
np.random.seed(None) 

# get the initial state of the RNG 
st0 = np.random.get_state() 

# draw some random numbers 
print(np.random.randint(0, 100, 10)) 
# [ 8 76 76 33 77 26 3 1 68 21] 

# set the state back to what it was originally 
np.random.set_state(st0) 

# draw again 
print(np.random.randint(0, 100, 10)) 
# [ 8 76 76 33 77 26 3 1 68 21] 

Answer 2

Проверьте первый элемент массива, возвращенного np.random.get_state(), кажется, именно случайное семя мне.