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