root / Pi / RNG / MWC.c @ 284
Historique | Voir | Annoter | Télécharger (331 octet)
| 1 | 7 | equemene | #include <stdio.h> |
|---|---|---|---|
| 2 | 7 | equemene | #include <math.h> |
| 3 | 7 | equemene | |
| 4 | 7 | equemene | #define znew (z=36969*(z&65535)+(z>>16)) |
| 5 | 7 | equemene | #define wnew (w=18000*(w&65535)+(w>>16)) |
| 6 | 7 | equemene | #define MWC ((znew<<16)+wnew) |
| 7 | 7 | equemene | #define MWCfp (MWC + 1.0f) * 2.328306435454494e-10f |
| 8 | 7 | equemene | |
| 9 | 7 | equemene | int main(void) { |
| 10 | 7 | equemene | int i,z,w;
|
| 11 | 7 | equemene | |
| 12 | 7 | equemene | z=37;
|
| 13 | 7 | equemene | w=91;
|
| 14 | 7 | equemene | |
| 15 | 7 | equemene | for (i=1;i<100;i++) { |
| 16 | 7 | equemene | |
| 17 | 7 | equemene | printf("%i %i %i %i\n",i,znew,wnew,MWC);
|
| 18 | 7 | equemene | |
| 19 | 7 | equemene | } |
| 20 | 7 | equemene | |
| 21 | 7 | equemene | } |