OpenMP параллельная работа для выделенного массива

Question

Я хочу сделать некоторые элементарные вычисления на массивах в 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 

Я хочу сейчас, чтобы заменить статические массивы х и р с размещаемыми массивами, так что я тип:

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

но что запуск программы в серийном без ошибок и компилятор не учитывает строку «! $ omp parallel workhare».

Какие варианты следует использовать для распараллеливания в этом случае? Я пробовал с omp parallel do с петлями, но он намного медленнее.

Я компиляции моего кода с 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.