I’m looking for a reliable library/module to sample a poisson distribution in Fortran. Thanks for any suggestions!
The interface for the function might be
function sample_poisson(mu) result(k)
real(dp), intent(in) :: mu
integer :: k
! ...
end function