OpenMPの共有ベクター

Question

私は使用しているプログラムをparallizeしようとしていますが、次の質問があります。 複数のスレッドが同じベクトルでもベクトルの異なる要素を読み書きする必要がある場合、パフォーマンスが低下しますか？私の気持ちは、私のプログラムがそれを視覚化することでほとんど遅くならないという理由があります。次のコードを取る：

#include <vector> 

int main(){ 

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

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

#pragma omp parallel for \ 
    private(x) \ 
    shared(numbers, results) 
     for(int j = 0; j < 10; j++){ 

      x = 2 * numbers[j] + 5; 
#pragma omp critical // do I need this ? 
      { 
       results[j] = x;  
      } 
     } 

    return 0; 

} 

明らかに実際のプログラムでは、はるかに高価な操作を行いますが、この例では唯一の私の質問を説明し ばなりません。それで、forループは高速かつ完全に並列処理できますか、あるいは、異なるスレッドはお互いに待ち合わせる必要がありますか？ベクトルの異なる要素をすべて読み込んでいますが、一度に1つのスレッドしかベクトル番号にアクセスできないためです。

書き込み操作で同じ質問：すべてのスレッドがベクトル結果の異なる要素に書き込むため、クリティカルなプラグマが必要ですか、問題ありませんか？ 私は得ることができるあらゆる援助に満足していますし、これを行う良い方法があるかどうかを知ることも良いでしょう（ベクトルをまったく使用しないかもしれませんが、単純な配列やポインタなど？） 私はベクトルaren場合によってはスレッドセーフで、ポインタを使用することをお勧めします：OpenMP and STL vector

ありがとうございました！

Answer 1

複数のスレッド内のベクトルの問題の大部分は、サイズ変更が必要な場合、ベクターの内容全体をメモリ内の新しい場所（より大きな割り当てられたチャンク）にコピーすると思いますこれに並列にアクセスすると、削除されたオブジェクトを読み込もうとしました。不足については

あなたは、配列のサイズを変更していない場合は、私が同時読み取りとのトラブルがあったことがなかったしているベクターに書き込む（明らかに限り、私は二度同じ要素を書いていないよと）

（実際にそのクリティカルセクションの外で行われる量に応じて）

クリティカルセクションステートメントを削除することができます上記を念頭に置いて）。

Answer 2

同じ要素が同時に変更されることはないため、重要なセクティノのためにスピードアップする必要はありません。クリティカルセクションの部分を削除しても問題ありません。

メモリアクセスが線形でない場合（これは例のとおりです）、スレッドはキャッシュ（同じキャッシュラインの要素を書き込む）と戦う可能性があるため、スケジュール戦略でもプレイできます。あなたの場合のように要素の数が与えられてループに分岐がない場合（したがって、それらはほぼ同じ速度で実行されます）、デフォルトのIIRCであるstaticが最適です。

（ところで、あなたがprivate(x)を避けるために、ループ内xを宣言することができますし、sharedディレクティブがIIRC（私はそれを使用したことがない）暗黙的に指定されます。）