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

回复此楼

» 猜你喜欢

博士去二本高校当辅导员，如何调整心态？已经有8人回复
27细胞生物学申博已经有5人回复
面上有专家说实验设备不是我们单位的已经有5人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有20人回复
江苏省自然基金什么时候出结果已经有6人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有11人回复
闲聊已经有3人回复
如何查询2026年获批项目名称？已经有3人回复
两块石头已经有10人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有15人回复

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

查看全部散金贴

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