21 dp_h[] = { 0.0, 5.84960938e-01,}, /* 0x3f15c000 */
22 dp_l[] = { 0.0, 1.56322085e-06,}, /* 0x35d1cfdc */
26 thrd    =  3.33333343e-01, /* 0x3eaaaaab */
29 two24	=  16777216.0,	/* 0x4b800000 */
33 L1  =  6.0000002384e-01, /* 0x3f19999a */
34 L2  =  4.2857143283e-01, /* 0x3edb6db7 */
35 L3  =  3.3333334327e-01, /* 0x3eaaaaab */
36 L4  =  2.7272811532e-01, /* 0x3e8ba305 */
37 L5  =  2.3066075146e-01, /* 0x3e6c3255 */
38 L6  =  2.0697501302e-01, /* 0x3e53f142 */
39 P1   =  1.6666667163e-01, /* 0x3e2aaaab */
40 P2   = -2.7777778450e-03, /* 0xbb360b61 */
41 P3   =  6.6137559770e-05, /* 0x388ab355 */
42 P4   = -1.6533901999e-06, /* 0xb5ddea0e */
43 P5   =  4.1381369442e-08, /* 0x3331bb4c */
44 lg2  =  6.9314718246e-01, /* 0x3f317218 */
45 lg2_h  =  6.93145752e-01, /* 0x3f317200 */
46 lg2_l  =  1.42860654e-06, /* 0x35bfbe8c */
48 cp    =  9.6179670095e-01, /* 0x3f76384f =2/(3ln2) */
49 cp_h  =  9.6191406250e-01, /* 0x3f764000 =12b cp */
50 cp_l  = -1.1736857402e-04, /* 0xb8f623c6 =tail of cp_h */
51 ivln2    =  1.4426950216e+00, /* 0x3fb8aa3b =1/ln2 */
52 ivln2_h  =  1.4426879883e+00, /* 0x3fb8aa00 =16b 1/ln2*/
53 ivln2_l  =  7.0526075433e-06; /* 0x36eca570 =1/ln2 tail*/
65 	ix = hx&0x7fffffff;  iy = hy&0x7fffffff;  in powf()
67     /* y==zero: x**0 = 1 */  in powf()
68 	if(iy==0) return one;  in powf()
71 	if (hx==0x3f800000) return one;  in powf()
74 	if(ix > 0x7f800000 ||  in powf()
75 	   iy > 0x7f800000)  in powf()
78     /* determine if y is an odd int when x < 0  in powf()
79      * yisint = 0	... y is not an integer  in powf()
83 	yisint  = 0;  in powf()
84 	if(hx<0) {  in powf()
85 	    if(iy>=0x4b800000) yisint = 2; /* even integer y */  in powf()
86 	    else if(iy>=0x3f800000) {  in powf()
87 		k = (iy>>23)-0x7f;	   /* exponent */  in powf()
94 	if (iy==0x7f800000) {	/* y is +-inf */  in powf()
95 	    if (ix==0x3f800000)  in powf()
97 	    else if (ix > 0x3f800000)/* (|x|>1)**+-inf = inf,0 */  in powf()
98 	        return (hy>=0)? y: zero;  in powf()
99 	    else			/* (|x|<1)**-,+inf = inf,0 */  in powf()
100 	        return (hy<0)?-y: zero;  in powf()
102 	if(iy==0x3f800000) {	/* y is  +-1 */  in powf()
103 	    if(hy<0) return one/x; else return x;  in powf()
105 	if(hy==0x40000000) return x*x; /* y is  2 */  in powf()
106 	if(hy==0x3f000000) {	/* y is  0.5 */  in powf()
107 	    if(hx>=0)	/* x >= +0 */  in powf()
113 	if(ix==0x7f800000||ix==0||ix==0x3f800000){  in powf()
114 	    z = ax;			/*x is +-0,+-inf,+-1*/  in powf()
115 	    if(hy<0) z = one/z;	/* z = (1/|x|) */  in powf()
116 	    if(hx<0) {  in powf()
117 		if(((ix-0x3f800000)|yisint)==0) {  in powf()
120 		    z = -z;		/* (x<0)**odd = -(|x|**odd) */  in powf()
127     /* (x<0)**(non-int) is NaN */  in powf()
128 	if((n|yisint)==0) return (x-x)/(x-x);  in powf()
131 	if((n|(yisint-1))==0) sn = -one;/* (-ve)**(odd int) */  in powf()
134 	if(iy>0x4d000000) { /* if |y| > 2**27 */  in powf()
136 	    if(ix<0x3f7ffff6) return (hy<0)? sn*huge*huge:sn*tiny*tiny;  in powf()
137 	    if(ix>0x3f800007) return (hy>0)? sn*huge*huge:sn*tiny*tiny;  in powf()
146 	    SET_FLOAT_WORD(t1,is&0xfffff000);  in powf()
150 	    n = 0;  in powf()
152 	    if(ix<0x00800000)  in powf()
154 	    n  += ((ix)>>23)-0x7f;  in powf()
155 	    j  = ix&0x007fffff;  in powf()
157 	    ix = j|0x3f800000;		/* normalize ix */  in powf()
158 	    if(j<=0x1cc471) k=0;	/* |x|<sqrt(3/2) */  in powf()
159 	    else if(j<0x5db3d7) k=1;	/* |x|<sqrt(3)   */  in powf()
160 	    else {k=0;n+=1;ix -= 0x00800000;}  in powf()
164 	    u = ax-bp[k];		/* bp[0]=1.0, bp[1]=1.5 */  in powf()
169 	    SET_FLOAT_WORD(s_h,is&0xfffff000);  in powf()
171 	    is = ((ix>>1)&0xfffff000)|0x20000000;  in powf()
172 	    SET_FLOAT_WORD(t_h,is+0x00400000+(k<<21));  in powf()
182 	    SET_FLOAT_WORD(t_h,is&0xfffff000);  in powf()
190 	    SET_FLOAT_WORD(p_h,is&0xfffff000);  in powf()
198 	    SET_FLOAT_WORD(t1,is&0xfffff000);  in powf()
204 	SET_FLOAT_WORD(y1,is&0xfffff000);  in powf()
209 	if (j>0x43000000)				/* if z > 128 */  in powf()
211 	else if (j==0x43000000) {			/* if z == 128 */  in powf()
214 	else if ((j&0x7fffffff)>0x43160000)		/* z <= -150 */  in powf()
216 	else if (j==0xc3160000){			/* z == -150 */  in powf()
222 	i = j&0x7fffffff;  in powf()
223 	k = (i>>23)-0x7f;  in powf()
224 	n = 0;  in powf()
225 	if(i>0x3f000000) {		/* if |z| > 0.5, set n = [z+0.5] */  in powf()
226 	    n = j+(0x00800000>>(k+1));  in powf()
227 	    k = ((n&0x7fffffff)>>23)-0x7f;	/* new k for n */  in powf()
228 	    SET_FLOAT_WORD(t,n&~(0x007fffff>>k));  in powf()
229 	    n = ((n&0x007fffff)|0x00800000)>>(23-k);  in powf()
230 	    if(j<0) n = -n;  in powf()
235 	SET_FLOAT_WORD(t,is&0xffff8000);  in powf()
246      * sign bit of z is 0.  in powf()
247      * sign bit of j will indicate sign of 0x7f-biased exponent.  in powf()
250 	if((j>>23)<=0) z = scalbnf(z,n);	/* subnormal output */  in powf()