添加Openmp指令时出现奇怪的gfortran编译错误

Question

我有一个名为pot.f的旧版fortran源文件，我需要将OpenMP应用于并行，如下所示，但是我会出现有关意外结束状态等错误消息。但是当我$OMP通过添加其他注释掉行时！在第一列中，没有错误。

对我来说真的很奇怪。谁能告诉我出了什么问题？

subroutine pot_osc(rvp,R_pot,e_pot,pe_pot,ftmp,gtmp,vtmp,natoms)
   implicit none
   include 'sizes.h'
   include 'constants.h'
   include 'omp_lib.h'

   double precision  ftmp(maxatoms,3),gtmp(3),R_pot(maxatoms,3)

   !$OMP PARALLEL WORKSHARE  SHARED(gtmp,ftmp)
   !$OMP PARALLEL NUM_THREADS(16)
      gtmp = 0d0
      ftmp = 0d0
   !$OMP END PARALLEL WORKSHARE 
   return
end

subroutine pot_asym(rvp,vtmp)
   implicit none
   include 'constants.h'  
   return
end

错误讯息：

 end
   1

Error: Unexpected END statement at (1)


  subroutine pot_asym(rvp,vtmp)
  1
Error: Unclassifiable statement at (1)

Answer 1

您parallel在第二个OpenMP指令中开始第二部分，该指令不以终止end parallel。因此，OpenMP指令应显示为

!$OMP PARALLEL WORKSHARE SHARED(gtmp,ftmp) NUM_THREADS(16)
       gtmp = 0d0
       ftmp = 0d0
!$OMP END PARALLEL WORKSHARE

或者如果您想保留换行符

!$OMP PARALLEL WORKSHARE SHARED(gtmp,ftmp) &
!$OMP NUM_THREADS(16)
       gtmp = 0d0
       ftmp = 0d0
!$OMP END PARALLEL WORKSHARE

过去，正是这种初始化遇到了一些问题。似乎使用gfortran主线程进行编译时可以完成所有工作。更糟糕的是，借助“初次接触原理”，整个阵列位于与第一个线程关联的内存中。在我们的CCNUMA机器上，这导致了巨大的速度下降。

为了解决这个问题，我使用了显式循环进行初始化：

!$OMP PARALLEL DO SHARED(gtmp,ftmp) NUM_THREADS(16)
       do i=1,maxatoms
         ftmp(i,:) = 0d0
       enddo
!$OMP END PARALLEL DO
!      No need to do three elements in parallel
       gtmp = 0d0

我不知道他们是否解决了这个问题，但是从那时起，我对共享内存中的数组使用这种初始化方式。