我正在编写一个运行Monte Carlo模拟的程序;具体来说,我使用Metropolis算法。该程序需要生成可能达数十亿的“随机”数。我知道Mersenne twister在Monte Carlo模拟中非常流行,但我想确保以最佳方式初始化发生器。
目前,我正在使用以下方法计算32位种子:
mt19937_64 prng; //pseudo random number generator
unsigned long seed; //store seed so that every run can follow the same sequence
unsigned char seed_count; //to help keep seeds from repeating because of temporal proximity
unsigned long genSeed() {
return ( static_cast<unsigned long>(time(NULL)) << 16 )
| ( (static_cast<unsigned long>(clock()) & 0xFF) << 8 )
| ( (static_cast<unsigned long>(seed_count++) & 0xFF) );
}
//...
seed = genSeed();
prng.seed(seed);
我有一种感觉,有更好的方法来确保不重复的新种子,而且我相当确定mt19937_64可以用超过32位的数值进行初始化。是否有人有什么建议吗?