def _logpdf(self, x, a, c):
        return np.log(abs(c)) + special.xlogy(c*a - 1, x) - x**c - special.gammaln(a)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/gengamma_gen._logpdf/original