scipy.special.gammaln

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

y=gammaln(z) returns the base e logarithm of the absolute value of the gamma function of z: ln(|gamma(z)|)