def _logpdf(self, x, beta):
        return np.log(beta) - x**beta - special.gammaln(1. / beta)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/halfgennorm_gen._logpdf/swapped operands of `-`