def _isf(self, q, a, c):
        val1 = special.gammaincinv(a, q)
        val2 = special.gammainccinv(a, q)
        return np.where(0 > c, val2, val1)**(1.0/c)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/gengamma_gen._isf/swapped operands of `>`