Question

对于以下每个代码段，请使用OpenMP pragma使循环并行，或解释为什么代码段不适合并行执行。

a. for (i = 0; i < sqrt(x); i++)
 a[i] = 2.3 * i;
 if (i < 10)
 b[i] = a[i];
}

b. flag = 0;
 for (i = 0; i < n && !flag; i++)
 a[i] = 2.3 * i;
 if (a[i] < b[i])
 flag = 1;
}

c. for (i = 0; i < n && !flag; i++)
 a[i] = foo(i);
d. for (i = 0; i < n && !flag; i++) {
 a[i] = foo(i);
 if (a[i] < b[i])
 a[i] = b[i];
}

e. for (i = 0; i < n && !flag; i++) {
 a[i] = foo(i);
 if (a[i] < b[i])
 break;
}

f. dotp = 0;
for (i = 0; i < n; i++)
 dotp += a[i] * b[i];
g. for (i = k; i < 2 * k; i++)
 a[i] = a[i] + a[i – k];
h. for (i = k; i < n; i++) {
 a[i] = c * a[i – k];

对于上述问题的任何帮助都会非常受欢迎。无论如何......

Answer 1

我不会做你的硬件，但我会给出一个提示。在使用OpenMp for循环时，您应该对变量的范围保持警惕。例如：

#pragma omp parallel for
for(int x=0; x < width; x++)
{
    for(int y=0; y < height; y++)
    {
        finalImage[x][y] = RenderPixel(x,y, &sceneData);
    }
}

没问题，因为x和y是私有变量。

怎么样？

int x,y;
#pragma omp parallel for
for(x=0; x < width; x++)
{
    for(y=0; y < height; y++)
    {
        finalImage[x][y] = RenderPixel(x,y, &sceneData);
    }
}

在这里，我们在for循环之外定义了x和y。现在考虑y。每个线程都会在没有任何同步的情况下访问/写入它，因此会发生数据争用，这很可能导致逻辑错误。

阅读更多here并祝你好运。

OPenMP：何时可以并行化循环

1 个答案: