^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 1) // SPDX-License-Identifier: GPL-2.0-or-later
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 2) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 3) * Linux/PA-RISC Project (http://www.parisc-linux.org/)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 4) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 5) * Floating-point emulation code
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 6) * Copyright (C) 2001 Hewlett-Packard (Paul Bame) <bame@debian.org>
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 7) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 8) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 9) * BEGIN_DESC
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 10) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 11) * File:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 12) * @(#) pa/spmath/dfmpy.c $Revision: 1.1 $
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 13) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 14) * Purpose:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 15) * Double Precision Floating-point Multiply
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 16) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 17) * External Interfaces:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 18) * dbl_fmpy(srcptr1,srcptr2,dstptr,status)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 19) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 20) * Internal Interfaces:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 21) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 22) * Theory:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 23) * <<please update with a overview of the operation of this file>>
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 24) *
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 25) * END_DESC
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 26) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 27)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 28)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 29) #include "float.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 30) #include "dbl_float.h"
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 31)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 32) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 33) * Double Precision Floating-point Multiply
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 34) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 35)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 36) int
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 37) dbl_fmpy(
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 38) dbl_floating_point *srcptr1,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 39) dbl_floating_point *srcptr2,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 40) dbl_floating_point *dstptr,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 41) unsigned int *status)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 42) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 43) register unsigned int opnd1p1, opnd1p2, opnd2p1, opnd2p2;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 44) register unsigned int opnd3p1, opnd3p2, resultp1, resultp2;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 45) register int dest_exponent, count;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 46) register boolean inexact = FALSE, guardbit = FALSE, stickybit = FALSE;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 47) boolean is_tiny;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 48)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 49) Dbl_copyfromptr(srcptr1,opnd1p1,opnd1p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 50) Dbl_copyfromptr(srcptr2,opnd2p1,opnd2p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 51)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 52) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 53) * set sign bit of result
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 54) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 55) if (Dbl_sign(opnd1p1) ^ Dbl_sign(opnd2p1))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 56) Dbl_setnegativezerop1(resultp1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 57) else Dbl_setzerop1(resultp1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 58) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 59) * check first operand for NaN's or infinity
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 60) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 61) if (Dbl_isinfinity_exponent(opnd1p1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 62) if (Dbl_iszero_mantissa(opnd1p1,opnd1p2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 63) if (Dbl_isnotnan(opnd2p1,opnd2p2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 64) if (Dbl_iszero_exponentmantissa(opnd2p1,opnd2p2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 65) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 66) * invalid since operands are infinity
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 67) * and zero
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 68) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 69) if (Is_invalidtrap_enabled())
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 70) return(INVALIDEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 71) Set_invalidflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 72) Dbl_makequietnan(resultp1,resultp2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 73) Dbl_copytoptr(resultp1,resultp2,dstptr);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 74) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 75) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 76) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 77) * return infinity
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 78) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 79) Dbl_setinfinity_exponentmantissa(resultp1,resultp2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 80) Dbl_copytoptr(resultp1,resultp2,dstptr);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 81) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 82) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 83) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 84) else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 85) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 86) * is NaN; signaling or quiet?
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 87) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 88) if (Dbl_isone_signaling(opnd1p1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 89) /* trap if INVALIDTRAP enabled */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 90) if (Is_invalidtrap_enabled())
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 91) return(INVALIDEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 92) /* make NaN quiet */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 93) Set_invalidflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 94) Dbl_set_quiet(opnd1p1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 95) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 96) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 97) * is second operand a signaling NaN?
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 98) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 99) else if (Dbl_is_signalingnan(opnd2p1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 100) /* trap if INVALIDTRAP enabled */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 101) if (Is_invalidtrap_enabled())
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 102) return(INVALIDEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 103) /* make NaN quiet */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 104) Set_invalidflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 105) Dbl_set_quiet(opnd2p1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 106) Dbl_copytoptr(opnd2p1,opnd2p2,dstptr);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 107) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 108) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 109) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 110) * return quiet NaN
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 111) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 112) Dbl_copytoptr(opnd1p1,opnd1p2,dstptr);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 113) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 114) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 115) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 116) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 117) * check second operand for NaN's or infinity
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 118) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 119) if (Dbl_isinfinity_exponent(opnd2p1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 120) if (Dbl_iszero_mantissa(opnd2p1,opnd2p2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 121) if (Dbl_iszero_exponentmantissa(opnd1p1,opnd1p2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 122) /* invalid since operands are zero & infinity */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 123) if (Is_invalidtrap_enabled())
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 124) return(INVALIDEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 125) Set_invalidflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 126) Dbl_makequietnan(opnd2p1,opnd2p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 127) Dbl_copytoptr(opnd2p1,opnd2p2,dstptr);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 128) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 129) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 130) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 131) * return infinity
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 132) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 133) Dbl_setinfinity_exponentmantissa(resultp1,resultp2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 134) Dbl_copytoptr(resultp1,resultp2,dstptr);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 135) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 136) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 137) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 138) * is NaN; signaling or quiet?
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 139) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 140) if (Dbl_isone_signaling(opnd2p1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 141) /* trap if INVALIDTRAP enabled */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 142) if (Is_invalidtrap_enabled()) return(INVALIDEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 143) /* make NaN quiet */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 144) Set_invalidflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 145) Dbl_set_quiet(opnd2p1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 146) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 147) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 148) * return quiet NaN
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 149) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 150) Dbl_copytoptr(opnd2p1,opnd2p2,dstptr);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 151) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 152) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 153) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 154) * Generate exponent
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 155) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 156) dest_exponent = Dbl_exponent(opnd1p1) + Dbl_exponent(opnd2p1) -DBL_BIAS;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 157)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 158) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 159) * Generate mantissa
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 160) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 161) if (Dbl_isnotzero_exponent(opnd1p1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 162) /* set hidden bit */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 163) Dbl_clear_signexponent_set_hidden(opnd1p1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 164) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 165) else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 166) /* check for zero */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 167) if (Dbl_iszero_mantissa(opnd1p1,opnd1p2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 168) Dbl_setzero_exponentmantissa(resultp1,resultp2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 169) Dbl_copytoptr(resultp1,resultp2,dstptr);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 170) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 171) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 172) /* is denormalized, adjust exponent */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 173) Dbl_clear_signexponent(opnd1p1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 174) Dbl_leftshiftby1(opnd1p1,opnd1p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 175) Dbl_normalize(opnd1p1,opnd1p2,dest_exponent);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 176) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 177) /* opnd2 needs to have hidden bit set with msb in hidden bit */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 178) if (Dbl_isnotzero_exponent(opnd2p1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 179) Dbl_clear_signexponent_set_hidden(opnd2p1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 180) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 181) else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 182) /* check for zero */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 183) if (Dbl_iszero_mantissa(opnd2p1,opnd2p2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 184) Dbl_setzero_exponentmantissa(resultp1,resultp2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 185) Dbl_copytoptr(resultp1,resultp2,dstptr);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 186) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 187) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 188) /* is denormalized; want to normalize */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 189) Dbl_clear_signexponent(opnd2p1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 190) Dbl_leftshiftby1(opnd2p1,opnd2p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 191) Dbl_normalize(opnd2p1,opnd2p2,dest_exponent);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 192) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 193)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 194) /* Multiply two source mantissas together */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 195)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 196) /* make room for guard bits */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 197) Dbl_leftshiftby7(opnd2p1,opnd2p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 198) Dbl_setzero(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 199) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 200) * Four bits at a time are inspected in each loop, and a
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 201) * simple shift and add multiply algorithm is used.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 202) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 203) for (count=1;count<=DBL_P;count+=4) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 204) stickybit |= Dlow4p2(opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 205) Dbl_rightshiftby4(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 206) if (Dbit28p2(opnd1p2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 207) /* Twoword_add should be an ADDC followed by an ADD. */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 208) Twoword_add(opnd3p1, opnd3p2, opnd2p1<<3 | opnd2p2>>29,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 209) opnd2p2<<3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 210) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 211) if (Dbit29p2(opnd1p2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 212) Twoword_add(opnd3p1, opnd3p2, opnd2p1<<2 | opnd2p2>>30,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 213) opnd2p2<<2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 214) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 215) if (Dbit30p2(opnd1p2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 216) Twoword_add(opnd3p1, opnd3p2, opnd2p1<<1 | opnd2p2>>31,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 217) opnd2p2<<1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 218) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 219) if (Dbit31p2(opnd1p2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 220) Twoword_add(opnd3p1, opnd3p2, opnd2p1, opnd2p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 221) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 222) Dbl_rightshiftby4(opnd1p1,opnd1p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 223) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 224) if (Dbit3p1(opnd3p1)==0) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 225) Dbl_leftshiftby1(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 226) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 227) else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 228) /* result mantissa >= 2. */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 229) dest_exponent++;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 230) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 231) /* check for denormalized result */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 232) while (Dbit3p1(opnd3p1)==0) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 233) Dbl_leftshiftby1(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 234) dest_exponent--;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 235) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 236) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 237) * check for guard, sticky and inexact bits
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 238) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 239) stickybit |= Dallp2(opnd3p2) << 25;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 240) guardbit = (Dallp2(opnd3p2) << 24) >> 31;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 241) inexact = guardbit | stickybit;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 242)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 243) /* align result mantissa */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 244) Dbl_rightshiftby8(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 245)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 246) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 247) * round result
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 248) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 249) if (inexact && (dest_exponent>0 || Is_underflowtrap_enabled())) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 250) Dbl_clear_signexponent(opnd3p1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 251) switch (Rounding_mode()) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 252) case ROUNDPLUS:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 253) if (Dbl_iszero_sign(resultp1))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 254) Dbl_increment(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 255) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 256) case ROUNDMINUS:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 257) if (Dbl_isone_sign(resultp1))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 258) Dbl_increment(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 259) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 260) case ROUNDNEAREST:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 261) if (guardbit) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 262) if (stickybit || Dbl_isone_lowmantissap2(opnd3p2))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 263) Dbl_increment(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 264) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 265) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 266) if (Dbl_isone_hidden(opnd3p1)) dest_exponent++;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 267) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 268) Dbl_set_mantissa(resultp1,resultp2,opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 269)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 270) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 271) * Test for overflow
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 272) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 273) if (dest_exponent >= DBL_INFINITY_EXPONENT) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 274) /* trap if OVERFLOWTRAP enabled */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 275) if (Is_overflowtrap_enabled()) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 276) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 277) * Adjust bias of result
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 278) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 279) Dbl_setwrapped_exponent(resultp1,dest_exponent,ovfl);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 280) Dbl_copytoptr(resultp1,resultp2,dstptr);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 281) if (inexact)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 282) if (Is_inexacttrap_enabled())
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 283) return (OVERFLOWEXCEPTION | INEXACTEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 284) else Set_inexactflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 285) return (OVERFLOWEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 286) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 287) inexact = TRUE;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 288) Set_overflowflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 289) /* set result to infinity or largest number */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 290) Dbl_setoverflow(resultp1,resultp2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 291) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 292) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 293) * Test for underflow
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 294) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 295) else if (dest_exponent <= 0) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 296) /* trap if UNDERFLOWTRAP enabled */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 297) if (Is_underflowtrap_enabled()) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 298) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 299) * Adjust bias of result
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 300) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 301) Dbl_setwrapped_exponent(resultp1,dest_exponent,unfl);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 302) Dbl_copytoptr(resultp1,resultp2,dstptr);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 303) if (inexact)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 304) if (Is_inexacttrap_enabled())
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 305) return (UNDERFLOWEXCEPTION | INEXACTEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 306) else Set_inexactflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 307) return (UNDERFLOWEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 308) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 309)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 310) /* Determine if should set underflow flag */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 311) is_tiny = TRUE;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 312) if (dest_exponent == 0 && inexact) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 313) switch (Rounding_mode()) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 314) case ROUNDPLUS:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 315) if (Dbl_iszero_sign(resultp1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 316) Dbl_increment(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 317) if (Dbl_isone_hiddenoverflow(opnd3p1))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 318) is_tiny = FALSE;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 319) Dbl_decrement(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 320) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 321) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 322) case ROUNDMINUS:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 323) if (Dbl_isone_sign(resultp1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 324) Dbl_increment(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 325) if (Dbl_isone_hiddenoverflow(opnd3p1))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 326) is_tiny = FALSE;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 327) Dbl_decrement(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 328) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 329) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 330) case ROUNDNEAREST:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 331) if (guardbit && (stickybit ||
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 332) Dbl_isone_lowmantissap2(opnd3p2))) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 333) Dbl_increment(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 334) if (Dbl_isone_hiddenoverflow(opnd3p1))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 335) is_tiny = FALSE;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 336) Dbl_decrement(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 337) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 338) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 339) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 340) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 341)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 342) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 343) * denormalize result or set to signed zero
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 344) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 345) stickybit = inexact;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 346) Dbl_denormalize(opnd3p1,opnd3p2,dest_exponent,guardbit,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 347) stickybit,inexact);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 348)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 349) /* return zero or smallest number */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 350) if (inexact) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 351) switch (Rounding_mode()) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 352) case ROUNDPLUS:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 353) if (Dbl_iszero_sign(resultp1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 354) Dbl_increment(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 355) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 356) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 357) case ROUNDMINUS:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 358) if (Dbl_isone_sign(resultp1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 359) Dbl_increment(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 360) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 361) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 362) case ROUNDNEAREST:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 363) if (guardbit && (stickybit ||
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 364) Dbl_isone_lowmantissap2(opnd3p2))) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 365) Dbl_increment(opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 366) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 367) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 368) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 369) if (is_tiny) Set_underflowflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 370) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 371) Dbl_set_exponentmantissa(resultp1,resultp2,opnd3p1,opnd3p2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 372) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 373) else Dbl_set_exponent(resultp1,dest_exponent);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 374) /* check for inexact */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 375) Dbl_copytoptr(resultp1,resultp2,dstptr);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 376) if (inexact) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 377) if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 378) else Set_inexactflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 379) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 380) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 381) }
Orange Pi5 kernel
Deprecated Linux kernel 5.10.110 for OrangePi 5/5B/5+ boards
3 Commits
0 Branches
0 Tags