Index: fusl/src/math/__cosl.c |
diff --git a/fusl/src/math/__cosl.c b/fusl/src/math/__cosl.c |
index fa522ddd7ae263b02a0dfb00b04486aa4ca90af5..b9aee53f13778374cb1d6e5436d3baf04e024d42 100644 |
--- a/fusl/src/math/__cosl.c |
+++ b/fusl/src/math/__cosl.c |
@@ -12,7 +12,6 @@ |
* ==================================================== |
*/ |
- |
#include "libm.h" |
#if (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384 |
@@ -43,16 +42,16 @@ |
* almost for free from the complications needed to search for the best |
* higher coefficients. |
*/ |
-static const long double |
-C1 = 0.0416666666666666666136L; /* 0xaaaaaaaaaaaaaa9b.0p-68 */ |
-static const double |
-C2 = -0.0013888888888888874, /* -0x16c16c16c16c10.0p-62 */ |
-C3 = 0.000024801587301571716, /* 0x1a01a01a018e22.0p-68 */ |
-C4 = -0.00000027557319215507120, /* -0x127e4fb7602f22.0p-74 */ |
-C5 = 0.0000000020876754400407278, /* 0x11eed8caaeccf1.0p-81 */ |
-C6 = -1.1470297442401303e-11, /* -0x19393412bd1529.0p-89 */ |
-C7 = 4.7383039476436467e-14; /* 0x1aac9d9af5c43e.0p-97 */ |
-#define POLY(z) (z*(C1+z*(C2+z*(C3+z*(C4+z*(C5+z*(C6+z*C7))))))) |
+static const long double C1 = |
+ 0.0416666666666666666136L; /* 0xaaaaaaaaaaaaaa9b.0p-68 */ |
+static const double C2 = -0.0013888888888888874, /* -0x16c16c16c16c10.0p-62 */ |
+ C3 = 0.000024801587301571716, /* 0x1a01a01a018e22.0p-68 */ |
+ C4 = -0.00000027557319215507120, /* -0x127e4fb7602f22.0p-74 */ |
+ C5 = 0.0000000020876754400407278, /* 0x11eed8caaeccf1.0p-81 */ |
+ C6 = -1.1470297442401303e-11, /* -0x19393412bd1529.0p-89 */ |
+ C7 = 4.7383039476436467e-14; /* 0x1aac9d9af5c43e.0p-97 */ |
+#define POLY(z) \ |
+ (z * (C1 + z * (C2 + z * (C3 + z * (C4 + z * (C5 + z * (C6 + z * C7))))))) |
#elif LDBL_MANT_DIG == 113 |
/* |
* ld128 version of __cos.c. See __cos.c for most comments. |
@@ -66,31 +65,37 @@ C7 = 4.7383039476436467e-14; /* 0x1aac9d9af5c43e.0p-97 */ |
* that is 1 ulp below 0.5, but we want it to be precisely 0.5. See |
* above for more details. |
*/ |
-static const long double |
-C1 = 0.04166666666666666666666666666666658424671L, |
-C2 = -0.001388888888888888888888888888863490893732L, |
-C3 = 0.00002480158730158730158730158600795304914210L, |
-C4 = -0.2755731922398589065255474947078934284324e-6L, |
-C5 = 0.2087675698786809897659225313136400793948e-8L, |
-C6 = -0.1147074559772972315817149986812031204775e-10L, |
-C7 = 0.4779477332386808976875457937252120293400e-13L; |
-static const double |
-C8 = -0.1561920696721507929516718307820958119868e-15, |
-C9 = 0.4110317413744594971475941557607804508039e-18, |
-C10 = -0.8896592467191938803288521958313920156409e-21, |
-C11 = 0.1601061435794535138244346256065192782581e-23; |
-#define POLY(z) (z*(C1+z*(C2+z*(C3+z*(C4+z*(C5+z*(C6+z*(C7+ \ |
- z*(C8+z*(C9+z*(C10+z*C11))))))))))) |
+static const long double C1 = 0.04166666666666666666666666666666658424671L, |
+ C2 = -0.001388888888888888888888888888863490893732L, |
+ C3 = 0.00002480158730158730158730158600795304914210L, |
+ C4 = -0.2755731922398589065255474947078934284324e-6L, |
+ C5 = 0.2087675698786809897659225313136400793948e-8L, |
+ C6 = -0.1147074559772972315817149986812031204775e-10L, |
+ C7 = 0.4779477332386808976875457937252120293400e-13L; |
+static const double C8 = -0.1561920696721507929516718307820958119868e-15, |
+ C9 = 0.4110317413744594971475941557607804508039e-18, |
+ C10 = -0.8896592467191938803288521958313920156409e-21, |
+ C11 = 0.1601061435794535138244346256065192782581e-23; |
+#define POLY(z) \ |
+ (z * \ |
+ (C1 + \ |
+ z * (C2 + \ |
+ z * (C3 + \ |
+ z * (C4 + \ |
+ z * (C5 + \ |
+ z * (C6 + \ |
+ z * (C7 + \ |
+ z * (C8 + \ |
+ z * (C9 + z * (C10 + z * C11))))))))))) |
#endif |
-long double __cosl(long double x, long double y) |
-{ |
- long double hz,z,r,w; |
+long double __cosl(long double x, long double y) { |
+ long double hz, z, r, w; |
- z = x*x; |
- r = POLY(z); |
- hz = 0.5*z; |
- w = 1.0-hz; |
- return w + (((1.0-w)-hz) + (z*r-x*y)); |
+ z = x * x; |
+ r = POLY(z); |
+ hz = 0.5 * z; |
+ w = 1.0 - hz; |
+ return w + (((1.0 - w) - hz) + (z * r - x * y)); |
} |
#endif |