24小时热门版块排行榜

>论坛更新日志 (3202)
>虫友互识 (703)
>休闲灌水 (109)
>导师招生 (104)
>博后之家 (92)
>基金申请 (92)
>论文投稿 (87)
>文献求助 (77)
>教师之家 (47)
>硕博家园 (43)
>考博 (25)
>药学 (21)
>招聘信息布告栏 (15)
>医学 (11)
>公派出国 (11)
>考研 (11)

返回列表

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

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

回复此楼

» 猜你喜欢

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

查看全部散金贴

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