OpenMP 2の共有ベクトル

Question

私の前の質問 Shared vectors in OpenMP 異なるスレッドが の異なるスレッドがベクトルの異なる要素にアクセスする限り、共有ベクトルで異なるスレッドを読み書きできるようにすることができました。 異なるスレッドが、次のようなベクトルのすべての（たまに同じである）要素を読み取る必要がある場合はどうなりますか？一度に一つのスレッドがベクトルの任意の要素を読むことができるか、そうでないので

#include <vector> 

int main(){ 

vector<double> numbers; 
vector<double> results(10); 
double x; 

//write 25 values in vector numbers 
for (int i =0; i<25; i++){ 
    numbers.push_back(cos(i)); 
} 

#pragma omp parallel for default(none) \ 
shared(numbers, results) \ 
private(x) 
    for (int j = 0; j < 10; j++){ 
     for(int k = 0; k < 25; k++){ 
      x += 2 * numbers[j] * numbers[k] + 5 * numbers[j * k/25]; 
     } 
     results[j] = x;  
    } 

    return 0; 

} 

この並列化は遅いでしょうか？私はこの問題をfirstprivate(numbers)という節で解決できますか？

すべてのスレッドが独自のベクトルを得るようにベクトルの配列を作成することは意味がありますか？例えば

：

vector<double> numbersx[**-number of threads-**]; 

Answer 1

読取素子を共有するために多くの仕事がないこと複数のスレッドからのベクトルは問題ではありません。あなたのコードには同期がないので、同時にアクセスされます。

ベクトルのサイズが大きければ、キャッシュの問題もありませんが、ベクトルが大きければキャッシュアクセスパターンによって速度が低下することがあります。その場合、数字データのコピーを別にすると、パフォーマンスが向上します。

Answer 2

より良いアプローチ：それはちょっとによる事実に遅くなります

#include <vector> 

int main(){ 

vector<double> numbers; 
vector<double> results(10); 

//write 25 values in vector numbers 
for (int i =0; i<25; i++){ 
    numbers.push_back(cos(i)); 
} 

#pragma omp parallel for 
    for (int j = 0; j < 10; j++){ 
     double x = 0; // make x local var 
     for(int k = 0; k < 25; k++){ 
      x += 2 * numbers[j] * numbers[k] + 5 * numbers[j * k/25]; 
     } 
     results[j] = x; // no race here  
    } 

    return 0; 

} 

同じの