rand function is slow with a parfor loop

3 ビュー (過去 30 日間)
James 2011 年 11 月 18 日
I'm trying to multithread a stochastic simulation using the parallel computing toolbox. I've profiled the program extensively to ensure that it is optimized. However, when I'm using PCT with 4 local workers, I find that the rand function consumes an inordinate amount of computational time far beyond what it uses in a single threaded program. Here is a representation of the code in question.
if(~(rand <= Pr))
I am stumped as to why this would be 10-15x slower than the single threaded case, unless rand is a synchronized function. I'm using a MacPro with 16 GB RAM and 2 quad core processors (MATLAB R2008b) for reference.

回答 (3 件)

Alex 2011 年 11 月 18 日
All computer based random number generators need a 'seed'. I don't know about Matlab, but this seed, if not provided directly by the user, is based off of the current system time. (by using the same 'seed', you can then get repeatable results, even in a 'random' environment).
I think, 2008b does not have the updated random generator methods.
http://www.mathworks.com/help/techdoc/math/bsn94u0-1.html This link talks about some of the older changes that have been made to the rand function.

Walter Roberson
Walter Roberson 2011 年 11 月 18 日
2008b does have the updated random number generator methods.
MATLAB uses a constant seed for each session, but a different seed for each parallel worker. This is described more in Peter Perkin's posting at http://www.mathworks.com/matlabcentral/newsreader/view_thread/299841

James 2011 年 11 月 18 日
Each worker uses a provided seed, not the default. Even if they used the worker ID as a seed, that wouldn't explain why calling rand is much slower in a parfor loop than a for loop.
  1 件のコメント
Konrad Malkowski
Konrad Malkowski 2011 年 11 月 30 日
James could you provide a more complete example that reproduces the issue?


Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by