if c == 0:
        return _gumbel_r_link(x, logsf, phat[1:], i - 1)