24小时热门版块排行榜

返回列表

ruthxu

金虫 (正式写手)

应助: 2 (幼儿园)
金币: 1109.2
帖子: 319
在线: 138小时
虫号: 1159022

[交流] Gromacs 中的随机数发生器

各位虫友，大家在用groamcs 的langevin dynamic simulation （sd 算法）时，有没有注意到随机数发生器的原理？
原来一直没有关注过随机数发生器，最近编程用到才想起Gromacs中的随机数发生器，但原理现在也不是很懂。代码中（Gmx_random.c)：
gmx_rng_uniform_uint32(gmx_rng_t rng)
{
  unsigned int y;

  if(rng->mti==624)
gmx_rng_update(rng);
  y=rng->mt[rng->mti++];

  y ^= (y >> 11);
  y ^= (y << 7) & 0x9d2c5680UL;
  y ^= (y << 15) & 0xefc60000UL;
  y ^= (y >> 18);

  return y;
}

看了一些资料，似乎是twister中的tempering，但不能确定，有没有虫友对随机数比较懂，给介绍一下呢？
谢谢！

回复此楼