def _logpdf(self, x, a, c):
        negxc = -x**c
        exm1c = -special.expm1(negxc)
        logp = (log(a) + log(c) + special.xlogy(a - 1.0, exm1c) +
                negxc + special.xlogy(1.0 - c, x))
        return logp
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/exponweib_gen._logpdf/swapped operands of `-`