openmp OpenMP reductions Approximation of PI hand-crafting the #pragma omp reduction
Example
h = 1.0 / n;
#pragma omp parallel private(x) shared(n, h)
{
double thread_area = 0; // Private / local variable
#pragma omp for
for (i = 1; i <= n; i++)
{
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
}
pi = h * area;
The threads are spawned in the #pragma omp parallel. Each thread will have an independent/private thread_area that stores its partial addition. The following loop is distributed among threads using #pragma omp for. In this loop, each thread calculates its own thread_area and after this loop, the code sequentially aggregates the area atomically through #pragma omp atomic.
