Annotation of OpenXM_contrib/gmp/gmp.h, Revision 1.1
1.1 ! maekawa 1: /* gmp.h -- Definitions for GNU multiple precision functions.
! 2:
! 3: Copyright (C) 1991, 1993, 1994, 1995, 1996 Free Software Foundation, Inc.
! 4:
! 5: This file is part of the GNU MP Library.
! 6:
! 7: The GNU MP 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 GNU MP 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 GNU MP 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: #ifndef __GMP_H__
! 23:
! 24: #ifndef __GNU_MP__
! 25: #define __GNU_MP__ 2
! 26: #define __need_size_t
! 27: #include <stddef.h>
! 28: #undef __need_size_t
! 29:
! 30: #if defined (__STDC__) || defined (__cplusplus)
! 31: #define __gmp_const const
! 32: #else
! 33: #define __gmp_const
! 34: #endif
! 35:
! 36: #if defined (__GNUC__)
! 37: #define __gmp_inline __inline__
! 38: #else
! 39: #define __gmp_inline
! 40: #endif
! 41:
! 42: #ifndef _EXTERN_INLINE
! 43: #ifdef __GNUC__
! 44: #define _EXTERN_INLINE extern __inline__
! 45: #else
! 46: #define _EXTERN_INLINE static
! 47: #endif
! 48: #endif
! 49:
! 50: #ifdef _SHORT_LIMB
! 51: typedef unsigned int mp_limb_t;
! 52: typedef int mp_limb_signed_t;
! 53: #else
! 54: #ifdef _LONG_LONG_LIMB
! 55: typedef unsigned long long int mp_limb_t;
! 56: typedef long long int mp_limb_signed_t;
! 57: #else
! 58: typedef unsigned long int mp_limb_t;
! 59: typedef long int mp_limb_signed_t;
! 60: #endif
! 61: #endif
! 62:
! 63: typedef mp_limb_t * mp_ptr;
! 64: typedef __gmp_const mp_limb_t * mp_srcptr;
! 65: typedef long int mp_size_t;
! 66: typedef long int mp_exp_t;
! 67:
! 68: #ifndef __MP_SMALL__
! 69: typedef struct
! 70: {
! 71: int _mp_alloc; /* Number of *limbs* allocated and pointed
! 72: to by the D field. */
! 73: int _mp_size; /* abs(SIZE) is the number of limbs
! 74: the last field points to. If SIZE
! 75: is negative this is a negative
! 76: number. */
! 77: mp_limb_t *_mp_d; /* Pointer to the limbs. */
! 78: } __mpz_struct;
! 79: #else
! 80: typedef struct
! 81: {
! 82: short int _mp_alloc; /* Number of *limbs* allocated and pointed
! 83: to by the D field. */
! 84: short int _mp_size; /* abs(SIZE) is the number of limbs
! 85: the last field points to. If SIZE
! 86: is negative this is a negative
! 87: number. */
! 88: mp_limb_t *_mp_d; /* Pointer to the limbs. */
! 89: } __mpz_struct;
! 90: #endif
! 91: #endif /* __GNU_MP__ */
! 92:
! 93: /* User-visible types. */
! 94: typedef __mpz_struct MP_INT;
! 95: typedef __mpz_struct mpz_t[1];
! 96:
! 97: /* Structure for rational numbers. Zero is represented as 0/any, i.e.
! 98: the denominator is ignored. Negative numbers have the sign in
! 99: the numerator. */
! 100: typedef struct
! 101: {
! 102: __mpz_struct _mp_num;
! 103: __mpz_struct _mp_den;
! 104: #if 0
! 105: int _mp_num_alloc; /* Number of limbs allocated
! 106: for the numerator. */
! 107: int _mp_num_size; /* The absolute value of this field is the
! 108: length of the numerator; the sign is the
! 109: sign of the entire rational number. */
! 110: mp_ptr _mp_num; /* Pointer to the numerator limbs. */
! 111: int _mp_den_alloc; /* Number of limbs allocated
! 112: for the denominator. */
! 113: int _mp_den_size; /* Length of the denominator. (This field
! 114: should always be positive.) */
! 115: mp_ptr _mp_den; /* Pointer to the denominator limbs. */
! 116: #endif
! 117: } __mpq_struct;
! 118:
! 119: typedef __mpq_struct MP_RAT;
! 120: typedef __mpq_struct mpq_t[1];
! 121:
! 122: typedef struct
! 123: {
! 124: int _mp_prec; /* Max precision, in number of `mp_limb_t's.
! 125: Set by mpf_init and modified by
! 126: mpf_set_prec. The area pointed to
! 127: by the `d' field contains `prec' + 1
! 128: limbs. */
! 129: int _mp_size; /* abs(SIZE) is the number of limbs
! 130: the last field points to. If SIZE
! 131: is negative this is a negative
! 132: number. */
! 133: mp_exp_t _mp_exp; /* Exponent, in the base of `mp_limb_t'. */
! 134: mp_limb_t *_mp_d; /* Pointer to the limbs. */
! 135: } __mpf_struct;
! 136:
! 137: /* typedef __mpf_struct MP_FLOAT; */
! 138: typedef __mpf_struct mpf_t[1];
! 139:
! 140: /* Types for function declarations in gmp files. */
! 141: /* ??? Should not pollute user name space with these ??? */
! 142: typedef __gmp_const __mpz_struct *mpz_srcptr;
! 143: typedef __mpz_struct *mpz_ptr;
! 144: typedef __gmp_const __mpf_struct *mpf_srcptr;
! 145: typedef __mpf_struct *mpf_ptr;
! 146: typedef __gmp_const __mpq_struct *mpq_srcptr;
! 147: typedef __mpq_struct *mpq_ptr;
! 148:
! 149: #ifndef _PROTO
! 150: #if defined (__STDC__) || defined (__cplusplus)
! 151: #define _PROTO(x) x
! 152: #else
! 153: #define _PROTO(x) ()
! 154: #endif
! 155: #endif
! 156:
! 157: #ifndef __MPN
! 158: #if defined (__STDC__) || defined (__cplusplus)
! 159: #define __MPN(x) __mpn_##x
! 160: #else
! 161: #define __MPN(x) __mpn_/**/x
! 162: #endif
! 163: #endif
! 164:
! 165: #if defined (FILE) || defined (H_STDIO) || defined (_H_STDIO) \
! 166: || defined (_STDIO_H) || defined (_STDIO_H_) || defined (__STDIO_H__) \
! 167: || defined (_STDIO_INCLUDED)
! 168: #define _GMP_H_HAVE_FILE 1
! 169: #endif
! 170:
! 171: void mp_set_memory_functions _PROTO ((void *(*) (size_t),
! 172: void *(*) (void *, size_t, size_t),
! 173: void (*) (void *, size_t)));
! 174: extern __gmp_const int mp_bits_per_limb;
! 175:
! 176: /**************** Integer (i.e. Z) routines. ****************/
! 177:
! 178: #if defined (__cplusplus)
! 179: extern "C" {
! 180: #endif
! 181: void *_mpz_realloc _PROTO ((mpz_ptr, mp_size_t));
! 182:
! 183: void mpz_abs _PROTO ((mpz_ptr, mpz_srcptr));
! 184: void mpz_add _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 185: void mpz_add_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 186: void mpz_and _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 187: void mpz_array_init _PROTO ((mpz_ptr, mp_size_t, mp_size_t));
! 188: void mpz_cdiv_q _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 189: unsigned long int mpz_cdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 190: void mpz_cdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
! 191: unsigned long int mpz_cdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
! 192: void mpz_cdiv_r _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 193: unsigned long int mpz_cdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 194: unsigned long int mpz_cdiv_ui _PROTO ((mpz_srcptr, unsigned long int));
! 195: void mpz_clear _PROTO ((mpz_ptr));
! 196: void mpz_clrbit _PROTO ((mpz_ptr, unsigned long int));
! 197: int mpz_cmp _PROTO ((mpz_srcptr, mpz_srcptr));
! 198: int mpz_cmp_si _PROTO ((mpz_srcptr, signed long int));
! 199: int mpz_cmp_ui _PROTO ((mpz_srcptr, unsigned long int));
! 200: void mpz_com _PROTO ((mpz_ptr, mpz_srcptr));
! 201: void mpz_divexact _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 202: void mpz_fac_ui _PROTO ((mpz_ptr, unsigned long int));
! 203: void mpz_fdiv_q _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 204: void mpz_fdiv_q_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 205: unsigned long int mpz_fdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 206: void mpz_fdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
! 207: unsigned long int mpz_fdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
! 208: void mpz_fdiv_r _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 209: void mpz_fdiv_r_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 210: unsigned long int mpz_fdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 211: unsigned long int mpz_fdiv_ui _PROTO ((mpz_srcptr, unsigned long int));
! 212: void mpz_gcd _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 213: unsigned long int mpz_gcd_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 214: void mpz_gcdext _PROTO ((mpz_ptr, mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
! 215: double mpz_get_d _PROTO ((mpz_srcptr));
! 216: /* signed */ long int mpz_get_si _PROTO ((mpz_srcptr));
! 217: char *mpz_get_str _PROTO ((char *, int, mpz_srcptr));
! 218: unsigned long int mpz_get_ui _PROTO ((mpz_srcptr));
! 219: mp_limb_t mpz_getlimbn _PROTO ((mpz_srcptr, mp_size_t));
! 220: unsigned long int mpz_hamdist _PROTO ((mpz_srcptr, mpz_srcptr));
! 221: void mpz_init _PROTO ((mpz_ptr));
! 222: #ifdef _GMP_H_HAVE_FILE
! 223: size_t mpz_inp_binary _PROTO ((mpz_ptr, FILE *));
! 224: size_t mpz_inp_raw _PROTO ((mpz_ptr, FILE *));
! 225: size_t mpz_inp_str _PROTO ((mpz_ptr, FILE *, int));
! 226: #endif
! 227: void mpz_init_set _PROTO ((mpz_ptr, mpz_srcptr));
! 228: void mpz_init_set_d _PROTO ((mpz_ptr, double));
! 229: void mpz_init_set_si _PROTO ((mpz_ptr, signed long int));
! 230: int mpz_init_set_str _PROTO ((mpz_ptr, const char *, int));
! 231: void mpz_init_set_ui _PROTO ((mpz_ptr, unsigned long int));
! 232: int mpz_invert _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 233: void mpz_ior _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 234: int mpz_jacobi _PROTO ((mpz_srcptr, mpz_srcptr));
! 235: int mpz_legendre _PROTO ((mpz_srcptr, mpz_srcptr));
! 236: void mpz_mod _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 237: void mpz_mul _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 238: void mpz_mul_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 239: void mpz_mul_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 240: void mpz_neg _PROTO ((mpz_ptr, mpz_srcptr));
! 241: #ifdef _GMP_H_HAVE_FILE
! 242: size_t mpz_out_binary _PROTO ((FILE *, mpz_srcptr));
! 243: size_t mpz_out_raw _PROTO ((FILE *, mpz_srcptr));
! 244: size_t mpz_out_str _PROTO ((FILE *, int, mpz_srcptr));
! 245: #endif
! 246: int mpz_perfect_square_p _PROTO ((mpz_srcptr));
! 247: unsigned long int mpz_popcount _PROTO ((mpz_srcptr));
! 248: void mpz_pow_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 249: void mpz_powm _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr, mpz_srcptr));
! 250: void mpz_powm_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int, mpz_srcptr));
! 251: int mpz_probab_prime_p _PROTO ((mpz_srcptr, int));
! 252: void mpz_random _PROTO ((mpz_ptr, mp_size_t));
! 253: void mpz_random2 _PROTO ((mpz_ptr, mp_size_t));
! 254: unsigned long int mpz_scan0 _PROTO ((mpz_srcptr, unsigned long int));
! 255: unsigned long int mpz_scan1 _PROTO ((mpz_srcptr, unsigned long int));
! 256: void mpz_set _PROTO ((mpz_ptr, mpz_srcptr));
! 257: void mpz_set_d _PROTO ((mpz_ptr, double));
! 258: void mpz_set_f _PROTO ((mpz_ptr, mpf_srcptr));
! 259: void mpz_set_q _PROTO ((mpz_ptr, mpq_srcptr));
! 260: void mpz_set_si _PROTO ((mpz_ptr, signed long int));
! 261: int mpz_set_str _PROTO ((mpz_ptr, const char *, int));
! 262: void mpz_set_ui _PROTO ((mpz_ptr, unsigned long int));
! 263: void mpz_setbit _PROTO ((mpz_ptr, unsigned long int));
! 264: size_t mpz_size _PROTO ((mpz_srcptr));
! 265: size_t mpz_sizeinbase _PROTO ((mpz_srcptr, int));
! 266: void mpz_sqrt _PROTO ((mpz_ptr, mpz_srcptr));
! 267: void mpz_sqrtrem _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr));
! 268: void mpz_sub _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 269: void mpz_sub_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 270: void mpz_tdiv_q _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 271: void mpz_tdiv_q_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 272: void mpz_tdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 273: void mpz_tdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
! 274: void mpz_tdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
! 275: void mpz_tdiv_r _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
! 276: void mpz_tdiv_r_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 277: void mpz_tdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
! 278: void mpz_ui_pow_ui _PROTO ((mpz_ptr, unsigned long int, unsigned long int));
! 279:
! 280: /**************** Rational (i.e. Q) routines. ****************/
! 281:
! 282: void mpq_init _PROTO ((mpq_ptr));
! 283: void mpq_clear _PROTO ((mpq_ptr));
! 284: void mpq_set _PROTO ((mpq_ptr, mpq_srcptr));
! 285: void mpq_set_ui _PROTO ((mpq_ptr, unsigned long int, unsigned long int));
! 286: void mpq_set_si _PROTO ((mpq_ptr, signed long int, unsigned long int));
! 287: void mpq_set_z _PROTO ((mpq_ptr, mpz_srcptr));
! 288: void mpq_add _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
! 289: void mpq_sub _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
! 290: void mpq_mul _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
! 291: void mpq_div _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
! 292: void mpq_neg _PROTO ((mpq_ptr, mpq_srcptr));
! 293: int mpq_cmp _PROTO ((mpq_srcptr, mpq_srcptr));
! 294: int mpq_cmp_ui _PROTO ((mpq_srcptr, unsigned long int, unsigned long int));
! 295: int mpq_equal _PROTO ((mpq_srcptr, mpq_srcptr));
! 296: void mpq_inv _PROTO ((mpq_ptr, mpq_srcptr));
! 297: void mpq_set_num _PROTO ((mpq_ptr, mpz_srcptr));
! 298: void mpq_set_den _PROTO ((mpq_ptr, mpz_srcptr));
! 299: void mpq_get_num _PROTO ((mpz_ptr, mpq_srcptr));
! 300: void mpq_get_den _PROTO ((mpz_ptr, mpq_srcptr));
! 301: double mpq_get_d _PROTO ((mpq_srcptr));
! 302: void mpq_canonicalize _PROTO ((mpq_ptr));
! 303:
! 304: /**************** Float (i.e. F) routines. ****************/
! 305:
! 306: void mpf_abs _PROTO ((mpf_ptr, mpf_srcptr));
! 307: void mpf_add _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
! 308: void mpf_add_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
! 309: void mpf_clear _PROTO ((mpf_ptr));
! 310: int mpf_cmp _PROTO ((mpf_srcptr, mpf_srcptr));
! 311: int mpf_cmp_si _PROTO ((mpf_srcptr, signed long int));
! 312: int mpf_cmp_ui _PROTO ((mpf_srcptr, unsigned long int));
! 313: void mpf_div _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
! 314: void mpf_div_2exp _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
! 315: void mpf_div_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
! 316: void mpf_dump _PROTO ((mpf_srcptr));
! 317: int mpf_eq _PROTO ((mpf_srcptr, mpf_srcptr, unsigned long int));
! 318: double mpf_get_d _PROTO ((mpf_srcptr));
! 319: unsigned long int mpf_get_prec _PROTO ((mpf_srcptr));
! 320: char *mpf_get_str _PROTO ((char *, mp_exp_t *, int, size_t, mpf_srcptr));
! 321: void mpf_init _PROTO ((mpf_ptr));
! 322: void mpf_init2 _PROTO ((mpf_ptr, unsigned long int));
! 323: #ifdef _GMP_H_HAVE_FILE
! 324: size_t mpf_inp_str _PROTO ((mpf_ptr, FILE *, int));
! 325: #endif
! 326: void mpf_init_set _PROTO ((mpf_ptr, mpf_srcptr));
! 327: void mpf_init_set_d _PROTO ((mpf_ptr, double));
! 328: void mpf_init_set_si _PROTO ((mpf_ptr, signed long int));
! 329: int mpf_init_set_str _PROTO ((mpf_ptr, const char *, int));
! 330: void mpf_init_set_ui _PROTO ((mpf_ptr, unsigned long int));
! 331: void mpf_mul _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
! 332: void mpf_mul_2exp _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
! 333: void mpf_mul_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
! 334: void mpf_neg _PROTO ((mpf_ptr, mpf_srcptr));
! 335: #ifdef _GMP_H_HAVE_FILE
! 336: size_t mpf_out_str _PROTO ((FILE *, int, size_t, mpf_srcptr));
! 337: #endif
! 338: void mpf_random2 _PROTO ((mpf_ptr, mp_size_t, mp_exp_t));
! 339: void mpf_reldiff _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
! 340: void mpf_set _PROTO ((mpf_ptr, mpf_srcptr));
! 341: void mpf_set_d _PROTO ((mpf_ptr, double));
! 342: void mpf_set_default_prec _PROTO ((unsigned long int));
! 343: void mpf_set_prec _PROTO ((mpf_ptr, unsigned long int));
! 344: void mpf_set_prec_raw _PROTO ((mpf_ptr, unsigned long int));
! 345: void mpf_set_q _PROTO ((mpf_ptr, mpq_srcptr));
! 346: void mpf_set_si _PROTO ((mpf_ptr, signed long int));
! 347: int mpf_set_str _PROTO ((mpf_ptr, const char *, int));
! 348: void mpf_set_ui _PROTO ((mpf_ptr, unsigned long int));
! 349: void mpf_set_z _PROTO ((mpf_ptr, mpz_srcptr));
! 350: size_t mpf_size _PROTO ((mpf_srcptr));
! 351: void mpf_sqrt _PROTO ((mpf_ptr, mpf_srcptr));
! 352: void mpf_sqrt_ui _PROTO ((mpf_ptr, unsigned long int));
! 353: void mpf_sub _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
! 354: void mpf_sub_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
! 355: void mpf_ui_div _PROTO ((mpf_ptr, unsigned long int, mpf_srcptr));
! 356: void mpf_ui_sub _PROTO ((mpf_ptr, unsigned long int, mpf_srcptr));
! 357: #if defined (__cplusplus)
! 358: }
! 359: #endif
! 360: /************ Low level positive-integer (i.e. N) routines. ************/
! 361:
! 362: /* This is ugly, but we need to make usr calls reach the prefixed function. */
! 363: #define mpn_add __MPN(add)
! 364: #define mpn_add_1 __MPN(add_1)
! 365: #define mpn_add_n __MPN(add_n)
! 366: #define mpn_addmul_1 __MPN(addmul_1)
! 367: #define mpn_bdivmod __MPN(bdivmod)
! 368: #define mpn_cmp __MPN(cmp)
! 369: #define mpn_divmod_1 __MPN(divmod_1)
! 370: #define mpn_divrem __MPN(divrem)
! 371: #define mpn_divrem_1 __MPN(divrem_1)
! 372: #define mpn_dump __MPN(dump)
! 373: #define mpn_gcd __MPN(gcd)
! 374: #define mpn_gcd_1 __MPN(gcd_1)
! 375: #define mpn_gcdext __MPN(gcdext)
! 376: #define mpn_get_str __MPN(get_str)
! 377: #define mpn_hamdist __MPN(hamdist)
! 378: #define mpn_lshift __MPN(lshift)
! 379: #define mpn_mod_1 __MPN(mod_1)
! 380: #define mpn_mul __MPN(mul)
! 381: #define mpn_mul_1 __MPN(mul_1)
! 382: #define mpn_mul_n __MPN(mul_n)
! 383: #define mpn_perfect_square_p __MPN(perfect_square_p)
! 384: #define mpn_popcount __MPN(popcount)
! 385: #define mpn_preinv_mod_1 __MPN(preinv_mod_1)
! 386: #define mpn_random2 __MPN(random2)
! 387: #define mpn_rshift __MPN(rshift)
! 388: #define mpn_scan0 __MPN(scan0)
! 389: #define mpn_scan1 __MPN(scan1)
! 390: #define mpn_set_str __MPN(set_str)
! 391: #define mpn_sqrtrem __MPN(sqrtrem)
! 392: #define mpn_sub __MPN(sub)
! 393: #define mpn_sub_1 __MPN(sub_1)
! 394: #define mpn_sub_n __MPN(sub_n)
! 395: #define mpn_submul_1 __MPN(submul_1)
! 396: #define mpn_udiv_w_sdiv __MPN(udiv_w_sdiv)
! 397:
! 398: #if defined (__cplusplus)
! 399: extern "C" {
! 400: #endif
! 401: mp_limb_t mpn_add _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr,mp_size_t));
! 402: mp_limb_t mpn_add_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
! 403: mp_limb_t mpn_add_n _PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t));
! 404: mp_limb_t mpn_addmul_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
! 405: mp_limb_t mpn_bdivmod _PROTO ((mp_ptr, mp_ptr, mp_size_t, mp_srcptr, mp_size_t, unsigned long int));
! 406: int mpn_cmp _PROTO ((mp_srcptr, mp_srcptr, mp_size_t));
! 407: mp_limb_t mpn_divmod_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
! 408: mp_limb_t mpn_divrem _PROTO ((mp_ptr, mp_size_t, mp_ptr, mp_size_t, mp_srcptr, mp_size_t));
! 409: mp_limb_t mpn_divrem_1 _PROTO ((mp_ptr, mp_size_t, mp_srcptr, mp_size_t, mp_limb_t));
! 410: void mpn_dump _PROTO ((mp_srcptr, mp_size_t));
! 411: mp_size_t mpn_gcd _PROTO ((mp_ptr, mp_ptr, mp_size_t, mp_ptr, mp_size_t));
! 412: mp_limb_t mpn_gcd_1 _PROTO ((mp_srcptr, mp_size_t, mp_limb_t));
! 413: mp_size_t mpn_gcdext _PROTO ((mp_ptr, mp_ptr, mp_ptr, mp_size_t, mp_ptr, mp_size_t));
! 414: size_t mpn_get_str _PROTO ((unsigned char *, int, mp_ptr, mp_size_t));
! 415: unsigned long int mpn_hamdist _PROTO ((mp_srcptr, mp_srcptr, mp_size_t));
! 416: mp_limb_t mpn_lshift _PROTO ((mp_ptr, mp_srcptr, mp_size_t, unsigned int));
! 417: mp_limb_t mpn_mod_1 _PROTO ((mp_srcptr, mp_size_t, mp_limb_t));
! 418: mp_limb_t mpn_mul _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t));
! 419: mp_limb_t mpn_mul_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
! 420: void mpn_mul_n _PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t));
! 421: int mpn_perfect_square_p _PROTO ((mp_srcptr, mp_size_t));
! 422: unsigned long int mpn_popcount _PROTO ((mp_srcptr, mp_size_t));
! 423: mp_limb_t mpn_preinv_mod_1 _PROTO ((mp_srcptr, mp_size_t, mp_limb_t, mp_limb_t));
! 424: void mpn_random2 _PROTO ((mp_ptr, mp_size_t));
! 425: mp_limb_t mpn_rshift _PROTO ((mp_ptr, mp_srcptr, mp_size_t, unsigned int));
! 426: unsigned long int mpn_scan0 _PROTO ((mp_srcptr, unsigned long int));
! 427: unsigned long int mpn_scan1 _PROTO ((mp_srcptr, unsigned long int));
! 428: mp_size_t mpn_set_str _PROTO ((mp_ptr, const unsigned char *, size_t, int));
! 429: mp_size_t mpn_sqrtrem _PROTO ((mp_ptr, mp_ptr, mp_srcptr, mp_size_t));
! 430: mp_limb_t mpn_sub _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr,mp_size_t));
! 431: mp_limb_t mpn_sub_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
! 432: mp_limb_t mpn_sub_n _PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t));
! 433: mp_limb_t mpn_submul_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
! 434: #if defined (__cplusplus)
! 435: }
! 436: #endif
! 437:
! 438: #if defined (__GNUC__) || defined (_FORCE_INLINES)
! 439: _EXTERN_INLINE mp_limb_t
! 440: #if defined (__STDC__) || defined (__cplusplus)
! 441: mpn_add_1 (register mp_ptr res_ptr,
! 442: register mp_srcptr s1_ptr,
! 443: register mp_size_t s1_size,
! 444: register mp_limb_t s2_limb)
! 445: #else
! 446: mpn_add_1 (res_ptr, s1_ptr, s1_size, s2_limb)
! 447: register mp_ptr res_ptr;
! 448: register mp_srcptr s1_ptr;
! 449: register mp_size_t s1_size;
! 450: register mp_limb_t s2_limb;
! 451: #endif
! 452: {
! 453: register mp_limb_t x;
! 454:
! 455: x = *s1_ptr++;
! 456: s2_limb = x + s2_limb;
! 457: *res_ptr++ = s2_limb;
! 458: if (s2_limb < x)
! 459: {
! 460: while (--s1_size != 0)
! 461: {
! 462: x = *s1_ptr++ + 1;
! 463: *res_ptr++ = x;
! 464: if (x != 0)
! 465: goto fin;
! 466: }
! 467:
! 468: return 1;
! 469: }
! 470:
! 471: fin:
! 472: if (res_ptr != s1_ptr)
! 473: {
! 474: mp_size_t i;
! 475: for (i = 0; i < s1_size - 1; i++)
! 476: res_ptr[i] = s1_ptr[i];
! 477: }
! 478: return 0;
! 479: }
! 480:
! 481: _EXTERN_INLINE mp_limb_t
! 482: #if defined (__STDC__) || defined (__cplusplus)
! 483: mpn_add (register mp_ptr res_ptr,
! 484: register mp_srcptr s1_ptr,
! 485: register mp_size_t s1_size,
! 486: register mp_srcptr s2_ptr,
! 487: register mp_size_t s2_size)
! 488: #else
! 489: mpn_add (res_ptr, s1_ptr, s1_size, s2_ptr, s2_size)
! 490: register mp_ptr res_ptr;
! 491: register mp_srcptr s1_ptr;
! 492: register mp_size_t s1_size;
! 493: register mp_srcptr s2_ptr;
! 494: register mp_size_t s2_size;
! 495: #endif
! 496: {
! 497: mp_limb_t cy_limb = 0;
! 498:
! 499: if (s2_size != 0)
! 500: cy_limb = mpn_add_n (res_ptr, s1_ptr, s2_ptr, s2_size);
! 501:
! 502: if (s1_size - s2_size != 0)
! 503: cy_limb = mpn_add_1 (res_ptr + s2_size,
! 504: s1_ptr + s2_size,
! 505: s1_size - s2_size,
! 506: cy_limb);
! 507: return cy_limb;
! 508: }
! 509:
! 510: _EXTERN_INLINE mp_limb_t
! 511: #if defined (__STDC__) || defined (__cplusplus)
! 512: mpn_sub_1 (register mp_ptr res_ptr,
! 513: register mp_srcptr s1_ptr,
! 514: register mp_size_t s1_size,
! 515: register mp_limb_t s2_limb)
! 516: #else
! 517: mpn_sub_1 (res_ptr, s1_ptr, s1_size, s2_limb)
! 518: register mp_ptr res_ptr;
! 519: register mp_srcptr s1_ptr;
! 520: register mp_size_t s1_size;
! 521: register mp_limb_t s2_limb;
! 522: #endif
! 523: {
! 524: register mp_limb_t x;
! 525:
! 526: x = *s1_ptr++;
! 527: s2_limb = x - s2_limb;
! 528: *res_ptr++ = s2_limb;
! 529: if (s2_limb > x)
! 530: {
! 531: while (--s1_size != 0)
! 532: {
! 533: x = *s1_ptr++;
! 534: *res_ptr++ = x - 1;
! 535: if (x != 0)
! 536: goto fin;
! 537: }
! 538:
! 539: return 1;
! 540: }
! 541:
! 542: fin:
! 543: if (res_ptr != s1_ptr)
! 544: {
! 545: mp_size_t i;
! 546: for (i = 0; i < s1_size - 1; i++)
! 547: res_ptr[i] = s1_ptr[i];
! 548: }
! 549: return 0;
! 550: }
! 551:
! 552: _EXTERN_INLINE mp_limb_t
! 553: #if defined (__STDC__) || defined (__cplusplus)
! 554: mpn_sub (register mp_ptr res_ptr,
! 555: register mp_srcptr s1_ptr,
! 556: register mp_size_t s1_size,
! 557: register mp_srcptr s2_ptr,
! 558: register mp_size_t s2_size)
! 559: #else
! 560: mpn_sub (res_ptr, s1_ptr, s1_size, s2_ptr, s2_size)
! 561: register mp_ptr res_ptr;
! 562: register mp_srcptr s1_ptr;
! 563: register mp_size_t s1_size;
! 564: register mp_srcptr s2_ptr;
! 565: register mp_size_t s2_size;
! 566: #endif
! 567: {
! 568: mp_limb_t cy_limb = 0;
! 569:
! 570: if (s2_size != 0)
! 571: cy_limb = mpn_sub_n (res_ptr, s1_ptr, s2_ptr, s2_size);
! 572:
! 573: if (s1_size - s2_size != 0)
! 574: cy_limb = mpn_sub_1 (res_ptr + s2_size,
! 575: s1_ptr + s2_size,
! 576: s1_size - s2_size,
! 577: cy_limb);
! 578: return cy_limb;
! 579: }
! 580: #endif /* __GNUC__ */
! 581:
! 582: /* Allow faster testing for negative, zero, and positive. */
! 583: #define mpz_sgn(Z) ((Z)->_mp_size < 0 ? -1 : (Z)->_mp_size > 0)
! 584: #define mpf_sgn(F) ((F)->_mp_size < 0 ? -1 : (F)->_mp_size > 0)
! 585: #define mpq_sgn(Q) ((Q)->_mp_num._mp_size < 0 ? -1 : (Q)->_mp_num._mp_size > 0)
! 586:
! 587: /* Allow direct user access to numerator and denominator of a mpq_t object. */
! 588: #define mpq_numref(Q) (&((Q)->_mp_num))
! 589: #define mpq_denref(Q) (&((Q)->_mp_den))
! 590:
! 591: /* When using GCC, optimize certain common comparisons. */
! 592: #if defined (__GNUC__)
! 593: #define mpz_cmp_ui(Z,UI) \
! 594: (__builtin_constant_p (UI) && (UI) == 0 \
! 595: ? mpz_sgn (Z) : mpz_cmp_ui (Z,UI))
! 596: #define mpz_cmp_si(Z,UI) \
! 597: (__builtin_constant_p (UI) && (UI) == 0 ? mpz_sgn (Z) \
! 598: : __builtin_constant_p (UI) && (UI) > 0 ? mpz_cmp_ui (Z,UI) \
! 599: : mpz_cmp_si (Z,UI))
! 600: #define mpq_cmp_ui(Q,NUI,DUI) \
! 601: (__builtin_constant_p (NUI) && (NUI) == 0 \
! 602: ? mpq_sgn (Q) : mpq_cmp_ui (Q,NUI,DUI))
! 603: #endif
! 604:
! 605: #define mpn_divmod(qp,np,nsize,dp,dsize) mpn_divrem (qp,0,np,nsize,dp,dsize)
! 606: #if 0
! 607: #define mpn_divmod_1(qp,np,nsize,dlimb) mpn_divrem_1 (qp,0,np,nsize,dlimb)
! 608: #endif
! 609:
! 610: /* Compatibility with GMP 1. */
! 611: #define mpz_mdiv mpz_fdiv_q
! 612: #define mpz_mdivmod mpz_fdiv_qr
! 613: #define mpz_mmod mpz_fdiv_r
! 614: #define mpz_mdiv_ui mpz_fdiv_q_ui
! 615: #define mpz_mdivmod_ui(q,r,n,d) \
! 616: ((r == 0) ? mpz_fdiv_q_ui (q,n,d) : mpz_fdiv_qr_ui (q,r,n,d))
! 617: #define mpz_mmod_ui(r,n,d) \
! 618: ((r == 0) ? mpz_fdiv_ui (n,d) : mpz_fdiv_r_ui (r,n,d))
! 619:
! 620: /* Useful synonyms, but not quite compatible with GMP 1. */
! 621: #define mpz_div mpz_fdiv_q
! 622: #define mpz_divmod mpz_fdiv_qr
! 623: #define mpz_div_ui mpz_fdiv_q_ui
! 624: #define mpz_divmod_ui mpz_fdiv_qr_ui
! 625: #define mpz_mod_ui mpz_fdiv_r_ui
! 626: #define mpz_div_2exp mpz_fdiv_q_2exp
! 627: #define mpz_mod_2exp mpz_fdiv_r_2exp
! 628:
! 629: #define __GNU_MP_VERSION 2
! 630: #define __GNU_MP_VERSION_MINOR 0
! 631: #define __GMP_H__
! 632: #endif /* __GMP_H__ */
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to gmp.h CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / gmp