Index: fusl/src/math/asinf.c |
diff --git a/fusl/src/math/asinf.c b/fusl/src/math/asinf.c |
index bcd304a3425991689cbf14dd1cd575a43a8b1117..3d5f637be79b9b9999dea7939bd2fcf3f04b3d66 100644 |
--- a/fusl/src/math/asinf.c |
+++ b/fusl/src/math/asinf.c |
@@ -14,48 +14,43 @@ |
*/ |
#include "libm.h" |
-static const double |
-pio2 = 1.570796326794896558e+00; |
+static const double pio2 = 1.570796326794896558e+00; |
static const float |
-/* coefficients for R(x^2) */ |
-pS0 = 1.6666586697e-01, |
-pS1 = -4.2743422091e-02, |
-pS2 = -8.6563630030e-03, |
-qS1 = -7.0662963390e-01; |
+ /* coefficients for R(x^2) */ |
+ pS0 = 1.6666586697e-01, |
+ pS1 = -4.2743422091e-02, pS2 = -8.6563630030e-03, qS1 = -7.0662963390e-01; |
-static float R(float z) |
-{ |
- float_t p, q; |
- p = z*(pS0+z*(pS1+z*pS2)); |
- q = 1.0f+z*qS1; |
- return p/q; |
+static float R(float z) { |
+ float_t p, q; |
+ p = z * (pS0 + z * (pS1 + z * pS2)); |
+ q = 1.0f + z * qS1; |
+ return p / q; |
} |
-float asinf(float x) |
-{ |
- double s; |
- float z; |
- uint32_t hx,ix; |
+float asinf(float x) { |
+ double s; |
+ float z; |
+ uint32_t hx, ix; |
- GET_FLOAT_WORD(hx, x); |
- ix = hx & 0x7fffffff; |
- if (ix >= 0x3f800000) { /* |x| >= 1 */ |
- if (ix == 0x3f800000) /* |x| == 1 */ |
- return x*pio2 + 0x1p-120f; /* asin(+-1) = +-pi/2 with inexact */ |
- return 0/(x-x); /* asin(|x|>1) is NaN */ |
- } |
- if (ix < 0x3f000000) { /* |x| < 0.5 */ |
- /* if 0x1p-126 <= |x| < 0x1p-12, avoid raising underflow */ |
- if (ix < 0x39800000 && ix >= 0x00800000) |
- return x; |
- return x + x*R(x*x); |
- } |
- /* 1 > |x| >= 0.5 */ |
- z = (1 - fabsf(x))*0.5f; |
- s = sqrt(z); |
- x = pio2 - 2*(s+s*R(z)); |
- if (hx >> 31) |
- return -x; |
- return x; |
+ GET_FLOAT_WORD(hx, x); |
+ ix = hx & 0x7fffffff; |
+ if (ix >= 0x3f800000) { /* |x| >= 1 */ |
+ if (ix == 0x3f800000) /* |x| == 1 */ |
+ return x * pio2 + 0x1p-120f; /* asin(+-1) = +-pi/2 with inexact */ |
+ return 0 / (x - x); /* asin(|x|>1) is NaN */ |
+ } |
+ if (ix < 0x3f000000) { /* |x| < 0.5 */ |
+ /* if 0x1p-126 <= |x| < 0x1p-12, avoid raising underflow */ |
+ if (ix < 0x39800000 && ix >= 0x00800000) |
+ return x; |
+ return x + x * R(x * x); |
+ } |
+ /* 1 > |x| >= 0.5 */ |
+ z = (1 - fabsf(x)) * 0.5f; |
+ s = sqrt(z); |
+ x = pio2 - 2 * (s + s * R(z)); |
+ if (hx >> 31) |
+ return -x; |
+ return x; |
} |