Changeset 3180

Timestamp:

02/14/08 03:20:22 (10 months ago)

Author:

Don Clugston

Message:

More fixes for quadruple reals. In particular, ieeeMean was horribly broken.

Files:

• trunk/tango/math/IEEE.d (modified) (19 diffs)

trunk/tango/math/IEEE.d

@@ -229 +229 @@
-            INVALID_MASK   = 0xF80
-        }
+    } else { //SPARC -- not yet implemented
+        enum : int {
+            INEXACT_MASK   = 0,
+            UNDERFLOW_MASK = 0,
+            OVERFLOW_MASK  = 0,
+            DIVBYZERO_MASK = 0,
+            INVALID_MASK   = 0
+        }
-    }
-private:
@@ -444 +452 @@
-        vu[EXPONENTPOS_SHORT] = cast(ushort)((0x8000 & vu[EXPONENTPOS_SHORT]) | 0x3FFE);
-    }
-
-
+
-        if (ex) { // If exponent is non-zero
-            if (ex == EXPMASK) {   // infinity or NaN
@@ -470 +477 @@
-            // multi-byte left shift
-            vl[1]<<=1;
-1
+l
-            *vl <<= 1;
-        } while ((vl[1]&0x0000_8000_0000_0000)== 0);
@@ -476 +483 @@
-        vu[EXPONENTPOS_SHORT] = cast(ushort)((0x8000 & vu[EXPONENTPOS_SHORT]) | 0x3FFE);
-    }
-
-
-
+
-    if (ex) { // If exponent is non-zero
-        if (ex == EXPMASK) {   // infinity or NaN
@@ -510 +515 @@
-        ve[EXPONENTPOS_SHORT] = sgn;
-    }
+  }else { //static if(real.mant_dig==106) // doubledouble
+        assert(0, "Unsupported");
-  }
-    return value;
@@ -737 +744 @@
-        }
-    } else {
-BUG
+NOTE
-        return tango.stdc.math.scalbnl(x, n);
-    }
@@ -903 +910 @@
-    ushort *p = cast(ushort *)&d;
-    version(LittleEndian) {
-in
+shor
-    } else {
-in
+shor
-    }
-    return e!=0 && e != 0x7FF0;
@@ -1009 +1016 @@
-    uint *p = cast(uint *)&f;
-
-0000) == 0 && *p & 0x007F
+_0000) == 0 && *p & 0x007F_
-}
-
@@ -1028 +1035 @@
-    uint *p = cast(uint *)&d;
-
-
+
+
+
+
+
-}
-
@@ -1077 +1088 @@
-        return ((*cast(ulong *)&x) & 0x7FFF_FFFF_FFFF_FFFF) == 0;
-    } else static if (real.mant_dig == 113) { // quadruple   
-e
-    return (ps[MANTISSA_LSB]==0 || (ps[MANTISSA_MSB]& 0x7FFF_FFFF_FFFF_FFFF) == 0)
+x
+return (ps[MANTISSA_LSB] | (ps[MANTISSA_MSB]& 0x7FFF_FFFF_FFFF_FFFF)) == 0
-    } else { // real80
-        ushort* pe = cast(ushort *)&x;
@@ -1469 +1480 @@
-    // they could have 0 or 1 bits in common.
-
-64
+80 or quadruple
-    int bitsdiff = ( ((pa[EXPONENTPOS_SHORT]&0x7FFF) + (pb[EXPONENTPOS_SHORT]&0x7FFF)-1)>>1) - pd[EXPONENTPOS_SHORT];
-
@@ -1486 +1497 @@
-    // Avoid out-by-1 errors when factor is almost 2.
-    return (bitsdiff == 0) ? (pa[EXPONENTPOS_SHORT] == pb[EXPONENTPOS_SHORT]) : 0;
-
-
-
+
-    int bitsdiff = (( ((pa[EXPONENTPOS_SHORT]&0x7FF0) + (pb[EXPONENTPOS_SHORT]&0x7FF0)-0x10)>>1) 
-                 - (pd[EXPONENTPOS_SHORT]&0x7FF0))>>4;
@@ -1507 +1516 @@
-    if (bitsdiff == 0 && !((pa[EXPONENTPOS_SHORT] ^ pb[EXPONENTPOS_SHORT])&EXPMASK)) return 1;
-    else return 0;
-
-
+
+
+
-}
-
@@ -1644 +1654 @@
-
-    T u;
-x87, 80-bit reals
+real80
-        // There's slight additional complexity because they are actually
-        // 79-bit reals...
@@ -1653 +1663 @@
-        ushort *ye = cast(ushort *)&y;
-        ulong *yl = cast(ulong *)&y;
-
+ (Bonus: this prevents overflows)
-        ulong m = ((*xl) & 0x7FFF_FFFF_FFFF_FFFF) + ((*yl) & 0x7FFF_FFFF_FFFF_FFFF);
-
@@ -1674 +1684 @@
-        ulong *xl = cast(ulong *)&x;
-        ulong *yl = cast(ulong *)&y;
-
-
-
+
+
+
+
+
+
+
+
+
+
-    } else static if (T.mant_dig == double.mant_dig) {
-        ulong *ul = cast(ulong *)&u;
@@ -1684 +1701 @@
-        m |= ((*xl) & 0x8000_0000_0000_0000);
-        *ul = m;
-
+ 
-        uint *ul = cast(uint *)&u;
-        uint *xl = cast(uint *)&x;
@@ -1691 +1708 @@
-        m |= ((*xl) & 0x8000_0000);
-        *ul = m;
+    } else {
+        assert(0, "Not implemented");
-    }
-    return u;
@@ -1732 +1751 @@
-real NaN(ulong payload)
-{
-
+
+
+
-      ulong v = 2; // no implied bit. quiet bit = 1
-    } else {
-      ulong v = 3; // implied bit = 1, quiet bit = 1
-    }
-