9 static int gauss_dev(struct gauss_state *gs)
18 vr = gs->stddev * (r / (FRAND32_MAX + 1.0));
20 return vr - gs->stddev / 2;
23 unsigned long long gauss_next(struct gauss_state *gs)
25 unsigned long long sum = 0;
28 for (i = 0; i < GAUSS_ITERS; i++)
29 sum += __rand(&gs->r) % (gs->nranges + 1);
31 sum = (sum + GAUSS_ITERS - 1) / GAUSS_ITERS;
34 int dev = gauss_dev(gs);
36 while (dev + sum >= gs->nranges)
41 return __hash_u64(sum) % gs->nranges;
44 void gauss_init(struct gauss_state *gs, unsigned long nranges, double dev,
47 memset(gs, 0, sizeof(*gs));
48 init_rand_seed(&gs->r, seed, 0);
49 gs->nranges = nranges;
52 gs->stddev = ceil((double) (nranges * 100.0) / dev);
53 if (gs->stddev > nranges / 2)
54 gs->stddev = nranges / 2;