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