使用基於編譯指示原子的約簡來逼近 PI

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 atomic
  area += (4.0 / (1.0 + x*x));
}
pi = h * area;

在此示例中,每個執行緒執行迭代計數的子集,並且它們以原子方式累積到共享變數區域中,這確保不會丟失更新。我們可以在這裡使用#pragma atomic,因為給定的操作(+ =)可以原子方式完成,與使用#pragma omp critical 相比,這簡化了可讀性。