gonzui

t2ex

bsd_source/lib/libc/src_bsd/complex/s_csqrt.c - 3.8KB - 144 lines

81:                 else {
82:                         r = fabs (x);
83:                         r = sqrt (r);
93:         if (x == 0.0) {
94:                 r = fabs (y);
95:                 r = sqrt (0.5*r);
102:         /* Rescale to avoid internal overflow or underflow.  */
103:         if ((fabs(x) > 4.0) || (fabs(y) > 4.0)) {
104:                 x *= 0.25;
121:                 t = sqrt(0.5 * r + 0.5 * x);
122:                 r = scale * fabs((0.5 * y) / t);
123:                 t *= scale;
126:                 r = sqrt( 0.5 * r - 0.5 * x );
127:                 t = scale * fabs( (0.5 * y) / r );
128:                 r *= scale;