C - 并行化循环 omp

Question

我有一个问题:我必须使用 OMP 并行化这段代码。

存在数据依赖问题，不知道如何解决。有什么建议么？

for (n = 2; n < N+1; n++) {
  dz = *(dynamic_d + n-1)*z;
  *(dynamic_A + n) = *(dynamic_A + n-1) + dz * (*(dynamic_A + n-2));
  *(dynamic_B + n) = *(dynamic_B + n-1) + dz * (*(dynamic_B + n-2));
}

Answer 1

由于依赖关系，您无法并行化循环迭代，但您可以使用部分拆分 dynamic_A 与 dynamic_B 的计算:

#pragma omp parallel sections
{
    #pragma omp section
    {
        // NOTE: Declare n and dz locally so that it is private!
        for (int n = 2; n < N+1; n++) {
            my_type dz = dynamic_d[n-1] * z;
            dynamic_A[n] = dynamic_A[n-1] + dz * dynamic_A[n-2];
        }
    }
    #pragma omp section
    {
        for (int n = 2; n < N+1; n++) {
            my_type dz = dynamic_d[n-1] * z;
            dynamic_B[n] = dynamic_B[n-1] + dz * dynamic_B[n-2];
        }
    }
}

请使用数组索引，而不是邪恶的指针算术引用异常。