{

g_rand[set] = g_mult * g_rand[set] + g_const;

return (unsigned int)((g_rand[set] >> 32) & PAF_RAND_MAX);

{

double f = exp(-lambda);

for (int i = 1; i <= x; i++) {

f *= lambda / double(i);

}

return f;

{

double f = exp(-lambda);

double c = f;

for (int i = 1; i <= x; i++) {

f *= lambda / double(i);

c += f;

}

return c;

{

if (p <= 0) {

return 0;

}

if (p >= 1) {

// passing this 1 will cause an infinite loop, try this instead:

p = 1-1e-9;

}

double f = exp(-lambda);

int i = 0;

for (double c = f; c < p; c += f) {

i++;

f *= lambda / double(i);

}

return i;