Index: fusl/src/math/sincosl.c |
diff --git a/fusl/src/math/sincosl.c b/fusl/src/math/sincosl.c |
index d3ac1c4c8c29e11067d55a3a99fa82eb888c5dd5..a4f6c2823547607ba142cdef500b1aa94345403a 100644 |
--- a/fusl/src/math/sincosl.c |
+++ b/fusl/src/math/sincosl.c |
@@ -2,59 +2,58 @@ |
#include "libm.h" |
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 |
-void sincosl(long double x, long double *sin, long double *cos) |
-{ |
- double sind, cosd; |
- sincos(x, &sind, &cosd); |
- *sin = sind; |
- *cos = cosd; |
+void sincosl(long double x, long double* sin, long double* cos) { |
+ double sind, cosd; |
+ sincos(x, &sind, &cosd); |
+ *sin = sind; |
+ *cos = cosd; |
} |
#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384 |
-void sincosl(long double x, long double *sin, long double *cos) |
-{ |
- union ldshape u = {x}; |
- unsigned n; |
- long double y[2], s, c; |
+void sincosl(long double x, long double* sin, long double* cos) { |
+ union ldshape u = {x}; |
+ unsigned n; |
+ long double y[2], s, c; |
- u.i.se &= 0x7fff; |
- if (u.i.se == 0x7fff) { |
- *sin = *cos = x - x; |
- return; |
- } |
- if (u.f < M_PI_4) { |
- if (u.i.se < 0x3fff - LDBL_MANT_DIG) { |
- /* raise underflow if subnormal */ |
- if (u.i.se == 0) FORCE_EVAL(x*0x1p-120f); |
- *sin = x; |
- /* raise inexact if x!=0 */ |
- *cos = 1.0 + x; |
- return; |
- } |
- *sin = __sinl(x, 0, 0); |
- *cos = __cosl(x, 0); |
- return; |
- } |
- n = __rem_pio2l(x, y); |
- s = __sinl(y[0], y[1], 1); |
- c = __cosl(y[0], y[1]); |
- switch (n & 3) { |
- case 0: |
- *sin = s; |
- *cos = c; |
- break; |
- case 1: |
- *sin = c; |
- *cos = -s; |
- break; |
- case 2: |
- *sin = -s; |
- *cos = -c; |
- break; |
- case 3: |
- default: |
- *sin = -c; |
- *cos = s; |
- break; |
- } |
+ u.i.se &= 0x7fff; |
+ if (u.i.se == 0x7fff) { |
+ *sin = *cos = x - x; |
+ return; |
+ } |
+ if (u.f < M_PI_4) { |
+ if (u.i.se < 0x3fff - LDBL_MANT_DIG) { |
+ /* raise underflow if subnormal */ |
+ if (u.i.se == 0) |
+ FORCE_EVAL(x * 0x1p-120f); |
+ *sin = x; |
+ /* raise inexact if x!=0 */ |
+ *cos = 1.0 + x; |
+ return; |
+ } |
+ *sin = __sinl(x, 0, 0); |
+ *cos = __cosl(x, 0); |
+ return; |
+ } |
+ n = __rem_pio2l(x, y); |
+ s = __sinl(y[0], y[1], 1); |
+ c = __cosl(y[0], y[1]); |
+ switch (n & 3) { |
+ case 0: |
+ *sin = s; |
+ *cos = c; |
+ break; |
+ case 1: |
+ *sin = c; |
+ *cos = -s; |
+ break; |
+ case 2: |
+ *sin = -s; |
+ *cos = -c; |
+ break; |
+ case 3: |
+ default: |
+ *sin = -c; |
+ *cos = s; |
+ break; |
+ } |
} |
#endif |