gonzui

t2ex

bsd_source/lib/libc/src_bsd/complex/s_csqrtf.c - 3.4KB - 131 lines

87:         if (x == 0.0f) {
88:                 r = fabsf(y);
89:                 r = sqrtf(0.5f*r);
97:         /* Rescale to avoid internal overflow or underflow.  */
98:         if ((fabsf(x) > 4.0f) || (fabsf(y) > 4.0f)) {
99:                 x *= 0.25f;
116:                 t = sqrtf( 0.5f * r + 0.5f * x );
117:                 r = scale * fabsf((0.5f * y) / t);
118:                 t *= scale;
121:                 r = sqrtf(0.5f * r - 0.5f * x);
122:                 t = scale * fabsf((0.5f * y) / r);
123:                 r *= scale;