1 files changed, 165 insertions, 0 deletions

diff --git a/arch/mips/math-emu/ieee754int.h b/arch/mips/math-emu/ieee754int.h
new file mode 100644
index 000000000000..4a5a81d6b893
--- /dev/null
+++ b/arch/mips/math-emu/ieee754int.h
@@ -0,0 +1,165 @@
+/*
+ * IEEE754 floating point
+ * common internal header file
+ */
+/*
+ * MIPS floating point support
+ * Copyright (C) 1994-2000 Algorithmics Ltd.
+ * http://www.algor.co.uk
+ *
+ * ########################################################################
+ *
+ *  This program is free software; you can distribute it and/or modify it
+ *  under the terms of the GNU General Public License (Version 2) as
+ *  published by the Free Software Foundation.
+ *
+ *  This program is distributed in the hope it will be useful, but WITHOUT
+ *  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ *  FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+ *  for more details.
+ *
+ *  You should have received a copy of the GNU General Public License along
+ *  with this program; if not, write to the Free Software Foundation, Inc.,
+ *  59 Temple Place - Suite 330, Boston MA 02111-1307, USA.
+ *
+ * ########################################################################
+ */
+
+
+#include "ieee754.h"
+
+#define DP_EBIAS        1023
+#define DP_EMIN         (-1022)
+#define DP_EMAX         1023
+#define DP_MBITS        52
+
+#define SP_EBIAS        127
+#define SP_EMIN         (-126)
+#define SP_EMAX         127
+#define SP_MBITS        23
+
+#define DP_MBIT(x)      ((u64)1 << (x))
+#define DP_HIDDEN_BIT   DP_MBIT(DP_MBITS)
+#define DP_SIGN_BIT     DP_MBIT(63)
+
+#define SP_MBIT(x)      ((u32)1 << (x))
+#define SP_HIDDEN_BIT   SP_MBIT(SP_MBITS)
+#define SP_SIGN_BIT     SP_MBIT(31)
+
+
+#define SPSIGN(sp)      (sp.parts.sign)
+#define SPBEXP(sp)      (sp.parts.bexp)
+#define SPMANT(sp)      (sp.parts.mant)
+
+#define DPSIGN(dp)      (dp.parts.sign)
+#define DPBEXP(dp)      (dp.parts.bexp)
+#define DPMANT(dp)      (dp.parts.mant)
+
+#define CLPAIR(x,y)     ((x)*6+(y))
+
+#define CLEARCX \
+  (ieee754_csr.cx = 0)
+
+#define SETCX(x) \
+  (ieee754_csr.cx |= (x),ieee754_csr.sx |= (x))
+
+#define SETANDTESTCX(x) \
+  (SETCX(x),ieee754_csr.mx & (x))
+
+#define TSTX()  \
+        (ieee754_csr.cx & ieee754_csr.mx)
+
+
+#define COMPXSP \
+  unsigned xm; int xe; int xs; int xc
+
+#define COMPYSP \
+  unsigned ym; int ye; int ys; int yc
+
+#define EXPLODESP(v,vc,vs,ve,vm) \
+{\
+    vs = SPSIGN(v);\
+    ve = SPBEXP(v);\
+    vm = SPMANT(v);\
+    if(ve == SP_EMAX+1+SP_EBIAS){\
+        if(vm == 0)\
+          vc = IEEE754_CLASS_INF;\
+        else if(vm & SP_MBIT(SP_MBITS-1)) \
+          vc = IEEE754_CLASS_SNAN;\
+        else \
+          vc = IEEE754_CLASS_QNAN;\
+    } else if(ve == SP_EMIN-1+SP_EBIAS) {\
+        if(vm) {\
+            ve = SP_EMIN;\
+            vc = IEEE754_CLASS_DNORM;\
+        } else\
+          vc = IEEE754_CLASS_ZERO;\
+    } else {\
+        ve -= SP_EBIAS;\
+        vm |= SP_HIDDEN_BIT;\
+        vc = IEEE754_CLASS_NORM;\
+    }\
+}
+#define EXPLODEXSP EXPLODESP(x,xc,xs,xe,xm)
+#define EXPLODEYSP EXPLODESP(y,yc,ys,ye,ym)
+
+
+#define COMPXDP \
+u64 xm; int xe; int xs; int xc
+
+#define COMPYDP \
+u64 ym; int ye; int ys; int yc
+
+#define EXPLODEDP(v,vc,vs,ve,vm) \
+{\
+    vm = DPMANT(v);\
+    vs = DPSIGN(v);\
+    ve = DPBEXP(v);\
+    if(ve == DP_EMAX+1+DP_EBIAS){\
+        if(vm == 0)\
+          vc = IEEE754_CLASS_INF;\
+        else if(vm & DP_MBIT(DP_MBITS-1)) \
+          vc = IEEE754_CLASS_SNAN;\
+        else \
+          vc = IEEE754_CLASS_QNAN;\
+    } else if(ve == DP_EMIN-1+DP_EBIAS) {\
+        if(vm) {\
+            ve = DP_EMIN;\
+            vc = IEEE754_CLASS_DNORM;\
+        } else\
+          vc = IEEE754_CLASS_ZERO;\
+    } else {\
+        ve -= DP_EBIAS;\
+        vm |= DP_HIDDEN_BIT;\
+        vc = IEEE754_CLASS_NORM;\
+    }\
+}
+#define EXPLODEXDP EXPLODEDP(x,xc,xs,xe,xm)
+#define EXPLODEYDP EXPLODEDP(y,yc,ys,ye,ym)
+
+#define FLUSHDP(v,vc,vs,ve,vm) \
+        if(vc==IEEE754_CLASS_DNORM) {\
+            if(ieee754_csr.nod) {\
+                SETCX(IEEE754_INEXACT);\
+                vc = IEEE754_CLASS_ZERO;\
+                ve = DP_EMIN-1+DP_EBIAS;\
+                vm = 0;\
+                v = ieee754dp_zero(vs);\
+            }\
+        }
+
+#define FLUSHSP(v,vc,vs,ve,vm) \
+        if(vc==IEEE754_CLASS_DNORM) {\
+            if(ieee754_csr.nod) {\
+                SETCX(IEEE754_INEXACT);\
+                vc = IEEE754_CLASS_ZERO;\
+                ve = SP_EMIN-1+SP_EBIAS;\
+                vm = 0;\
+                v = ieee754sp_zero(vs);\
+            }\
+        }
+
+#define FLUSHXDP FLUSHDP(x,xc,xs,xe,xm)
+#define FLUSHYDP FLUSHDP(y,yc,ys,ye,ym)
+#define FLUSHXSP FLUSHSP(x,xc,xs,xe,xm)
+#define FLUSHYSP FLUSHSP(y,yc,ys,ye,ym)