Bugfix/feature 1053 overflow poisson rng

This isn't all your fault, but instead of C-style casts `(double*)0`, we want C++-style `static_cast<double*>(0)`.  Would you mind changing all these?