def _logpdf(self, x, b):
        xb = x**b
        f = 1 + log(b) + special.xlogy(b - 1.0, x) + xb - exp(xb)
        return f
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/exponpow_gen._logpdf/original