mirror of
https://github.com/foxsen/archbase.git
synced 2026-02-04 10:54:05 +08:00
18 lines
371 B
C
18 lines
371 B
C
#include <stdio.h>
|
|
#include <omp.h>
|
|
int main(){
|
|
int i;
|
|
int num_steps=1000000;
|
|
double x,pi,step,sum=0.0;
|
|
step = 1.0/(double) num_steps;
|
|
#pragma omp parallel for private(i, x), reduction(+:sum)
|
|
for(i=0;i<num_steps;i++)
|
|
{
|
|
x=(i+0.5)*step;
|
|
sum = sum+4.0/(1.0+x*x);
|
|
}
|
|
pi = step*sum;
|
|
printf("pi %1f\n", pi);
|
|
return 0;
|
|
}
|