OpenMP C++ - 如何並行化這個函數？

Question

我想這個並行化功能，但我是新開放的MP，我會很感激，如果有人可以幫助我：

void my_function(float** A,int nbNeurons,int nbOutput, float* p, float* amp){ 
    float t=0; 
    for(int r=0;r<nbNeurons;r++){ 
     t+=p[r]; 
    } 

    for(int i=0;i<nbOutput;i++){ 
     float coef=0; 
     for(int r=0;r<nbNeurons;r++){ 
     coef+=p[r]*A[r][i]; 
     } 
    amp[i]=coef/t; 
    } 
} 

我不知道如何正確地進行並行化因雙循環，目前，我只想到做一個： #pragma omp parallel for reduction(+:t)

但我認爲這不是通過openMp更快地獲得計算的最佳方式。

在此先感謝，

Answer 1

首先：我們需要了解上下文。你的分析器在哪裏告訴你最花時間？

一般來說，粗粒度平行化效果最好，@Alex說：平行於外部for循環。

void my_function(float** A,int nbNeurons,int nbOutput, float* p, float* amp) 
{ 
    float t=0; 
    for(int r=0;r<nbNeurons;r++) 
     t+=p[r]; 

#pragma parallel omp for 
    for(int i=0;i<nbOutput;i++){ 
     float coef=0; 
     for(int r=0;r<nbNeurons;r++){ 
      coef+=p[r]*A[r][i]; 
     } 
     amp[i]=coef/t; 
    } 
} 

根據實際體積，它可能是有趣的計算在後臺t和移動劃分出並行循環：

void my_function(float** A,int nbNeurons,int nbOutput, float* p, float* amp) 
{ 
    float t=0; 
#pragma omp parallel shared(amp) 
    { 
#pragma omp single nowait // only a single thread executes this 
     { 
      for(int r=0;r<nbNeurons;r++) 
       t+=p[r]; 
     } 

#pragma omp for 
     for(int i=0;i<nbOutput;i++){ 
      float coef=0; 
      for(int r=0;r<nbNeurons;r++){ 
       coef+=p[r]*A[r][i]; 
      } 
      amp[i]=coef; 
     } 

#pragma omp barrier 
#pragma omp master // only a single thread executes this 
     { 
      for(int i=0; i<nbOutput; i++){ 
       amp[i] /= t; 
      } 
     } 
    } 
} 

注意未經測試的代碼。 OMP有時候會有棘手的語義，所以我可能錯過了那裏的'共享'聲明。儘管如此，探測器不會很快通知您。