diff options
Diffstat (limited to 'libquadmath/math/acoshq.c')
| -rw-r--r-- | libquadmath/math/acoshq.c | 24 |
1 files changed, 12 insertions, 12 deletions
diff --git a/libquadmath/math/acoshq.c b/libquadmath/math/acoshq.c index 9845a8e..6e500a8 100644 --- a/libquadmath/math/acoshq.c +++ b/libquadmath/math/acoshq.c @@ -1,4 +1,4 @@ -/* acoshq.c -- __float128 version of e_acosh.c. +/* e_acoshl.c -- long double version of e_acosh.c. * Conversion to long double by Jakub Jelinek, jj@ultra.linux.cz. */ @@ -16,11 +16,11 @@ /* acoshq(x) * Method : * Based on - * acoshl(x) = logl [ x + sqrtl(x*x-1) ] + * acoshl(x) = logq [ x + sqrtq(x*x-1) ] * we have - * acoshl(x) := logl(x)+ln2, if x is large; else - * acoshl(x) := logl(2x-1/(sqrtl(x*x-1)+x)) if x>2; else - * acoshl(x) := log1pl(t+sqrtl(2.0*t+t*t)); where t=x-1. + * acoshl(x) := logq(x)+ln2, if x is large; else + * acoshl(x) := logq(2x-1/(sqrtq(x*x-1)+x)) if x>2; else + * acoshl(x) := log1pq(t+sqrtq(2.0*t+t*t)); where t=x-1. * * Special cases: * acoshl(x) is NaN with signal if x<1. @@ -30,11 +30,11 @@ #include "quadmath-imp.h" static const __float128 -one = 1.0Q, +one = 1.0, ln2 = 0.6931471805599453094172321214581766Q; __float128 -acoshq (__float128 x) +acoshq(__float128 x) { __float128 t; uint64_t lx; @@ -44,16 +44,16 @@ acoshq (__float128 x) return (x-x)/(x-x); } else if(hx >=0x4035000000000000LL) { /* x > 2**54 */ if(hx >=0x7fff000000000000LL) { /* x is inf of NaN */ - return x+x; + return x+x; } else - return logq(x)+ln2; /* acoshl(huge)=logl(2x) */ + return logq(x)+ln2; /* acoshl(huge)=logq(2x) */ } else if(((hx-0x3fff000000000000LL)|lx)==0) { - return 0.0Q; /* acosh(1) = 0 */ + return 0; /* acosh(1) = 0 */ } else if (hx > 0x4000000000000000LL) { /* 2**28 > x > 2 */ t=x*x; - return logq(2.0Q*x-one/(x+sqrtq(t-one))); + return logq(2*x-one/(x+sqrtq(t-one))); } else { /* 1<x<2 */ t = x-one; - return log1pq(t+sqrtq(2.0Q*t+t*t)); + return log1pq(t+sqrtq(2*t+t*t)); } } |