Numerically stable contractions (expm1 and friends).
Switch branches/tags
Nothing to show
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Permalink
Failed to load latest commit information.
README.org
floating-point-contractions.asd
floating-point-contractions.lisp
package.lisp

README.org

This system provides some numerically stable, contracted versions of common but unstable floating-point operations.

Export(AKA)Computes
log1+log1plog (1 + x)
log1-log1mlog (1 - x)
log1+/x(log (1 + x))/x
exp-1expm1exp(x) - 1
exp-1/x(exp(x)-1)/x
expt-1powm1(x^y)-1
log1-explog1mexplog(1 - exp(x))
log1+explog1pexplog(1 + exp(x))
log2-explog2mexplog(2 - exp(x))
logexp-1logexpm1log(exp(a) -1)
hypotsqrt(x^2 + y^2)