pub fn gaussian_ziggurat(r: &mut Rng, sigma: f64) -> f64