181  2.6176170809902549338711E-20L,
182 -1.0126791927256478897086E-20L,
184  1.2207982955417546912101E-20L,
186  1.3164426894366316434230E-20L,
187 -1.8527916071632873716786E-20L,
188  1.8950325588932570796551E-20L,
189  1.5564775779538780478155E-20L,
191 -2.0208749253662532228949E-20L,
192  1.4966292219224761844552E-20L,
195 -1.2327176863327626135542E-20L,
220 #define F W
223 #define G W
226 #define H W
234 static _Thread_local long double w, W, Wa, Wb, ya, yb, u;  variable
298 w = floorl(y);  in powl()
301 if( w == y )  in powl()
310 	yb = 0.5L * fabsl(w);  in powl()
366 	i = w;  in powl()
367 	w = floorl(x);  in powl()
368 	if( (w == x) && (fabsl(y) < 32768.0) )  in powl()
370 		w = powil( x, (int) y );  in powl()
371 		return( w );  in powl()
415 w = x * ( z * __polevll( x, P, 3 ) / __p1evll( x, Q, 3 ) );  in powl()
416 w = w - ldexpl( z, -1 );   /*  w - 0.5 * z  */  in powl()
421 z = LOG2EA * w;  in powl()
422 z += w;  in powl()
427 w = -i;  in powl()
428 w = ldexpl( w, -LNXT );	/* divide by NXT */  in powl()
429 w += e;  in powl()
430 /* Now base 2 log of x is w + z. */  in powl()
440 /* (w+z)(ya+yb)  in powl()
441  * = w*ya + w*yb + z*y  in powl()
443 F = z * y  +  w * yb;  in powl()
447 G = Fa + w * ya;  in powl()
453 w = ldexpl( Ga+Ha, LNXT );  in powl()
456 if( w > MEXP )  in powl()
459 if( w < MNEXP )  in powl()
462 e = w;  in powl()
487 w = douba( e );  in powl()
488 z = w * z;      /*    2**-e * ( 1 + (2**Hb-1) )    */  in powl()
489 z = z + w;  in powl()
498 	w = ldexpl( y, -1 );  in powl()
499 	w = floorl(w);  in powl()
500 	w = ldexpl( w, 1 );  in powl()
501 	if( w != y )  in powl()