35     float hfsq, f, s, z, R, w, t1, t2, dk;  in logf()  local
57     f = x - (float)1.0;  in logf()
59         if (f == zero) {  in logf()
67         R = f * f * ((float)0.5 - (float)0.33333333333333333 * f);  in logf()
69             return f - R;  in logf()
72             return dk * ln2_hi - ((R - dk * ln2_lo) - f);  in logf()
75     s = f / ((float)2.0 + f);  in logf()
86         hfsq = (float)0.5 * f * f;  in logf()
88             return f - (hfsq - s * (hfsq + R));  in logf()
90             return dk * ln2_hi - ((hfsq - (s * (hfsq + R) + dk * ln2_lo)) - f);  in logf()
93             return f - s * (f - R);  in logf()
95             return dk * ln2_hi - ((s * (f - R) - dk * ln2_lo) - f);  in logf()