scipy.special.log1p

scipy.special.log1p(x[, out]) = <ufunc 'log1p'>

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