^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/sfdiv.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) * Single Precision Floating-point Divide
^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) * sgl_fdiv(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 "sgl_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) * Single Precision Floating-point Divide
^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) sgl_fdiv (sgl_floating_point * srcptr1, sgl_floating_point * srcptr2,
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 38) sgl_floating_point * dstptr, unsigned int *status)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 39) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 40) register unsigned int opnd1, opnd2, opnd3, result;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 41) register int dest_exponent, count;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 42) register boolean inexact = FALSE, guardbit = FALSE, stickybit = FALSE;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 43) boolean is_tiny;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 44)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 45) opnd1 = *srcptr1;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 46) opnd2 = *srcptr2;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 47) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 48) * set sign bit of result
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 49) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 50) if (Sgl_sign(opnd1) ^ Sgl_sign(opnd2)) Sgl_setnegativezero(result);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 51) else Sgl_setzero(result);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 52) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 53) * check first operand for NaN's or infinity
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 54) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 55) if (Sgl_isinfinity_exponent(opnd1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 56) if (Sgl_iszero_mantissa(opnd1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 57) if (Sgl_isnotnan(opnd2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 58) if (Sgl_isinfinity(opnd2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 59) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 60) * invalid since both operands
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 61) * are infinity
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 62) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 63) if (Is_invalidtrap_enabled())
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 64) return(INVALIDEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 65) Set_invalidflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 66) Sgl_makequietnan(result);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 67) *dstptr = result;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 68) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 69) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 70) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 71) * return infinity
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 72) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 73) Sgl_setinfinity_exponentmantissa(result);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 74) *dstptr = result;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 75) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 76) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 77) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 78) else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 79) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 80) * is NaN; signaling or quiet?
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 81) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 82) if (Sgl_isone_signaling(opnd1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 83) /* trap if INVALIDTRAP enabled */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 84) if (Is_invalidtrap_enabled())
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 85) return(INVALIDEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 86) /* make NaN quiet */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 87) Set_invalidflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 88) Sgl_set_quiet(opnd1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 89) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 90) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 91) * is second operand a signaling NaN?
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 92) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 93) else if (Sgl_is_signalingnan(opnd2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 94) /* trap if INVALIDTRAP enabled */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 95) if (Is_invalidtrap_enabled())
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 96) return(INVALIDEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 97) /* make NaN quiet */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 98) Set_invalidflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 99) Sgl_set_quiet(opnd2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 100) *dstptr = opnd2;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 101) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 102) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 103) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 104) * return quiet NaN
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 105) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 106) *dstptr = opnd1;
^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) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 111) * check second operand for NaN's or infinity
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 112) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 113) if (Sgl_isinfinity_exponent(opnd2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 114) if (Sgl_iszero_mantissa(opnd2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 115) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 116) * return zero
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 117) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 118) Sgl_setzero_exponentmantissa(result);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 119) *dstptr = result;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 120) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 121) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 122) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 123) * is NaN; signaling or quiet?
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 124) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 125) if (Sgl_isone_signaling(opnd2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 126) /* trap if INVALIDTRAP enabled */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 127) if (Is_invalidtrap_enabled()) return(INVALIDEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 128) /* make NaN quiet */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 129) Set_invalidflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 130) Sgl_set_quiet(opnd2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 131) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 132) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 133) * return quiet NaN
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 134) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 135) *dstptr = opnd2;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 136) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 137) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 138) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 139) * check for division by zero
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 140) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 141) if (Sgl_iszero_exponentmantissa(opnd2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 142) if (Sgl_iszero_exponentmantissa(opnd1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 143) /* invalid since both operands are zero */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 144) if (Is_invalidtrap_enabled()) return(INVALIDEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 145) Set_invalidflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 146) Sgl_makequietnan(result);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 147) *dstptr = result;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 148) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 149) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 150) if (Is_divisionbyzerotrap_enabled())
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 151) return(DIVISIONBYZEROEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 152) Set_divisionbyzeroflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 153) Sgl_setinfinity_exponentmantissa(result);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 154) *dstptr = result;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 155) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 156) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 157) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 158) * Generate exponent
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 159) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 160) dest_exponent = Sgl_exponent(opnd1) - Sgl_exponent(opnd2) + SGL_BIAS;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 161)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 162) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 163) * Generate mantissa
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 164) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 165) if (Sgl_isnotzero_exponent(opnd1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 166) /* set hidden bit */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 167) Sgl_clear_signexponent_set_hidden(opnd1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 168) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 169) else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 170) /* check for zero */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 171) if (Sgl_iszero_mantissa(opnd1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 172) Sgl_setzero_exponentmantissa(result);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 173) *dstptr = result;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 174) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 175) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 176) /* is denormalized; want to normalize */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 177) Sgl_clear_signexponent(opnd1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 178) Sgl_leftshiftby1(opnd1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 179) Sgl_normalize(opnd1,dest_exponent);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 180) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 181) /* opnd2 needs to have hidden bit set with msb in hidden bit */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 182) if (Sgl_isnotzero_exponent(opnd2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 183) Sgl_clear_signexponent_set_hidden(opnd2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 184) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 185) else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 186) /* is denormalized; want to normalize */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 187) Sgl_clear_signexponent(opnd2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 188) Sgl_leftshiftby1(opnd2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 189) while(Sgl_iszero_hiddenhigh7mantissa(opnd2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 190) Sgl_leftshiftby8(opnd2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 191) dest_exponent += 8;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 192) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 193) if(Sgl_iszero_hiddenhigh3mantissa(opnd2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 194) Sgl_leftshiftby4(opnd2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 195) dest_exponent += 4;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 196) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 197) while(Sgl_iszero_hidden(opnd2)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 198) Sgl_leftshiftby1(opnd2);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 199) dest_exponent += 1;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 200) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 201) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 202)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 203) /* Divide the source mantissas */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 204)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 205) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 206) * A non_restoring divide algorithm is used.
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 207) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 208) Sgl_subtract(opnd1,opnd2,opnd1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 209) Sgl_setzero(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 210) for (count=1;count<=SGL_P && Sgl_all(opnd1);count++) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 211) Sgl_leftshiftby1(opnd1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 212) Sgl_leftshiftby1(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 213) if (Sgl_iszero_sign(opnd1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 214) Sgl_setone_lowmantissa(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 215) Sgl_subtract(opnd1,opnd2,opnd1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 216) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 217) else Sgl_addition(opnd1,opnd2,opnd1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 218) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 219) if (count <= SGL_P) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 220) Sgl_leftshiftby1(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 221) Sgl_setone_lowmantissa(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 222) Sgl_leftshift(opnd3,SGL_P-count);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 223) if (Sgl_iszero_hidden(opnd3)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 224) Sgl_leftshiftby1(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 225) dest_exponent--;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 226) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 227) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 228) else {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 229) if (Sgl_iszero_hidden(opnd3)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 230) /* need to get one more bit of result */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 231) Sgl_leftshiftby1(opnd1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 232) Sgl_leftshiftby1(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 233) if (Sgl_iszero_sign(opnd1)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 234) Sgl_setone_lowmantissa(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 235) Sgl_subtract(opnd1,opnd2,opnd1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 236) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 237) else Sgl_addition(opnd1,opnd2,opnd1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 238) dest_exponent--;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 239) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 240) if (Sgl_iszero_sign(opnd1)) guardbit = TRUE;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 241) stickybit = Sgl_all(opnd1);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 242) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 243) inexact = guardbit | stickybit;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 244)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 245) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 246) * round result
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 247) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 248) if (inexact && (dest_exponent > 0 || Is_underflowtrap_enabled())) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 249) Sgl_clear_signexponent(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 250) switch (Rounding_mode()) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 251) case ROUNDPLUS:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 252) if (Sgl_iszero_sign(result))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 253) Sgl_increment_mantissa(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 254) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 255) case ROUNDMINUS:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 256) if (Sgl_isone_sign(result))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 257) Sgl_increment_mantissa(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 258) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 259) case ROUNDNEAREST:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 260) if (guardbit) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 261) if (stickybit || Sgl_isone_lowmantissa(opnd3))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 262) Sgl_increment_mantissa(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 263) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 264) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 265) if (Sgl_isone_hidden(opnd3)) dest_exponent++;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 266) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 267) Sgl_set_mantissa(result,opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 268)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 269) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 270) * Test for overflow
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 271) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 272) if (dest_exponent >= SGL_INFINITY_EXPONENT) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 273) /* trap if OVERFLOWTRAP enabled */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 274) if (Is_overflowtrap_enabled()) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 275) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 276) * Adjust bias of result
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 277) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 278) Sgl_setwrapped_exponent(result,dest_exponent,ovfl);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 279) *dstptr = result;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 280) if (inexact)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 281) if (Is_inexacttrap_enabled())
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 282) return(OVERFLOWEXCEPTION | INEXACTEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 283) else Set_inexactflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 284) return(OVERFLOWEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 285) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 286) Set_overflowflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 287) /* set result to infinity or largest number */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 288) Sgl_setoverflow(result);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 289) inexact = TRUE;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 290) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 291) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 292) * Test for underflow
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 293) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 294) else if (dest_exponent <= 0) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 295) /* trap if UNDERFLOWTRAP enabled */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 296) if (Is_underflowtrap_enabled()) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 297) /*
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 298) * Adjust bias of result
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 299) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 300) Sgl_setwrapped_exponent(result,dest_exponent,unfl);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 301) *dstptr = result;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 302) if (inexact)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 303) if (Is_inexacttrap_enabled())
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 304) return(UNDERFLOWEXCEPTION | INEXACTEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 305) else Set_inexactflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 306) return(UNDERFLOWEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 307) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 308)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 309) /* Determine if should set underflow flag */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 310) is_tiny = TRUE;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 311) if (dest_exponent == 0 && inexact) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 312) switch (Rounding_mode()) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 313) case ROUNDPLUS:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 314) if (Sgl_iszero_sign(result)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 315) Sgl_increment(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 316) if (Sgl_isone_hiddenoverflow(opnd3))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 317) is_tiny = FALSE;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 318) Sgl_decrement(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 319) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 320) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 321) case ROUNDMINUS:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 322) if (Sgl_isone_sign(result)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 323) Sgl_increment(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 324) if (Sgl_isone_hiddenoverflow(opnd3))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 325) is_tiny = FALSE;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 326) Sgl_decrement(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 327) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 328) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 329) case ROUNDNEAREST:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 330) if (guardbit && (stickybit ||
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 331) Sgl_isone_lowmantissa(opnd3))) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 332) Sgl_increment(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 333) if (Sgl_isone_hiddenoverflow(opnd3))
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 334) is_tiny = FALSE;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 335) Sgl_decrement(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 336) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 337) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 338) }
^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) * denormalize result or set to signed zero
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 343) */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 344) stickybit = inexact;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 345) Sgl_denormalize(opnd3,dest_exponent,guardbit,stickybit,inexact);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 346)
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 347) /* return rounded number */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 348) if (inexact) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 349) switch (Rounding_mode()) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 350) case ROUNDPLUS:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 351) if (Sgl_iszero_sign(result)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 352) Sgl_increment(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 353) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 354) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 355) case ROUNDMINUS:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 356) if (Sgl_isone_sign(result)) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 357) Sgl_increment(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 358) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 359) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 360) case ROUNDNEAREST:
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 361) if (guardbit && (stickybit ||
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 362) Sgl_isone_lowmantissa(opnd3))) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 363) Sgl_increment(opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 364) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 365) break;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 366) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 367) if (is_tiny) Set_underflowflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 368) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 369) Sgl_set_exponentmantissa(result,opnd3);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 370) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 371) else Sgl_set_exponent(result,dest_exponent);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 372) *dstptr = result;
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 373) /* check for inexact */
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 374) if (inexact) {
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 375) if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 376) else Set_inexactflag();
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 377) }
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 378) return(NOEXCEPTION);
^8f3ce5b39 (kx 2023-10-28 12:00:06 +0300 379) }
Orange Pi5 kernel
Deprecated Linux kernel 5.10.110 for OrangePi 5/5B/5+ boards
3 Commits
0 Branches
0 Tags