pymc.math.log1mexp

pymc.math.log1mexp(x, *, negative_input=False)

Return log(1 - exp(-x)).

This function is numerically more stable than the naive approach.

For details, see https://cran.r-project.org/web/packages/Rmpfr/vignettes/log1mexp-note.pdf

References

[Machler2012]

Martin Mächler (2012). “Accurately computing log(1-exp(- mid a mid)) Assessed by the Rmpfr package”