diff options
Diffstat (limited to 'vendor/gmp-6.3.0/mpf/fits_s.h')
-rw-r--r-- | vendor/gmp-6.3.0/mpf/fits_s.h | 71 |
1 files changed, 71 insertions, 0 deletions
diff --git a/vendor/gmp-6.3.0/mpf/fits_s.h b/vendor/gmp-6.3.0/mpf/fits_s.h new file mode 100644 index 0000000..80e74be --- /dev/null +++ b/vendor/gmp-6.3.0/mpf/fits_s.h @@ -0,0 +1,71 @@ +/* mpf_fits_s*_p -- test whether an mpf fits a C signed type. + +Copyright 2001, 2002, 2014 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-impl.h" + + +/* Notice this is equivalent to mpz_set_f + mpz_fits_s*_p. */ + +int +FUNCTION (mpf_srcptr f) __GMP_NOTHROW +{ + mp_size_t fs, fn; + mp_srcptr fp; + mp_exp_t exp; + mp_limb_t fl; + + exp = EXP(f); + if (exp < 1) + return 1; /* -1 < f < 1 truncates to zero, so fits */ + + fs = SIZ (f); + fp = PTR(f); + fn = ABS (fs); + + if (exp == 1) + { + fl = fp[fn-1]; + } +#if GMP_NAIL_BITS != 0 + else if (exp == 2 && MAXIMUM > GMP_NUMB_MAX) + { + fl = fp[fn-1]; + if ((fl >> GMP_NAIL_BITS) != 0) + return 0; + fl = (fl << GMP_NUMB_BITS); + if (fn >= 2) + fl |= fp[fn-2]; + } +#endif + else + return 0; + + return fl <= (fs >= 0 ? (mp_limb_t) MAXIMUM : NEG_CAST (mp_limb_t, MINIMUM)); +} |