OpenMPの不均等な負荷ループ

Question

用せず、私は、次の

したがって

while(counter < MAX) { 
    #pragma omp parallel reduction(+:counter) 
    { 
    // do monte carlo stuff 
    // if a certain condition is met, counter is incremented 

    } 
} 

ように見えるOpenMPのコードを持って、アイデアは、平行部があれば、カウンタが一定以下であるように使用可能なスレッドによって実行されることをあります値。シナリオによっては（私はここでMCのことをするのでランダムなので）、計算は他のものより長くかかるかもしれないので、並列の最後の暗黙的な障壁のために明白になる作業者の間に不均衡がありますセクション。

#pragma omp parallelのように見えるのは、これを回避する方法（nowaitディレクティブ/ダイナミックスケジューリング）があるかもしれませんが、forループの上位反復番号がわからないため、これを使用できません。

どのようなアイデアやデザインパターンそのような状況に対処するには？

よろしくお願いいたします。

Answer 1

すべてを1つの並列セクションで実行し、counterにアトミックにアクセスします。

int counter = 0; 
#pragma omp parallel 
while(1) { 
    int local_counter; 
    #pragma omp atomic read 
    local_counter = counter; 
    if (local_counter >= MAX) { 
      break; 
    } 
    // do monte carlo stuff 
    // if a certain condition is met, counter is incremented 
    if (certain_condition) { 
     #pragma omp atomic update 
     counter++; 
    } 
} 

原子アクセスのため、while条件で直接チェックすることはできません。 このコードはオーバーシュートします。つまり、パラレルセクションの後にcounter > MAXが可能です。 counterは多くのスレッドによって共有され、読み書きされることに注意してください。