summaryrefslogtreecommitdiff
path: root/src/math/sincosl.c
diff options
context:
space:
mode:
authorSzabolcs Nagy <nsz@port70.net>2013-09-03 18:50:58 +0000
committerSzabolcs Nagy <nsz@port70.net>2013-09-05 11:30:08 +0000
commitea9bb95a5b36c0a3d2ed8fb03808745b406c2633 (patch)
tree5c0d395bfff168d97c7055826a30b822eba7293b /src/math/sincosl.c
parentbcd797a5ba4631c031919dad832d670e564212e9 (diff)
downloadmusl-ea9bb95a5b36c0a3d2ed8fb03808745b406c2633.tar.gz
musl-ea9bb95a5b36c0a3d2ed8fb03808745b406c2633.tar.bz2
musl-ea9bb95a5b36c0a3d2ed8fb03808745b406c2633.tar.xz
musl-ea9bb95a5b36c0a3d2ed8fb03808745b406c2633.zip
math: long double trigonometric cleanup (cosl, sinl, sincosl, tanl)
ld128 support was added to internal kernel functions (__cosl, __sinl, __tanl, __rem_pio2l) from freebsd (not tested, but should be a good start for when ld128 arch arrives) __rem_pio2l had some code cleanup, the freebsd ld128 code seems to gather the results of a large reduction with precision loss (fixed the bug but a todo comment was added for later investigation) the old copyright was removed from the non-kernel wrapper functions (cosl, sinl, sincosl, tanl) since these are trivial and the interesting parts and comments had been already rewritten.
Diffstat (limited to 'src/math/sincosl.c')
-rw-r--r--src/math/sincosl.c19
1 files changed, 6 insertions, 13 deletions
diff --git a/src/math/sincosl.c b/src/math/sincosl.c
index 5db69bd6..2c600801 100644
--- a/src/math/sincosl.c
+++ b/src/math/sincosl.c
@@ -9,25 +9,19 @@ void sincosl(long double x, long double *sin, long double *cos)
#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 IEEEl2bits u;
+ union ldshape u = {x};
unsigned n;
long double y[2], s, c;
- u.e = x;
- u.bits.sign = 0;
-
- /* x = nan or inf */
- if (u.bits.exp == 0x7fff) {
+ u.i.se &= 0x7fff;
+ if (u.i.se == 0x7fff) {
*sin = *cos = x - x;
return;
}
-
- /* |x| < (double)pi/4 */
- if (u.e < M_PI_4) {
- /* |x| < 0x1p-64 */
- if (u.bits.exp < 0x3fff - 64) {
+ if (u.f < M_PI_4) {
+ if (u.i.se < 0x3fff - LDBL_MANT_DIG) {
/* raise underflow if subnormal */
- if (u.bits.exp == 0) FORCE_EVAL(x*0x1p-120f);
+ if (u.i.se == 0) FORCE_EVAL(x*0x1p-120f);
*sin = x;
/* raise inexact if x!=0 */
*cos = 1.0 + x;
@@ -37,7 +31,6 @@ void sincosl(long double x, long double *sin, long double *cos)
*cos = __cosl(x, 0);
return;
}
-
n = __rem_pio2l(x, y);
s = __sinl(y[0], y[1], 1);
c = __cosl(y[0], y[1]);