使用 pragma omp 減少子句逼近 PI

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

在此示例中，每個執行緒執行迭代計數的子集。每個執行緒都有其本地私有區域副本，並且在並行區域的末尾它們都應用加法運算（+），以便為區域生成最終值。