long rk_poisson(rk_state *state, double lam) { if (lam >= 10) { return rk_poisson_ptrs(state, lam); } else if (lam == 0) { return 0; } else { return rk_poisson_mult(state, lam); } }
long poissonsample(double lam, mocapy::RandomGen* rg) { if (lam >= 10) { return rk_poisson_ptrs(lam, rg); } else if (lam == 0) { return 0; } else { return rk_poisson_mult(lam, rg); } }