823906a3efcf66897eac8a4c89052d9153bca49e,pymc3/math.py,,log1mexp_numpy,#Any#,240

Before Change

    For details, see
    https://cran.r-project.org/web/packages/Rmpfr/vignettes/log1mexp-note.pdf
    
    return np.where(x < 0.6931471805599453, np.log(-np.expm1(-x)), np.log1p(-np.exp(-x)))


def flatten_list(tensors):
    return tt.concatenate([var.ravel() for var in tensors])

After Change

    https://cran.r-project.org/web/packages/Rmpfr/vignettes/log1mexp-note.pdf
    
    x = np.asarray(x)
    out = np.empty_like(x)
    mask = x < 0.6931471805599453  // log(2)
    out[mask] = np.log(-np.expm1(-x[mask]))
    mask = ~mask
    out[mask] = np.log1p(-np.exp(-x[mask]))
    return out


def flatten_list(tensors):
    return tt.concatenate([var.ravel() for var in tensors])

Instances