ieee754dp.c (2209bcb1310ffa9ee1af12573f1413581c712b15) ieee754dp.c (47fa0c0251413db66a9018fbac6f6266201195ae)
1/* IEEE754 floating point arithmetic
2 * double precision: common utilities
3 */
4/*
5 * MIPS floating point support
6 * Copyright (C) 1994-2000 Algorithmics Ltd.
7 *
8 * ########################################################################

--- 137 unchanged lines hidden (view full) ---

146 SETCX(IEEE754_UNDERFLOW);
147 SETCX(IEEE754_INEXACT);
148
149 switch(ieee754_csr.rm) {
150 case IEEE754_RN:
151 case IEEE754_RZ:
152 return ieee754dp_zero(sn);
153 case IEEE754_RU: /* toward +Infinity */
1/* IEEE754 floating point arithmetic
2 * double precision: common utilities
3 */
4/*
5 * MIPS floating point support
6 * Copyright (C) 1994-2000 Algorithmics Ltd.
7 *
8 * ########################################################################

--- 137 unchanged lines hidden (view full) ---

146 SETCX(IEEE754_UNDERFLOW);
147 SETCX(IEEE754_INEXACT);
148
149 switch(ieee754_csr.rm) {
150 case IEEE754_RN:
151 case IEEE754_RZ:
152 return ieee754dp_zero(sn);
153 case IEEE754_RU: /* toward +Infinity */
154 if(sn == 0)
154 if (sn == 0)
155 return ieee754dp_min(0);
156 else
157 return ieee754dp_zero(1);
158 case IEEE754_RD: /* toward -Infinity */
155 return ieee754dp_min(0);
156 else
157 return ieee754dp_zero(1);
158 case IEEE754_RD: /* toward -Infinity */
159 if(sn == 0)
159 if (sn == 0)
160 return ieee754dp_zero(0);
161 else
162 return ieee754dp_min(1);
163 }
164 }
165
166 if (xe == DP_EMIN - 1
167 && get_rounding(sn, xm) >> (DP_MBITS + 1 + 3))

--- 77 unchanged lines hidden ---
160 return ieee754dp_zero(0);
161 else
162 return ieee754dp_min(1);
163 }
164 }
165
166 if (xe == DP_EMIN - 1
167 && get_rounding(sn, xm) >> (DP_MBITS + 1 + 3))

--- 77 unchanged lines hidden ---