def _logpmf(self, k, mu):
        Pk = gamln(k + 1) - mu - special.xlogy(k, mu)
        return Pk
dataset/ETHPy150Open scipy/scipy/scipy/stats/_discrete_distns.py/poisson_gen._logpmf/swapped operands of `-`