在混合MPI/OpenMP中进行MPI调用的线程

我在复制的混合MPI/OpenMP代码中发现了一个问题以下面引用的代码中最简单的形式。我正在使用2个线程每个MPI等级。这两个线程随后被用于OpenMP"Section"要进行多次计算，其中之一是对两个不同的向量a和B进行"mpi_allreduce"调用，这两个向量的结果存储在W和WW中。问题是每次我运行程序时我最终得到了不同的输出。我认为MPI调用重叠并且缩减的阵列W和WW即使在它们有不同的名字，但我不确定。关于如何克服的任何评论这个问题是受欢迎的。

详细信息：MPI线程级别在代码中初始化为MPI_THRAD_MULTIPLE但我也尝试过串行和漏斗（同样的问题）。

我编译了mpiifort-openmp allreduce_omp_mpi.f90和运行我使用的：

导出OMP_NUM_THREADS=2mpirun-np3/a.out

      PROGRAM HELLO
      use mpi 
      use omp_lib
      IMPLICIT NONE
      INTEGER nthreads, tid 
      Integer Provided,mpi_err,myid,nproc
      CHARACTER(MPI_MAX_PROCESSOR_NAME):: hostname
      INTEGER :: nhostchars
      integer :: i
      real*8 :: A(1000), B(1000), W(1000),WW(1000)
      provided=0
      !Initialize MPI context
      call mpi_init_thread(MPI_THREAD_MULTIPLE,provided,mpi_err)
      CALL mpi_comm_rank(mpi_comm_world,myid,mpi_err)
      CALL mpi_comm_size(mpi_comm_world,nproc,mpi_err)
      CALL mpi_get_processor_name(hostname,nhostchars,mpi_err)
      !Initialize arrays
      A=1.0
      B=2.0
      !Check if MPI_THREAD_MULTIPLE is available
      if (provided >= MPI_THREAD_MULTIPLE) then
      write(6,*) ' mpi_thread_multiple provided',myid
      else
      write(6,*) ' not mpi_thread_multiple provided',myid
      endif
!$OMP PARALLEL PRIVATE(nthreads, tid) NUM_THREADS(2)
!$omp sections
!$omp section
       call mpi_allreduce(A,W,1000,mpi_double_precision,mpi_sum,mpi_comm_world,mpi_err)
!$omp section
       call mpi_allreduce(B,WW,1000,mpi_double_precision,mpi_sum,mpi_comm_world,mpi_err)
!$omp end sections
!$OMP END PARALLEL
       write(6,*) 'W',(w(i),i=1,10)
       write(6,*) 'WW',(ww(i),i=1,10)
      CALL mpi_finalize(mpi_err)
      END