log1pl
long double log1pl(long double);
*rx = log1pl(Am1 + sqrtl(Am1*(A+1)));
*rx = log1pl((y-1) + sqrtl((y-1)*(y+1)));
rx = log1pl(4*ax / sum_squares(ax-1, ay)) / 4;
__weak_alias(log1pl, _log1pl)
RETURNI(log1pl(t+sqrtl(2.0*t+t*t)));
t = 0.5*log1pl(t+t*x/(one-x));
t = 0.5*log1pl((x+x)/(one-x));
w =log1pl(fabsl(x)+t/(one+sqrtl(one+t)));
RETURNI(CMPLXL(log1pl(ay2l + t + sh) / 2, v));
RETURNI(CMPLXL(log1pl(ay * ay) / 2, v));
__weak_alias(log1pl, _log1pl)
__weak_alias(log1pl, _log1pl)
CHECKL_NAN(i, log1pl, log1pf_invalid[i]);
CHECKL_NAN(i, log1pl, log1p_invalid[i]);
CHECKL_NAN(i, log1pl, log1pl_invalid[i]);
CHECKL_EQ(0, log1pl, -1., -HUGE_VALL);
CHECKL_EQ(01, log1pl, -LDBL_MIN, -LDBL_MIN);
CHECKL_EQ(1, log1pl, -0., 0);
CHECKL_EQ(2, log1pl, +0., 0);
CHECKL_EQ(3, log1pl, 1, logl(2));
CHECKL_EQ(0, log1pl, -LDBL_DENORM_MIN, -LDBL_DENORM_MIN);
ATF_CHECK_MSG(fabsl((log1pl(em1) - 1)/1) < 2*LDBL_EPSILON,
"log1pl(e)=%La=%.34Lg", log1pl(em1), log1pl(em1));
ATF_CHECK_MSG(fabsl((log1pl(e2m1) - 2)/2) < 2*LDBL_EPSILON,
"log1pl(e^2)=%La=%.34Lg", log1pl(em1), log1pl(em1));
CHECKL_EQ(0, log1pl, INFINITY, INFINITY);