rand function is slow with a parfor loop
2 次查看(过去 30 天)
显示 更早的评论
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.
0 个评论
回答(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.
0 个评论
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
0 个评论
James
2011-11-18
1 个评论
Konrad Malkowski
2011-11-30
James could you provide a more complete example that reproduces the issue?
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Parallel for-Loops (parfor) 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!