OpenMPが変数を正しく宣言しています

Question

私はOpenMPと並列化を全く新しくしています。この問題は非常にシンプルですが、その原因や解決策を探す方法がわかりません。私のコードは次のようになります。また、FindPrimeを呼び出して、その数が素数であるかどうかをチェックするisPrime主な機能はあり

int FindPrime(int size) { 
    long long next_p = 2; 
    #pragma omp parallel \ 
     shared(next_p,size) 
    long long tid = omp_get_thread_num(); 
    long long tnum = omp_get_num_threads(); 
    for(int long long i = tid; i < size; i=i+tnum) { 
     long long j; 
     j=next_p+1; 
     while (!IsPrime(j)) { 
      j++; 
     } 
     next_p = j; 
    } 
    return next_p; 
} 

。

私はグラムで、このコードをコンパイルしてみてください++とフラグ-fopenmpと-lpthread、私は結果の行に

for(int long long i = tid; i < size; i=i+tnum){ 

を指し

error: 'tid' was not declared in this scope 

を得ることは、それだけで文句と私にも奇妙です約tidであり、tnumではありません。これもOpenMP関数からのものです。

Answer 1

あなたが直面する問題は、あなたが書いたparallel指示の範囲に由来します。あなたはその範囲を定義するブロックを作成するために{と}を使用していないので、この1は、明らかに、すなわち

long long tid = omp_get_thread_num(); 

だからあなたtid変数は今だけ1行の範囲を持って、それを次の一行に制限されています（その宣言）、したがって、ループをforループで使用する場合は宣言しません。問題を解決するには

、あなただけのように、並列ブロックをマークする{と}を追加する必要があります。

int FindPrime(int size) { 
    long long next_p = 2; 
    #pragma omp parallel shared(next_p,size) 
    { // start of the parallel block 
     long long tid = omp_get_thread_num(); 
     long long tnum = omp_get_num_threads(); 
     for(long long i = tid; i < size; i=i+tnum) { 
      long long j; 
      j=next_p+1; 
      while (!IsPrime(j)) { 
       j++; 
      } 
      next_p = j; 
     } 
    } // end of the parallel block 
    return next_p; 
} 

この方法で、あなたのコードは並列化されます。しかし、私はそれが何をすべきかはわかりませんが、私はそれがあなたが期待していることをしないと確信しています...実際には、スレッド間のnext_pの更新に関する競争条件があります。あなたのアルゴリズムを考え直すべきだと思う。

最後に、forループを並列化する方法は、#pragma omp for schedule(static, 1)とまったく同じです。このアプローチに固執したい場合は、おそらくそれを使用することを検討する必要があります。