24小时热门版块排行榜

返回列表

当前主题已经存档。

buddy_84

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 792
散金: 300
红花: 1
帖子: 166
在线: 111小时
虫号: 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 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[教师之家] 上海大学实验技术岗位非升即走 +10	嘻嘻哈哈乐呵呵 2026-05-15	11/550	2026-05-19 10:03 by 嘻嘻哈哈乐呵呵
[基金申请] 同样的基金本子，换个专家直接从C变A！ +3	国自然国社科中 2026-05-19	3/150	2026-05-19 08:50 by Equinoxhua
[考博] 博士申请 +5	星…… 2026-05-18	6/300	2026-05-18 23:49 by 糊糊涂涂好
[基金申请] 国自然上会要求 +5	无名者登山 2026-05-18	9/450	2026-05-18 17:50 by BlakeReary
[硕博家园] 考博自荐 +5	科研狗111 2026-05-13	6/300	2026-05-18 11:22 by 糊糊涂涂好
[基金申请] 青C资助名额大幅增加！ +12	西葫芦炒鸡蛋 2026-05-13	16/800	2026-05-18 10:02 by Equinoxhua
[文学芳草园] 半夜喝咖啡 +3	myrtle 2026-05-15	5/250	2026-05-18 01:03 by 小沈2018
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4	xx7gd5zq4e 2026-05-15	6/300	2026-05-17 19:36 by Equinoxhua
[基金申请] 这年头没有找到涵评专家，还有中面上的可能吗 +13	dd921ww 2026-05-12	16/800	2026-05-17 12:38 by 于轩
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5	cjf4bx70cj 2026-05-14	7/350	2026-05-17 06:55 by 11n4dfd8yn
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3	k37jurhrau 2026-05-16	3/150	2026-05-17 01:37 by ue3ir18jc3
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4	k37jurhrau 2026-05-16	4/200	2026-05-17 01:35 by ue3ir18jc3
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3	ky2p12rrjj 2026-05-15	4/200	2026-05-17 00:50 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3	ky2p12rrjj 2026-05-15	3/150	2026-05-17 00:45 by ue3ir18jc3
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4	x0mp7owy2b 2026-05-15	4/200	2026-05-17 00:25 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4	x0mp7owy2b 2026-05-15	4/200	2026-05-16 17:45 by j6b2pdz07o
[硕博家园] 申请博士 +3	呃?呃 2026-05-15	3/150	2026-05-16 11:01 by a4742549
[文学芳草园] 风把牡丹吹跑了 +5	myrtle 2026-05-12	9/450	2026-05-15 15:27 by myrtle
[考博] 26应届毕业生考博求助 +3	wo一定上岸 2026-05-13	3/150	2026-05-14 21:47 by 明海天涯
[考博] 材料类只有一篇综述能申博么 +4	乐逍遥谷 2026-05-13	4/200	2026-05-14 12:05 by zhyzzh