OpenMP ungleiche Last ohne Schleife

Question

I einen OpenMP Code, die

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

    } 
} 

daher wie folgen

sieht, ist die Idee, dass der parallele Abschnitt durch die verfügbaren Threads ausgeführt wird, solange der Zähler unter einem gewissen ist Wert. Abhängig vom Szenario (ich mache MC-Sachen hier, also ist es zufällig), könnten die Berechnungen länger dauern als andere, so dass es hier ein Ungleichgewicht zwischen den Arbeitern gibt, was sich aufgrund der impliziten Barriere am Ende der Parallele zeigt Sektion.

Es scheint, als hätte #pragma omp parallel für Möglichkeiten, dies zu umgehen (d. H. Nowait Direktive/dynamische Planung), aber ich kann dies nicht verwenden, da ich keine obere Iterationsnummer für die for-Schleife kenne.

Irgendwelche Ideen/Entwurfsmuster, wie man mit solch einer Situation umgeht?

Mit freundlichen Grüßen!

Answer 1

Führen Sie alles in einem einzigen parallelen Abschnitt und greifen Sie auf die counter atomar.

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++; 
    } 
} 

Sie können nicht direkt in der while-Bedingung, wegen der atomaren Zugriff überprüfen. Beachten Sie, dass dieser Code überschwingt, d. H. counter > MAX ist nach dem parallelen Abschnitt möglich. Denken Sie daran, dass counter freigegeben ist und von vielen Threads gelesen/aktualisiert wird.