是否可以在程序的多个文件中使用 OpenMP 指令？

Question

我有一个由多个 .c 文件和多个 .h 文件组成的 C 程序。我希望在主函数中有一个 #pragma omp parallel 指令（以便所有线程仅创建一次），然后执行其他 OpenMP 操作，例如 #pragma omp for 在其他文件中。

但是，我似乎无法做到这一点。编译主文件时，它抱怨 #pragma omp parallel 的 private() 和 shared() 位中提到的一些变量指令不存在（在该文件中它们不存在 - 因为它们在另一个文件中），并且在编译另一个文件时它抱怨我有一个 #pragma omp for 没有封闭的#pragma omp parallel。

代码在文件之间很好地分割，我不想将它们全部放回到一个文件中。有什么办法解决这个问题吗？

Answer 1

你绝对可以：

outer.c:

#include <stdio.h>
#include <stdlib.h>
#include <omp.h>

void inner(int n);

int main(int argc, char **argv) {
    int n=15;

    #pragma omp parallel shared(n) num_threads(4)
    {
        inner(n);
    }

    return 0;
}

inner.c:

#include <omp.h>
#include <stdio.h>

void inner(int n) {
    int thread = omp_get_thread_num();
    printf("%3d: got %d\n", thread, n);

    #pragma omp for
    for  (int i=0;i<n;i++) {
        int newthread=omp_get_thread_num();
        printf("%3d: doing iter %d.\n",newthread,i);
    }
}

和running:

$ make
gcc -fopenmp -std=c99   -c -o outer.o outer.c
gcc -fopenmp -std=c99   -c -o inner.o inner.c
gcc -o nested outer.o inner.o -fopenmp -std=c99 -lgomp    
$ ./nested 
  3: got 15
  3: doing iter 12.
  3: doing iter 13.
  3: doing iter 14.
  0: got 15
  0: doing iter 0.
  0: doing iter 1.
  0: doing iter 2.
  0: doing iter 3.
  1: got 15
  1: doing iter 4.
  1: doing iter 5.
  1: doing iter 6.
  1: doing iter 7.
  2: got 15
  2: doing iter 8.
  2: doing iter 9.
  2: doing iter 10.
  2: doing iter 11.

但是不行，你不能在一个例程中设置另一个例程中变量的共享属性——它们只是不在范围内。您无法设置它们的共享，就像您无法设置它们的值一样。

一旦你启动了（比如说）inner，那里的一切都是私人的；您必须将任何共享的内容作为共享传递。

只是为了澄清有关“所有内容都是私有的”的一点：上面的内容与 thread、i 和 newthread 没有任何不同，

    int n=15;

    #pragma omp parallel shared(n) num_threads(4)
    {
        int thread = omp_get_thread_num();
        printf("%3d: got %d\n", thread, n);

        #pragma omp for
        for  (int i=0;i<n;i++) {
            int newthread=omp_get_thread_num();
            printf("%3d: doing iter %d.\n",newthread,i);
        }
    }

因为 thread、i 和 newthread 是在并行块内定义的 - 无论是否在函数内部 - 它们都必须是私有的。

Answer 2

您是否尝试过“extern XY”使其他文件中的变量已知？

难道不能对声明相关变量的所有文件使用通用标头吗？

它们是局部变量吗（那么你无论如何也不能这样做）？