Annotation of OpenXM_contrib/gmp/mpfr/mpfr.h, Revision 1.1
1.1 ! maekawa 1: /* mpfr.h -- Include file for mpfr.
! 2:
! 3: Copyright (C) 1999 PolKA project, Inria Lorraine and Loria
! 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 Library General Public License as published by
! 9: the Free Software Foundation; either version 2 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 Library General Public
! 15: License for more details.
! 16:
! 17: You should have received a copy of the GNU Library 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 <stdio.h>
! 23:
! 24: /* Cygnus does not know about *rand48 functions */
! 25: #ifdef __CYGWIN32__
! 26: #define mrand48 rand
! 27: #define drand48 rand
! 28: #define lrand48 rand
! 29: #define srand48 srand
! 30: #endif
! 31:
! 32: /* Definition of rounding modes */
! 33:
! 34: #define GMP_RNDN 0
! 35: #define GMP_RNDZ 1
! 36: #define GMP_RNDU 2
! 37: #define GMP_RNDD 3
! 38:
! 39: /* Definitions of types and their semantics */
! 40:
! 41: typedef struct {
! 42: unsigned long int _mp_prec; /* WARNING : for the mpfr type, the precision */
! 43: /* should be understood as the number of BITS,*/
! 44: /* not the number of mp_limb_t's. This means */
! 45: /* that the corresponding number of allocated
! 46: limbs is >= ceil(_mp_prec/BITS_PER_MP_LIMB) */
! 47: mp_size_t _mp_size; /* abs(_mp_size) is the number of allocated
! 48: limbs the field _mp_d points to.
! 49: The sign is that of _mp_size.
! 50: The number 0 is such that _mp_d[k-1]=0
! 51: where k = ceil(_mp_prec/BITS_PER_MP_LIMB) */
! 52: mp_exp_t _mp_exp;
! 53: mp_limb_t *_mp_d;
! 54: }
! 55: __mpfr_struct;
! 56:
! 57: /*
! 58: The number represented is
! 59:
! 60: sign(_mp_size)*(_mp_d[k-1]/B+_mp_d[k-2]/B^2+...+_mp_d[0]/B^k)*2^_mp_exp
! 61:
! 62: where k=ceil(_mp_prec/BITS_PER_MP_LIMB) and B=2^BITS_PER_MP_LIMB.
! 63:
! 64: For the msb (most significant bit) normalized representation, we must have
! 65: _mp_d[k-1]>=B/2, unless the number is zero (in that case its sign is still
! 66: given by sign(_mp_size)).
! 67:
! 68: We must also have the last k*BITS_PER_MP_LIMB-_mp_prec bits set to zero.
! 69: */
! 70:
! 71: typedef __mpfr_struct mpfr_t[1];
! 72: typedef __mpfr_struct *mpfr_ptr;
! 73: typedef __gmp_const __mpfr_struct *mpfr_srcptr;
! 74:
! 75:
! 76: /* Prototypes */
! 77:
! 78: #ifndef _PROTO
! 79: #if defined (__STDC__) || defined (__cplusplus)
! 80: #define _PROTO(x) x
! 81: #else
! 82: #define _PROTO(x) ()
! 83: #endif
! 84: #endif
! 85:
! 86: /* bit 31 of _mp_size is used for sign,
! 87: bit 30 of _mp_size is used for Nan flag,
! 88: remaining bits are used to store the number of allocated limbs */
! 89: #define FLAG_NAN(x) (((x)->_mp_size >> 30)&1)
! 90: #define SET_NAN(x) ((x)->_mp_size |= (1<<30))
! 91: #define ABSSIZE(x) ((x)->_mp_size & ((1<<30)-1))
! 92: #define SIZE(x) ((x)->_mp_size)
! 93: #define EXP(x) ((x)->_mp_exp)
! 94: #define MANT(x) ((x)->_mp_d)
! 95: #define SIGN(x) (((x)->_mp_size >> 31) ? -1 : 1)
! 96: #define ISNONNEG(x) (SIGN(x)>=0)
! 97: #define ISNEG(x) (SIGN(x)==-1)
! 98: #define CHANGE_SIGN(x) (SIZE(x) = SIZE(x) ^ (1<<31))
! 99: #define PREC(x) ((x)->_mp_prec)
! 100: #define NOTZERO(x) (MANT(x)[(PREC(x)-1)/BITS_PER_MP_LIMB])
! 101: #define SET_ZERO(x) (MANT(x)[(PREC(x)-1)/BITS_PER_MP_LIMB] = 0)
! 102:
! 103: /* reallocates the mantissa of x to q bits and sets the precision to q */
! 104: #define _mpfr_realloc(x, q) { \
! 105: (x)->_mp_d = (mp_ptr) (*_mp_reallocate_func) \
! 106: ((x)->_mp_d, (x)->_mp_prec>>3 + 1, (q)>>3 + 1); \
! 107: (x)->_mp_prec = q; }
! 108:
! 109: void mpfr_init2 _PROTO ((mpfr_ptr, unsigned long int));
! 110: int mpfr_round_raw _PROTO ((mp_limb_t *, mp_limb_t *, unsigned long, char,
! 111: unsigned long, char));
! 112: int mpfr_round_raw2 _PROTO ((mp_limb_t *, unsigned long, char, char,
! 113: unsigned long));
! 114: void mpfr_round _PROTO ((mpfr_ptr, char, unsigned long));
! 115: int mpfr_can_round _PROTO ((mpfr_ptr, unsigned long, unsigned char,
! 116: unsigned char, unsigned long));
! 117: int mpfr_can_round_raw _PROTO ((mp_limb_t *, unsigned long, int,
! 118: unsigned long,
! 119: unsigned char, unsigned char, unsigned long));
! 120: void mpfr_set_d _PROTO ((mpfr_ptr, double, unsigned char));
! 121: int mpfr_set_z _PROTO ((mpfr_ptr, mpz_srcptr, unsigned char));
! 122: double mpfr_get_d _PROTO ((mpfr_srcptr));
! 123: double mpfr_get_d2 _PROTO ((mpfr_srcptr, long));
! 124: void mpfr_set_f _PROTO ((mpfr_ptr, mpf_srcptr, char));
! 125: void mpfr_set_si _PROTO ((mpfr_ptr, long, unsigned char));
! 126: void mpfr_set_ui _PROTO ((mpfr_ptr, unsigned long, unsigned char));
! 127: void mpfr_print_raw _PROTO ((mpfr_srcptr));
! 128: void mpfr_random _PROTO ((mpfr_ptr));
! 129: void mpfr_clear _PROTO ((mpfr_ptr));
! 130: void mpfr_set_str_raw _PROTO ((mpfr_ptr, char *));
! 131: void mpfr_get_str_raw _PROTO ((char *, mpfr_srcptr));
! 132: char* mpfr_get_str _PROTO ((char *, mp_exp_t *, int, size_t, mpfr_srcptr, unsigned char));
! 133: size_t mpfr_out_str _PROTO ((FILE *, int, size_t, mpfr_srcptr, unsigned char));
! 134: void mpfr_mul _PROTO ((mpfr_ptr, mpfr_srcptr, mpfr_srcptr, unsigned char));
! 135: void mpfr_pow_ui _PROTO ((mpfr_ptr, mpfr_srcptr, unsigned int, unsigned char));
! 136: void mpfr_ui_pow_ui _PROTO ((mpfr_ptr, unsigned int, unsigned int, unsigned char));
! 137: mp_limb_t mpn_divrem_n _PROTO ((mp_limb_t *, mp_limb_t *, mp_limb_t *, mp_size_t));
! 138: mp_size_t kara_sqrtrem _PROTO ((mp_limb_t *, mp_limb_t *, mp_limb_t *, mp_size_t));
! 139: void mpfr_div _PROTO ((mpfr_ptr, mpfr_srcptr, mpfr_srcptr, unsigned char));
! 140: void mpfr_agm _PROTO ((mpfr_ptr, mpfr_srcptr, mpfr_srcptr, unsigned char));
! 141: int mpfr_sqrt _PROTO ((mpfr_ptr, mpfr_srcptr, unsigned char));
! 142: void mpfr_add _PROTO ((mpfr_ptr, mpfr_srcptr, mpfr_srcptr, unsigned char));
! 143: int mpfr_add_one_ulp _PROTO ((mpfr_ptr));
! 144: void mpfr_sub _PROTO ((mpfr_ptr, mpfr_srcptr, mpfr_srcptr, unsigned char));
! 145: void mpfr_set4 _PROTO ((mpfr_ptr, mpfr_srcptr, unsigned char, int));
! 146: void mpfr_pi _PROTO ((mpfr_ptr, unsigned char));
! 147: void mpfr_log2 _PROTO ((mpfr_ptr, unsigned char));
! 148: int mpfr_log _PROTO ((mpfr_ptr, mpfr_srcptr, unsigned char));
! 149: int mpfr_exp _PROTO ((mpfr_ptr, mpfr_srcptr, unsigned char));
! 150: int mpfr_zeta _PROTO ((mpfr_ptr, mpfr_srcptr, unsigned char));
! 151: void mpfr_mul_ui _PROTO((mpfr_ptr, mpfr_srcptr, unsigned long, unsigned char));
! 152: void mpfr_set_machine_rnd_mode _PROTO ((unsigned char));
! 153: int mpfr_cmp3 _PROTO ((mpfr_srcptr, mpfr_srcptr, long int));
! 154: int mpfr_cmp_ui_2exp _PROTO ((mpfr_srcptr, unsigned long int, int));
! 155: int mpfr_cmp_si_2exp _PROTO ((mpfr_srcptr, long int, int));
! 156: int mpfr_cmp2 _PROTO ((mpfr_srcptr, mpfr_srcptr));
! 157: void mpfr_mul_2exp _PROTO((mpfr_ptr, mpfr_srcptr, unsigned long int,unsigned char));
! 158: void mpfr_div_2exp _PROTO((mpfr_ptr, mpfr_srcptr, unsigned long int,unsigned char));
! 159: void mpfr_set_prec _PROTO((mpfr_ptr, unsigned long int));
! 160: void mpfr_set_default_prec _PROTO((unsigned long int));
! 161: extern mp_size_t __gmp_default_fp_bit_precision;
! 162: extern char __gmp_default_rounding_mode;
! 163: char * mpfr_print_rnd_mode _PROTO((unsigned char));
! 164: void mpfr_neg _PROTO((mpfr_ptr, mpfr_srcptr, unsigned char));
! 165: int mpfr_sub_one_ulp _PROTO((mpfr_ptr x));
! 166: int mpfr_div_ui _PROTO((mpfr_ptr y, mpfr_srcptr x, unsigned long u, unsigned char rnd_mode));
! 167: unsigned long int mpfr_get_prec _PROTO((mpfr_t x));
! 168:
! 169: #define mpfr_init(x) mpfr_init2(x, __gmp_default_fp_bit_precision)
! 170: #define mpfr_cmp_ui(b,i) mpfr_cmp_ui_2exp(b,i,0)
! 171: #define mpfr_cmp_si(b,i) mpfr_cmp_si_2exp(b,i,0)
! 172: #define mpfr_set(a,b,r) mpfr_set4(a,b,r,SIGN(b))
! 173: #define mpfr_cmp(b,c) mpfr_cmp3(b,c,1)
! 174:
! 175: #if (BITS_PER_MP_LIMB==32)
! 176: #define MPFR_LIMBS_PER_DOUBLE 2
! 177: #elif (BITS_PER_MP_LIMB==64)
! 178: #define MPFR_LIMBS_PER_DOUBLE 1
! 179: #endif
! 180:
! 181: /* gmp-2.0.2 had only one threshold for both multiplication and squaring */
! 182: #ifndef KARATSUBA_MUL_THRESHOLD
! 183: #ifdef KARATSUBA_THRESHOLD
! 184: #define KARATSUBA_MUL_THRESHOLD KARATSUBA_THRESHOLD
! 185: #else
! 186: #define KARATSUBA_MUL_THRESHOLD 16
! 187: #endif
! 188: #endif
! 189:
! 190: #define mpfr_init_set_si(x, i, p, rnd) \
! 191: mpfr_init2((x), (p)); mpfr_set_si((x), (i), (rnd));
! 192: #define mpfr_init_set_ui(x, i, p, rnd) \
! 193: mpfr_init2((x), (p)); mpfr_set_ui((x), (i), (rnd));
! 194: #define mpfr_init_set_d(x, d, p, rnd) \
! 195: mpfr_init2((x), (p)); mpfr_set_d((x), (d), (rnd));
! 196: #define mpfr_init_set(x, y, p, rnd) \
! 197: mpfr_init2((x), (p)); mpfr_set((x), (y), (rnd));
! 198: #define mpfr_init_set_f(x, y, p, rnd) \
! 199: mpfr_init2((x), (p)); mpfr_set_f((x), (y), (rnd));
! 200: #define mpfr_init_set_str(x, y, p, rnd) \
! 201: mpfr_init2((x), (p)); mpfr_set_str((x), (y), (rnd));
! 202: #define mpfr_init_set_str_raw(x, y, p, rnd) \
! 203: mpfr_init2((x), (p)); mpfr_set_str_raw((x), (y), (rnd));
! 204:
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>