openmp并行化内部循环

Question

我有一个 Fortran 代码。代码包含两个循环。我想并行化 ONLY 内部循环。这是我的代码：

!$OMP PARALLEL PRIVATE(i,j)
do i=1, N 
  !$OMP PARALLEL DO
  do j=1, M
    ! do some calculations
  end do
  !$OMP END PARALLEL DO
end do
!$OMP END PARALLEL

并行化是否正确？我不确定，是否需要将!$OMP PARALLEL PRIVATE(i,j)放在开头？我应该省略它并在第二个循环之前只声明 PRIVATE(i) 吗？无论如何，我有点困惑，请解释什么是正确的行为。

Answer 1

为什么i 是私有的？不是所有线程都需要吗？它在整个内部循环期间不得更改（因为它是外部循环的循环计数器）。如果它被声明为私有，除非使用firstprivate，否则它在并行部分的开头是未定义的。

Fortran 中的 OpenMP 注意到 j 是并行循环的循环计数器，因此它是隐式私有的。所以没有必要明确声明它。

接下来你应该避免嵌套的 OpenMP 部分（即第二个!$OMP PARALLEL）

因为喜欢直白，所以写成

!$OMP PARALLEL PRIVATE(j) SHARED(i)
do i=1, N 
  !$OMP DO
  do j=1, M
    ! do some calculations
  end do
  !$OMP END DO
end do
!$OMP END PARALLEL