pub fn exponential<R: RngCore + ?Sized>(rng: &mut R, mean: f64) -> f64
Draw an exponential variate with the given mean via inverse-CDF.
mean
-mean * ln(1 - u), computed with ln_1p(-u) for accuracy.
-mean * ln(1 - u)
ln_1p(-u)