52  *		for x = +Inf				=> 0;
54  *		for x = +0				=> +Inf + divide-by-zero;
55  *		for x = -0				=> -Inf + divide-by-zero.
67  *		m0 = 0.5 * (1 + k/64) for m = [0.5,         0.5+127/256), k = [0, 63];
68  *		m0 = 1.0 * (0 + k/64) for m = [0.5+127/256, 1.0+127/128), k = [64, 127];
97 	if (n_n == 0)						\
125 		n_n = 0;  in __vrsqrt()
132 			if (hx >= 0x7ff00000)		/* X = NaN or Inf	*/  in __vrsqrt()
140 			if (hx < 0x00100000)		/* X = denormal, zero or negative	*/  in __vrsqrt()
143 				ax = hx & 0x7fffffff;  in __vrsqrt()
147 				if ((ax | lx) == 0)	/* |X| = zero	*/  in __vrsqrt()
151 				else if (hx >= 0)	/* X = denormal	*/  in __vrsqrt()
161 					sqrt_exp0 = (0x817 - (hx0 >> 21)) << 20;  in __vrsqrt()
162 					ind0 = (((hx0 >> 10) & 0x7f8) + 8) & -16;  in __vrsqrt()
164 					resh0 = (hx0 & 0x001fffff) | 0x3fe00000;  in __vrsqrt()
165 					res_ch0 = (resh0 + 0x00002000) & 0x7fffc000;  in __vrsqrt()
168 					LO(&res_c0) = 0;  in __vrsqrt()
170 					dexp_hi0 = ((double*)((char*)__vlibm_TBL_rsqrt + ind0))[0];  in __vrsqrt()
179 					LO(&dsqrt_exp0) = 0;  in __vrsqrt()
192 		if (n_n > 0)  in __vrsqrt()
195 	if (n > 0)  in __vrsqrt()
199 		if (hx >= 0x7ff00000)		/* X = NaN or Inf	*/  in __vrsqrt()
204 		else if (hx < 0x00100000)	/* X = denormal, zero or negative	*/  in __vrsqrt()
206 			ax = hx & 0x7fffffff;  in __vrsqrt()
210 			if ((ax | lx) == 0)	/* |X| = zero	*/  in __vrsqrt()
214 			else if (hx >= 0)	/* X = denormal	*/  in __vrsqrt()
224 				sqrt_exp0 = (0x817 - (hx0 >> 21)) << 20;  in __vrsqrt()
225 				ind0 = (((hx0 >> 10) & 0x7f8) + 8) & -16;  in __vrsqrt()
227 				resh0 = (hx0 & 0x001fffff) | 0x3fe00000;  in __vrsqrt()
228 				res_ch0 = (resh0 + 0x00002000) & 0x7fffc000;  in __vrsqrt()
231 				LO(&res_c0) = 0;  in __vrsqrt()
233 				dexp_hi0 = ((double*)((char*)__vlibm_TBL_rsqrt + ind0))[0];  in __vrsqrt()
242 				LO(&dsqrt_exp0) = 0;  in __vrsqrt()
259 			sqrt_exp0 = (0x5fe - (hx >> 21)) << 20;  in __vrsqrt()
260 			ind0 = (((hx >> 10) & 0x7f8) + 8) & -16;  in __vrsqrt()
262 			resh0 = (hx & 0x001fffff) | 0x3fe00000;  in __vrsqrt()
263 			res_ch0 = (resh0 + 0x00002000) & 0x7fffc000;  in __vrsqrt()
267 			LO(&res_c0) = 0;  in __vrsqrt()
269 			dexp_hi0 = ((double*)((char*)__vlibm_TBL_rsqrt + ind0))[0];  in __vrsqrt()
278 			LO(&dsqrt_exp0) = 0;  in __vrsqrt()
302 	LO(&dsqrt_exp0) = 0;  in __vrsqrt_n()
303 	LO(&dsqrt_exp1) = 0;  in __vrsqrt_n()
304 	LO(&dsqrt_exp2) = 0;  in __vrsqrt_n()
305 	LO(&res_c0) = 0;  in __vrsqrt_n()
306 	LO(&res_c1) = 0;  in __vrsqrt_n()
307 	LO(&res_c2) = 0;  in __vrsqrt_n()
323 		sqrt_exp0 = (0x5fe - (hx0 >> 21)) << 20;  in __vrsqrt_n()
324 		sqrt_exp1 = (0x5fe - (hx1 >> 21)) << 20;  in __vrsqrt_n()
325 		sqrt_exp2 = (0x5fe - (hx2 >> 21)) << 20;  in __vrsqrt_n()
326 		ind0 = (((hx0 >> 10) & 0x7f8) + 8) & -16;  in __vrsqrt_n()
327 		ind1 = (((hx1 >> 10) & 0x7f8) + 8) & -16;  in __vrsqrt_n()
328 		ind2 = (((hx2 >> 10) & 0x7f8) + 8) & -16;  in __vrsqrt_n()
330 		resh0 = (hx0 & 0x001fffff) | 0x3fe00000;  in __vrsqrt_n()
331 		resh1 = (hx1 & 0x001fffff) | 0x3fe00000;  in __vrsqrt_n()
332 		resh2 = (hx2 & 0x001fffff) | 0x3fe00000;  in __vrsqrt_n()
333 		res_ch0 = (resh0 + 0x00002000) & 0x7fffc000;  in __vrsqrt_n()
334 		res_ch1 = (resh1 + 0x00002000) & 0x7fffc000;  in __vrsqrt_n()
335 		res_ch2 = (resh2 + 0x00002000) & 0x7fffc000;  in __vrsqrt_n()
343 		dexp_hi0 = ((double*)((char*)__vlibm_TBL_rsqrt + ind0))[0];  in __vrsqrt_n()
344 		dexp_hi1 = ((double*)((char*)__vlibm_TBL_rsqrt + ind1))[0];  in __vrsqrt_n()
345 		dexp_hi2 = ((double*)((char*)__vlibm_TBL_rsqrt + ind2))[0];  in __vrsqrt_n()
380 	for(; n > 0 ; n--)  in __vrsqrt_n()
384 		sqrt_exp0 = (0x5fe - (hx0 >> 21)) << 20;  in __vrsqrt_n()
385 		ind0 = (((hx0 >> 10) & 0x7f8) + 8) & -16;  in __vrsqrt_n()
387 		resh0 = (hx0 & 0x001fffff) | 0x3fe00000;  in __vrsqrt_n()
388 		res_ch0 = (resh0 + 0x00002000) & 0x7fffc000;  in __vrsqrt_n()
392 		LO(&res_c0) = 0;  in __vrsqrt_n()
396 		dexp_hi0 = ((double*)((char*)__vlibm_TBL_rsqrt + ind0))[0];  in __vrsqrt_n()
405 		LO(&dsqrt_exp0) = 0;  in __vrsqrt_n()