91 	double hfsq,f,s,z,R,w,t1,t2,dk;  in log()  local
111 	f = x-1.0;  in log()
113 	    if(f==zero) {  in log()
121 	    R = f*f*(0.5-0.33333333333333333*f);  in log()
122 	    if(k==0) return f-R; else {dk=(double)k;  in log()
123 	    	     return dk*ln2_hi-((R-dk*ln2_lo)-f);}  in log()
125  	s = f/(2.0+f);   in log()
136 	    hfsq=0.5*f*f;  in log()
137 	    if(k==0) return f-(hfsq-s*(hfsq+R)); else  in log()
138 		     return dk*ln2_hi-((hfsq-(s*(hfsq+R)+dk*ln2_lo))-f);  in log()
140 	    if(k==0) return f-s*(f-R); else  in log()
141 		     return dk*ln2_hi-((s*(f-R)-dk*ln2_lo)-f);  in log()