fn poisson_sample(rng: &mut Xoshiro256PlusPlus, lambda: f64) -> u32
Knuth’s algorithm for Poisson sampling (λ < 30).