24小时热门版块排行榜

>论坛更新日志 (4562)
>休闲灌水 (1108)
>导师招生 (1023)
>虫友互识 (438)
>文献求助 (368)
>博后之家 (274)
>招聘信息布告栏 (199)
>论文投稿 (160)
>考博 (129)
>硕博家园 (119)
>找工作 (108)
>第一性原理 (93)
>生物科学 (38)
>绿色求助(高悬赏) (37)
>基金申请 (36)
>考研 (32)

返回列表

田山东

捐助贵宾 (著名写手)

应助: 5 (幼儿园)
金币: 77
帖子: 1577
在线: 138.1小时
虫号: 800307

[交流] 这样居然也能产生随机数

本来有个程序，要重复调用一个数组，但是数组对角元为-0.5---0.5之间的随机数，没想到下面的语句也可以实现：
N=20
Gj=0.8
CALL RANDOM_SEED()
!!第一次生成随机数
   CALL RANDOM_NUMBER(Gj)
do j=1,N
   H0(j,j)=(Gj(j)-0.5)*0.2*t
   enddo
CALL SLRD(N,ENG,H0,HL,HR,SL,ETA,EPS)
CALL RANDOM_NUMBER(Gj)
  !!第二次生成随机数
   do j=1,N
   H0(j,j)=(Gj(j)-0.5)*0.2*t
   enddo

CALL SLRD(N,ENG,H0,HR,HL,SR,ETA,EPS)
!!第三次生成随机数
CALL RANDOM_NUMBER(Gj)
   do j=1,N
   H0(j,j)=(Gj(j)-0.5)*0.2*t
   enddo
经过验证，三次生成的随机数不同。而且，数次运行程序，每次生成的随机数也是互不相同的。看来这样很合理了。
发个帖，作个记号吧。

回复此楼

» 猜你喜欢

面上没中，邀请各位路过的虫友分析一下分数已经有15人回复
博士去二本高校当辅导员，如何调整心态？已经有10人回复
面上有专家说实验设备不是我们单位的已经有6人回复
27细胞生物学申博已经有5人回复
江苏省自然基金什么时候出结果已经有6人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有11人回复
闲聊已经有3人回复
两块石头已经有10人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有15人回复
初秋的晨风已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2012-12-25 15:29:57

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

615780236

禁虫 (知名作家)

★
田山东(金币+1): 谢谢参与

本帖内容被屏蔽

2楼2012-12-25 21:43:23

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

田山东

捐助贵宾 (著名写手)

应助: 5 (幼儿园)
金币: 77
帖子: 1577
在线: 138.1小时
虫号: 800307

引用回帖:

2楼: Originally posted by 615780236 at 2012-12-25 21:43:23
恭喜楼主
N=20
Gj=0.8
CALL RANDOM_SEED()
!!第一次生成随机数
   CALL RANDOM_NUMBER(Gj)
      do j=1,N
   H0(j,j)=(Gj(j)-0.5)*0.2*t
   enddo
      CALL SLRD(N,ENG,H0,HL,HR,SL,ET ...

谢谢参与。

回复此楼

3楼2012-12-26 10:59:55

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

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

返回列表