Question

我遇到一个问题，我的以下代码使用openmp运行速度较慢：

chunk = nx/nthreads;
int i, j;
for(int t = 0; t < n; t++){
     #pragma omp parallel for default(shared) private(i, j) schedule(static,chunk) 
     for(i = 1; i < nx/2+1; i++){
        for(j = 1; j < nx-1; j++){
            T_c[i][j] =0.25*(T_p[i-1][j] +T_p[i+1][j]+T_p[i][j-1]+T_p[i][j+1]);
            T_c[nx-i+1][j] = T_c[i][j];
        }
    }
    copyT(T_p, T_c, nx);
}
print2file(T_c, nx, file);

问题是当我运行多个线程时，计算时间会更长。

Answer 1

首先，在外循环的每次迭代中重新启动并行区域，从而增加巨大开销。

其次，一半的线程只是坐在那里什么都不做，因为你的块大小是它应该的两倍 - 它是nx/nthreads而并行循环的迭代次数是{{1}因此总共有nx/2个块。除了你试图实现的是复制(nx/2)/(nx/nthreads) = nthreads/2的行为。

schedule(static)

如果您修改#pragma omp parallel for (int t = 0; t < n; t++) { #pragma omp for schedule(static) for (int i = 1; i < nx/2+1; i++) { for (int j = 1; j < nx-1; j++) { T_c[i][j] = 0.25*(T_p[i-1][j]+T_p[i+1][j]+T_p[i][j-1]+T_p[i][j+1]); T_c[nx-i-1][j] = T_c[i][j]; } } #pragma omp single copyT(T_p, T_c, nx); } print2file(T_c, nx, file);也使用并行copyT，则应删除for构造。您不需要single，因为这是默认设置。您不要声明并行循环default(shared)的循环变量 - 即使此变量来自外部作用域（因此在区域中隐式共享），OpenMP也会自动将其设置为私有。只需在循环控件中声明所有循环变量，它就会自动运行并应用默认的共享规则。

第二个半月，你的内循环中可能存在（可能）错误。第二个分配声明应为：

private

（或T_c[nx-i-1][j] = T_c[i][j];如果你没有在下方留下光环），否则T_c[nx-i][j]等于i，那么你将访问1之外的T_c[nx][...] T_c的界限。

第三，一般性提示：不是将一个数组复制到另一个数组，而是使用指向这些数组的指针，并在每次迭代结束时交换两个指针。

Answer 2

我发现至少有三个问题可能导致您发布的代码段性能不佳：

块大小太小，无法在线程间分配任何增益。
在循环内打开和关闭parallel区域可能会影响性能。
两个最里面的循环看起来是独立的，你只能并行化其中一个（失去了利用更广泛的迭代空间的可能性）。

您可以在下面找到我将对代码进行的一些修改：

// Moving the omp parallel you open/close the parallel 
// region only one time, not n times
#pragma omp parallel default(shared)
for(int t = 0; t < n; t++){
     // With collapse you parallelize over an iteration space that is 
     // composed of (nx/2+1)*(nx-1) elements not only (nx/2+1)
     #pragma omp for collapse(2) schedule(static)
     for(int i = 1; i < nx/2+1; i++){
        for(int j = 1; j < nx-1; j++){
            T_c[i][j] =0.25*(T_p[i-1][j] +T_p[i+1][j]+T_p[i][j-1]+T_p[i][j+1]);
            T_c[nx-i+1][j] = T_c[i][j];
        }
    }
    // As the iteration space is very small and the work done 
    // at each iteration is not much, static schedule will likely be the best option
    // as it is the one that adds the least overhead for scheduling
    copyT(T_p, T_c, nx);
}
print2file(T_c, nx, file);

openmp慢了多个线程，想不通

2 个答案: