def _sf(self, x, a):
        fac = 0.5*special.gammainc(a, abs(x))
        return where(0 > x, 0.5-fac, 0.5+fac)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/dgamma_gen._sf/swapped operands of `>`