add accurate exp(x) - 1 via expm1()