并行OpenMP缩减vs.函数定义

Parallel OpenMP reduction vs. function definition?

本文关键字:函数 定义 vs 缩减 OpenMP 并行      更新时间:2023-10-16

我正在使用OpenMP,但问题是我正在声明/定义一个函数如下:

void compute_image(double pixel[nb], double &sum)
{
        #pragma omp parallel for reduction(+:sum)
    for (int j=0;j<640;j++)
    {
    if ...
    sum=sum+pixel[0];
    ....
    }
....
}

我现在意识到的是:

Error   2   error C3030: 'sum' : variable in 'reduction' clause/directive cannot have reference type    C:Users...test.cpp    930

事实上,我无法摆脱OpenMP。有解决方案吗?

您可以将sum=sum+pixel[0]放在#pragma omp atomic#pragma omp critical行下,而不是简化。

另一个选项可能是在omp部分之前有一个double local_sum = sum;,在local_sum上减少,然后在for循环之后有sum = local_sum;