log1pf
float log1pf(float);
*rx = log1pf(Am1 + sqrtf(Am1 * (A + 1)));
*rx = log1pf((y - 1) + sqrtf((y - 1) * (y + 1)));
rx = log1pf(4 * ax / sum_squares(ax - 1, ay)) / 4;
return( z*log1pf(x+x) );
__weak_alias(log1pf, _log1pf)
log1pf(float x)
return log1pf(t+sqrtf((float)2.0*t+t*t));
t = (float)0.5*log1pf(t+t*x/(one-x));
t = (float)0.5*log1pf((x+x)/(one-x));
w =log1pf(fabsf(x)+t/(one+__ieee754_sqrtf(one+t)));
__weak_alias(log1pf, _log1pf)
log1pf(float x)
CHECK_NAN(i, log1pf, log1pf_invalid[i]);
CHECK_EQ(0, log1pf, -1., -HUGE_VALF);
CHECK_EQ(0, log1pf, -FLT_MIN, -FLT_MIN);
CHECK_EQ(1, log1pf, -0., 0);
CHECK_EQ(2, log1pf, +0., 0);
CHECK_EQ(3, log1pf, 1, logf(2));
CHECK_EQ(0, log1pf, -FLT_DENORM_MIN, -FLT_DENORM_MIN);
ATF_CHECK_MSG(fabsf((log1pf(em1) - 1)/1) < 2*FLT_EPSILON,
"log1pf(e)=%a=%.8g", log1pf(em1), log1pf(em1));
ATF_CHECK_MSG(fabsf((log1pf(e2m1) - 2)/2) < 2*FLT_EPSILON,
"log1pf(e^2)=%a=%.8g", log1pf(em1), log1pf(em1));
CHECK_EQ(0, log1pf, INFINITY, INFINITY);