如何正确并行化嵌套的 for 循环

Question

我正在使用 OpenMP 来并行化嵌套的标量 for 循环：

double P[N][N];
double x=0.0,y=0.0;

for (int i=0; i<N; i++)
{
    for (int j=0; j<N; j++)
    {
        P[i][j]=someLongFunction(x,y);
        y+=1;
    }
    x+=1;
}

在这个循环中，重要的是矩阵 P 在标量和并行版本中必须相同：

我所有可能的尝试都没有成功...

Answer 1

这里的问题是您添加了迭代到迭代的依赖项：

x+=1;
y+=1;

因此，就目前的代码而言，它是不可并行的。尝试这样做会导致错误的结果。 （你可能已经看到了）

幸运的是，在您的情况下，您可以直接计算它们而无需引入此依赖项：

for (int i=0; i<N; i++)
{
    for (int j=0; j<N; j++)
    {
        P[i][j]=someLongFunction((double)i, (double)N*i + j);
    }
}

现在您可以尝试在此上面添加一个 OpenMP 编译指示，看看它是否有效：

#pragma omp parallel for
for (int i=0; i<N; i++)
{
    for (int j=0; j<N; j++)
    {
        P[i][j]=someLongFunction((double)i, (double)N*i + j);
    }
}