def _entropy(self, a, c):
        val = special.psi(a)
        return a*(1-val) + 1.0/c*val + special.gammaln(a) - np.log(abs(c))
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/gengamma_gen._entropy/original