Poisson deviate generator in noise.f90 module sometimes hangs (infinite loop ?)

[marcdegraef]

The current implementation of the Poisson random number generator sometime has issues and is generally rather slow; should be replaced by a more efficient routine, for instance the code from https://jblevins.org/mirror/amiller/random.f90