Index: fusl/src/math/acosh.c |
diff --git a/fusl/src/math/acosh.c b/fusl/src/math/acosh.c |
index badbf9081e644ff9df80bd2ae86cfc32e32a888b..a0299e5171deb204183b3f4d720470c62e77a2ae 100644 |
--- a/fusl/src/math/acosh.c |
+++ b/fusl/src/math/acosh.c |
@@ -1,24 +1,26 @@ |
#include "libm.h" |
-#if FLT_EVAL_METHOD==2 |
+#if FLT_EVAL_METHOD == 2 |
#undef sqrt |
#define sqrt sqrtl |
#endif |
/* acosh(x) = log(x + sqrt(x*x-1)) */ |
-double acosh(double x) |
-{ |
- union {double f; uint64_t i;} u = {.f = x}; |
- unsigned e = u.i >> 52 & 0x7ff; |
+double acosh(double x) { |
+ union { |
+ double f; |
+ uint64_t i; |
+ } u = {.f = x}; |
+ unsigned e = u.i >> 52 & 0x7ff; |
- /* x < 1 domain error is handled in the called functions */ |
+ /* x < 1 domain error is handled in the called functions */ |
- if (e < 0x3ff + 1) |
- /* |x| < 2, up to 2ulp error in [1,1.125] */ |
- return log1p(x-1 + sqrt((x-1)*(x-1)+2*(x-1))); |
- if (e < 0x3ff + 26) |
- /* |x| < 0x1p26 */ |
- return log(2*x - 1/(x+sqrt(x*x-1))); |
- /* |x| >= 0x1p26 or nan */ |
- return log(x) + 0.693147180559945309417232121458176568; |
+ if (e < 0x3ff + 1) |
+ /* |x| < 2, up to 2ulp error in [1,1.125] */ |
+ return log1p(x - 1 + sqrt((x - 1) * (x - 1) + 2 * (x - 1))); |
+ if (e < 0x3ff + 26) |
+ /* |x| < 0x1p26 */ |
+ return log(2 * x - 1 / (x + sqrt(x * x - 1))); |
+ /* |x| >= 0x1p26 or nan */ |
+ return log(x) + 0.693147180559945309417232121458176568; |
} |