#include "WELL512a.h" #include #define MASK32 0xffffffffU #define JMAX 16 static unsigned int B[JMAX]; static void init (unsigned int *A) { int i; A[0] = 123456789; for (i = 1; i < JMAX; i++) A[i] = (663608941 * A[i - 1]) & MASK32; } int main (void) { int i; double u, som = 0; int n = 1000000; init (B); InitWELLRNG512a (B); for (i = 0; i < n; i++) { u = WELLRNG512a(); som += u; // printf ("%f\n", u); } printf ("somme = %f\n", som/n); return 0; }