def _logpdf(self, x, c, d):
        return log(c) + log(d) + special.xlogy(1-c, x) + special.xlog1py(-d-1, x**c)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/burr12_gen._logpdf/swapped operands of `-`