scipy.special.log1p

scipy.special.log1p(x[, out])

y=log1p(x) calculates log(1+x) for use when x is near zero.