PI 手工制作 pragma omp 减少的近似

int i;
int n = 1000000;
double area = 0;
double h = 1.0 / n;

#pragma omp parallel shared(n, h)
{
  double thread_area = 0;                      // Private / local variable

  #pragma omp for
  for (i = 1; i <= n; i++)
  {
    double x = h * (i - 0.5);
    thread_area += (4.0 / (1.0 + x*x));
  }

  #pragma omp atomic                       // Applies the reduction manually
  area += thread_area;                     // All threads aggregate into area
}
double pi = h * area;

线程在#pragma omp parallel 中生成。每个线程都有一个独立/私有的 thread_area，用于存储其部分添加内容。以下循环使用#pragma omp for 在线程之间分配。在这个循环中，每个线程计算自己的 thread_area，在此循环之后，代码按原子顺序聚合该区域