fabsf
fabsl, fabs, fabsf, x)
} else if (x >= FLT_EPSILON * fabsf(y - 1)) {
} else if (x >= FLT_EPSILON * fabsf(y - 1)) {
ax = fabsf(x);
ay = fabsf(y);
ax = fabsf(x);
ay = fabsf(y);
rx = fabsf(cimagf(w));
return (CMPLXF(fabsf(ry), rx));
return (CMPLXF(fabsf(ry), copysignf(rx, cimagf(z))));
ax = fabsf(x);
ay = fabsf(y);
ax = fabsf(x);
ay = fabsf(y);
#define isinf(x) (fabsf(x) == INFINITY)
w = one-fabsf(x);
else z=atanf(fabsf(y/x)); /* safe to do y/x */
t = expm1f(fabsf(x));
t = expf(fabsf(x));
if (ix < 0x42b17217) return half*expf(fabsf(x));
return __ldexp_expf(fabsf(x), -1);
a = fabsf(a);
b = fabsf(b);
w = fabsl(x+0.0L)-fabsf(y+0);
x = fabsf(x);
y = fabsf(x);
tmp = tmp*logf(fabsf(v*tmp));
if (fabsf(z) >= fabsf(w))
x = fabsf(x);
return -logf(fabsf(x));
nadj = logf(pi/fabsf(t*x));
ax = fabsf(x);
x = fabsf(x);
p = fabsf(p);
t = expm1f(fabsf(x));
if (ix < 0x42b17217) return h*expf(fabsf(x));
return h*2.0F*__ldexp_expf(fabsf(x), -1);
float fabsf(float) __pure2;
w = logf(fabsf(x))+ln2;
t = fabsf(x);
w =log1pf(fabsf(x)+t/(one+sqrtf(one+t)));
x = fabsf(x);
h = expf(fabsf(x)) / 2;
z = __ldexp_cexpf(CMPLXF(fabsf(x), y), -1);
ax = fabsf(x);
ay = fabsf(y);
h = expf(fabsf(x)) / 2;
z = __ldexp_cexpf(CMPLXF(fabsf(x), y), -1);
return (CMPLXF(fabsf(b - b), copysignf(a, b)));
return (CMPLXF(fabsf(b) / (2 * t), copysignf(t, b)));
float exp_mx = expf(-fabsf(x));
s = fabsf(x)-one;
x = fabsf(x);
s = fabsf(x)-one;
x = fabsf(x);
float ax = fabsf(x);
float ay = fabsf(y);
float ax = fabsf(x);
float ay = fabsf(y);
if(ix==0) return (float)-1.0/fabsf(x);
y = fabsf(y);
t = expm1f(two*fabsf(x));
t = expm1f(-two*fabsf(x));
if (fabsf(i) > 0x1.0p-4f)