pub fn ln_exp_1p(x: f32) -> f32