我应该使用哪个 openMP 编译指示来阻塞 for 循环？

Question

我的算法（求解泊松方程）是完全可并行的——前提是所有线程在每次迭代结束时同步。

Function f, fNext;
init(f);
#pragma omp parallel
for(int step=0; step<maxITER; step++) {
#pragma omp for
   for(int i=0; i<N; i++) {
      for(int j=0; j<N; j++) {
         fNext(i,j) = someOperator( f(i,j) );
      }
   }
   f = fNext;
}//Threads must synchronize here

#pragma omp for 在继续下一次迭代之前是否确保线程同步？

Answer 1

是的。来自 OpenMP 规范（例如，v 3.1，但这已在从一开始），在“工作共享结构：”下

循环结构末尾有一个隐式屏障，除非
指定了 nowait 子句。

也就是说，在 for 循环结束时，除非您执行类似 #pragma omp for nowait 的操作，否则会存在隐含的屏障，因此没有线程会执行 f=fNext直到所有线程都完成 for 循环。