openmp critique
Suite à cette question , pour le code ci-dessous (à partir de l'exemple de documentation MS OpenMP )
// omp_critical.cpp
// compile with: /openmp
#include <omp.h>
#include <stdio.h>
#include <stdlib.h>
#define SIZE 10
int main()
{
int i;
int max;
int a[SIZE];
for (i = 0; i < SIZE; i++)
{
a[i] = rand();
printf_s("%d\n", a[i]);
}
max = a[0];
#pragma omp parallel for num_threads(4)
for (i = 1; i < SIZE; i++)
{
if (a[i] > max)
{
#pragma omp critical
{
// compare a[i] and max again because max
// could have been changed by another thread after
// the comparison outside the critical section
if (a[i] > max)
max = a[i];
}
}
}
printf_s("max = %d\n", max);
}
Puis-je retirer l'extérieur si tester et faire
max = a[0];
#pragma omp parallel for num_threads(4)
for (i = 1; i < SIZE; i++)
{
#pragma omp critical
{
// compare a[i] and max again because max
// could have been changed by another thread after
// the comparison outside the critical section
if (a[i] > max)
max = a[i];
}
}
Réponses
Vous pouvez, mais cela aboutit effectivement à une exécution séquentielle. Les threads attendent constamment d'entrer dans la section critique de sorte qu'un seul thread exécute le corps de la boucle à la fois. Par conséquent, vous obtenez les mêmes performances (peut-être même pire en raison de la surcharge de synchronisation) qu'une boucle série simple.
L'exemple de la documentation MS se synchronise uniquement si une nouvelle valeur maximale a été rencontrée. Cela permet de traiter toutes les valeurs inférieures jusqu'à ce point en parallèle.
Comme suggéré dans les commentaires, utilisez une construction de réduction.