24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

astringent

铜虫 (著名写手)

应助: 16 (小学生)
金币: 2210.7
帖子: 1180
在线: 354.8小时
虫号: 900839

[交流] 【求助】为什么两种写法得到的值不一样？

为什么把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

回复此楼

» 猜你喜欢

招博士已经有4人回复
救命帖已经有7人回复
限项规定已经有6人回复
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向已经有3人回复
英文综述是否需要润色及查重已经有5人回复
为什么nbs上溴没有产物点出现呢已经有9人回复
最失望的一年已经有18人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-02-13 06:57:02

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

snoopyzhao

至尊木虫 (职业作家)

★ ★ ★
astringent(金币+10): 2011-02-13 20:36:39
余泽成(金币+3): 谢谢参与应助！ 2011-02-20 18:22:58

你这个程序片断有太多的东西没有交待，所以不能从程序片断上来解读

在下面的这段代码中，你加了一个判断才能进行 neworig 的赋值，但在方法一中，没有这样的判断，因此，可能会导致 neworig 的值在两个程序中不一样，那么可能会直接影响最后的结果……

CODE:

             if (h(ntel)/=0) then

               if (t0.eq.0) then

                   neworig = .true.

               elseif (ntel-time0(t0).eq.it0)then

                   neworig = .true.

               endif

              endif