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