| 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
|
|
|