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

回复此楼

» 猜你喜欢

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

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

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

返回列表