Create variable number of random streams (unique, deterministic, independent)

Question

AJ 2022-4-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1696515-create-variable-number-of-random-streams-unique-deterministic-independent

评论： AJ 2022-4-14

For our simulations, we wish to have repeatable (pseudo-)random number generation.

The problem with using a "global " random stream is that it is possibly shared with other MATLAB processes, which can corrput (make non-repeatable) the random generation.

Now, this is solved in part by the 'NumStreams' argument to RandStream.create, e.g.:

[s1,s2,s3,s4] = RandStream.create('mrg32k3a','NumStreams',4)

This creates four unique, deterministic, independent random number streams.

The problem: The left hand side is fixed at four elements.

The question: How do I generalize this to handle a variable number of random number streams? For example:

[???] = RandStream.create('mrg32k3a','NumStreams',K)

I can think of some ad hoc approaches (e.g. using one stream to generate seeds for the remaining streams) but I worry about true independence of the streams.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Steven Lord 2022-4-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1696515-create-variable-number-of-random-streams-unique-deterministic-independent#answer_942820

See the "Substreams" and "Restore State of Random Number Generator to Reproduce Output" sections on this documentation page for more information.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

AJ 2022-4-14

Thanks, glad to know it was right there the whole time. : )

And thanks, Steven, for your plethora of contributions over the years; you've been a constant force of good for as long as I can remember using Answers, and the forum that preceded it.

请先登录，再进行评论。

Answer 2

AJ 2022-4-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1696515-create-variable-number-of-random-streams-unique-deterministic-independent#answer_942785

在 MATLAB Online 中打开

With some experimentation, the use of cell arrays is a one solution:

K = 4;
myStreams = cell(1,K);
[myStreams{:}] = RandStream.create('mrg32k3a','NumStreams',K)

I would have preferred a standard array (rather than a cell array), but "Concatenation of RandStream objects is not allowed. Use a cell array to contain multiple objects."