我想在Fortran 90中对数组进行元素计算，同时将我的代码与openmp并行化 . 我现在有以下代码：
program test
implicit none
integer,parameter :: n=50
integer :: i
integer(8) :: t1,t2,freq
real(8) :: seq(n),r(n,n,n,n)
real(8),dimension(n,n,n,n) :: x
call system_clock(COUNT_RATE=freq)
seq=[(i,i=1,n)]
x=spread(spread(spread(seq,2,n),3,n),4,n)
call system_clock(t1)
!$omp parallel workshare
! do some array calculation
r=atan(exp(-x))
!$omp end parallel workshare
call system_clock(t2)
print*, sum(r)
print '(f6.3)',(t2-t1)/real(freq)
end program test
我现在想用可分配的数组替换静态数组x和r，所以我输入：
real(8),dimension(:,:,:,:),allocatable :: x,r
allocate(x(n,n,n,n))
allocate(r(n,n,n,n))
但程序串行运行没有错误，编译器没有考虑“！$ omp parallel workshare”这一行 .
在这种情况下，我应该使用哪些选项进行并行化？我已尝试使用带循环的 omp parallel do ，但速度要慢得多 .
我在windows上用gfortran 5.1.0编译我的代码：
gfortran -ffree-form test.f -o main.exe -O3 -fopenmp -fno-automatic
1 回答
我之前在
gfortran遇到过这个问题 . 解决方案是以下列形式指定数组：这是reference .