c - openmp for 循环

Question

我有以下代码:

for(i=0; i<num1; i++)
{
 j=i%num2;
 a[j]=do_compute(j);
}

假设 num2 << num1。如果我想使用omp for，是否需要使用j作为私有(private)？即如下:

#pragma omp parallel private(i,j)
{
 #pragma omp for
 for(i=0; i<num1; i++)
 {
  j=i%num2;
  a[j]=do_compute(j);
 }
}

或者这样就足够了:

#pragma omp parallel
{
 #pragma omp for
 for(i=0; i<num1; i++)
 {
  j=i%num2;
  a[j]=do_compute(j);
 }
}

谢谢。

Answer 1

如果 j 仅在循环范围内使用，我会执行以下操作:

[...]
unsigned int i;
#pragma omp parallel for
for(i = 0; i < num1; ++i)
{
    unsigned int j = i % num2;
    a[j] = do_compute(j);
}
[...]

不想在这里挑剔，但考虑到 do_compute 是确定性的并且没有副作用，您最好只使用

[...]
unsigned int i;
#pragma omp parallel for
for(i = 0; i < num2; ++i)
    a[i] = do_compute(i);
[...]