gonzui

t2ex

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

82:                         r = fabs (x);
83:                         r = sqrt (r);
84:                         if (x < 0.0) {
94:                 r = fabs (y);
95:                 r = sqrt (0.5*r);
96:                 if (y > 0)
120:         if (x > 0) {
121:                 t = sqrt(0.5 * r + 0.5 * x);
122:                 r = scale * fabs((0.5 * y) / t);
125:         else {
126:                 r = sqrt( 0.5 * r - 0.5 * x );
127:                 t = scale * fabs( (0.5 * y) / r );