如何在 OpenMP 中找到多个二维数组的最大值

Question

我想并行化这段代码：

        #pragma omp parallel for 
        for (i=0; i<=imax+1; i++) { // combined loops
            for (j=1; j<=jmax+1; j++) {
                umax = max(fabs(u[i][j]), umax);
                vmax = max(fabs(v[j][i]), vmax);
            }
        }

但是 umax 和 vmax 存在竞争条件。有没有办法在不使用 #pragma omp critical?

序列化代码的情况下解决这个问题？

Answer 1

是的，您只需要使用 OpenMP 减少条款。

    #pragma omp parallel for reduction(max:umax, vmax)
    for (i=0; i<=imax+1; i++) { // combined loops
        for (j=1; j<=jmax+1; j++) {
            umax = max(fabs(u[i][j]), umax);
            vmax = max(fabs(v[j][i]), vmax);
        }
    }