support ld80 pseudo-denormal invalid bit patterns; treat them as nan

[oweals/musl.git] / src / math / __tanl.c

diff --git a/src/math/__tanl.c b/src/math/__tanl.c
index e39e9df44b5b79cd9713cc49e9bf8b6dab6d8b09..50ba21409444d2735a612e3d4d065fbf535bb96a 100644 (file)
--- a/src/math/__tanl.c
+++ b/src/math/__tanl.c
@@ -51,8 +51,7 @@ long double __tanl(long double x, long double y, int iy) {
        int i;
 
        iy = iy == 1 ? -1 : 1;        /* XXX recover original interface */
-       // FIXME: this is wrong, use copysign, signbit or union bithack
-       osign = x >= 0 ? 1.0 : -1.0;  /* XXX slow, probably wrong for -0 */
+       osign = copysignl(1.0, x);
        if (fabsl(x) >= 0.67434) {
                if (x < 0) {
                        x = -x;