40  *	1.  (anything) ** 0  is 1
43  *	4.  NAN ** (anything except 0) is NAN
45  *	6.  +-(|x| > 1) **  -INF is +0
46  *	7.  +-(|x| < 1) **  +INF is +0
49  *	10. +0 ** (+anything except 0, NAN)               is +0
50  *	11. -0 ** (+anything except 0, NAN, odd integer)  is +0
51  *	12. +0 ** (-anything except 0, NAN)               is +INF
52  *	13. -0 ** (-anything except 0, NAN, odd integer)  is +INF
53  *	14. -0 ** (odd integer) = -( +0 ** (odd integer) )
54  *	15. +INF ** (+anything except 0,NAN) is +INF
55  *	16. +INF ** (-anything except 0,NAN) is +0
56  *	17. -INF ** (anything)  = -0 ** (-anything)
58  *	19. (-anything except 0 and inf) ** (non-integer) is NAN
114    0 <= x <= 0.5
157   ix = hx & 0x7fffffff;  in powl()
161   iy = hy & 0x7fffffff;  in powl()
164   /* y==zero: x**0 = 1 */  in powl()
165   if ((iy | q.parts32.mswlo | q.parts32.lswhi | q.parts32.lswlo) == 0)  in powl()
171   if (x == -1.0L && iy == 0x7fff0000  in powl()
172       && (q.parts32.mswlo | q.parts32.lswhi | q.parts32.lswlo) == 0)  in powl()
176   if ((ix > 0x7fff0000)  in powl()
177       || ((ix == 0x7fff0000)  in powl()
178 	  && ((p.parts32.mswlo | p.parts32.lswhi | p.parts32.lswlo) != 0))  in powl()
179       || (iy > 0x7fff0000)  in powl()
180       || ((iy == 0x7fff0000)  in powl()
181 	  && ((q.parts32.mswlo | q.parts32.lswhi | q.parts32.lswlo) != 0)))  in powl()
184   /* determine if y is an odd int when x < 0  in powl()
185    * yisint = 0       ... y is not an integer  in powl()
189   yisint = 0;  in powl()
190   if (hx < 0)  in powl()
192       if (iy >= 0x40700000)	/* 2^113 */  in powl()
194       else if (iy >= 0x3fff0000)	/* 1.0 */  in powl()
208   if ((q.parts32.mswlo | q.parts32.lswhi | q.parts32.lswlo) == 0)  in powl()
210       if (iy == 0x7fff0000)	/* y is +-inf */  in powl()
212 	  if (((ix - 0x3fff0000) | p.parts32.mswlo | p.parts32.lswhi |  in powl()
213 	    p.parts32.lswlo) == 0)  in powl()
215 	  else if (ix >= 0x3fff0000)	/* (|x|>1)**+-inf = inf,0 */  in powl()
216 	    return (hy >= 0) ? y : zero;  in powl()
217 	  else			/* (|x|<1)**-,+inf = inf,0 */  in powl()
218 	    return (hy < 0) ? -y : zero;  in powl()
220       if (iy == 0x3fff0000)  in powl()
222 	  if (hy < 0)  in powl()
227       if (hy == 0x40000000)  in powl()
229       if (hy == 0x3ffe0000)  in powl()
231 	  if (hx >= 0)		/* x >= +0 */  in powl()
238   if ((p.parts32.mswlo | p.parts32.lswhi | p.parts32.lswlo) == 0)  in powl()
240       if (ix == 0x7fff0000 || ix == 0 || ix == 0x3fff0000)  in powl()
242 	  z = ax;		/*x is +-0,+-inf,+-1 */  in powl()
243 	  if (hy < 0)  in powl()
245 	  if (hx < 0)  in powl()
247 	      if (((ix - 0x3fff0000) | yisint) == 0)  in powl()
252 		z = -z;		/* (x<0)**odd = -(|x|**odd) */  in powl()
258   /* (x<0)**(non-int) is NaN */  in powl()
259   if (((((u_int32_t) hx >> 31) - 1) | yisint) == 0)  in powl()
265   if (iy > 0x401d654b)  in powl()
268       if (iy > 0x407d654b)  in powl()
270 	  if (ix <= 0x3ffeffff)  in powl()
271 	    return (hy < 0) ? huge * huge : tiny * tiny;  in powl()
272 	  if (ix >= 0x3fff0000)  in powl()
273 	    return (hy > 0) ? huge * huge : tiny * tiny;  in powl()
276       if (ix < 0x3ffeffff)  in powl()
277 	return (hy < 0) ? huge * huge : tiny * tiny;  in powl()
278       if (ix > 0x3fff0000)  in powl()
279 	return (hy > 0) ? huge * huge : tiny * tiny;  in powl()
282   n = 0;  in powl()
284   if (ix < 0x00010000)  in powl()
291   n += ((ix) >> 16) - 0x3fff;  in powl()
292   j = ix & 0x0000ffff;  in powl()
294   ix = j | 0x3fff0000;		/* normalize ix */  in powl()
295   if (j <= 0x3988)  in powl()
296     k = 0;			/* |x|<sqrt(3/2) */  in powl()
297   else if (j < 0xbb67)  in powl()
301       k = 0;  in powl()
303       ix -= 0x00010000;  in powl()
311   u = ax - bp[k];		/* bp[0]=1.0, bp[1]=1.5 */  in powl()
317   o.parts32.lswlo = 0;  in powl()
318   o.parts32.lswhi &= 0xf8000000;  in powl()
323   o.parts32.lswlo = 0;  in powl()
324   o.parts32.lswhi &= 0xf8000000;  in powl()
330   u = LN[0] + s2 * (LN[1] + s2 * (LN[2] + s2 * (LN[3] + s2 * LN[4])));  in powl()
331   v = LD[0] + s2 * (LD[1] + s2 * (LD[2] + s2 * (LD[3] + s2 * (LD[4] + s2))));  in powl()
337   o.parts32.lswlo = 0;  in powl()
338   o.parts32.lswhi &= 0xf8000000;  in powl()
347   o.parts32.lswlo = 0;  in powl()
348   o.parts32.lswhi &= 0xf8000000;  in powl()
357   o.parts32.lswlo = 0;  in powl()
358   o.parts32.lswhi &= 0xf8000000;  in powl()
364   if (((((u_int32_t) hx >> 31) - 1) | (yisint - 1)) == 0)  in powl()
370   o.parts32.lswlo = 0;  in powl()
371   o.parts32.lswhi &= 0xf8000000;  in powl()
378   if (j >= 0x400d0000) /* z >= 16384 */  in powl()
381       if (((j - 0x400d0000) | o.parts32.mswlo | o.parts32.lswhi |  in powl()
382 	o.parts32.lswlo) != 0)  in powl()
390   else if ((j & 0x7fffffff) >= 0x400d01b9)	/* z <= -16495 */  in powl()
393       if (((j - 0xc00d01bc) | o.parts32.mswlo | o.parts32.lswhi |  in powl()
395 	  != 0)  in powl()
404   i = j & 0x7fffffff;  in powl()
405   k = (i >> 16) - 0x3fff;  in powl()
406   n = 0;  in powl()
407   if (i > 0x3ffe0000)  in powl()
415   o.parts32.lswlo = 0;  in powl()
416   o.parts32.lswhi &= 0xf8000000;  in powl()
424   u = PN[0] + t * (PN[1] + t * (PN[2] + t * (PN[3] + t * PN[4])));  in powl()
425   v = PD[0] + t * (PD[1] + t * (PD[2] + t * (PD[3] + t)));  in powl()
432   if ((j >> 16) <= 0)  in powl()