| /* | |
| * Linux/PA-RISC Project (http://www.parisc-linux.org/) | |
| * | |
| * Floating-point emulation code | |
| * Copyright (C) 2001 Hewlett-Packard (Paul Bame) <bame@debian.org> | |
| * | |
| * This program is free software; you can redistribute it and/or modify | |
| * it under the terms of the GNU General Public License as published by | |
| * the Free Software Foundation; either version 2, or (at your option) | |
| * any later version. | |
| * | |
| * This program is distributed in the hope that it will be useful, | |
| * but WITHOUT ANY WARRANTY; without even the implied warranty of | |
| * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | |
| * GNU General Public License for more details. | |
| * | |
| * You should have received a copy of the GNU General Public License | |
| * along with this program; if not, write to the Free Software | |
| * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA | |
| */ | |
| /* | |
| * BEGIN_DESC | |
| * | |
| * File: | |
| * @(#) pa/spmath/fcnvfx.c $Revision: 1.1 $ | |
| * | |
| * Purpose: | |
| * Single Floating-point to Single Fixed-point | |
| * Single Floating-point to Double Fixed-point | |
| * Double Floating-point to Single Fixed-point | |
| * Double Floating-point to Double Fixed-point | |
| * | |
| * External Interfaces: | |
| * dbl_to_dbl_fcnvfx(srcptr,nullptr,dstptr,status) | |
| * dbl_to_sgl_fcnvfx(srcptr,nullptr,dstptr,status) | |
| * sgl_to_dbl_fcnvfx(srcptr,nullptr,dstptr,status) | |
| * sgl_to_sgl_fcnvfx(srcptr,nullptr,dstptr,status) | |
| * | |
| * Internal Interfaces: | |
| * | |
| * Theory: | |
| * <<please update with a overview of the operation of this file>> | |
| * | |
| * END_DESC | |
| */ | |
| #include "float.h" | |
| #include "sgl_float.h" | |
| #include "dbl_float.h" | |
| #include "cnv_float.h" | |
| /* | |
| * Single Floating-point to Single Fixed-point | |
| */ | |
| /*ARGSUSED*/ | |
| int | |
| sgl_to_sgl_fcnvfx( | |
| sgl_floating_point *srcptr, | |
| sgl_floating_point *nullptr, | |
| int *dstptr, | |
| sgl_floating_point *status) | |
| { | |
| register unsigned int src, temp; | |
| register int src_exponent, result; | |
| register boolean inexact = FALSE; | |
| src = *srcptr; | |
| src_exponent = Sgl_exponent(src) - SGL_BIAS; | |
| /* | |
| * Test for overflow | |
| */ | |
| if (src_exponent > SGL_FX_MAX_EXP) { | |
| /* check for MININT */ | |
| if ((src_exponent > SGL_FX_MAX_EXP + 1) || | |
| Sgl_isnotzero_mantissa(src) || Sgl_iszero_sign(src)) { | |
| if (Sgl_iszero_sign(src)) result = 0x7fffffff; | |
| else result = 0x80000000; | |
| if (Is_invalidtrap_enabled()) { | |
| return(INVALIDEXCEPTION); | |
| } | |
| Set_invalidflag(); | |
| *dstptr = result; | |
| return(NOEXCEPTION); | |
| } | |
| } | |
| /* | |
| * Generate result | |
| */ | |
| if (src_exponent >= 0) { | |
| temp = src; | |
| Sgl_clear_signexponent_set_hidden(temp); | |
| Int_from_sgl_mantissa(temp,src_exponent); | |
| if (Sgl_isone_sign(src)) result = -Sgl_all(temp); | |
| else result = Sgl_all(temp); | |
| /* check for inexact */ | |
| if (Sgl_isinexact_to_fix(src,src_exponent)) { | |
| inexact = TRUE; | |
| /* round result */ | |
| switch (Rounding_mode()) { | |
| case ROUNDPLUS: | |
| if (Sgl_iszero_sign(src)) result++; | |
| break; | |
| case ROUNDMINUS: | |
| if (Sgl_isone_sign(src)) result--; | |
| break; | |
| case ROUNDNEAREST: | |
| if (Sgl_isone_roundbit(src,src_exponent)) { | |
| if (Sgl_isone_stickybit(src,src_exponent) | |
| || (Sgl_isone_lowmantissa(temp))) | |
| if (Sgl_iszero_sign(src)) result++; | |
| else result--; | |
| } | |
| } | |
| } | |
| } | |
| else { | |
| result = 0; | |
| /* check for inexact */ | |
| if (Sgl_isnotzero_exponentmantissa(src)) { | |
| inexact = TRUE; | |
| /* round result */ | |
| switch (Rounding_mode()) { | |
| case ROUNDPLUS: | |
| if (Sgl_iszero_sign(src)) result++; | |
| break; | |
| case ROUNDMINUS: | |
| if (Sgl_isone_sign(src)) result--; | |
| break; | |
| case ROUNDNEAREST: | |
| if (src_exponent == -1) | |
| if (Sgl_isnotzero_mantissa(src)) | |
| if (Sgl_iszero_sign(src)) result++; | |
| else result--; | |
| } | |
| } | |
| } | |
| *dstptr = result; | |
| if (inexact) { | |
| if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); | |
| else Set_inexactflag(); | |
| } | |
| return(NOEXCEPTION); | |
| } | |
| /* | |
| * Single Floating-point to Double Fixed-point | |
| */ | |
| /*ARGSUSED*/ | |
| int | |
| sgl_to_dbl_fcnvfx( | |
| sgl_floating_point *srcptr, | |
| unsigned int *nullptr, | |
| dbl_integer *dstptr, | |
| unsigned int *status) | |
| { | |
| register int src_exponent, resultp1; | |
| register unsigned int src, temp, resultp2; | |
| register boolean inexact = FALSE; | |
| src = *srcptr; | |
| src_exponent = Sgl_exponent(src) - SGL_BIAS; | |
| /* | |
| * Test for overflow | |
| */ | |
| if (src_exponent > DBL_FX_MAX_EXP) { | |
| /* check for MININT */ | |
| if ((src_exponent > DBL_FX_MAX_EXP + 1) || | |
| Sgl_isnotzero_mantissa(src) || Sgl_iszero_sign(src)) { | |
| if (Sgl_iszero_sign(src)) { | |
| resultp1 = 0x7fffffff; | |
| resultp2 = 0xffffffff; | |
| } | |
| else { | |
| resultp1 = 0x80000000; | |
| resultp2 = 0; | |
| } | |
| if (Is_invalidtrap_enabled()) { | |
| return(INVALIDEXCEPTION); | |
| } | |
| Set_invalidflag(); | |
| Dint_copytoptr(resultp1,resultp2,dstptr); | |
| return(NOEXCEPTION); | |
| } | |
| Dint_set_minint(resultp1,resultp2); | |
| Dint_copytoptr(resultp1,resultp2,dstptr); | |
| return(NOEXCEPTION); | |
| } | |
| /* | |
| * Generate result | |
| */ | |
| if (src_exponent >= 0) { | |
| temp = src; | |
| Sgl_clear_signexponent_set_hidden(temp); | |
| Dint_from_sgl_mantissa(temp,src_exponent,resultp1,resultp2); | |
| if (Sgl_isone_sign(src)) { | |
| Dint_setone_sign(resultp1,resultp2); | |
| } | |
| /* check for inexact */ | |
| if (Sgl_isinexact_to_fix(src,src_exponent)) { | |
| inexact = TRUE; | |
| /* round result */ | |
| switch (Rounding_mode()) { | |
| case ROUNDPLUS: | |
| if (Sgl_iszero_sign(src)) { | |
| Dint_increment(resultp1,resultp2); | |
| } | |
| break; | |
| case ROUNDMINUS: | |
| if (Sgl_isone_sign(src)) { | |
| Dint_decrement(resultp1,resultp2); | |
| } | |
| break; | |
| case ROUNDNEAREST: | |
| if (Sgl_isone_roundbit(src,src_exponent)) | |
| if (Sgl_isone_stickybit(src,src_exponent) || | |
| (Dint_isone_lowp2(resultp2))) | |
| if (Sgl_iszero_sign(src)) { | |
| Dint_increment(resultp1,resultp2); | |
| } | |
| else { | |
| Dint_decrement(resultp1,resultp2); | |
| } | |
| } | |
| } | |
| } | |
| else { | |
| Dint_setzero(resultp1,resultp2); | |
| /* check for inexact */ | |
| if (Sgl_isnotzero_exponentmantissa(src)) { | |
| inexact = TRUE; | |
| /* round result */ | |
| switch (Rounding_mode()) { | |
| case ROUNDPLUS: | |
| if (Sgl_iszero_sign(src)) { | |
| Dint_increment(resultp1,resultp2); | |
| } | |
| break; | |
| case ROUNDMINUS: | |
| if (Sgl_isone_sign(src)) { | |
| Dint_decrement(resultp1,resultp2); | |
| } | |
| break; | |
| case ROUNDNEAREST: | |
| if (src_exponent == -1) | |
| if (Sgl_isnotzero_mantissa(src)) | |
| if (Sgl_iszero_sign(src)) { | |
| Dint_increment(resultp1,resultp2); | |
| } | |
| else { | |
| Dint_decrement(resultp1,resultp2); | |
| } | |
| } | |
| } | |
| } | |
| Dint_copytoptr(resultp1,resultp2,dstptr); | |
| if (inexact) { | |
| if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); | |
| else Set_inexactflag(); | |
| } | |
| return(NOEXCEPTION); | |
| } | |
| /* | |
| * Double Floating-point to Single Fixed-point | |
| */ | |
| /*ARGSUSED*/ | |
| int | |
| dbl_to_sgl_fcnvfx( | |
| dbl_floating_point *srcptr, | |
| unsigned int *nullptr, | |
| int *dstptr, | |
| unsigned int *status) | |
| { | |
| register unsigned int srcp1,srcp2, tempp1,tempp2; | |
| register int src_exponent, result; | |
| register boolean inexact = FALSE; | |
| Dbl_copyfromptr(srcptr,srcp1,srcp2); | |
| src_exponent = Dbl_exponent(srcp1) - DBL_BIAS; | |
| /* | |
| * Test for overflow | |
| */ | |
| if (src_exponent > SGL_FX_MAX_EXP) { | |
| /* check for MININT */ | |
| if (Dbl_isoverflow_to_int(src_exponent,srcp1,srcp2)) { | |
| if (Dbl_iszero_sign(srcp1)) result = 0x7fffffff; | |
| else result = 0x80000000; | |
| if (Is_invalidtrap_enabled()) { | |
| return(INVALIDEXCEPTION); | |
| } | |
| Set_invalidflag(); | |
| *dstptr = result; | |
| return(NOEXCEPTION); | |
| } | |
| } | |
| /* | |
| * Generate result | |
| */ | |
| if (src_exponent >= 0) { | |
| tempp1 = srcp1; | |
| tempp2 = srcp2; | |
| Dbl_clear_signexponent_set_hidden(tempp1); | |
| Int_from_dbl_mantissa(tempp1,tempp2,src_exponent); | |
| if (Dbl_isone_sign(srcp1) && (src_exponent <= SGL_FX_MAX_EXP)) | |
| result = -Dbl_allp1(tempp1); | |
| else result = Dbl_allp1(tempp1); | |
| /* check for inexact */ | |
| if (Dbl_isinexact_to_fix(srcp1,srcp2,src_exponent)) { | |
| inexact = TRUE; | |
| /* round result */ | |
| switch (Rounding_mode()) { | |
| case ROUNDPLUS: | |
| if (Dbl_iszero_sign(srcp1)) result++; | |
| break; | |
| case ROUNDMINUS: | |
| if (Dbl_isone_sign(srcp1)) result--; | |
| break; | |
| case ROUNDNEAREST: | |
| if (Dbl_isone_roundbit(srcp1,srcp2,src_exponent)) | |
| if (Dbl_isone_stickybit(srcp1,srcp2,src_exponent) || | |
| (Dbl_isone_lowmantissap1(tempp1))) | |
| if (Dbl_iszero_sign(srcp1)) result++; | |
| else result--; | |
| } | |
| /* check for overflow */ | |
| if ((Dbl_iszero_sign(srcp1) && result < 0) || | |
| (Dbl_isone_sign(srcp1) && result > 0)) { | |
| if (Dbl_iszero_sign(srcp1)) result = 0x7fffffff; | |
| else result = 0x80000000; | |
| if (Is_invalidtrap_enabled()) { | |
| return(INVALIDEXCEPTION); | |
| } | |
| Set_invalidflag(); | |
| *dstptr = result; | |
| return(NOEXCEPTION); | |
| } | |
| } | |
| } | |
| else { | |
| result = 0; | |
| /* check for inexact */ | |
| if (Dbl_isnotzero_exponentmantissa(srcp1,srcp2)) { | |
| inexact = TRUE; | |
| /* round result */ | |
| switch (Rounding_mode()) { | |
| case ROUNDPLUS: | |
| if (Dbl_iszero_sign(srcp1)) result++; | |
| break; | |
| case ROUNDMINUS: | |
| if (Dbl_isone_sign(srcp1)) result--; | |
| break; | |
| case ROUNDNEAREST: | |
| if (src_exponent == -1) | |
| if (Dbl_isnotzero_mantissa(srcp1,srcp2)) | |
| if (Dbl_iszero_sign(srcp1)) result++; | |
| else result--; | |
| } | |
| } | |
| } | |
| *dstptr = result; | |
| if (inexact) { | |
| if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); | |
| else Set_inexactflag(); | |
| } | |
| return(NOEXCEPTION); | |
| } | |
| /* | |
| * Double Floating-point to Double Fixed-point | |
| */ | |
| /*ARGSUSED*/ | |
| int | |
| dbl_to_dbl_fcnvfx( | |
| dbl_floating_point *srcptr, | |
| unsigned int *nullptr, | |
| dbl_integer *dstptr, | |
| unsigned int *status) | |
| { | |
| register int src_exponent, resultp1; | |
| register unsigned int srcp1, srcp2, tempp1, tempp2, resultp2; | |
| register boolean inexact = FALSE; | |
| Dbl_copyfromptr(srcptr,srcp1,srcp2); | |
| src_exponent = Dbl_exponent(srcp1) - DBL_BIAS; | |
| /* | |
| * Test for overflow | |
| */ | |
| if (src_exponent > DBL_FX_MAX_EXP) { | |
| /* check for MININT */ | |
| if ((src_exponent > DBL_FX_MAX_EXP + 1) || | |
| Dbl_isnotzero_mantissa(srcp1,srcp2) || Dbl_iszero_sign(srcp1)) { | |
| if (Dbl_iszero_sign(srcp1)) { | |
| resultp1 = 0x7fffffff; | |
| resultp2 = 0xffffffff; | |
| } | |
| else { | |
| resultp1 = 0x80000000; | |
| resultp2 = 0; | |
| } | |
| if (Is_invalidtrap_enabled()) { | |
| return(INVALIDEXCEPTION); | |
| } | |
| Set_invalidflag(); | |
| Dint_copytoptr(resultp1,resultp2,dstptr); | |
| return(NOEXCEPTION); | |
| } | |
| } | |
| /* | |
| * Generate result | |
| */ | |
| if (src_exponent >= 0) { | |
| tempp1 = srcp1; | |
| tempp2 = srcp2; | |
| Dbl_clear_signexponent_set_hidden(tempp1); | |
| Dint_from_dbl_mantissa(tempp1,tempp2,src_exponent,resultp1, | |
| resultp2); | |
| if (Dbl_isone_sign(srcp1)) { | |
| Dint_setone_sign(resultp1,resultp2); | |
| } | |
| /* check for inexact */ | |
| if (Dbl_isinexact_to_fix(srcp1,srcp2,src_exponent)) { | |
| inexact = TRUE; | |
| /* round result */ | |
| switch (Rounding_mode()) { | |
| case ROUNDPLUS: | |
| if (Dbl_iszero_sign(srcp1)) { | |
| Dint_increment(resultp1,resultp2); | |
| } | |
| break; | |
| case ROUNDMINUS: | |
| if (Dbl_isone_sign(srcp1)) { | |
| Dint_decrement(resultp1,resultp2); | |
| } | |
| break; | |
| case ROUNDNEAREST: | |
| if (Dbl_isone_roundbit(srcp1,srcp2,src_exponent)) | |
| if (Dbl_isone_stickybit(srcp1,srcp2,src_exponent) || | |
| (Dint_isone_lowp2(resultp2))) | |
| if (Dbl_iszero_sign(srcp1)) { | |
| Dint_increment(resultp1,resultp2); | |
| } | |
| else { | |
| Dint_decrement(resultp1,resultp2); | |
| } | |
| } | |
| } | |
| } | |
| else { | |
| Dint_setzero(resultp1,resultp2); | |
| /* check for inexact */ | |
| if (Dbl_isnotzero_exponentmantissa(srcp1,srcp2)) { | |
| inexact = TRUE; | |
| /* round result */ | |
| switch (Rounding_mode()) { | |
| case ROUNDPLUS: | |
| if (Dbl_iszero_sign(srcp1)) { | |
| Dint_increment(resultp1,resultp2); | |
| } | |
| break; | |
| case ROUNDMINUS: | |
| if (Dbl_isone_sign(srcp1)) { | |
| Dint_decrement(resultp1,resultp2); | |
| } | |
| break; | |
| case ROUNDNEAREST: | |
| if (src_exponent == -1) | |
| if (Dbl_isnotzero_mantissa(srcp1,srcp2)) | |
| if (Dbl_iszero_sign(srcp1)) { | |
| Dint_increment(resultp1,resultp2); | |
| } | |
| else { | |
| Dint_decrement(resultp1,resultp2); | |
| } | |
| } | |
| } | |
| } | |
| Dint_copytoptr(resultp1,resultp2,dstptr); | |
| if (inexact) { | |
| if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); | |
| else Set_inexactflag(); | |
| } | |
| return(NOEXCEPTION); | |
| } |