log1p(x) := log(1+x)

expm1(x) := exp(x) - 1