native/libfdlibm/e_sqrt.c

113         unsigned r,t1,s1,ix1,q1;
117         ix1 = __LO(x);          /* low word of x */
126             if(((ix0&(~sign))|ix1)==0) return x;/* sqrt(+-0) = +-0 */
135                 ix0 |= (ix1>>11); ix1 <<= 21;
139             ix0 |= (ix1>>(32-i));
140             ix1 <<= i;
145             ix0 += ix0 + ((ix1&sign)>>31);
146             ix1 += ix1;
151         ix0 += ix0 + ((ix1&sign)>>31);
152         ix1 += ix1;
163             ix0 += ix0 + ((ix1&sign)>>31);
164             ix1 += ix1;
172             if((t<ix0)||((t==ix0)&&(t1<=ix1))) {
176                 if (ix1 < t1) ix0 -= 1;
177                 ix1 -= t1;
180             ix0 += ix0 + ((ix1&sign)>>31);
181             ix1 += ix1;
186         if((ix0|ix1)!=0) {
199         ix1 =  q1>>1;
200         if ((q&1)==1) ix1 |= sign;
203         __LO(z) = ix1;