diff options
Diffstat (limited to 'arch/mips/math-emu/dp_tlong.c')
-rw-r--r-- | arch/mips/math-emu/dp_tlong.c | 68 |
1 files changed, 26 insertions, 42 deletions
diff --git a/arch/mips/math-emu/dp_tlong.c b/arch/mips/math-emu/dp_tlong.c index 133ce2ba0012..9cdc145b75e0 100644 --- a/arch/mips/math-emu/dp_tlong.c +++ b/arch/mips/math-emu/dp_tlong.c | |||
@@ -5,8 +5,6 @@ | |||
5 | * MIPS floating point support | 5 | * MIPS floating point support |
6 | * Copyright (C) 1994-2000 Algorithmics Ltd. | 6 | * Copyright (C) 1994-2000 Algorithmics Ltd. |
7 | * | 7 | * |
8 | * ######################################################################## | ||
9 | * | ||
10 | * This program is free software; you can distribute it and/or modify it | 8 | * This program is free software; you can distribute it and/or modify it |
11 | * under the terms of the GNU General Public License (Version 2) as | 9 | * under the terms of the GNU General Public License (Version 2) as |
12 | * published by the Free Software Foundation. | 10 | * published by the Free Software Foundation. |
@@ -18,19 +16,21 @@ | |||
18 | * | 16 | * |
19 | * You should have received a copy of the GNU General Public License along | 17 | * You should have received a copy of the GNU General Public License along |
20 | * with this program; if not, write to the Free Software Foundation, Inc., | 18 | * with this program; if not, write to the Free Software Foundation, Inc., |
21 | * 59 Temple Place - Suite 330, Boston MA 02111-1307, USA. | 19 | * 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA. |
22 | * | ||
23 | * ######################################################################## | ||
24 | */ | 20 | */ |
25 | 21 | ||
26 | |||
27 | #include "ieee754dp.h" | 22 | #include "ieee754dp.h" |
28 | 23 | ||
29 | s64 ieee754dp_tlong(ieee754dp x) | 24 | s64 ieee754dp_tlong(union ieee754dp x) |
30 | { | 25 | { |
26 | u64 residue; | ||
27 | int round; | ||
28 | int sticky; | ||
29 | int odd; | ||
30 | |||
31 | COMPXDP; | 31 | COMPXDP; |
32 | 32 | ||
33 | CLEARCX; | 33 | ieee754_clearcx(); |
34 | 34 | ||
35 | EXPLODEXDP; | 35 | EXPLODEXDP; |
36 | FLUSHXDP; | 36 | FLUSHXDP; |
@@ -39,10 +39,12 @@ s64 ieee754dp_tlong(ieee754dp x) | |||
39 | case IEEE754_CLASS_SNAN: | 39 | case IEEE754_CLASS_SNAN: |
40 | case IEEE754_CLASS_QNAN: | 40 | case IEEE754_CLASS_QNAN: |
41 | case IEEE754_CLASS_INF: | 41 | case IEEE754_CLASS_INF: |
42 | SETCX(IEEE754_INVALID_OPERATION); | 42 | ieee754_setcx(IEEE754_INVALID_OPERATION); |
43 | return ieee754di_xcpt(ieee754di_indef(), "dp_tlong", x); | 43 | return ieee754di_indef(); |
44 | |||
44 | case IEEE754_CLASS_ZERO: | 45 | case IEEE754_CLASS_ZERO: |
45 | return 0; | 46 | return 0; |
47 | |||
46 | case IEEE754_CLASS_DNORM: | 48 | case IEEE754_CLASS_DNORM: |
47 | case IEEE754_CLASS_NORM: | 49 | case IEEE754_CLASS_NORM: |
48 | break; | 50 | break; |
@@ -53,18 +55,13 @@ s64 ieee754dp_tlong(ieee754dp x) | |||
53 | return -0x8000000000000000LL; | 55 | return -0x8000000000000000LL; |
54 | /* Set invalid. We will only use overflow for floating | 56 | /* Set invalid. We will only use overflow for floating |
55 | point overflow */ | 57 | point overflow */ |
56 | SETCX(IEEE754_INVALID_OPERATION); | 58 | ieee754_setcx(IEEE754_INVALID_OPERATION); |
57 | return ieee754di_xcpt(ieee754di_indef(), "dp_tlong", x); | 59 | return ieee754di_indef(); |
58 | } | 60 | } |
59 | /* oh gawd */ | 61 | /* oh gawd */ |
60 | if (xe > DP_MBITS) { | 62 | if (xe > DP_FBITS) { |
61 | xm <<= xe - DP_MBITS; | 63 | xm <<= xe - DP_FBITS; |
62 | } else if (xe < DP_MBITS) { | 64 | } else if (xe < DP_FBITS) { |
63 | u64 residue; | ||
64 | int round; | ||
65 | int sticky; | ||
66 | int odd; | ||
67 | |||
68 | if (xe < -1) { | 65 | if (xe < -1) { |
69 | residue = xm; | 66 | residue = xm; |
70 | round = 0; | 67 | round = 0; |
@@ -75,51 +72,38 @@ s64 ieee754dp_tlong(ieee754dp x) | |||
75 | * so we do it in two steps. Be aware that xe | 72 | * so we do it in two steps. Be aware that xe |
76 | * may be -1 */ | 73 | * may be -1 */ |
77 | residue = xm << (xe + 1); | 74 | residue = xm << (xe + 1); |
78 | residue <<= 63 - DP_MBITS; | 75 | residue <<= 63 - DP_FBITS; |
79 | round = (residue >> 63) != 0; | 76 | round = (residue >> 63) != 0; |
80 | sticky = (residue << 1) != 0; | 77 | sticky = (residue << 1) != 0; |
81 | xm >>= DP_MBITS - xe; | 78 | xm >>= DP_FBITS - xe; |
82 | } | 79 | } |
83 | odd = (xm & 0x1) != 0x0; | 80 | odd = (xm & 0x1) != 0x0; |
84 | switch (ieee754_csr.rm) { | 81 | switch (ieee754_csr.rm) { |
85 | case IEEE754_RN: | 82 | case FPU_CSR_RN: |
86 | if (round && (sticky || odd)) | 83 | if (round && (sticky || odd)) |
87 | xm++; | 84 | xm++; |
88 | break; | 85 | break; |
89 | case IEEE754_RZ: | 86 | case FPU_CSR_RZ: |
90 | break; | 87 | break; |
91 | case IEEE754_RU: /* toward +Infinity */ | 88 | case FPU_CSR_RU: /* toward +Infinity */ |
92 | if ((round || sticky) && !xs) | 89 | if ((round || sticky) && !xs) |
93 | xm++; | 90 | xm++; |
94 | break; | 91 | break; |
95 | case IEEE754_RD: /* toward -Infinity */ | 92 | case FPU_CSR_RD: /* toward -Infinity */ |
96 | if ((round || sticky) && xs) | 93 | if ((round || sticky) && xs) |
97 | xm++; | 94 | xm++; |
98 | break; | 95 | break; |
99 | } | 96 | } |
100 | if ((xm >> 63) != 0) { | 97 | if ((xm >> 63) != 0) { |
101 | /* This can happen after rounding */ | 98 | /* This can happen after rounding */ |
102 | SETCX(IEEE754_INVALID_OPERATION); | 99 | ieee754_setcx(IEEE754_INVALID_OPERATION); |
103 | return ieee754di_xcpt(ieee754di_indef(), "dp_tlong", x); | 100 | return ieee754di_indef(); |
104 | } | 101 | } |
105 | if (round || sticky) | 102 | if (round || sticky) |
106 | SETCX(IEEE754_INEXACT); | 103 | ieee754_setcx(IEEE754_INEXACT); |
107 | } | 104 | } |
108 | if (xs) | 105 | if (xs) |
109 | return -xm; | 106 | return -xm; |
110 | else | 107 | else |
111 | return xm; | 108 | return xm; |
112 | } | 109 | } |
113 | |||
114 | |||
115 | u64 ieee754dp_tulong(ieee754dp x) | ||
116 | { | ||
117 | ieee754dp hb = ieee754dp_1e63(); | ||
118 | |||
119 | /* what if x < 0 ?? */ | ||
120 | if (ieee754dp_lt(x, hb)) | ||
121 | return (u64) ieee754dp_tlong(x); | ||
122 | |||
123 | return (u64) ieee754dp_tlong(ieee754dp_sub(x, hb)) | | ||
124 | (1ULL << 63); | ||
125 | } | ||