OpenMP Programm ist langsamer als sequenzielle ein

Question

Wenn ich den folgenden Code versuchen

double start = omp_get_wtime(); 

long i; 

#pragma omp parallel for 
    for (i = 0; i <= 1000000000; i++) { 
     double x = rand(); 
    } 

    double end = omp_get_wtime(); 

    printf("%f\n", end - start); 

Ausführungszeit ca. 168 Sekunden, während die sequentielle Version verbringt nur 20 Sekunden.

Ich bin immer noch ein Neuling in der parallelen Programmierung. Wie könnte ich eine parallele Version bekommen, die schneller ist als die sequentielle Version?

Answer 1

Der Zufallszahlengenerator rand(3) verwendet globale Statusvariablen (versteckt in der (g) libc-Implementierung). Der Zugriff auf sie aus mehreren Threads führt zu Cache-Problemen und ist auch nicht threadsicher. Sie sollten den rand_r(3) Anruf mit seed Parameter privat auf das Gewinde verwenden:

long i; 
unsigned seed; 

#pragma omp parallel private(seed) 
{ 
    // Initialise the random number generator with different seed in each thread 
    // The following constants are chosen arbitrarily... use something more sensible 
    seed = 25234 + 17*omp_get_thread_num(); 
    #pragma omp for 
    for (i = 0; i <= 1000000000; i++) { 
     double x = rand_r(&seed); 
    } 
} 

Beachten Sie, dass diese unterschiedlichen Strom von Zufallszahlen, wenn parallel ausgeführt produzieren wird, als wenn in Serien ausgeführt. Ich würde auch erand48(3) als eine bessere (Pseudo-) Zufallszahlenquelle empfehlen.