| 查看: 466 | 回复: 2 | |||
[交流]
【求助】为什么两种写法得到的值不一样？
|
|
为什么把if语句放到不同的位置，结果就不一样了呢？请大家帮忙分析一下，多谢。 写法一（判断语句为h0(t)/=0）： t0=0 do ntel=1, npoints neworig = .false. if (t0.eq.0) then neworig = .true. elseif (ntel-time0(t0).eq.it0)then neworig = .true. endif if (neworig) then t0 = t0 + 1 ! update number of t=0 time0(t0) = ntel ! store the time of t=0 h0(t0) = h(ntel) endif do t=1, t0 ! update tcf for t=0 delt=ntel-time0(t)+1 ! actual time minus t=0 if(h0(t)/=0.and.(delt.lt.npoints)) then ntime(delt) = ntime(delt)+1d0 tcf(delt)=tcf(delt)+h(ntel)*h0(t) tcf0(delt)=tcf0(delt)+h0(t)*h0(t) endif enddo enddo 写法二（判断语句为if (h(ntel)/=0) then） t0=0 do ntel=1, npoints neworig = .false. if (h(ntel)/=0) then if (t0.eq.0) then neworig = .true. elseif (ntel-time0(t0).eq.it0)then neworig = .true. endif endif if (neworig) then t0 = t0 + 1 ! update number of t=0 time0(t0) = ntel ! store the time of t=0 h0(t0) = h(ntel) endif do t=1, t0 ! update tcf for t=0 delt=ntel-time0(t)+1 ! actual time minus t=0 if（delt.lt.npoints) then ntime(delt) = ntime(delt)+1d0 tcf(delt)=tcf(delt)+h(ntel)*h0(t) tcf0(delt)=tcf0(delt)+h0(t)*h0(t) endif enddo enddo |
回复此楼» 猜你喜欢
情人节自我反思：在爱情中有过遗憾吗？
已经有10人回复
-
今年春晚有几个节目很不错，点赞！
已经有7人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
山东征女友，坐标济南
+1/171
-
大连海事大学轮机学院博士名额1个
+1/170
-
坐标广州，征女友
+2/138
-
天津科技大学海洋与环境学院殷焕顺团队招博士生1名---分析化学领域
+1/82
-
贺电中定位于“积极作用”，是不是对基金委工作不够满意？
+1/73
-
一个陌生女人的来信
+1/60
-
海法大学线上开放日
+1/33
-
国家级青年人才课题组招收2026级硕士研究生
+1/30
-
海南大学海洋技术与装备学院-科研助理招聘（可读博）--膜分离水处理方向
+1/29
-
国家级青年人才课题组招收2026级硕士研究生
+1/28
-
海南大学海洋技术与装备学院-科研助理招聘（可读博）膜分离水处理方向
+1/27
-
清华大学深圳国际研究生院招聘-博士后（长期有效）
+1/27
-
代朋友发 88公务员诚征男友
+1/20
-
2026年天津科技大学“新能源催化与膜材料团队”研究生招生
+1/17
-
英国布里斯托大学诚招博士生，博士后和联合培养生
+1/17
-
英国南安普顿大学禅铎课题组诚招气候动力方向博士后
+1/7
-
湖南大学-分析检测技术和生物柔性传感器-招收1名博士研究生 (2026年，第二批)
+1/6
-
中国科学技术大学 工程科学学院 国家级人才团队 诚聘博士后
+1/4
-
墨尔本大学(QS13)急招CSC博士(补齐全奖)/访问学者/博士后 (材料/生物医学/器官芯片等)
+1/3
-
澳科大招收2026年秋季药物递送/生物材料方向硕士研究生（3月5日18:00报名截止）
+1/2
2楼2011-02-13 10:18:59
3楼2011-02-13 10:24:20