22 ln2_hi =   6.9313812256e-01,	/* 0x3f317180 */
23 ln2_lo =   9.0580006145e-06,	/* 0x3717f7d1 */
24 two25 =    3.355443200e+07,	/* 0x4c000000 */
43 	ax = hx&0x7fffffff;  in log1pf()
46 	if (hx < 0x3ed413d0) {			/* 1+x < sqrt(2)+  */  in log1pf()
47 	    if(ax>=0x3f800000) {		/* x <= -1.0 */  in log1pf()
51 	    if(ax<0x38000000) {			/* |x| < 2**-15 */  in log1pf()
53 	            &&ax<0x33800000) 		/* |x| < 2**-24 */  in log1pf()
58 	    if(hx>0||hx<=((int32_t)0xbe95f619)) {  in log1pf()
59 		k=0;f=x;hu=1;}		/* sqrt(2)/2- <= 1+x < sqrt(2)+ */  in log1pf()
61 	if (hx >= 0x7f800000) return x+x;  in log1pf()
62 	if(k!=0) {  in log1pf()
63 	    if(hx<0x5a000000) {  in log1pf()
68 	        c  = (k>0)? (float)1.0-(u-x):x-(u-(float)1.0);  in log1pf()
74 		c  = 0;  in log1pf()
76 	    hu &= 0x007fffff;  in log1pf()
80 	     * strict bounds than the one here so that the k==0 case is  in log1pf()
82 	     * using the correction term but don't use it if k==0.  in log1pf()
84 	    if(hu<0x3504f4) {			/* u < sqrt(2) */  in log1pf()
85 	        SET_FLOAT_WORD(u,hu|0x3f800000);/* normalize u */  in log1pf()
88 		SET_FLOAT_WORD(u,hu|0x3f000000);	/* normalize u/2 */  in log1pf()
89 	        hu = (0x00800000-hu)>>2;  in log1pf()
94 	if(hu==0) {	/* |f| < 2**-20 */  in log1pf()
96 		if(k==0) {  in log1pf()
104 	    if(k==0) return f-R; else  in log1pf()
110 	if(k==0) return f-(hfsq-s*(hfsq+R)); else  in log1pf()