def _pdf(self, x, c):
        cx = c*x
        logex2 = where((c == 0)*(x == x), 0.0, special.log1p(-cx))
        logpex2 = where((c == 0)*(x == x), -x, logex2/c)
        pex2 = exp(logpex2)
        # Handle special cases
        logpdf = where((cx == 1) | (cx == -inf), -inf, -pex2+logex2-logpex2)
        putmask(logpdf, (c == 1) & (x == 1), 0.0)
        return exp(logpdf)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/genextreme_gen._pdf/swapped operands of `-`