X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Flog10f.c;h=e10749b5bbeedb9f482f82961c58bcba7daf010a;hb=697acde67e0da4d73b46445ed536fe9923d515c7;hp=4175cce289de7304bfceb81cf4661b21259184d8;hpb=b69f695acedd4ce2798ef9ea28d834ceccc789bd;p=musl diff --git a/src/math/log10f.c b/src/math/log10f.c index 4175cce2..e10749b5 100644 --- a/src/math/log10f.c +++ b/src/math/log10f.c @@ -23,8 +23,6 @@ ivln10lo = -3.1689971365e-05, /* 0xb804ead9 */ log10_2hi = 3.0102920532e-01, /* 0x3e9a2080 */ log10_2lo = 7.9034151668e-07; /* 0x355427db */ -static const float zero = 0.0; - float log10f(float x) { float f,hfsq,hi,lo,r,y; @@ -35,9 +33,9 @@ float log10f(float x) k = 0; if (hx < 0x00800000) { /* x < 2**-126 */ if ((hx&0x7fffffff) == 0) - return -two25/zero; /* log(+-0)=-inf */ + return -two25/0.0f; /* log(+-0)=-inf */ if (hx < 0) - return (x-x)/zero; /* log(-#) = NaN */ + return (x-x)/0.0f; /* log(-#) = NaN */ /* subnormal number, scale up x */ k -= 25; x *= two25; @@ -46,15 +44,15 @@ float log10f(float x) if (hx >= 0x7f800000) return x+x; if (hx == 0x3f800000) - return zero; /* log(1) = +0 */ + return 0.0f; /* log(1) = +0 */ k += (hx>>23) - 127; hx &= 0x007fffff; i = (hx+(0x4afb0d))&0x800000; SET_FLOAT_WORD(x, hx|(i^0x3f800000)); /* normalize x or x/2 */ k += i>>23; y = (float)k; - f = x - (float)1.0; - hfsq = (float)0.5*f*f; + f = x - 1.0f; + hfsq = 0.5f * f * f; r = __log1pf(f); // FIXME