Index: fusl/src/math/sinhf.c |
diff --git a/fusl/src/math/sinhf.c b/fusl/src/math/sinhf.c |
index 6ad19ea2b0c8d2d639d49feef32106a80b7f0235..705ba77a9b507ae82fc5a73a8287442cd9b6bacc 100644 |
--- a/fusl/src/math/sinhf.c |
+++ b/fusl/src/math/sinhf.c |
@@ -1,31 +1,33 @@ |
#include "libm.h" |
-float sinhf(float x) |
-{ |
- union {float f; uint32_t i;} u = {.f = x}; |
- uint32_t w; |
- float t, h, absx; |
+float sinhf(float x) { |
+ union { |
+ float f; |
+ uint32_t i; |
+ } u = {.f = x}; |
+ uint32_t w; |
+ float t, h, absx; |
- h = 0.5; |
- if (u.i >> 31) |
- h = -h; |
- /* |x| */ |
- u.i &= 0x7fffffff; |
- absx = u.f; |
- w = u.i; |
+ h = 0.5; |
+ if (u.i >> 31) |
+ h = -h; |
+ /* |x| */ |
+ u.i &= 0x7fffffff; |
+ absx = u.f; |
+ w = u.i; |
- /* |x| < log(FLT_MAX) */ |
- if (w < 0x42b17217) { |
- t = expm1f(absx); |
- if (w < 0x3f800000) { |
- if (w < 0x3f800000 - (12<<23)) |
- return x; |
- return h*(2*t - t*t/(t+1)); |
- } |
- return h*(t + t/(t+1)); |
- } |
+ /* |x| < log(FLT_MAX) */ |
+ if (w < 0x42b17217) { |
+ t = expm1f(absx); |
+ if (w < 0x3f800000) { |
+ if (w < 0x3f800000 - (12 << 23)) |
+ return x; |
+ return h * (2 * t - t * t / (t + 1)); |
+ } |
+ return h * (t + t / (t + 1)); |
+ } |
- /* |x| > logf(FLT_MAX) or nan */ |
- t = 2*h*__expo2f(absx); |
- return t; |
+ /* |x| > logf(FLT_MAX) or nan */ |
+ t = 2 * h * __expo2f(absx); |
+ return t; |
} |