Lines Matching refs:b
45 double a, b, rx, ry, scale, t; in csqrt() local
48 b = cimag(z); in csqrt()
52 return (CMPLX(0, b)); in csqrt()
53 if (isinf(b)) in csqrt()
54 return (CMPLX(INFINITY, b)); in csqrt()
56 t = (b - b) / (b - b); /* raise invalid if b is not a NaN */ in csqrt()
67 return (CMPLX(fabs(b - b), copysign(a, b))); in csqrt()
69 return (CMPLX(a, copysign(b - b, b))); in csqrt()
71 if (isnan(b)) { in csqrt()
73 return (CMPLX(b + 0.0L + t, b + 0.0L + t)); /* NaN + NaN i */ in csqrt()
77 if (fabs(a) >= THRESH || fabs(b) >= THRESH) { in csqrt()
85 if (fabs(b) >= 0x1p-1020) in csqrt()
86 b *= 0.25; in csqrt()
93 if (fabs(a) < 0x1p-1022 && fabs(b) < 0x1p-1022) { in csqrt()
95 b *= 0x1p54; in csqrt()
101 t = sqrt((a + hypot(a, b)) * 0.5); in csqrt()
103 ry = scale * b / (2 * t); in csqrt()
105 t = sqrt((-a + hypot(a, b)) * 0.5); in csqrt()
106 rx = scale * fabs(b) / (2 * t); in csqrt()
107 ry = copysign(scale * t, b); in csqrt()