24小时热门版块排行榜    

查看: 360  |  回复: 1

ruthxu

金虫 (正式写手)


[交流] 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,但不能确定,有没有虫友对随机数比较懂,给介绍一下呢?
谢谢!
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
嘉特2楼
2012-12-18 21:00   回复  
ruthxu(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 ruthxu 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见