diff options
author | Peter Breitenlohner <peb@mppmu.mpg.de> | 2014-06-18 11:12:01 +0000 |
---|---|---|
committer | Peter Breitenlohner <peb@mppmu.mpg.de> | 2014-06-18 11:12:01 +0000 |
commit | 0c2f0333a4e539325da80cdf3fb182cd17134665 (patch) | |
tree | de5daa66b5595292a33c25110aa58ff8bff75543 /Build/source/libs/gmp/gmp-6.0.0/mpf/cmp_si.c | |
parent | db612d5045ef55da386a2c3350d4811d1da1a858 (diff) |
Add libs/gmp and libs/mpfr in preparation for MetaPost 2.0
git-svn-id: svn://tug.org/texlive/trunk@34274 c570f23f-e606-0410-a88d-b1316a301751
Diffstat (limited to 'Build/source/libs/gmp/gmp-6.0.0/mpf/cmp_si.c')
-rw-r--r-- | Build/source/libs/gmp/gmp-6.0.0/mpf/cmp_si.c | 118 |
1 files changed, 118 insertions, 0 deletions
diff --git a/Build/source/libs/gmp/gmp-6.0.0/mpf/cmp_si.c b/Build/source/libs/gmp/gmp-6.0.0/mpf/cmp_si.c new file mode 100644 index 00000000000..eaa8b87da9b --- /dev/null +++ b/Build/source/libs/gmp/gmp-6.0.0/mpf/cmp_si.c @@ -0,0 +1,118 @@ +/* mpf_cmp_si -- Compare a float with a signed integer. + +Copyright 1993-1995, 1999-2002, 2004, 2012 Free Software Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library is free software; you can redistribute it and/or modify +it under the terms of either: + + * the GNU Lesser General Public License as published by the Free + Software Foundation; either version 3 of the License, or (at your + option) any later version. + +or + + * the GNU General Public License as published by the Free Software + Foundation; either version 2 of the License, or (at your option) any + later version. + +or both in parallel, as here. + +The GNU MP Library 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 copies of the GNU General Public License and the +GNU Lesser General Public License along with the GNU MP Library. If not, +see https://www.gnu.org/licenses/. */ + +#include "gmp.h" +#include "gmp-impl.h" + +int +mpf_cmp_si (mpf_srcptr u, long int vval) __GMP_NOTHROW +{ + mp_srcptr up; + mp_size_t usize; + mp_exp_t uexp; + mp_limb_t ulimb; + int usign; + unsigned long abs_vval; + + uexp = u->_mp_exp; + usize = u->_mp_size; + + /* 1. Are the signs different? */ + if ((usize < 0) == (vval < 0)) /* don't use xor, type size may differ */ + { + /* U and V are both non-negative or both negative. */ + if (usize == 0) + /* vval >= 0 */ + return -(vval != 0); + if (vval == 0) + /* usize >= 0 */ + return usize != 0; + /* Fall out. */ + } + else + { + /* Either U or V is negative, but not both. */ + return usize >= 0 ? 1 : -1; + } + + /* U and V have the same sign and are both non-zero. */ + + usign = usize >= 0 ? 1 : -1; + usize = ABS (usize); + abs_vval = ABS_CAST (unsigned long, vval); + + /* 2. Are the exponents different (V's exponent == 1)? */ +#if GMP_NAIL_BITS != 0 + if (uexp > 1 + (abs_vval > GMP_NUMB_MAX)) + return usign; + if (uexp < 1 + (abs_vval > GMP_NUMB_MAX)) + return -usign; +#else + if (uexp > 1) + return usign; + if (uexp < 1) + return -usign; +#endif + + up = u->_mp_d; + + ulimb = up[usize - 1]; +#if GMP_NAIL_BITS != 0 + if (usize >= 2 && uexp == 2) + { + if ((ulimb >> GMP_NAIL_BITS) != 0) + return usign; + ulimb = (ulimb << GMP_NUMB_BITS) | up[usize - 2]; + usize--; + } +#endif + usize--; + + /* 3. Compare the most significant mantissa limb with V. */ + if (ulimb > abs_vval) + return usign; + else if (ulimb < abs_vval) + return -usign; + + /* Ignore zeroes at the low end of U. */ + while (*up == 0) + { + up++; + usize--; + } + + /* 4. Now, if the number of limbs are different, we have a difference + since we have made sure the trailing limbs are not zero. */ + if (usize > 0) + return usign; + + /* Wow, we got zero even if we tried hard to avoid it. */ + return 0; +} |