如何在 OpenMP 中找到多个 2D arrays 的最大值

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);
        }
    }