Révision 79 Pi/C/Hybrid/Pi_Hybrid.c
| Pi_Hybrid.c (revision 79) | ||
|---|---|---|
| 27 | 27 |
#define CONG (jcong=69069*jcong+1234567) |
| 28 | 28 |
#define KISS ((MWC^CONG)+SHR3) |
| 29 | 29 |
|
| 30 |
#define ITERATIONS 1000000000 |
|
| 31 |
|
|
| 30 | 32 |
#define MWCfp MWC * 2.328306435454494e-10f |
| 31 | 33 |
#define KISSfp KISS * 2.328306435454494e-10f |
| 34 |
#define SHR3fp SHR3 * 2.328306435454494e-10f |
|
| 35 |
#define CONGfp CONG * 2.328306435454494e-10f |
|
| 32 | 36 |
|
| 33 |
#define ITERATIONS 1000000000 |
|
| 34 |
|
|
| 35 | 37 |
#define PROCESS 1 |
| 36 | 38 |
|
| 37 | 39 |
#ifdef LONG |
| ... | ... | |
| 55 | 57 |
|
| 56 | 58 |
LENGTH MainLoopGlobal(LENGTH iterations,unsigned int seed_w,unsigned int seed_z) |
| 57 | 59 |
{
|
| 60 |
|
|
| 61 |
#if defined TCONG |
|
| 62 |
unsigned int jcong=seed_z; |
|
| 63 |
#elif defined TSHR3 |
|
| 64 |
unsigned int jsr=seed_w; |
|
| 65 |
#elif defined TMWC |
|
| 58 | 66 |
unsigned int z=seed_z; |
| 59 | 67 |
unsigned int w=seed_w; |
| 60 |
|
|
| 68 |
#elif defined TKISS |
|
| 69 |
unsigned int jcong=seed_z; |
|
| 70 |
unsigned int jsr=seed_w; |
|
| 71 |
unsigned int z=seed_z; |
|
| 72 |
unsigned int w=seed_w; |
|
| 73 |
#endif |
|
| 74 |
|
|
| 61 | 75 |
LENGTH total=0; |
| 62 | 76 |
|
| 63 | 77 |
for (LENGTH i=0;i<iterations;i++) {
|
| 64 | 78 |
|
| 65 |
float x=MWCfp ; |
|
| 66 |
float y=MWCfp ; |
|
| 79 |
#if defined TINT32 |
|
| 80 |
#define THEONE 1073741824 |
|
| 81 |
#if defined TCONG |
|
| 82 |
unsigned int x=CONG>>17 ; |
|
| 83 |
unsigned int y=CONG>>17 ; |
|
| 84 |
#elif defined TSHR3 |
|
| 85 |
unsigned int x=SHR3>>17 ; |
|
| 86 |
unsigned int y=SHR3>>17 ; |
|
| 87 |
#elif defined TMWC |
|
| 88 |
unsigned int x=MWC>>17 ; |
|
| 89 |
unsigned int y=MWC>>17 ; |
|
| 90 |
#elif defined TKISS |
|
| 91 |
unsigned int x=KISS>>17 ; |
|
| 92 |
unsigned int y=KISS>>17 ; |
|
| 93 |
#endif |
|
| 94 |
#elif defined TINT64 |
|
| 95 |
#define THEONE 4611686018427387904 |
|
| 96 |
#if defined TCONG |
|
| 97 |
unsigned long x=(unsigned long)(CONG>>1) ; |
|
| 98 |
unsigned long y=(unsigned long)(CONG>>1) ; |
|
| 99 |
#elif defined TSHR3 |
|
| 100 |
unsigned long x=(unsigned long)(SHR3>>1) ; |
|
| 101 |
unsigned long y=(unsigned long)(SHR3>>1) ; |
|
| 102 |
#elif defined TMWC |
|
| 103 |
unsigned long x=(unsigned long)(MWC>>1) ; |
|
| 104 |
unsigned long y=(unsigned long)(MWC>>1) ; |
|
| 105 |
#elif defined TKISS |
|
| 106 |
unsigned long x=(unsigned long)(KISS>>1) ; |
|
| 107 |
unsigned long y=(unsigned long)(KISS>>1) ; |
|
| 108 |
#endif |
|
| 109 |
#elif defined TFP32 |
|
| 110 |
#define THEONE 1.0f |
|
| 111 |
#if defined TCONG |
|
| 112 |
float x=CONGfp ; |
|
| 113 |
float y=CONGfp ; |
|
| 114 |
#elif defined TSHR3 |
|
| 115 |
float x=SHR3fp ; |
|
| 116 |
float y=SHR3fp ; |
|
| 117 |
#elif defined TMWC |
|
| 118 |
float x=MWCfp ; |
|
| 119 |
float y=MWCfp ; |
|
| 120 |
#elif defined TKISS |
|
| 121 |
float x=KISSfp ; |
|
| 122 |
float y=KISSfp ; |
|
| 123 |
#endif |
|
| 124 |
#elif defined TFP64 |
|
| 125 |
#define THEONE 1.0f |
|
| 126 |
#if defined TCONG |
|
| 127 |
double x=(double)CONGfp ; |
|
| 128 |
double y=(double)CONGfp ; |
|
| 129 |
#elif defined TSHR3 |
|
| 130 |
double x=(double)SHR3fp ; |
|
| 131 |
double y=(double)SHR3fp ; |
|
| 132 |
#elif defined TMWC |
|
| 133 |
double x=(double)MWCfp ; |
|
| 134 |
double y=(double)MWCfp ; |
|
| 135 |
#elif defined TKISS |
|
| 136 |
double x=(double)KISSfp ; |
|
| 137 |
double y=(double)KISSfp ; |
|
| 138 |
#endif |
|
| 139 |
#endif |
|
| 67 | 140 |
|
| 68 | 141 |
// Matching test |
| 69 |
int inside=((x*x+y*y) < 1.0f) ? 1:0;
|
|
| 142 |
unsigned long inside=((x*x+y*y) < THEONE) ? 1:0;
|
|
| 70 | 143 |
total+=inside; |
| 71 | 144 |
} |
| 72 | 145 |
|
Formats disponibles : Unified diff