高难度的三维以上数值如何快速存储问题
本帖最后由 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零柒玖玖（企鹅）
计算量太小，讨论效率没有意义。程序启动时间占比太大。 这个代码就是错的，根本不能运行，那你怎么对比的时间呢。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 li913 发表于 2021-2-4 15:51
计算量太小，讨论效率没有意义。程序启动时间占比太大。
已经更改了。 li913 发表于 2021-2-4 16:00
这个代码就是错的，根本不能运行，那你怎么对比的时间呢。MODULE mod
REAL,DIME ...
我已经重新更改了 其一，多重循环在优化编译后，一般会变成单重循环的机器语言。建义关闭优化参数，重新编译，看看有无区别。
其二，尽量多用Fortran的Array功能， 如上面给出例子中的数组赋值，或者数组Where语句等。少用循环。 风平老涡 发表于 2021-2-6 06:23
其一，多重循环在优化编译后，一般会变成单重循环的机器语言。建义关闭优化参数，重新编译，看看有无区别。 ...
已经尝试过，没用
页:
[1]