| 1 |
1 |
//
|
| 2 |
2 |
// Estimation of Pi using Monte Carlo exploration process
|
|
3 |
// Cecill v2 Emmanuel QUEMENER <emmanuel.quemener@gmail.com>
|
| 3 |
4 |
// gcc -std=c99 -O3 -o Pi Pi.c -lm
|
| 4 |
5 |
//
|
| 5 |
6 |
|
| ... | ... | |
| 19 |
20 |
|
| 20 |
21 |
#define MWCfp MWC * 2.328306435454494e-10f
|
| 21 |
22 |
#define KISSfp KISS * 2.328306435454494e-10f
|
|
23 |
#define SHR3fp SHR3 * 2.328306435454494e-10f
|
|
24 |
#define CONGfp CONG * 2.328306435454494e-10f
|
| 22 |
25 |
|
| 23 |
26 |
#define ITERATIONS 1000000000
|
| 24 |
27 |
|
| ... | ... | |
| 32 |
35 |
|
| 33 |
36 |
LENGTH MainLoopGlobal(LENGTH iterations,unsigned int seed_w,unsigned int seed_z)
|
| 34 |
37 |
{
|
|
38 |
#if defined TCONG
|
|
39 |
unsigned int jcong=seed_z;
|
|
40 |
#elif defined TSHR3
|
|
41 |
unsigned int jsr=seed_w;
|
|
42 |
#elif defined TMWC
|
| 35 |
43 |
unsigned int z=seed_z;
|
| 36 |
44 |
unsigned int w=seed_w;
|
|
45 |
#elif defined TKISS
|
|
46 |
unsigned int jcong=seed_z;
|
|
47 |
unsigned int jsr=seed_w;
|
|
48 |
unsigned int z=seed_z;
|
|
49 |
unsigned int w=seed_w;
|
|
50 |
#endif
|
| 37 |
51 |
|
| 38 |
52 |
LENGTH total=0;
|
| 39 |
|
|
|
53 |
|
| 40 |
54 |
for (LENGTH i=0;i<iterations;i++) {
|
| 41 |
55 |
|
| 42 |
|
float x=MWCfp ;
|
| 43 |
|
float y=MWCfp ;
|
|
56 |
#if defined TINT32
|
|
57 |
#define THEONE 1073741824
|
|
58 |
#if defined TCONG
|
|
59 |
unsigned int x=CONG>>17 ;
|
|
60 |
unsigned int y=CONG>>17 ;
|
|
61 |
#elif defined TSHR3
|
|
62 |
unsigned int x=SHR3>>17 ;
|
|
63 |
unsigned int y=SHR3>>17 ;
|
|
64 |
#elif defined TMWC
|
|
65 |
unsigned int x=MWC>>17 ;
|
|
66 |
unsigned int y=MWC>>17 ;
|
|
67 |
#elif defined TKISS
|
|
68 |
unsigned int x=KISS>>17 ;
|
|
69 |
unsigned int y=KISS>>17 ;
|
|
70 |
#endif
|
|
71 |
#elif defined TINT64
|
|
72 |
#define THEONE 4611686018427387904
|
|
73 |
#if defined TCONG
|
|
74 |
unsigned long x=(unsigned long)(CONG>>1) ;
|
|
75 |
unsigned long y=(unsigned long)(CONG>>1) ;
|
|
76 |
#elif defined TSHR3
|
|
77 |
unsigned long x=(unsigned long)(SHR3>>1) ;
|
|
78 |
unsigned long y=(unsigned long)(SHR3>>1) ;
|
|
79 |
#elif defined TMWC
|
|
80 |
unsigned long x=(unsigned long)(MWC>>1) ;
|
|
81 |
unsigned long y=(unsigned long)(MWC>>1) ;
|
|
82 |
#elif defined TKISS
|
|
83 |
unsigned long x=(unsigned long)(KISS>>1) ;
|
|
84 |
unsigned long y=(unsigned long)(KISS>>1) ;
|
|
85 |
#endif
|
|
86 |
#elif defined TFP32
|
|
87 |
#define THEONE 1.0f
|
|
88 |
#if defined TCONG
|
|
89 |
float x=CONGfp ;
|
|
90 |
float y=CONGfp ;
|
|
91 |
#elif defined TSHR3
|
|
92 |
float x=SHR3fp ;
|
|
93 |
float y=SHR3fp ;
|
|
94 |
#elif defined TMWC
|
|
95 |
float x=MWCfp ;
|
|
96 |
float y=MWCfp ;
|
|
97 |
#elif defined TKISS
|
|
98 |
float x=KISSfp ;
|
|
99 |
float y=KISSfp ;
|
|
100 |
#endif
|
|
101 |
#elif defined TFP64
|
|
102 |
#define THEONE 1.0f
|
|
103 |
#if defined TCONG
|
|
104 |
double x=(double)CONGfp ;
|
|
105 |
double y=(double)CONGfp ;
|
|
106 |
#elif defined TSHR3
|
|
107 |
double x=(double)SHR3fp ;
|
|
108 |
double y=(double)SHR3fp ;
|
|
109 |
#elif defined TMWC
|
|
110 |
double x=(double)MWCfp ;
|
|
111 |
double y=(double)MWCfp ;
|
|
112 |
#elif defined TKISS
|
|
113 |
double x=(double)KISSfp ;
|
|
114 |
double y=(double)KISSfp ;
|
|
115 |
#endif
|
|
116 |
#endif
|
| 44 |
117 |
|
| 45 |
|
// Matching test
|
| 46 |
|
int inside=((x*x+y*y) < 1.0f) ? 1:0;
|
|
118 |
unsigned long inside=((x*x+y*y) < THEONE) ? 1:0;
|
| 47 |
119 |
total+=inside;
|
| 48 |
120 |
}
|
| 49 |
121 |
|
| ... | ... | |
| 53 |
125 |
LENGTH splitter(LENGTH iterations,unsigned int process,
|
| 54 |
126 |
unsigned int seed_w,unsigned int seed_z)
|
| 55 |
127 |
{
|
| 56 |
|
LENGTH inside[1024],insides=0;
|
|
128 |
// Trashy static allocation...
|
|
129 |
LENGTH inside[4096],insides=0;
|
| 57 |
130 |
|
| 58 |
131 |
#pragma omp parallel for
|
| 59 |
132 |
for (int i=0 ; i<process; i++) {
|