/* * Linux/PA-RISC Project (http://www.parisc-linux.org/) * * Floating-point emulation code * Copyright (C) 2001 Hewlett-Packard (Paul Bame) * * 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); }