並行OpenMP縮減與功能定義？

Question

我正在使用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。 有什么辦法嗎？

Answer 1

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

另一種選擇是讓double local_sum = sum; 在omp節之前，減少local_sum，然后sum = local_sum; 在for循環之后。