if a < 0:
        with np.errstate(divide='ignore'):
            nla = _norm_logcdf(a)
            nlb = _norm_logcdf(b)
            tab = np.log1p(-np.exp(nla - nlb))