log_expit function

Numerically accurately compute f(x) = log(x / (1/x)).