pub fn log_1plusx(x: f64) -> f64
This routine computes \log(1 + x) for x > -1 using an algorithm that is accurate for small x.