並行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行下。

另一個選項可能是在omp部分之前有一個double local_sum = sum;，在local_sum上減少，然後在for循環之後有sum = local_sum;。