openmp Approximation de PI en utilisant des réductions basées sur #pragma omp critical


Exemple

double area;
double h = 1.0 / n;
#pragma omp parallel for shared(n, h, area) 
for (i = 1; i <= n; i++)
{
  double x = h * (i - 0.5);
  #pragma omp critical
  {
    area += (4.0 / (1.0 + x*x));
  }
}
double pi = h * area;

Dans cet exemple, chaque thread exécute un sous-ensemble du nombre d'itérations et il s'accumule de manière atomique dans la zone de la variable partagée, ce qui garantit qu'il n'y a pas de mises à jour perdues.