Annotation of OpenXM_contrib/gmp/mpfr/exceptions.c, Revision 1.1
1.1 ! ohara 1: /* Exception flags and utilities.
! 2:
! 3: Copyright 2001 Free Software Foundation.
! 4:
! 5: This file is part of the MPFR Library.
! 6:
! 7: The MPFR Library is free software; you can redistribute it and/or modify
! 8: it under the terms of the GNU Lesser General Public License as published by
! 9: the Free Software Foundation; either version 2.1 of the License, or (at your
! 10: option) any later version.
! 11:
! 12: The MPFR Library is distributed in the hope that it will be useful, but
! 13: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
! 14: or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
! 15: License for more details.
! 16:
! 17: You should have received a copy of the GNU Lesser General Public License
! 18: along with the MPFR Library; see the file COPYING.LIB. If not, write to
! 19: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
! 20: MA 02111-1307, USA. */
! 21:
! 22: #include "gmp.h"
! 23: #include "gmp-impl.h"
! 24: #include "mpfr.h"
! 25: #include "mpfr-impl.h"
! 26:
! 27: unsigned int __mpfr_flags = 0;
! 28:
! 29: mp_exp_t __mpfr_emin = MPFR_EMIN_DEFAULT;
! 30: mp_exp_t __mpfr_emax = MPFR_EMAX_DEFAULT;
! 31:
! 32: #undef mpfr_get_emin
! 33:
! 34: mp_exp_t
! 35: mpfr_get_emin (void)
! 36: {
! 37: return __mpfr_emin;
! 38: }
! 39:
! 40: #undef mpfr_set_emin
! 41:
! 42: int
! 43: mpfr_set_emin (mp_exp_t exponent)
! 44: {
! 45: if (exponent >= MPFR_EMIN_MIN && exponent <= MPFR_EMIN_MAX)
! 46: {
! 47: __mpfr_emin = exponent;
! 48: return 0;
! 49: }
! 50: else
! 51: {
! 52: return 1;
! 53: }
! 54: }
! 55:
! 56: #undef mpfr_get_emax
! 57:
! 58: mp_exp_t
! 59: mpfr_get_emax (void)
! 60: {
! 61: return __mpfr_emax;
! 62: }
! 63:
! 64: #undef mpfr_set_emax
! 65:
! 66: int
! 67: mpfr_set_emax (mp_exp_t exponent)
! 68: {
! 69: if (exponent >= MPFR_EMAX_MIN && exponent <= MPFR_EMAX_MAX)
! 70: {
! 71: __mpfr_emax = exponent;
! 72: return 0;
! 73: }
! 74: else
! 75: {
! 76: return 1;
! 77: }
! 78: }
! 79:
! 80: #undef mpfr_clear_flags
! 81:
! 82: void
! 83: mpfr_clear_flags (void)
! 84: {
! 85: __mpfr_flags = 0;
! 86: }
! 87:
! 88: #undef mpfr_clear_underflow
! 89:
! 90: void
! 91: mpfr_clear_underflow (void)
! 92: {
! 93: __mpfr_flags &= MPFR_FLAGS_ALL ^ MPFR_FLAGS_UNDERFLOW;
! 94: }
! 95:
! 96: #undef mpfr_clear_overflow
! 97:
! 98: void
! 99: mpfr_clear_overflow (void)
! 100: {
! 101: __mpfr_flags &= MPFR_FLAGS_ALL ^ MPFR_FLAGS_OVERFLOW;
! 102: }
! 103:
! 104: #undef mpfr_clear_nanflag
! 105:
! 106: void
! 107: mpfr_clear_nanflag (void)
! 108: {
! 109: __mpfr_flags &= MPFR_FLAGS_ALL ^ MPFR_FLAGS_NAN;
! 110: }
! 111:
! 112: #undef mpfr_clear_inexflag
! 113:
! 114: void
! 115: mpfr_clear_inexflag (void)
! 116: {
! 117: __mpfr_flags &= MPFR_FLAGS_ALL ^ MPFR_FLAGS_INEXACT;
! 118: }
! 119:
! 120: #undef mpfr_check_range
! 121:
! 122: int
! 123: mpfr_check_range (mpfr_ptr x, mp_rnd_t rnd_mode)
! 124: {
! 125: if (MPFR_IS_FP(x) && MPFR_NOTZERO(x))
! 126: { /* x is a non-zero FP */
! 127: mp_exp_t exp = MPFR_EXP(x);
! 128: if (exp < __mpfr_emin)
! 129: return mpfr_set_underflow(x, rnd_mode, MPFR_SIGN(x));
! 130: if (exp > __mpfr_emax)
! 131: return mpfr_set_overflow(x, rnd_mode, MPFR_SIGN(x));
! 132: }
! 133: return 0;
! 134: }
! 135:
! 136: #undef mpfr_underflow_p
! 137:
! 138: int
! 139: mpfr_underflow_p (void)
! 140: {
! 141: return __mpfr_flags & MPFR_FLAGS_UNDERFLOW;
! 142: }
! 143:
! 144: #undef mpfr_overflow_p
! 145:
! 146: int
! 147: mpfr_overflow_p (void)
! 148: {
! 149: return __mpfr_flags & MPFR_FLAGS_OVERFLOW;
! 150: }
! 151:
! 152: #undef mpfr_nanflag_p
! 153:
! 154: int
! 155: mpfr_nanflag_p (void)
! 156: {
! 157: return __mpfr_flags & MPFR_FLAGS_NAN;
! 158: }
! 159:
! 160: #undef mpfr_inexflag_p
! 161:
! 162: int
! 163: mpfr_inexflag_p (void)
! 164: {
! 165: return __mpfr_flags & MPFR_FLAGS_INEXACT;
! 166: }
! 167:
! 168: #undef mpfr_set_underflow
! 169:
! 170: int
! 171: mpfr_set_underflow (mpfr_ptr x, mp_rnd_t rnd_mode, int sign)
! 172: {
! 173: int inex;
! 174:
! 175: MPFR_CLEAR_FLAGS(x);
! 176: if ((rnd_mode == GMP_RNDU && sign > 0)
! 177: || (rnd_mode == GMP_RNDD && sign < 0))
! 178: {
! 179: mp_size_t xn;
! 180: mp_limb_t *xp;
! 181:
! 182: MPFR_EXP(x) = __mpfr_emin;
! 183: xn = (MPFR_PREC(x)-1)/BITS_PER_MP_LIMB;
! 184: xp = MPFR_MANT(x);
! 185: xp[xn] = GMP_LIMB_HIGHBIT;
! 186: MPN_ZERO(xp, xn);
! 187: inex = 1;
! 188: }
! 189: else
! 190: {
! 191: MPFR_SET_ZERO(x);
! 192: inex = -1;
! 193: }
! 194: if (MPFR_SIGN(x) != sign)
! 195: MPFR_CHANGE_SIGN(x);
! 196: __mpfr_flags |= MPFR_FLAGS_INEXACT | MPFR_FLAGS_UNDERFLOW;
! 197: return sign > 0 ? inex : -inex;
! 198: }
! 199:
! 200: #undef mpfr_set_overflow
! 201:
! 202: int
! 203: mpfr_set_overflow (mpfr_ptr x, mp_rnd_t rnd_mode, int sign)
! 204: {
! 205: int inex;
! 206:
! 207: MPFR_CLEAR_FLAGS(x);
! 208: if ((rnd_mode == GMP_RNDU && sign < 0)
! 209: || (rnd_mode == GMP_RNDD && sign > 0))
! 210: {
! 211: mp_size_t xn, i;
! 212: int sh;
! 213: mp_limb_t *xp;
! 214:
! 215: MPFR_EXP(x) = __mpfr_emax;
! 216: xn = 1 + (MPFR_PREC(x) - 1) / BITS_PER_MP_LIMB;
! 217: sh = xn * BITS_PER_MP_LIMB - MPFR_PREC(x);
! 218: xp = MPFR_MANT(x);
! 219: xp[0] = MP_LIMB_T_MAX << sh;
! 220: for (i = 1; i < xn; i++)
! 221: xp[i] = MP_LIMB_T_MAX;
! 222: inex = -1;
! 223: }
! 224: else
! 225: {
! 226: MPFR_SET_INF(x);
! 227: inex = 1;
! 228: }
! 229: if (MPFR_SIGN(x) != sign)
! 230: MPFR_CHANGE_SIGN(x);
! 231: __mpfr_flags |= MPFR_FLAGS_INEXACT | MPFR_FLAGS_OVERFLOW;
! 232: return sign > 0 ? inex : -inex;
! 233: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to exceptions.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / gmp / mpfr