24小时热门版块排行榜

返回列表

当前主题已经存档。

buddy_84

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 806
散金: 300
红花: 1
帖子: 164
在线: 110.1小时
虫号: 283631
注册: 2006-10-08
专业: 原子和分子物理

[交流] 随机数发生器（Fortran77编写的，测试可用）

real function ranw(idum)
   Parameter (Mbig=2**30-2, Xinvers=1./Mbig)
   data ibit/ 1/
   Integer IX(55)
   save
   if (ibit.ne.0) then
      ibit=0
cc
cc    fill up the vector ix with some random integers, which are
cc    not all even
cc
      if (idum.eq.0) pause 'use nonzero value of idum'
      idum=abs(mod(idum,Mbig))
      ibit=0
      Ix(1)=871871
      Do i=2,55
         Ix(i)=mod(Ix(i-1)+idum,Ix(i-1))
         Ix(i)=max(mod(Ix(i),Mbig),idum)
      end do
      j=24
      k=55
cc
cc    warm up the generator
cc
      do i=1,1258
         Ix(k)=mod(Ix(k)+Ix(j),Mbig)
         j=j-1
         if (j.eq.0) j=55
         k=k-1
         if (k.eq.0) k=55
      end do
   end if
cc
cc this is where execution usually starts:
cc
   Ix(k)=mod(Ix(k)+Ix(j),Mbig)
   j=j-1
   if (j.eq.0) j=55
   k=k-1
   if (k.eq.0) k=55
   ranw=Ix(k)*Xinvers
   end

回复此楼

» 猜你喜欢

青椒八年已不青，大家都被折磨成啥样了？已经有11人回复
限项规定已经有9人回复
免疫学博士有名额，速联系已经有4人回复
交叉科学部支持青年基金，对三无青椒是个机会吗？已经有5人回复
国家基金申请书模板内插入图片不可调整大小? 已经有6人回复
国家级人才课题组招收2026年入学博士已经有5人回复
Fe3O4@SiO2合成已经有6人回复
青年基金C终止已经有4人回复
26申博求博导推荐-遥感图像处理方向已经有4人回复
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向已经有3人回复

1楼 2009-03-28 09:32:32

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

相关版块跳转我要订阅楼主 buddy_84 的主题更新

返回列表