fcnvxf.c revision 1.1.2.2
1 1.1.2.2 jdolecek /* $NetBSD: fcnvxf.c,v 1.1.2.2 2002/06/23 17:37:13 jdolecek Exp $ */ 2 1.1.2.2 jdolecek 3 1.1.2.2 jdolecek /* $OpenBSD: fcnvxf.c,v 1.5 2001/03/29 03:58:18 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 #include "../spmath/dbl_float.h" 47 1.1.2.2 jdolecek #include "../spmath/cnv_float.h" 48 1.1.2.2 jdolecek 49 1.1.2.2 jdolecek /* 50 1.1.2.2 jdolecek * Convert single fixed-point to single floating-point format 51 1.1.2.2 jdolecek */ 52 1.1.2.2 jdolecek int 53 1.1.2.2 jdolecek sgl_to_sgl_fcnvxf(srcptr,dstptr,status) 54 1.1.2.2 jdolecek 55 1.1.2.2 jdolecek int *srcptr; 56 1.1.2.2 jdolecek sgl_floating_point *dstptr; 57 1.1.2.2 jdolecek unsigned int *status; 58 1.1.2.2 jdolecek { 59 1.1.2.2 jdolecek register int src, dst_exponent; 60 1.1.2.2 jdolecek register unsigned int result = 0; 61 1.1.2.2 jdolecek 62 1.1.2.2 jdolecek src = *srcptr; 63 1.1.2.2 jdolecek /* 64 1.1.2.2 jdolecek * set sign bit of result and get magnitude of source 65 1.1.2.2 jdolecek */ 66 1.1.2.2 jdolecek if (src < 0) { 67 1.1.2.2 jdolecek Sgl_setone_sign(result); 68 1.1.2.2 jdolecek Int_negate(src); 69 1.1.2.2 jdolecek } 70 1.1.2.2 jdolecek else { 71 1.1.2.2 jdolecek Sgl_setzero_sign(result); 72 1.1.2.2 jdolecek /* Check for zero */ 73 1.1.2.2 jdolecek if (src == 0) { 74 1.1.2.2 jdolecek Sgl_setzero(result); 75 1.1.2.2 jdolecek *dstptr = result; 76 1.1.2.2 jdolecek return(NOEXCEPTION); 77 1.1.2.2 jdolecek } 78 1.1.2.2 jdolecek } 79 1.1.2.2 jdolecek /* 80 1.1.2.2 jdolecek * Generate exponent and normalized mantissa 81 1.1.2.2 jdolecek */ 82 1.1.2.2 jdolecek dst_exponent = 16; /* initialize for normalization */ 83 1.1.2.2 jdolecek /* 84 1.1.2.2 jdolecek * Check word for most significant bit set. Returns 85 1.1.2.2 jdolecek * a value in dst_exponent indicating the bit position, 86 1.1.2.2 jdolecek * between -1 and 30. 87 1.1.2.2 jdolecek */ 88 1.1.2.2 jdolecek Find_ms_one_bit(src,dst_exponent); 89 1.1.2.2 jdolecek /* left justify source, with msb at bit position 1 */ 90 1.1.2.2 jdolecek if (dst_exponent >= 0) src <<= dst_exponent; 91 1.1.2.2 jdolecek else src = 1 << 30; 92 1.1.2.2 jdolecek Sgl_set_mantissa(result, src >> (SGL_EXP_LENGTH-1)); 93 1.1.2.2 jdolecek Sgl_set_exponent(result, 30+SGL_BIAS - dst_exponent); 94 1.1.2.2 jdolecek 95 1.1.2.2 jdolecek /* check for inexact */ 96 1.1.2.2 jdolecek if (Int_isinexact_to_sgl(src)) { 97 1.1.2.2 jdolecek switch (Rounding_mode()) { 98 1.1.2.2 jdolecek case ROUNDPLUS: 99 1.1.2.2 jdolecek if (Sgl_iszero_sign(result)) 100 1.1.2.2 jdolecek Sgl_increment(result); 101 1.1.2.2 jdolecek break; 102 1.1.2.2 jdolecek case ROUNDMINUS: 103 1.1.2.2 jdolecek if (Sgl_isone_sign(result)) 104 1.1.2.2 jdolecek Sgl_increment(result); 105 1.1.2.2 jdolecek break; 106 1.1.2.2 jdolecek case ROUNDNEAREST: 107 1.1.2.2 jdolecek Sgl_roundnearest_from_int(src,result); 108 1.1.2.2 jdolecek } 109 1.1.2.2 jdolecek if (Is_inexacttrap_enabled()) { 110 1.1.2.2 jdolecek *dstptr = result; 111 1.1.2.2 jdolecek return(INEXACTEXCEPTION); 112 1.1.2.2 jdolecek } 113 1.1.2.2 jdolecek else Set_inexactflag(); 114 1.1.2.2 jdolecek } 115 1.1.2.2 jdolecek *dstptr = result; 116 1.1.2.2 jdolecek return(NOEXCEPTION); 117 1.1.2.2 jdolecek } 118 1.1.2.2 jdolecek 119 1.1.2.2 jdolecek /* 120 1.1.2.2 jdolecek * Single Fixed-point to Double Floating-point 121 1.1.2.2 jdolecek */ 122 1.1.2.2 jdolecek int 123 1.1.2.2 jdolecek sgl_to_dbl_fcnvxf(srcptr,dstptr,status) 124 1.1.2.2 jdolecek 125 1.1.2.2 jdolecek int *srcptr; 126 1.1.2.2 jdolecek dbl_floating_point *dstptr; 127 1.1.2.2 jdolecek unsigned int *status; 128 1.1.2.2 jdolecek { 129 1.1.2.2 jdolecek register int src, dst_exponent; 130 1.1.2.2 jdolecek register unsigned int resultp1 = 0, resultp2 = 0; 131 1.1.2.2 jdolecek 132 1.1.2.2 jdolecek src = *srcptr; 133 1.1.2.2 jdolecek /* 134 1.1.2.2 jdolecek * set sign bit of result and get magnitude of source 135 1.1.2.2 jdolecek */ 136 1.1.2.2 jdolecek if (src < 0) { 137 1.1.2.2 jdolecek Dbl_setone_sign(resultp1); 138 1.1.2.2 jdolecek Int_negate(src); 139 1.1.2.2 jdolecek } 140 1.1.2.2 jdolecek else { 141 1.1.2.2 jdolecek Dbl_setzero_sign(resultp1); 142 1.1.2.2 jdolecek /* Check for zero */ 143 1.1.2.2 jdolecek if (src == 0) { 144 1.1.2.2 jdolecek Dbl_setzero(resultp1,resultp2); 145 1.1.2.2 jdolecek Dbl_copytoptr(resultp1,resultp2,dstptr); 146 1.1.2.2 jdolecek return(NOEXCEPTION); 147 1.1.2.2 jdolecek } 148 1.1.2.2 jdolecek } 149 1.1.2.2 jdolecek /* 150 1.1.2.2 jdolecek * Generate exponent and normalized mantissa 151 1.1.2.2 jdolecek */ 152 1.1.2.2 jdolecek dst_exponent = 16; /* initialize for normalization */ 153 1.1.2.2 jdolecek /* 154 1.1.2.2 jdolecek * Check word for most significant bit set. Returns 155 1.1.2.2 jdolecek * a value in dst_exponent indicating the bit position, 156 1.1.2.2 jdolecek * between -1 and 30. 157 1.1.2.2 jdolecek */ 158 1.1.2.2 jdolecek Find_ms_one_bit(src,dst_exponent); 159 1.1.2.2 jdolecek /* left justify source, with msb at bit position 1 */ 160 1.1.2.2 jdolecek if (dst_exponent >= 0) src <<= dst_exponent; 161 1.1.2.2 jdolecek else src = 1 << 30; 162 1.1.2.2 jdolecek Dbl_set_mantissap1(resultp1, (src >> (DBL_EXP_LENGTH - 1))); 163 1.1.2.2 jdolecek Dbl_set_mantissap2(resultp2, (src << (33-DBL_EXP_LENGTH))); 164 1.1.2.2 jdolecek Dbl_set_exponent(resultp1, (30+DBL_BIAS) - dst_exponent); 165 1.1.2.2 jdolecek Dbl_copytoptr(resultp1,resultp2,dstptr); 166 1.1.2.2 jdolecek return(NOEXCEPTION); 167 1.1.2.2 jdolecek } 168 1.1.2.2 jdolecek 169 1.1.2.2 jdolecek /* 170 1.1.2.2 jdolecek * Double Fixed-point to Single Floating-point 171 1.1.2.2 jdolecek */ 172 1.1.2.2 jdolecek int 173 1.1.2.2 jdolecek dbl_to_sgl_fcnvxf(srcptr,dstptr,status) 174 1.1.2.2 jdolecek 175 1.1.2.2 jdolecek dbl_integer *srcptr; 176 1.1.2.2 jdolecek sgl_floating_point *dstptr; 177 1.1.2.2 jdolecek unsigned int *status; 178 1.1.2.2 jdolecek { 179 1.1.2.2 jdolecek int dst_exponent, srcp1; 180 1.1.2.2 jdolecek unsigned int result = 0, srcp2; 181 1.1.2.2 jdolecek 182 1.1.2.2 jdolecek Dint_copyfromptr(srcptr,srcp1,srcp2); 183 1.1.2.2 jdolecek /* 184 1.1.2.2 jdolecek * set sign bit of result and get magnitude of source 185 1.1.2.2 jdolecek */ 186 1.1.2.2 jdolecek if (srcp1 < 0) { 187 1.1.2.2 jdolecek Sgl_setone_sign(result); 188 1.1.2.2 jdolecek Dint_negate(srcp1,srcp2); 189 1.1.2.2 jdolecek } 190 1.1.2.2 jdolecek else { 191 1.1.2.2 jdolecek Sgl_setzero_sign(result); 192 1.1.2.2 jdolecek /* Check for zero */ 193 1.1.2.2 jdolecek if (srcp1 == 0 && srcp2 == 0) { 194 1.1.2.2 jdolecek Sgl_setzero(result); 195 1.1.2.2 jdolecek *dstptr = result; 196 1.1.2.2 jdolecek return(NOEXCEPTION); 197 1.1.2.2 jdolecek } 198 1.1.2.2 jdolecek } 199 1.1.2.2 jdolecek /* 200 1.1.2.2 jdolecek * Generate exponent and normalized mantissa 201 1.1.2.2 jdolecek */ 202 1.1.2.2 jdolecek dst_exponent = 16; /* initialize for normalization */ 203 1.1.2.2 jdolecek if (srcp1 == 0) { 204 1.1.2.2 jdolecek /* 205 1.1.2.2 jdolecek * Check word for most significant bit set. Returns 206 1.1.2.2 jdolecek * a value in dst_exponent indicating the bit position, 207 1.1.2.2 jdolecek * between -1 and 30. 208 1.1.2.2 jdolecek */ 209 1.1.2.2 jdolecek Find_ms_one_bit(srcp2,dst_exponent); 210 1.1.2.2 jdolecek /* left justify source, with msb at bit position 1 */ 211 1.1.2.2 jdolecek if (dst_exponent >= 0) { 212 1.1.2.2 jdolecek srcp1 = srcp2 << dst_exponent; 213 1.1.2.2 jdolecek srcp2 = 0; 214 1.1.2.2 jdolecek } 215 1.1.2.2 jdolecek else { 216 1.1.2.2 jdolecek srcp1 = srcp2 >> 1; 217 1.1.2.2 jdolecek srcp2 <<= 31; 218 1.1.2.2 jdolecek } 219 1.1.2.2 jdolecek /* 220 1.1.2.2 jdolecek * since msb set is in second word, need to 221 1.1.2.2 jdolecek * adjust bit position count 222 1.1.2.2 jdolecek */ 223 1.1.2.2 jdolecek dst_exponent += 32; 224 1.1.2.2 jdolecek } 225 1.1.2.2 jdolecek else { 226 1.1.2.2 jdolecek /* 227 1.1.2.2 jdolecek * Check word for most significant bit set. Returns 228 1.1.2.2 jdolecek * a value in dst_exponent indicating the bit position, 229 1.1.2.2 jdolecek * between -1 and 30. 230 1.1.2.2 jdolecek * 231 1.1.2.2 jdolecek */ 232 1.1.2.2 jdolecek Find_ms_one_bit(srcp1,dst_exponent); 233 1.1.2.2 jdolecek /* left justify source, with msb at bit position 1 */ 234 1.1.2.2 jdolecek if (dst_exponent > 0) { 235 1.1.2.2 jdolecek Variable_shift_double(srcp1,srcp2,(32-dst_exponent), 236 1.1.2.2 jdolecek srcp1); 237 1.1.2.2 jdolecek srcp2 <<= dst_exponent; 238 1.1.2.2 jdolecek } 239 1.1.2.2 jdolecek /* 240 1.1.2.2 jdolecek * If dst_exponent = 0, we don't need to shift anything. 241 1.1.2.2 jdolecek * If dst_exponent = -1, src = - 2**63 so we won't need to 242 1.1.2.2 jdolecek * shift srcp2. 243 1.1.2.2 jdolecek */ 244 1.1.2.2 jdolecek else srcp1 >>= -(dst_exponent); 245 1.1.2.2 jdolecek } 246 1.1.2.2 jdolecek Sgl_set_mantissa(result, (srcp1 >> (SGL_EXP_LENGTH - 1))); 247 1.1.2.2 jdolecek Sgl_set_exponent(result, (62+SGL_BIAS) - dst_exponent); 248 1.1.2.2 jdolecek 249 1.1.2.2 jdolecek /* check for inexact */ 250 1.1.2.2 jdolecek if (Dint_isinexact_to_sgl(srcp1,srcp2)) { 251 1.1.2.2 jdolecek switch (Rounding_mode()) { 252 1.1.2.2 jdolecek case ROUNDPLUS: 253 1.1.2.2 jdolecek if (Sgl_iszero_sign(result)) 254 1.1.2.2 jdolecek Sgl_increment(result); 255 1.1.2.2 jdolecek break; 256 1.1.2.2 jdolecek case ROUNDMINUS: 257 1.1.2.2 jdolecek if (Sgl_isone_sign(result)) 258 1.1.2.2 jdolecek Sgl_increment(result); 259 1.1.2.2 jdolecek break; 260 1.1.2.2 jdolecek case ROUNDNEAREST: 261 1.1.2.2 jdolecek Sgl_roundnearest_from_dint(srcp1,srcp2,result); 262 1.1.2.2 jdolecek } 263 1.1.2.2 jdolecek if (Is_inexacttrap_enabled()) { 264 1.1.2.2 jdolecek *dstptr = result; 265 1.1.2.2 jdolecek return(INEXACTEXCEPTION); 266 1.1.2.2 jdolecek } 267 1.1.2.2 jdolecek else Set_inexactflag(); 268 1.1.2.2 jdolecek } 269 1.1.2.2 jdolecek *dstptr = result; 270 1.1.2.2 jdolecek return(NOEXCEPTION); 271 1.1.2.2 jdolecek } 272 1.1.2.2 jdolecek 273 1.1.2.2 jdolecek /* 274 1.1.2.2 jdolecek * Double Fixed-point to Double Floating-point 275 1.1.2.2 jdolecek */ 276 1.1.2.2 jdolecek int 277 1.1.2.2 jdolecek dbl_to_dbl_fcnvxf(srcptr,dstptr,status) 278 1.1.2.2 jdolecek 279 1.1.2.2 jdolecek dbl_integer *srcptr; 280 1.1.2.2 jdolecek dbl_floating_point *dstptr; 281 1.1.2.2 jdolecek unsigned int *status; 282 1.1.2.2 jdolecek { 283 1.1.2.2 jdolecek register int srcp1, dst_exponent; 284 1.1.2.2 jdolecek register unsigned int srcp2, resultp1 = 0, resultp2 = 0; 285 1.1.2.2 jdolecek 286 1.1.2.2 jdolecek Dint_copyfromptr(srcptr,srcp1,srcp2); 287 1.1.2.2 jdolecek /* 288 1.1.2.2 jdolecek * set sign bit of result and get magnitude of source 289 1.1.2.2 jdolecek */ 290 1.1.2.2 jdolecek if (srcp1 < 0) { 291 1.1.2.2 jdolecek Dbl_setone_sign(resultp1); 292 1.1.2.2 jdolecek Dint_negate(srcp1,srcp2); 293 1.1.2.2 jdolecek } 294 1.1.2.2 jdolecek else { 295 1.1.2.2 jdolecek Dbl_setzero_sign(resultp1); 296 1.1.2.2 jdolecek /* Check for zero */ 297 1.1.2.2 jdolecek if (srcp1 == 0 && srcp2 ==0) { 298 1.1.2.2 jdolecek Dbl_setzero(resultp1,resultp2); 299 1.1.2.2 jdolecek Dbl_copytoptr(resultp1,resultp2,dstptr); 300 1.1.2.2 jdolecek return(NOEXCEPTION); 301 1.1.2.2 jdolecek } 302 1.1.2.2 jdolecek } 303 1.1.2.2 jdolecek /* 304 1.1.2.2 jdolecek * Generate exponent and normalized mantissa 305 1.1.2.2 jdolecek */ 306 1.1.2.2 jdolecek dst_exponent = 16; /* initialize for normalization */ 307 1.1.2.2 jdolecek if (srcp1 == 0) { 308 1.1.2.2 jdolecek /* 309 1.1.2.2 jdolecek * Check word for most significant bit set. Returns 310 1.1.2.2 jdolecek * a value in dst_exponent indicating the bit position, 311 1.1.2.2 jdolecek * between -1 and 30. 312 1.1.2.2 jdolecek */ 313 1.1.2.2 jdolecek Find_ms_one_bit(srcp2,dst_exponent); 314 1.1.2.2 jdolecek /* left justify source, with msb at bit position 1 */ 315 1.1.2.2 jdolecek if (dst_exponent >= 0) { 316 1.1.2.2 jdolecek srcp1 = srcp2 << dst_exponent; 317 1.1.2.2 jdolecek srcp2 = 0; 318 1.1.2.2 jdolecek } 319 1.1.2.2 jdolecek else { 320 1.1.2.2 jdolecek srcp1 = srcp2 >> 1; 321 1.1.2.2 jdolecek srcp2 <<= 31; 322 1.1.2.2 jdolecek } 323 1.1.2.2 jdolecek /* 324 1.1.2.2 jdolecek * since msb set is in second word, need to 325 1.1.2.2 jdolecek * adjust bit position count 326 1.1.2.2 jdolecek */ 327 1.1.2.2 jdolecek dst_exponent += 32; 328 1.1.2.2 jdolecek } 329 1.1.2.2 jdolecek else { 330 1.1.2.2 jdolecek /* 331 1.1.2.2 jdolecek * Check word for most significant bit set. Returns 332 1.1.2.2 jdolecek * a value in dst_exponent indicating the bit position, 333 1.1.2.2 jdolecek * between -1 and 30. 334 1.1.2.2 jdolecek */ 335 1.1.2.2 jdolecek Find_ms_one_bit(srcp1,dst_exponent); 336 1.1.2.2 jdolecek /* left justify source, with msb at bit position 1 */ 337 1.1.2.2 jdolecek if (dst_exponent > 0) { 338 1.1.2.2 jdolecek Variable_shift_double(srcp1,srcp2,(32-dst_exponent), 339 1.1.2.2 jdolecek srcp1); 340 1.1.2.2 jdolecek srcp2 <<= dst_exponent; 341 1.1.2.2 jdolecek } 342 1.1.2.2 jdolecek /* 343 1.1.2.2 jdolecek * If dst_exponent = 0, we don't need to shift anything. 344 1.1.2.2 jdolecek * If dst_exponent = -1, src = - 2**63 so we won't need to 345 1.1.2.2 jdolecek * shift srcp2. 346 1.1.2.2 jdolecek */ 347 1.1.2.2 jdolecek else srcp1 >>= -(dst_exponent); 348 1.1.2.2 jdolecek } 349 1.1.2.2 jdolecek Dbl_set_mantissap1(resultp1, srcp1 >> (DBL_EXP_LENGTH-1)); 350 1.1.2.2 jdolecek Shiftdouble(srcp1,srcp2,DBL_EXP_LENGTH-1,resultp2); 351 1.1.2.2 jdolecek Dbl_set_exponent(resultp1, (62+DBL_BIAS) - dst_exponent); 352 1.1.2.2 jdolecek 353 1.1.2.2 jdolecek /* check for inexact */ 354 1.1.2.2 jdolecek if (Dint_isinexact_to_dbl(srcp2)) { 355 1.1.2.2 jdolecek switch (Rounding_mode()) { 356 1.1.2.2 jdolecek case ROUNDPLUS: 357 1.1.2.2 jdolecek if (Dbl_iszero_sign(resultp1)) { 358 1.1.2.2 jdolecek Dbl_increment(resultp1,resultp2); 359 1.1.2.2 jdolecek } 360 1.1.2.2 jdolecek break; 361 1.1.2.2 jdolecek case ROUNDMINUS: 362 1.1.2.2 jdolecek if (Dbl_isone_sign(resultp1)) { 363 1.1.2.2 jdolecek Dbl_increment(resultp1,resultp2); 364 1.1.2.2 jdolecek } 365 1.1.2.2 jdolecek break; 366 1.1.2.2 jdolecek case ROUNDNEAREST: 367 1.1.2.2 jdolecek Dbl_roundnearest_from_dint(srcp2,resultp1, 368 1.1.2.2 jdolecek resultp2); 369 1.1.2.2 jdolecek } 370 1.1.2.2 jdolecek if (Is_inexacttrap_enabled()) { 371 1.1.2.2 jdolecek Dbl_copytoptr(resultp1,resultp2,dstptr); 372 1.1.2.2 jdolecek return(INEXACTEXCEPTION); 373 1.1.2.2 jdolecek } 374 1.1.2.2 jdolecek else Set_inexactflag(); 375 1.1.2.2 jdolecek } 376 1.1.2.2 jdolecek Dbl_copytoptr(resultp1,resultp2,dstptr); 377 1.1.2.2 jdolecek return(NOEXCEPTION); 378 1.1.2.2 jdolecek } 379
Indexes created Thu Oct 30 16:09:53 GMT 2025