96 	int hx,ix;  local
99 	ix = hx&0x7fffffff;
100 	if(ix>=0x7ff00000) return one/(x*x);
102 	if(ix >= 0x40000000) {	/* |x| >= 2.0 */
107 		if(ix<0x7fe00000) {  /* make sure x+x not overflow */
116 		if(ix>0x48000000) z = (invsqrtpi*cc)/ieee_sqrt(x);
123 	if(ix<0x3f200000) {	/* |x| < 2**-13 */
125 	        if(ix<0x3e400000) return one;	/* |x|<2**-27 */
132 	if(ix < 0x3FF00000) {	/* |x| < 1.00 */
165 	int hx,ix,lx;  local
168         ix = 0x7fffffff&hx;
171 	if(ix>=0x7ff00000) return  one/(x+x*x); 
172         if((ix|lx)==0) return -one/zero;
174         if(ix >= 0x40000000) {  /* |x| >= 2.0 */
194                 if(ix<0x7fe00000) {  /* make sure x+x not overflow */
199                 if(ix>0x48000000) z = (invsqrtpi*ss)/ieee_sqrt(x);
206 	if(ix<=0x3e400000) {	/* x < 2**-27 */
333 	int ix;  local
334 	ix = 0x7fffffff&__HI(x);
335 	if(ix>=0x40200000)     {p = pR8; q= pS8;}
336 	else if(ix>=0x40122E8B){p = pR5; q= pS5;}
337 	else if(ix>=0x4006DB6D){p = pR3; q= pS3;}
338 	else if(ix>=0x40000000){p = pR2; q= pS2;}
468 	int ix;  local
469 	ix = 0x7fffffff&__HI(x);
470 	if(ix>=0x40200000)     {p = qR8; q= qS8;}
471 	else if(ix>=0x40122E8B){p = qR5; q= qS5;}
472 	else if(ix>=0x4006DB6D){p = qR3; q= qS3;}
473 	else if(ix>=0x40000000){p = qR2; q= qS2;}