gonzui

t2ex

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

74:                 if (x < 0.0f) {
75:                         w = 0.0f + sqrtf(-x) * I;
76:                         return (w);
81:                 else {
82:                         w = sqrtf(x) + y * I;
83:                         return (w);
88:                 r = fabsf(y);
89:                 r = sqrtf(0.5f*r);
90:                 if(y > 0)
115:         if (x > 0) {
116:                 t = sqrtf( 0.5f * r + 0.5f * x );
117:                 r = scale * fabsf((0.5f * y) / t);
120:         else {
121:                 r = sqrtf(0.5f * r - 0.5f * x);
122:                 t = scale * fabsf((0.5f * y) / r);