OpenMPを使用してGPUにチームを配布するには？

Question

私はNvidia Geforce GT 740MをOpenMPとclang-3.8コンパイラを使った並列プログラミングに利用しようとしています。

CPUで並列処理すると、結果が得られます。しかし、GPUで処理すると、私の結果はほとんど乱数になります。

したがって、私はスレッドチームを正しく配布していないと思っていましたし、データ競合があるかもしれません。私はfor-loopsを別々にやらなければならないと思うが、間違いがどこにあるのか分からない。

#include <stdio.h> 
#include <stdlib.h> 
#include <omp.h> 


int main(int argc, char* argv[]) 
    { 
    const int n =100; float a = 3.0f; float b = 2.0f; 
    float *x = (float *) malloc(n * sizeof(float)); 
    float *y = (float *) malloc(n * sizeof(float)); 

    int i; 
    int j; 
    int k; 

    double start; 
    double end; 
    start = omp_get_wtime(); 


    for (k=0; k<n; k++){ 
     x[k] = 2.0f; 
     y[k] = 3.0f; 
    } 


    #pragma omp target data map(to:x[0:n]) map(tofrom:y[0:n]) map(to:i) map(to:j) 
    { 

     #pragma omp target teams 
     #pragma omp distribute 
     for(i = 0; i < n; i++) { 

      #pragma omp parallel for 
      for (j = 0; j < n; j++){ 
       y[j] = a*x[j] + y[j]; 
     } 
    } 


} 


end = omp_get_wtime(); 

printf("Work took %f seconds.\n", end - start); 

free(x); free(y); 

return 0; 
} 

私のGPUのアーキテクチャには何かがあるかもしれません。そこで私はこれを追加している：あなたの助けに感謝をトピックにかなり新しい

イム、そう:)

Answer 1

はい、ここでのレースがあります。異なるチームが配列 'y'の同じ要素を読み書きしています。おそらく、あなたはこれのようなものが欲しいですか？

for(i = 0; i < n; i++) { 
    #pragma omp target teams distribute parallel for 
    for (j = 0; j < n; j++){ 
    y[j] = a*x[j] + y[j]; 
    } 
}