40 	float hfsq,f,c,s,z,R,u;  in log1pf()  local
60 		k=0;f=x;hu=1;}		/* sqrt(2)/2- <= 1+x < sqrt(2)+ */  in log1pf()
92 	    f = u-(float)1.0;  in log1pf()
94 	hfsq=(float)0.5*f*f;  in log1pf()
96 	    if(f==zero) if(k==0) return zero;  in log1pf()
98 	    R = hfsq*((float)1.0-(float)0.66666666666666666*f);  in log1pf()
99 	    if(k==0) return f-R; else  in log1pf()
100 	    	     return k*ln2_hi-((R-(k*ln2_lo+c))-f);  in log1pf()
102  	s = f/((float)2.0+f);  in log1pf()
105 	if(k==0) return f-(hfsq-s*(hfsq+R)); else  in log1pf()
106 		 return k*ln2_hi-((hfsq-(s*(hfsq+R)+(k*ln2_lo+c)))-f);  in log1pf()