| 查看: 671 | 回复: 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 |
回复此楼» 猜你喜欢
郑州大学申博
已经有6人回复
-
40+中年人失业后，还能找到好工作吗？
已经有18人回复
-
人工智能虚构史料、伪造文献频现，多家期刊声明抵制
已经有5人回复
-
CsPbI2Br全无机钙钛矿SEM图像异常分析
已经有8人回复
-
申请科研项目纯靠实力能成功吗？
已经有12人回复
2楼2011-02-13 10:18:59
3楼2011-02-13 10:24:20