Efficient Jump Ahead for 𝔽2-Linear Random Number Generators

H Haramoto, M Matsumoto… - INFORMS Journal …, 2008 - pubsonline.informs.org
H Haramoto, M Matsumoto, T Nishimura, F Panneton, P L'Ecuyer
INFORMS Journal on Computing, 2008pubsonline.informs.org
The fastest long-period random number generators currently available are based on linear
recurrences modulo 2. So far, software that provides multiple disjoint streams and
substreams has not been available for these generators because of the lack of efficient jump-
ahead facilities. In principle, it suffices to multiply the state (ak-bit vector) by an appropriate
k× k binary matrix to find the new state far ahead in the sequence. However, when k is large
(eg, for a generator such as the popular Mersenne twister, for which k= 19,937), this matrix …
The fastest long-period random number generators currently available are based on linear recurrences modulo 2. So far, software that provides multiple disjoint streams and substreams has not been available for these generators because of the lack of efficient jump-ahead facilities. In principle, it suffices to multiply the state (a k-bit vector) by an appropriate k × k binary matrix to find the new state far ahead in the sequence. However, when k is large (e.g., for a generator such as the popular Mersenne twister, for which k = 19,937), this matrix-vector multiplication is slow, and a large amount of memory is required to store the k × k matrix. In this paper, we provide a faster algorithm to jump ahead by a large number of steps in a linear recurrence modulo 2. The method uses much less than the k2 bits of memory required by the matrix method. It is based on polynomial calculus modulo the characteristic polynomial of the recurrence, and uses a sliding window algorithm for the multiplication.
INFORMS