pub fn variable_poisson<F>( lambda: F, max_lambda: f64, domain: &Domain, ) -> Array2<f64>where F: Fn(&Array1<f64>) -> f64 + Sync + Send,
Poisson process on a d-dimensional region with variable intensity, using a rejection sampling algorithm.