math: raise invalid flag in ilogb*.c on +-0, +-inf and nan
authorSzabolcs Nagy <nsz@port70.net>
Mon, 12 Nov 2012 22:58:18 +0000 (23:58 +0100)
committerSzabolcs Nagy <nsz@port70.net>
Mon, 12 Nov 2012 22:58:18 +0000 (23:58 +0100)
src/math/ilogb.c
src/math/ilogbf.c
src/math/ilogbl.c

index 0a3a6a46352cdcf1eac8d4d7ac183e72d9fdd51a..5a1819d85f4c31e1de8302095586eddea260d21a 100644 (file)
@@ -8,13 +8,17 @@ int ilogb(double x)
 
        if (!e) {
                u.bits <<= 12;
-               if (u.bits == 0)
+               if (u.bits == 0) {
+                       FORCE_EVAL(0/0.0f);
                        return FP_ILOGB0;
+               }
                /* subnormal x */
                for (e = -0x3ff; u.bits < (uint64_t)1<<63; e--, u.bits<<=1);
                return e;
        }
-       if (e == 0x7ff)
+       if (e == 0x7ff) {
+               FORCE_EVAL(0/0.0f);
                return u.bits<<12 ? FP_ILOGBNAN : INT_MAX;
+       }
        return e - 0x3ff;
 }
index 272cbdac0ccf36dc2bb9d2ede7116853cad1748f..42cd62e2dd9b37598c716b73825090fd74eb588a 100644 (file)
@@ -8,13 +8,17 @@ int ilogbf(float x)
 
        if (!e) {
                u.bits <<= 9;
-               if (u.bits == 0)
+               if (u.bits == 0) {
+                       FORCE_EVAL(0/0.0f);
                        return FP_ILOGB0;
+               }
                /* subnormal x */
                for (e = -0x7f; u.bits < (uint32_t)1<<31; e--, u.bits<<=1);
                return e;
        }
-       if (e == 0xff)
+       if (e == 0xff) {
+               FORCE_EVAL(0/0.0f);
                return u.bits<<9 ? FP_ILOGBNAN : INT_MAX;
+       }
        return e - 0x7f;
 }
index ed9ddcbc42582ac0fb2c0ab5923a0bccb6ccd050..1512934f18c97df352ebb7106a22212bb7241246 100644 (file)
@@ -14,15 +14,19 @@ int ilogbl(long double x)
        int e = u.bits.exp;
 
        if (!e) {
-               if (m == 0)
+               if (m == 0) {
+                       FORCE_EVAL(0/0.0f);
                        return FP_ILOGB0;
+               }
                /* subnormal x */
                for (e = -0x3fff+1; m < (uint64_t)1<<63; e--, m<<=1);
                return e;
        }
-       if (e == 0x7fff)
+       if (e == 0x7fff) {
+               FORCE_EVAL(0/0.0f);
                /* in ld80 msb is set in inf */
                return m & (uint64_t)-1>>1 ? FP_ILOGBNAN : INT_MAX;
+       }
        return e - 0x3fff;
 }
 #endif