diff options
Diffstat (limited to 'ANDROID_3.4.5/arch/parisc/math-emu/frnd.c')
-rw-r--r-- | ANDROID_3.4.5/arch/parisc/math-emu/frnd.c | 252 |
1 files changed, 0 insertions, 252 deletions
diff --git a/ANDROID_3.4.5/arch/parisc/math-emu/frnd.c b/ANDROID_3.4.5/arch/parisc/math-emu/frnd.c deleted file mode 100644 index 904b3844..00000000 --- a/ANDROID_3.4.5/arch/parisc/math-emu/frnd.c +++ /dev/null @@ -1,252 +0,0 @@ -/* - * 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 - * - * Purpose: - * Single Floating-point Round to Integer - * Double Floating-point Round to Integer - * Quad Floating-point Round to Integer (returns unimplemented) - * - * External Interfaces: - * dbl_frnd(srcptr,nullptr,dstptr,status) - * sgl_frnd(srcptr,nullptr,dstptr,status) - * - * END_DESC -*/ - - -#include "float.h" -#include "sgl_float.h" -#include "dbl_float.h" -#include "cnv_float.h" - -/* - * Single Floating-point Round to Integer - */ - -/*ARGSUSED*/ -int -sgl_frnd(sgl_floating_point *srcptr, - unsigned int *nullptr, - sgl_floating_point *dstptr, - unsigned int *status) -{ - register unsigned int src, result; - register int src_exponent; - register boolean inexact = FALSE; - - src = *srcptr; - /* - * check source operand for NaN or infinity - */ - if ((src_exponent = Sgl_exponent(src)) == SGL_INFINITY_EXPONENT) { - /* - * is signaling NaN? - */ - if (Sgl_isone_signaling(src)) { - /* trap if INVALIDTRAP enabled */ - if (Is_invalidtrap_enabled()) return(INVALIDEXCEPTION); - /* make NaN quiet */ - Set_invalidflag(); - Sgl_set_quiet(src); - } - /* - * return quiet NaN or infinity - */ - *dstptr = src; - return(NOEXCEPTION); - } - /* - * Need to round? - */ - if ((src_exponent -= SGL_BIAS) >= SGL_P - 1) { - *dstptr = src; - return(NOEXCEPTION); - } - /* - * Generate result - */ - if (src_exponent >= 0) { - Sgl_clear_exponent_set_hidden(src); - result = src; - Sgl_rightshift(result,(SGL_P-1) - (src_exponent)); - /* 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)) Sgl_increment(result); - break; - case ROUNDMINUS: - if (Sgl_isone_sign(src)) Sgl_increment(result); - break; - case ROUNDNEAREST: - if (Sgl_isone_roundbit(src,src_exponent)) - if (Sgl_isone_stickybit(src,src_exponent) - || (Sgl_isone_lowmantissa(result))) - Sgl_increment(result); - } - } - Sgl_leftshift(result,(SGL_P-1) - (src_exponent)); - if (Sgl_isone_hiddenoverflow(result)) - Sgl_set_exponent(result,src_exponent + (SGL_BIAS+1)); - else Sgl_set_exponent(result,src_exponent + SGL_BIAS); - } - else { - result = src; /* set sign */ - Sgl_setzero_exponentmantissa(result); - /* check for inexact */ - if (Sgl_isnotzero_exponentmantissa(src)) { - inexact = TRUE; - /* round result */ - switch (Rounding_mode()) { - case ROUNDPLUS: - if (Sgl_iszero_sign(src)) - Sgl_set_exponent(result,SGL_BIAS); - break; - case ROUNDMINUS: - if (Sgl_isone_sign(src)) - Sgl_set_exponent(result,SGL_BIAS); - break; - case ROUNDNEAREST: - if (src_exponent == -1) - if (Sgl_isnotzero_mantissa(src)) - Sgl_set_exponent(result,SGL_BIAS); - } - } - } - *dstptr = result; - if (inexact) { - if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); - else Set_inexactflag(); - } - return(NOEXCEPTION); -} - -/* - * Double Floating-point Round to Integer - */ - -/*ARGSUSED*/ -int -dbl_frnd( - dbl_floating_point *srcptr, - unsigned int *nullptr, - dbl_floating_point *dstptr, - unsigned int *status) -{ - register unsigned int srcp1, srcp2, resultp1, resultp2; - register int src_exponent; - register boolean inexact = FALSE; - - Dbl_copyfromptr(srcptr,srcp1,srcp2); - /* - * check source operand for NaN or infinity - */ - if ((src_exponent = Dbl_exponent(srcp1)) == DBL_INFINITY_EXPONENT) { - /* - * is signaling NaN? - */ - if (Dbl_isone_signaling(srcp1)) { - /* trap if INVALIDTRAP enabled */ - if (Is_invalidtrap_enabled()) return(INVALIDEXCEPTION); - /* make NaN quiet */ - Set_invalidflag(); - Dbl_set_quiet(srcp1); - } - /* - * return quiet NaN or infinity - */ - Dbl_copytoptr(srcp1,srcp2,dstptr); - return(NOEXCEPTION); - } - /* - * Need to round? - */ - if ((src_exponent -= DBL_BIAS) >= DBL_P - 1) { - Dbl_copytoptr(srcp1,srcp2,dstptr); - return(NOEXCEPTION); - } - /* - * Generate result - */ - if (src_exponent >= 0) { - Dbl_clear_exponent_set_hidden(srcp1); - resultp1 = srcp1; - resultp2 = srcp2; - Dbl_rightshift(resultp1,resultp2,(DBL_P-1) - (src_exponent)); - /* 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)) - Dbl_increment(resultp1,resultp2); - break; - case ROUNDMINUS: - if (Dbl_isone_sign(srcp1)) - Dbl_increment(resultp1,resultp2); - break; - case ROUNDNEAREST: - if (Dbl_isone_roundbit(srcp1,srcp2,src_exponent)) - if (Dbl_isone_stickybit(srcp1,srcp2,src_exponent) - || (Dbl_isone_lowmantissap2(resultp2))) - Dbl_increment(resultp1,resultp2); - } - } - Dbl_leftshift(resultp1,resultp2,(DBL_P-1) - (src_exponent)); - if (Dbl_isone_hiddenoverflow(resultp1)) - Dbl_set_exponent(resultp1,src_exponent + (DBL_BIAS+1)); - else Dbl_set_exponent(resultp1,src_exponent + DBL_BIAS); - } - else { - resultp1 = srcp1; /* set sign */ - Dbl_setzero_exponentmantissa(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)) - Dbl_set_exponent(resultp1,DBL_BIAS); - break; - case ROUNDMINUS: - if (Dbl_isone_sign(srcp1)) - Dbl_set_exponent(resultp1,DBL_BIAS); - break; - case ROUNDNEAREST: - if (src_exponent == -1) - if (Dbl_isnotzero_mantissa(srcp1,srcp2)) - Dbl_set_exponent(resultp1,DBL_BIAS); - } - } - } - Dbl_copytoptr(resultp1,resultp2,dstptr); - if (inexact) { - if (Is_inexacttrap_enabled()) return(INEXACTEXCEPTION); - else Set_inexactflag(); - } - return(NOEXCEPTION); -} |