使用基于编译指示原子的约简来逼近 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 相比，这简化了可读性。