并行 OpenMP 缩减与函数定义？

Posted 2023-02-14

【中文标题】并行 OpenMP 缩减与函数定义？

【英文标题】：Parallel OpenMP reduction vs. function definition?

【发布时间】：2013-05-14 15:32:05

【问题描述】：

我正在使用 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。 有什么解决办法吗？

【参考方案1】：

您可以将sum=sum+pixel[0] 放在#pragma omp atomic 或#pragma omp critical 行下，而不是减少。

另一种选择是在 omp 部分之前有一个 double local_sum = sum;，在 local_sum 上减少，然后在 for 循环之后有 sum = local_sum;。

【讨论】：

使用第二个选项（local_sum）可能会比使用“omp atomic”或“omp critical”执行得更好，因为减少子句让每个线程都有自己的local_sum本地版本。使用“comp critical”强制线程排队等待更新 sum，使用“omp atomic”使线程竞争相同的内存位置。