sfdiv.c revision 1.1.2.2
1 1.1.2.2 jdolecek /* $NetBSD: sfdiv.c,v 1.1.2.2 2002/06/23 17:37:15 jdolecek Exp $ */ 2 1.1.2.2 jdolecek 3 1.1.2.2 jdolecek /* $OpenBSD: sfdiv.c,v 1.4 2001/03/29 03:58:19 mickey Exp $ */ 4 1.1.2.2 jdolecek 5 1.1.2.2 jdolecek /* 6 1.1.2.2 jdolecek * Copyright 1996 1995 by Open Software Foundation, Inc. 7 1.1.2.2 jdolecek * All Rights Reserved 8 1.1.2.2 jdolecek * 9 1.1.2.2 jdolecek * Permission to use, copy, modify, and distribute this software and 10 1.1.2.2 jdolecek * its documentation for any purpose and without fee is hereby granted, 11 1.1.2.2 jdolecek * provided that the above copyright notice appears in all copies and 12 1.1.2.2 jdolecek * that both the copyright notice and this permission notice appear in 13 1.1.2.2 jdolecek * supporting documentation. 14 1.1.2.2 jdolecek * 15 1.1.2.2 jdolecek * OSF DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE 16 1.1.2.2 jdolecek * INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS 17 1.1.2.2 jdolecek * FOR A PARTICULAR PURPOSE. 18 1.1.2.2 jdolecek * 19 1.1.2.2 jdolecek * IN NO EVENT SHALL OSF BE LIABLE FOR ANY SPECIAL, INDIRECT, OR 20 1.1.2.2 jdolecek * CONSEQUENTIAL DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM 21 1.1.2.2 jdolecek * LOSS OF USE, DATA OR PROFITS, WHETHER IN ACTION OF CONTRACT, 22 1.1.2.2 jdolecek * NEGLIGENCE, OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION 23 1.1.2.2 jdolecek * WITH THE USE OR PERFORMANCE OF THIS SOFTWARE. 24 1.1.2.2 jdolecek * 25 1.1.2.2 jdolecek */ 26 1.1.2.2 jdolecek /* 27 1.1.2.2 jdolecek * pmk1.1 28 1.1.2.2 jdolecek */ 29 1.1.2.2 jdolecek /* 30 1.1.2.2 jdolecek * (c) Copyright 1986 HEWLETT-PACKARD COMPANY 31 1.1.2.2 jdolecek * 32 1.1.2.2 jdolecek * To anyone who acknowledges that this file is provided "AS IS" 33 1.1.2.2 jdolecek * without any express or implied warranty: 34 1.1.2.2 jdolecek * permission to use, copy, modify, and distribute this file 35 1.1.2.2 jdolecek * for any purpose is hereby granted without fee, provided that 36 1.1.2.2 jdolecek * the above copyright notice and this notice appears in all 37 1.1.2.2 jdolecek * copies, and that the name of Hewlett-Packard Company not be 38 1.1.2.2 jdolecek * used in advertising or publicity pertaining to distribution 39 1.1.2.2 jdolecek * of the software without specific, written prior permission. 40 1.1.2.2 jdolecek * Hewlett-Packard Company makes no representations about the 41 1.1.2.2 jdolecek * suitability of this software for any purpose. 42 1.1.2.2 jdolecek */ 43 1.1.2.2 jdolecek 44 1.1.2.2 jdolecek #include "../spmath/float.h" 45 1.1.2.2 jdolecek #include "../spmath/sgl_float.h" 46 1.1.2.2 jdolecek 47 1.1.2.2 jdolecek /* 48 1.1.2.2 jdolecek * Single Precision Floating-point Divide 49 1.1.2.2 jdolecek */ 50 1.1.2.2 jdolecek int 51 1.1.2.2 jdolecek sgl_fdiv(srcptr1,srcptr2,dstptr,status) 52 1.1.2.2 jdolecek 53 1.1.2.2 jdolecek sgl_floating_point *srcptr1, *srcptr2, *dstptr; 54 1.1.2.2 jdolecek unsigned int *status; 55 1.1.2.2 jdolecek { 56 1.1.2.2 jdolecek register unsigned int opnd1, opnd2, opnd3, result; 57 1.1.2.2 jdolecek register int dest_exponent, count; 58 1.1.2.2 jdolecek register int inexact = FALSE, guardbit = FALSE, stickybit = FALSE; 59 1.1.2.2 jdolecek int is_tiny; 60 1.1.2.2 jdolecek 61 1.1.2.2 jdolecek opnd1 = *srcptr1; 62 1.1.2.2 jdolecek opnd2 = *srcptr2; 63 1.1.2.2 jdolecek /* 64 1.1.2.2 jdolecek * set sign bit of result 65 1.1.2.2 jdolecek */ 66 1.1.2.2 jdolecek if (Sgl_sign(opnd1) ^ Sgl_sign(opnd2)) Sgl_setnegativezero(result); 67 1.1.2.2 jdolecek else Sgl_setzero(result); 68 1.1.2.2 jdolecek /* 69 1.1.2.2 jdolecek * check first operand for NaN's or infinity 70 1.1.2.2 jdolecek */ 71 1.1.2.2 jdolecek if (Sgl_isinfinity_exponent(opnd1)) { 72 1.1.2.2 jdolecek if (Sgl_iszero_mantissa(opnd1)) { 73 1.1.2.2 jdolecek if (Sgl_isnotnan(opnd2)) { 74 1.1.2.2 jdolecek if (Sgl_isinfinity(opnd2)) { 75 1.1.2.2 jdolecek /* 76 1.1.2.2 jdolecek * invalid since both operands 77 1.1.2.2 jdolecek * are infinity 78 1.1.2.2 jdolecek */ 79 1.1.2.2 jdolecek if (Is_invalidtrap_enabled()) 80 1.1.2.2 jdolecek return(INVALIDEXCEPTION); 81 1.1.2.2 jdolecek Set_invalidflag(); 82 1.1.2.2 jdolecek Sgl_makequietnan(result); 83 1.1.2.2 jdolecek *dstptr = result; 84 1.1.2.2 jdolecek return(NOEXCEPTION); 85 1.1.2.2 jdolecek } 86 1.1.2.2 jdolecek /* 87 1.1.2.2 jdolecek * return infinity 88 1.1.2.2 jdolecek */ 89 1.1.2.2 jdolecek Sgl_setinfinity_exponentmantissa(result); 90 1.1.2.2 jdolecek *dstptr = result; 91 1.1.2.2 jdolecek return(NOEXCEPTION); 92 1.1.2.2 jdolecek } 93 1.1.2.2 jdolecek } 94 1.1.2.2 jdolecek else { 95 1.1.2.2 jdolecek /* 96 1.1.2.2 jdolecek * is NaN; signaling or quiet? 97 1.1.2.2 jdolecek */ 98 1.1.2.2 jdolecek if (Sgl_isone_signaling(opnd1)) { 99 1.1.2.2 jdolecek /* trap if INVALIDTRAP enabled */ 100 1.1.2.2 jdolecek if (Is_invalidtrap_enabled()) 101 1.1.2.2 jdolecek return(INVALIDEXCEPTION); 102 1.1.2.2 jdolecek /* make NaN quiet */ 103 1.1.2.2 jdolecek Set_invalidflag(); 104 1.1.2.2 jdolecek Sgl_set_quiet(opnd1); 105 1.1.2.2 jdolecek } 106 1.1.2.2 jdolecek /* 107 1.1.2.2 jdolecek * is second operand a signaling NaN? 108 1.1.2.2 jdolecek */ 109 1.1.2.2 jdolecek else if (Sgl_is_signalingnan(opnd2)) { 110 1.1.2.2 jdolecek /* trap if INVALIDTRAP enabled */ 111 1.1.2.2 jdolecek if (Is_invalidtrap_enabled()) 112 1.1.2.2 jdolecek return(INVALIDEXCEPTION); 113 1.1.2.2 jdolecek /* make NaN quiet */ 114 1.1.2.2 jdolecek Set_invalidflag(); 115 1.1.2.2 jdolecek Sgl_set_quiet(opnd2); 116 1.1.2.2 jdolecek *dstptr = opnd2; 117 1.1.2.2 jdolecek return(NOEXCEPTION); 118 1.1.2.2 jdolecek } 119 1.1.2.2 jdolecek /* 120 1.1.2.2 jdolecek * return quiet NaN 121 1.1.2.2 jdolecek */ 122 1.1.2.2 jdolecek *dstptr = opnd1; 123 1.1.2.2 jdolecek return(NOEXCEPTION); 124 1.1.2.2 jdolecek } 125 1.1.2.2 jdolecek } 126 1.1.2.2 jdolecek /* 127 1.1.2.2 jdolecek * check second operand for NaN's or infinity 128 1.1.2.2 jdolecek */ 129 1.1.2.2 jdolecek if (Sgl_isinfinity_exponent(opnd2)) { 130 1.1.2.2 jdolecek if (Sgl_iszero_mantissa(opnd2)) { 131 1.1.2.2 jdolecek /* 132 1.1.2.2 jdolecek * return zero 133 1.1.2.2 jdolecek */ 134 1.1.2.2 jdolecek Sgl_setzero_exponentmantissa(result); 135 1.1.2.2 jdolecek *dstptr = result; 136 1.1.2.2 jdolecek return(NOEXCEPTION); 137 1.1.2.2 jdolecek } 138 1.1.2.2 jdolecek /* 139 1.1.2.2 jdolecek * is NaN; signaling or quiet? 140 1.1.2.2 jdolecek */ 141 1.1.2.2 jdolecek if (Sgl_isone_signaling(opnd2)) { 142 1.1.2.2 jdolecek /* trap if INVALIDTRAP enabled */ 143 1.1.2.2 jdolecek if (Is_invalidtrap_enabled()) return(INVALIDEXCEPTION); 144 1.1.2.2 jdolecek /* make NaN quiet */ 145 1.1.2.2 jdolecek Set_invalidflag(); 146 1.1.2.2 jdolecek Sgl_set_quiet(opnd2); 147 1.1.2.2 jdolecek } 148 1.1.2.2 jdolecek /* 149 1.1.2.2 jdolecek * return quiet NaN 150 1.1.2.2 jdolecek */ 151 1.1.2.2 jdolecek *dstptr = opnd2; 152 1.1.2.2 jdolecek return(NOEXCEPTION); 153 1.1.2.2 jdolecek } 154 1.1.2.2 jdolecek /* 155 1.1.2.2 jdolecek * check for division by zero 156 1.1.2.2 jdolecek */ 157 1.1.2.2 jdolecek if (Sgl_iszero_exponentmantissa(opnd2)) { 158 1.1.2.2 jdolecek if (Sgl_iszero_exponentmantissa(opnd1)) { 159 1.1.2.2 jdolecek /* invalid since both operands are zero */ 160 1.1.2.2 jdolecek if (Is_invalidtrap_enabled()) return(INVALIDEXCEPTION); 161 1.1.2.2 jdolecek Set_invalidflag(); 162 1.1.2.2 jdolecek Sgl_makequietnan(result); 163 1.1.2.2 jdolecek *dstptr = result; 164 1.1.2.2 jdolecek return(NOEXCEPTION); 165 1.1.2.2 jdolecek } 166 1.1.2.2 jdolecek if (Is_divisionbyzerotrap_enabled()) 167 1.1.2.2 jdolecek return(DIVISIONBYZEROEXCEPTION); 168 1.1.2.2 jdolecek Set_divisionbyzeroflag(); 169 1.1.2.2 jdolecek Sgl_setinfinity_exponentmantissa(result); 170 1.1.2.2 jdolecek *dstptr = result; 171 1.1.2.2 jdolecek return(NOEXCEPTION); 172 1.1.2.2 jdolecek } 173 1.1.2.2 jdolecek /* 174 1.1.2.2 jdolecek * Generate exponent 175 1.1.2.2 jdolecek */ 176 1.1.2.2 jdolecek dest_exponent = Sgl_exponent(opnd1) - Sgl_exponent(opnd2) + SGL_BIAS; 177 1.1.2.2 jdolecek 178 1.1.2.2 jdolecek /* 179 1.1.2.2 jdolecek * Generate mantissa 180 1.1.2.2 jdolecek */ 181 1.1.2.2 jdolecek if (Sgl_isnotzero_exponent(opnd1)) { 182 1.1.2.2 jdolecek /* set hidden bit */ 183 1.1.2.2 jdolecek Sgl_clear_signexponent_set_hidden(opnd1); 184 1.1.2.2 jdolecek } 185 1.1.2.2 jdolecek else { 186 1.1.2.2 jdolecek /* check for zero */ 187 1.1.2.2 jdolecek if (Sgl_iszero_mantissa(opnd1)) { 188 1.1.2.2 jdolecek Sgl_setzero_exponentmantissa(result); 189 1.1.2.2 jdolecek *dstptr = result; 190 1.1.2.2 jdolecek return(NOEXCEPTION); 191 1.1.2.2 jdolecek } 192 1.1.2.2 jdolecek /* is denormalized; want to normalize */ 193 1.1.2.2 jdolecek Sgl_clear_signexponent(opnd1); 194 1.1.2.2 jdolecek Sgl_leftshiftby1(opnd1); 195 1.1.2.2 jdolecek Sgl_normalize(opnd1,dest_exponent); 196 1.1.2.2 jdolecek } 197 1.1.2.2 jdolecek /* opnd2 needs to have hidden bit set with msb in hidden bit */ 198 1.1.2.2 jdolecek if (Sgl_isnotzero_exponent(opnd2)) { 199 1.1.2.2 jdolecek Sgl_clear_signexponent_set_hidden(opnd2); 200 1.1.2.2 jdolecek } 201 1.1.2.2 jdolecek else { 202 1.1.2.2 jdolecek /* is denormalized; want to normalize */ 203 1.1.2.2 jdolecek Sgl_clear_signexponent(opnd2); 204 1.1.2.2 jdolecek Sgl_leftshiftby1(opnd2); 205 1.1.2.2 jdolecek while(Sgl_iszero_hiddenhigh7mantissa(opnd2)) { 206 1.1.2.2 jdolecek Sgl_leftshiftby8(opnd2); 207 1.1.2.2 jdolecek dest_exponent += 8; 208 1.1.2.2 jdolecek } 209 1.1.2.2 jdolecek if(Sgl_iszero_hiddenhigh3mantissa(opnd2)) { 210 1.1.2.2 jdolecek Sgl_leftshiftby4(opnd2); 211 1.1.2.2 jdolecek dest_exponent += 4; 212 1.1.2.2 jdolecek } 213 1.1.2.2 jdolecek while(Sgl_iszero_hidden(opnd2)) { 214 1.1.2.2 jdolecek Sgl_leftshiftby1(opnd2); 215 1.1.2.2 jdolecek dest_exponent += 1; 216 1.1.2.2 jdolecek } 217 1.1.2.2 jdolecek } 218 1.1.2.2 jdolecek 219 1.1.2.2 jdolecek /* Divide the source mantissas */ 220 1.1.2.2 jdolecek 221 1.1.2.2 jdolecek /* 222 1.1.2.2 jdolecek * A non_restoring divide algorithm is used. 223 1.1.2.2 jdolecek */ 224 1.1.2.2 jdolecek Sgl_subtract(opnd1,opnd2,opnd1); 225 1.1.2.2 jdolecek Sgl_setzero(opnd3); 226 1.1.2.2 jdolecek for (count=1;count<=SGL_P && Sgl_all(opnd1);count++) { 227 1.1.2.2 jdolecek Sgl_leftshiftby1(opnd1); 228 1.1.2.2 jdolecek Sgl_leftshiftby1(opnd3); 229 1.1.2.2 jdolecek if (Sgl_iszero_sign(opnd1)) { 230 1.1.2.2 jdolecek Sgl_setone_lowmantissa(opnd3); 231 1.1.2.2 jdolecek Sgl_subtract(opnd1,opnd2,opnd1); 232 1.1.2.2 jdolecek } 233 1.1.2.2 jdolecek else Sgl_addition(opnd1,opnd2,opnd1); 234 1.1.2.2 jdolecek } 235 1.1.2.2 jdolecek if (count <= SGL_P) { 236 1.1.2.2 jdolecek Sgl_leftshiftby1(opnd3); 237 1.1.2.2 jdolecek Sgl_setone_lowmantissa(opnd3); 238 1.1.2.2 jdolecek Sgl_leftshift(opnd3,SGL_P-count); 239 1.1.2.2 jdolecek if (Sgl_iszero_hidden(opnd3)) { 240 1.1.2.2 jdolecek Sgl_leftshiftby1(opnd3); 241 1.1.2.2 jdolecek dest_exponent--; 242 1.1.2.2 jdolecek } 243 1.1.2.2 jdolecek } 244 1.1.2.2 jdolecek else { 245 1.1.2.2 jdolecek if (Sgl_iszero_hidden(opnd3)) { 246 1.1.2.2 jdolecek /* need to get one more bit of result */ 247 1.1.2.2 jdolecek Sgl_leftshiftby1(opnd1); 248 1.1.2.2 jdolecek Sgl_leftshiftby1(opnd3); 249 1.1.2.2 jdolecek if (Sgl_iszero_sign(opnd1)) { 250 1.1.2.2 jdolecek Sgl_setone_lowmantissa(opnd3); 251 1.1.2.2 jdolecek Sgl_subtract(opnd1,opnd2,opnd1); 252 1.1.2.2 jdolecek } 253 1.1.2.2 jdolecek else Sgl_addition(opnd1,opnd2,opnd1); 254 1.1.2.2 jdolecek dest_exponent--; 255 1.1.2.2 jdolecek } 256 1.1.2.2 jdolecek if (Sgl_iszero_sign(opnd1)) guardbit = TRUE; 257 1.1.2.2 jdolecek stickybit = Sgl_all(opnd1); 258 1.1.2.2 jdolecek } 259 1.1.2.2 jdolecek inexact = guardbit | stickybit; 260 1.1.2.2 jdolecek 261 1.1.2.2 jdolecek /* 262 1.1.2.2 jdolecek * round result 263 1.1.2.2 jdolecek */ 264 1.1.2.2 jdolecek if (inexact && (dest_exponent > 0 || Is_underflowtrap_enabled())) { 265 1.1.2.2 jdolecek Sgl_clear_signexponent(opnd3); 266 1.1.2.2 jdolecek switch (Rounding_mode()) { 267 1.1.2.2 jdolecek case ROUNDPLUS: 268 1.1.2.2 jdolecek if (Sgl_iszero_sign(result)) 269 1.1.2.2 jdolecek Sgl_increment_mantissa(opnd3); 270 1.1.2.2 jdolecek break; 271 1.1.2.2 jdolecek case ROUNDMINUS: 272 1.1.2.2 jdolecek if (Sgl_isone_sign(result)) 273 1.1.2.2 jdolecek Sgl_increment_mantissa(opnd3); 274 1.1.2.2 jdolecek break; 275 1.1.2.2 jdolecek case ROUNDNEAREST: 276 1.1.2.2 jdolecek if (guardbit && 277 1.1.2.2 jdolecek (stickybit || Sgl_isone_lowmantissa(opnd3))) 278 1.1.2.2 jdolecek Sgl_increment_mantissa(opnd3); 279 1.1.2.2 jdolecek } 280 1.1.2.2 jdolecek if (Sgl_isone_hidden(opnd3)) dest_exponent++; 281 1.1.2.2 jdolecek } 282 1.1.2.2 jdolecek Sgl_set_mantissa(result,opnd3); 283 1.1.2.2 jdolecek 284 1.1.2.2 jdolecek /* 285 1.1.2.2 jdolecek * Test for overflow 286 1.1.2.2 jdolecek */ 287 1.1.2.2 jdolecek if (dest_exponent >= SGL_INFINITY_EXPONENT) { 288 1.1.2.2 jdolecek /* trap if OVERFLOWTRAP enabled */ 289 1.1.2.2 jdolecek if (Is_overflowtrap_enabled()) { 290 1.1.2.2 jdolecek /* 291 1.1.2.2 jdolecek * Adjust bias of result 292 1.1.2.2 jdolecek */ 293 1.1.2.2 jdolecek Sgl_setwrapped_exponent(result,dest_exponent,ovfl); 294 1.1.2.2 jdolecek *dstptr = result; 295 1.1.2.2 jdolecek if (inexact) { 296 1.1.2.2 jdolecek if (Is_inexacttrap_enabled()) 297 1.1.2.2 jdolecek return(OVERFLOWEXCEPTION | INEXACTEXCEPTION); 298 1.1.2.2 jdolecek else Set_inexactflag(); 299 1.1.2.2 jdolecek } 300 1.1.2.2 jdolecek return(OVERFLOWEXCEPTION); 301 1.1.2.2 jdolecek } 302 1.1.2.2 jdolecek Set_overflowflag(); 303 1.1.2.2 jdolecek /* set result to infinity or largest number */ 304 1.1.2.2 jdolecek Sgl_setoverflow(result); 305 1.1.2.2 jdolecek inexact = TRUE; 306 1.1.2.2 jdolecek } 307 1.1.2.2 jdolecek /* 308 1.1.2.2 jdolecek * Test for underflow 309 1.1.2.2 jdolecek */ 310 1.1.2.2 jdolecek else if (dest_exponent <= 0) { 311 1.1.2.2 jdolecek /* trap if UNDERFLOWTRAP enabled */ 312 1.1.2.2 jdolecek if (Is_underflowtrap_enabled()) { 313 1.1.2.2 jdolecek /* 314 1.1.2.2 jdolecek * Adjust bias of result 315 1.1.2.2 jdolecek */ 316 1.1.2.2 jdolecek Sgl_setwrapped_exponent(result,dest_exponent,unfl); 317 1.1.2.2 jdolecek *dstptr = result; 318 1.1.2.2 jdolecek if (inexact) { 319 1.1.2.2 jdolecek if (Is_inexacttrap_enabled()) 320 1.1.2.2 jdolecek return(UNDERFLOWEXCEPTION | INEXACTEXCEPTION); 321 1.1.2.2 jdolecek else Set_inexactflag(); 322 1.1.2.2 jdolecek } 323 1.1.2.2 jdolecek return(UNDERFLOWEXCEPTION); 324 1.1.2.2 jdolecek } 325 1.1.2.2 jdolecek 326 1.1.2.2 jdolecek /* Determine if should set underflow flag */ 327 1.1.2.2 jdolecek is_tiny = TRUE; 328 1.1.2.2 jdolecek if (dest_exponent == 0 && inexact) { 329 1.1.2.2 jdolecek switch (Rounding_mode()) { 330 1.1.2.2 jdolecek case ROUNDPLUS: 331 1.1.2.2 jdolecek if (Sgl_iszero_sign(result)) { 332 1.1.2.2 jdolecek Sgl_increment(opnd3); 333 1.1.2.2 jdolecek if (Sgl_isone_hiddenoverflow(opnd3)) 334 1.1.2.2 jdolecek is_tiny = FALSE; 335 1.1.2.2 jdolecek Sgl_decrement(opnd3); 336 1.1.2.2 jdolecek } 337 1.1.2.2 jdolecek break; 338 1.1.2.2 jdolecek case ROUNDMINUS: 339 1.1.2.2 jdolecek if (Sgl_isone_sign(result)) { 340 1.1.2.2 jdolecek Sgl_increment(opnd3); 341 1.1.2.2 jdolecek if (Sgl_isone_hiddenoverflow(opnd3)) 342 1.1.2.2 jdolecek is_tiny = FALSE; 343 1.1.2.2 jdolecek Sgl_decrement(opnd3); 344 1.1.2.2 jdolecek } 345 1.1.2.2 jdolecek break; 346 1.1.2.2 jdolecek case ROUNDNEAREST: 347 1.1.2.2 jdolecek if (guardbit && (stickybit || 348 1.1.2.2 jdolecek Sgl_isone_lowmantissa(opnd3))) { 349 1.1.2.2 jdolecek Sgl_increment(opnd3); 350 1.1.2.2 jdolecek if (Sgl_isone_hiddenoverflow(opnd3)) 351 1.1.2.2 jdolecek is_tiny = FALSE; 352 1.1.2.2 jdolecek Sgl_decrement(opnd3); 353 1.1.2.2 jdolecek } 354 1.1.2.2 jdolecek break; 355 1.1.2.2 jdolecek } 356 1.1.2.2 jdolecek } 357 1.1.2.2 jdolecek 358 1.1.2.2 jdolecek /* 359 1.1.2.2 jdolecek * denormalize result or set to signed zero 360 1.1.2.2 jdolecek */ 361 1.1.2.2 jdolecek stickybit = inexact; 362 1.1.2.2 jdolecek Sgl_denormalize(opnd3,dest_exponent,guardbit,stickybit,inexact); 363 1.1.2.2 jdolecek 364 1.1.2.2 jdolecek /* return rounded number */ 365 1.1.2.2 jdolecek if (inexact) { 366 1.1.2.2 jdolecek switch (Rounding_mode()) { 367 1.1.2.2 jdolecek case ROUNDPLUS: 368 1.1.2.2 jdolecek if (Sgl_iszero_sign(result)) { 369 1.1.2.2 jdolecek Sgl_increment(opnd3); 370 1.1.2.2 jdolecek } 371 1.1.2.2 jdolecek break; 372 1.1.2.2 jdolecek case ROUNDMINUS: 373 1.1.2.2 jdolecek if (Sgl_isone_sign(result)) { 374 1.1.2.2 jdolecek Sgl_increment(opnd3); 375 1.1.2.2 jdolecek } 376 1.1.2.2 jdolecek break; 377 1.1.2.2 jdolecek case ROUNDNEAREST: 378 1.1.2.2 jdolecek if (guardbit && (stickybit || 379 1.1.2.2 jdolecek Sgl_isone_lowmantissa(opnd3))) { 380 1.1.2.2 jdolecek Sgl_increment(opnd3); 381 1.1.2.2 jdolecek } 382 1.1.2.2 jdolecek break; 383 1.1.2.2 jdolecek } 384 1.1.2.2 jdolecek if (is_tiny) Set_underflowflag(); 385 1.1.2.2 jdolecek } 386 1.1.2.2 jdolecek Sgl_set_exponentmantissa(result,opnd3); 387 1.1.2.2 jdolecek } 388 1.1.2.2 jdolecek else Sgl_set_exponent(result,dest_exponent); 389 1.1.2.2 jdolecek *dstptr = result; 390 1.1.2.2 jdolecek /* check for inexact */ 391 1.1.2.2 jdolecek if (inexact) { 392 1.1.2.2 jdolecek if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); 393 1.1.2.2 jdolecek else Set_inexactflag(); 394 1.1.2.2 jdolecek } 395 1.1.2.2 jdolecek return(NOEXCEPTION); 396 1.1.2.2 jdolecek } 397
Indexes created Fri Oct 17 23:09:53 GMT 2025