24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

tc1788

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 49.4
散金: 3
红花: 1
帖子: 253
在线: 113.4小时
虫号: 760748
注册: 2009-04-30
专业: 金属材料的合金相、相变及

[交流] 【求助】产生0~1之间随机数的函数

产生0~1之间随机数的函数大家帮帮忙好吗

回复此楼

» 猜你喜欢

1楼 2009-07-14 21:48:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

liuyalihd

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 117.7
散金: 577
帖子: 59
在线: 8.7小时
虫号: 531310
注册: 2008-03-23
专业: 计算机体系结构

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
woshilsh(金币+2,VIP+0): 9-6 16:18

给你段代码看看吧！

#define IADD 453806245
#define IMUL 314159269
#define MASK 2147483647
#define SCALE  0.4656612873e-9

int randSeedP = 17;

void InitRand (int randSeedI)
{
  struct timeval tv;

  if (randSeedI != 0) randSeedP = randSeedI;
  else {
gettimeofday (&tv, 0);
randSeedP = tv.tv_usec;
  }
}

real RandR ()
{
  randSeedP = (randSeedP * IMUL + IADD) & MASK;
  return (randSeedP * SCALE);
}

#if NDIM == 2

void VRand (VecR *p)
{
  real s;

  s = 2. * M_PI * RandR ();
  p->x = cos (s);
  p->y = sin (s);
}

#elif NDIM == 3

void VRand (VecR *p)
{
  real s, x, y;

  s = 2.;
  while (s > 1.) {
x = 2. * RandR () - 1.;
y = 2. * RandR () - 1.;
s = Sqr (x) + Sqr (y);
  }
  p->z = 1. - 2. * s;
  s = 2. * sqrt (1. - s);
  p->x = s * x;
  p->y = s * y;
}

#endif