X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Facoshf.c;h=0f7aae2abce9ddd299599c0019ccec876cbd81db;hb=fe0260400eebb613338a720c9568c10ab4f17225;hp=57ce5cb8297e0bb88caa375bb269be2ad5ee5cdb;hpb=93a50a26cd0f9efc59cc83daae7b2d916b327ab1;p=musl diff --git a/src/math/acoshf.c b/src/math/acoshf.c index 57ce5cb8..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(2.0f*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; + t = x-1.0f; return log1pf(t + sqrtf(2.0f*t+t*t)); } }