Question

我有遗留的fortran源文件名为pot.f，我需要将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

您在第二个OpenMP指令中启动第二个parallel部分，该部分未由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

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

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

1 个答案: