关于大数组openmp并行的问题

千手斑 · 发表于 2019-4-7 22:04:39

向各位大神求助，我现在想要并行做一些矩阵的计算，发现只要矩阵稍微一大就无法运行了，不知道有人遇到过这种情况吗，该如何解决呢？这个代码很简单，串行完全没问题，而并行的时候，矩阵规模很小时没有问题，规模一大比如现在的2000就无法运行了。
   program main
   use omp_lib
   integer::i,j
   integer,parameter::nec0=2000
   common/matr/ an(nec0,nec0)

   write(6,*) 'hello word!'
   !$OMP PARALLEL PRIVATE(i,j)
   !$OMP DO
   do i=1,nec0
      do j=1,nec0
         an(i,j)=float(i+j)
      enddo
   enddo
   !$OMP END DO
   !$OMP END PARALLEL

   open(3,file='an.txt',status='unknown')
   do i=1,nec0
      write(3,*) (an(i,j),j=1,nec0)
   enddo
   close(3)

   stop
   end

li913 · 发表于 2019-4-12 10:58:58

不要用小于10的通道号；屏幕输出用星号，不用6.

千手斑 · 发表于 2019-4-11 20:18:47

li913 发表于 2019-4-9 17:59
我这里运行没问题。你把整个项目打包上传。

这个在vs里运行没问题，在simply fortran里运行会出现上述问题，是哪里需要设置吗？

li913 · 发表于 2019-4-9 17:59:40

我这里运行没问题。你把整个项目打包上传。

千手斑 · 发表于 2019-4-8 12:43:37

li913 发表于 2019-4-8 12:06
给错误提示，对于找错来说，这是最重要的，比你的代码更重要。
猜测可能是堆栈问题。 ...

会弹出对话框，test.exe已停止工作

li913 · 发表于 2019-4-8 12:06:30

给错误提示，对于找错来说，这是最重要的，比你的代码更重要。
猜测可能是堆栈问题。

千手斑 · 发表于 2019-4-7 22:34:11

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

01program main
02use omp_lib
03integer::i,j
04integer,parameter::nec0=2000
05common/matr/ an(nec0,nec0)
06 
07write(6,*) 'hello word!'
08!$OMP PARALLEL PRIVATE(i,j)
09!$OMP DO
10do i=1,nec0
11  do j=1,nec0
12      an(i,j)=float(i+j)
13  enddo
14enddo
15!$OMP END DO
16!$OMP END PARALLEL
17 
18open(3,file='an.txt',status='unknown')
19do i=1,nec0
20  write(3,*) (an(i,j),j=1,nec0)
21enddo
22close(3)
23 
24stop
25end

[并行] 关于大数组openmp并行的问题