openmp - 如何使用 OpenMP 在中间的所有线程的同步（`barrier`）中创建`omp parallel for`

Question

我有两个功能，do_step_one(i)和do_step_two(i)，i从0到N-1。

目前，我有这个（顺序）代码：

for(unsigned int i=0; i<N; i++) {
     do_step_one(i);
}

for(unsigned int i=0; i<N; i++) {
     do_step_two(i);
}

do_step_one()和的每次调用do_step2()都可以按任何顺序并行执行，但任何都do_step_two()需要结束所有的do_step_one()开始（它使用do_step_one()结果）。

我尝试了以下方法：

#omp parallel for
for(unsigned int i=0; i<N; i++) {
    do_step_one(i);

#omp barrier

    do_step_two(i);
}

但是 gcc 抱怨

convolve_slices.c:21：警告：屏障区域可能不会紧密嵌套在工作共享、关键、有序、主或显式任务区域内。

我有什么误解？如何解决这个问题？

Answer 1

顺便说一句，如果要确保不重新创建线程，请将并行声明和 for 声明分开：

#pragma omp parallel
{
  #pragma omp for
  for(unsigned int i=0; i<N; i++){
    do_step_one(i);
  }
  //implicit barrier here
  #pragma omp for
  for(unsigned int i=0; i<N; i++){
    do_step_two(i);
  }
}

Answer 2

我在这段代码中看到的一个问题是代码不符合规范:)

如果您需要结束所有的 do_step_one()，则需要以下内容：

#pragma omp parallel for
for(unsigned int i=0; i<N; i++){
     do_step_one(i);
}

#pragma omp parallel for
for(unsigned int i=0; i<N; i++){
     do_step_two(i);
}

这样做的结果将是第一个 for 的并行性，然后是第二个 for 的并行性。