OpenMP并行FORTRAN嵌套循环问题

Transpose · 发表于 2023-4-26 11:17:21

本帖最后由 Transpose 于 2023-4-26 11:23 编辑

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

t1(nm,km(1)*km(2)*km(3),na,na)

这个数组占的内存大概是180G，应该是内存不足了
100*100*1*75*177*177=23496750000
23496750000/1e9=23.5
23.5*8=188GB

Transpose · 发表于 2023-4-26 16:42:22

感觉只能修改算法，减少内存的占用了，这么大的数组直接定义，不太合理
另外注意数组的列优先顺序，使用ifort的时候，极易产生临时数组，这也会导致内存不足

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

t1(i,j,:,:) = th*(i-j+ic)

比如这一句，可能会产生临时数组

你的前一个程序，大量使用matmul,在ifort里面也是不提倡的，ifort的matmul优化不好，而且大量的产生临时数组

Transpose · 发表于 2023-4-27 09:04:19

这种数组传递方式需要写接口

http://fcode.cn/guide-103-1.html