scipy.special
scipy.special.
expm1
exp(x) - 1 for use when x is near zero.
scipy.special.log1p
scipy.special.cosm1