割り当て可能な配列のOpenMP並列ワークシェアレット

Question

Fortran 90の配列に対して要素ごとの計算を行い、openmpでコードを並列化したい。エラーなし

real(8),dimension(:,:,:,:),allocatable :: x,r 
allocate(x(n,n,n,n)) 
allocate(r(n,n,n,n)) 

が、シリアルでプログラム実行することと：私は割付け配列と静的配列xとRを置き換えるために、今したいので、私は入力し

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 

：私は今、次のコードを持っていますコンパイラは "！$ omp parallel workshare"という行を考慮しません。

この場合、どのようなオプションを使用して並列化する必要がありますか？私はループでomp parallel doを試しましたが、はるかに遅いです。

私はWindows上でのgfortran 5.1.0で、私のコードをコンパイルしています：

gfortran -ffree-form test.f -o main.exe -O3 -fopenmp -fno-automatic 

Answer 1

私は前にgfortranでこの問題に遭遇してきました。溶液は、次の形で配列を指定することである：ここ

!$omp parallel workshare 
! do some array calculation 
r(:,:,:,:) = atan(exp(-x)) 
!$omp end parallel workshare 

はreferenceあります。