Fortran Coder

标题: 高难度的三维以上数值如何快速存储问题 [打印本页]

作者: sam295040799 时间: 2021-2-4 15:39
标题: 高难度的三维以上数值如何快速存储问题
本帖最后由 sam295040799 于 2021-2-6 09:55 编辑

最近发现的一个问题，一直困扰我，看看没有大神给我解惑，有偿。

   MODULE CDAS1

      REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: indt,ce
      real dzz

      ENDMODULE CDAS1

   PROGRAM cs

   USE CDAS1

   IMPLICIT NONE

      integer i,j,k,n,maxx,maxy,maxz
      real time1,time2

      dzz=0.6
      maxz=80;maxy=80;maxx=80
      allocate(indt(maxx,maxy,maxz),ce(maxx,maxy,maxz))
      indt=0;ce=0

         do k=1,maxz
         do j=1,maxy
         do i=1,maxx
         if(i<40) indt(i,j,k)=1
         enddo
         enddo
         enddo

   call CPU_TIME(time1)

   do n=1,50000
      do k=1,maxz
      do j=1,maxy
      do i=1,maxx
      if( indt(i,j,k)==1) cycle
      ce(i,j,k)=dzz
      enddo
      enddo
      enddo
   enddo
   call CPU_TIME(time2)
   print*,'time',time2-time1

   ENDPROGRAM cs

这段程序在release下执行的，出现下面情况
这段程序里 if( indt(i,j,k)==1) cycle是对所有的i，j，k进行筛选，不满足条件的就进行ce(i,j,k)=dzz，但是结果发现没有if( indt(i,j,k)==1) cycle，ce(i,j,k)=dzz程序执行次数是80*80*80次，而有if( indt(i,j,k)==1) cycle，ce(i,j,k)=dzz程序执行次数是40*80*80次，两者执行次数相差一倍，但是耗时很接近。
关于这个问题的解释就是，fortran循环是列循环，如果循环没有中断，循环速度很快，但是如果是跳跃循环，那这个循环速度就会变慢，再加上if(indt(i,j,k)==1)cycle，也耗时。但是没有找到好的办法让这个循环里带if（）cycle 的代码执行次数和时间成比例。如果在debug执行，情况又不一样了

联系方式：二9五零4零柒玖玖（企鹅）

作者: li913 时间: 2021-2-4 15:51
计算量太小，讨论效率没有意义。程序启动时间占比太大。

作者: li913 时间: 2021-2-4 16:00
这个代码就是错的，根本不能运行，那你怎么对比的时间呢。

[Fortran] 纯文本查看 复制代码

MODULE mod

REAL,DIMENSION(:,:,:),ALLOCATABLE :: indt,ce

ENDMODULE mod



PROGRAM ceshi

USE mod

IMPLICIT NONE

integer maxx,maxy,maxz

real dzz

dzz=0.6

maxz=80;maxy=80;maxx=80

allocate(indt(maxx,maxy,maxz),ce(maxx,maxy,maxz))

indt=0;ce=0

indt(:39,:,:)=1

ce(:39,:,:)=dzz



ENDPROGRAM ceshi

作者: sam295040799 时间: 2021-2-4 17:35

li913 发表于 2021-2-4 15:51
计算量太小，讨论效率没有意义。程序启动时间占比太大。

已经更改了。

作者: sam295040799 时间: 2021-2-4 17:36

li913 发表于 2021-2-4 16:00
这个代码就是错的，根本不能运行，那你怎么对比的时间呢。[mw_shl_code=fortran,true]MODULE mod
REAL,DIME ...

我已经重新更改了

作者: 风平老涡 时间: 2021-2-6 06:23
其一，多重循环在优化编译后，一般会变成单重循环的机器语言。建义关闭优化参数，重新编译，看看有无区别。
其二，尽量多用Fortran的Array功能，如上面给出例子中的数组赋值，或者数组Where语句等。少用循环。

作者: sam295040799 时间: 2021-2-6 09:54

风平老涡发表于 2021-2-6 06:23
其一，多重循环在优化编译后，一般会变成单重循环的机器语言。建义关闭优化参数，重新编译，看看有无区别。 ...

已经尝试过，没用

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)