/* Compute 1-exp(p) See <http://cran.r-project.org/web/packages/Rmpfr/vignettes/log1mexp-note.pdf> */ double panda_log1mexp( double p) { return (p > M_LN2) ? log1p(-exp(-p)) : log(-expm1(-p)); }
float expm1f(float x) { return expm1(x); }