Fortran Coder

标题: 大数组但不是每次调用，设成可分配还是存为common计算更快 [打印本页]
作者: 石子    时间: 2023-7-17 19:02
标题: 大数组但不是每次调用，设成可分配还是存为common计算更快
本帖最后由 石子 于 2023-7-17 19:07 编辑

各位好，我有若干个二维数组，每维的大小大约为10000~20000之间，采用omp并行计算。
计算的目标数组为A与B，计算的总时间步大约为十万左右。为了更新A和B，每隔50步要调用subroutine1，关系大概如下：
[Fortran] 纯文本查看 复制代码


do t=1,Nt

...

if(mod(t,50)==0) then

 call subroutine1(A,B,c1,c2,w1,w2,n1,n2)

end if

...

end do


subroutine1中的关系为：A,B属性为in, c1,c2属性为out，其余的w1,w2和n1,n2为inout，即读入上一次的值并进行更新。
目前我的处理办法是subroutine中所有的数组都存在公共数据块COMMON中。我想W1,W2和N1,N2可不可以在subroutine1中每次分配，然后计算，然后释放，这样会省很多内存，或许会快一些？但是这样有个问题就是，上一步算的没法存下来。

想请问各位，这个问题有没有什么更合适的处理方式？目标是让计算更高效、速度。感谢

作者: 风平老涡    时间: 2023-7-17 23:33
本帖最后由 风平老涡 于 2023-7-17 23:39 编辑

你的一个数组单精度最小0.4G，最大1.6G，双精度加倍。8个数组就是3.2G~12G。如果W和N数组每次重新分配，计算速度降低并且无法保存前一次的结果。能够改进的是算法，可以看看数组的对称性及稀疏性。
作者: 石子    时间: 2023-7-18 10:58
本帖最后由 石子 于 2023-7-18 11:20 编辑
风平老涡 发表于 2023-7-17 23:33
你的一个数组单精度最小0.4G，最大1.6G，双精度加倍。8个数组就是3.2G~12G。如果W和N数组每次重新分配，计 ...

好的，谢谢呀。



欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2