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，在此迴圈之後，程式碼按原子順序聚合該區域