int main(int ac, char *av[]) {
if (ac == 1) {
printf("Usage: randmtzig <n>\n");
return (-1);
}
int n = atoi(av[1]);
time_t t1;
dsfmt_gv_init_gen_rand(0);
double *p; posix_memalign((void **)&p, 16, n*sizeof(double));
uint32_t *u; posix_memalign((void **)&u, 16, 2*n*sizeof(uint32_t));
t1 = clock();
dsfmt_gv_fill_array_close_open(p, n);
printf("Uniform fill (n): %f\n", (clock() - t1) / (double) CLOCKS_PER_SEC);
t1 = clock();
for (int i = 0; i < n; i++) p[i] = dsfmt_gv_genrand_close_open();
printf("Uniform (n): %f\n", (clock() - t1) / (double) CLOCKS_PER_SEC);
t1 = clock();
for (int i = 0; i < 2*n; i++) u[i] = dsfmt_gv_genrand_uint32();
printf("Uniform 32-bit ints (2*n): %f\n", (clock() - t1) / (double) CLOCKS_PER_SEC);
memset((void *)p, 0, n*sizeof(double));
t1 = clock();
for (int i = 0; i < n; i++) p[i] = randmtzig_gv_randn();
printf("Normal (n): %f\n", (clock() - t1) / (double) CLOCKS_PER_SEC);
for (int i = 0; i < 10; i++) printf("%lf\n", p[i]);
return 0;
}
void randmtzig_gv_fill_randn (double *p, randmtzig_idx_type n)
{
randmtzig_idx_type i;
for (i = 0; i < n; i++)
p[i] = randmtzig_gv_randn();
}
