OpenMP double for loop

Question

openMPを使用してマルチスレッドを適用したいと思います。 ここに書いた簡単なコードです。

vector<Vector3f> a; 
int i, j; 
for (i = 0; i<10; i++) 
{ 
    Vector3f b; 
#pragma omp parallel for private(j) 

    for (j = 0; j < 3; j++) 
    { 
     b[j] = j; 
    } 
    a.push_back(b); 
} 

for (i = 0; i < 10; i++) 
{ 
    cout << a[i] << endl; 
} 

私は作品のリクにそれを変更したい：の#pragma行が削除された場合に

parallel for1 
{ 
    for2 
} 

または

for1 
{ 
    parallel for2 
} 

コードは動作します。私はそれを使用すると動作しません。どうしたの？

/////////追加

は、実際に私は ダブルループ質問のためのより複雑な例にOpenMPのを使用しています。

ここにも、私はMPを適用しないとき、それはうまく動作します。 しかし、私はそれを適用すると、 ベクトルpush_back行でエラーが発生します。

vector<Class> B; 
for 1 
{ 
    #pragma omp parallel for private(j) 
    parallel for j 
    { 
     Class A; 
     B.push_back(A); // error!!!!!!! 
    } 
} 

B.push_back（A）行を消去すると、MPを適用すると同様に動作します。

正確なエラーメッセージが見つかりませんでしたが、ベクトルに関する例外エラーのようです。デバッグは、スレッドセーフではありません

void _Reallocate(size_type _Count) 
    { // move to array of exactly _Count elements 
    pointer _Ptr = this->_Getal().allocate(_Count); 

    _TRY_BEGIN 
    _Umove(this->_Myfirst, this->_Mylast, _Ptr); 

Answer 1

std::vector::push_backで停止し、あなたが複数のスレッドからの競合状態に対する任意の保護なしでそれを呼び出すことはできません。

代わりに、サイズがすでに正しいようにベクトルを準備し、operator[]で要素を挿入します。

また、重要な領域に挿入を保護することができます。一度に一つのスレッド

#pragma omp critical 
B.push_back(A); 

この方法では、エラーを修正しますが、コードを遅くする挿入を行います。

一般的には、並列化には適切な方法でアプローチしないと思いますが、より明確でより代表的な問題記述がない限り、より良いアドバイスを与える方法はありません。