X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Facoshf.c;h=0f7aae2abce9ddd299599c0019ccec876cbd81db;hb=9f58d06007818c43b9504e959ef81ab5f113b374;hp=30a3a943142fcbcfc29cc7315ab2a890072c6460;hpb=b69f695acedd4ce2798ef9ea28d834ceccc789bd;p=musl diff --git a/src/math/acoshf.c b/src/math/acoshf.c index 30a3a943..0f7aae2a 100644 --- a/src/math/acoshf.c +++ b/src/math/acoshf.c @@ -16,7 +16,6 @@ #include "libm.h" static const float -one = 1.0, ln2 = 6.9314718246e-01; /* 0x3f317218 */ float acoshf(float x) @@ -32,12 +31,12 @@ float acoshf(float x) return x + x; return logf(x) + ln2; /* acosh(huge)=log(2x) */ } else if (hx == 0x3f800000) { - return 0.0; /* acosh(1) = 0 */ + return 0.0f; /* acosh(1) = 0 */ } else if (hx > 0x40000000) { /* 2**28 > x > 2 */ t = x*x; - return logf((float)2.0*x - one/(x+sqrtf(t-one))); + return logf(2.0f*x - 1.0f/(x+sqrtf(t-1.0f))); } else { /* 1 < x < 2 */ - t = x-one; - return log1pf(t + sqrtf((float)2.0*t+t*t)); + t = x-1.0f; + return log1pf(t + sqrtf(2.0f*t+t*t)); } }