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 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 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 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 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 (_STDIO_H_) || defined (__STDIO_H__) || defined (H_STDIO) 166 #define _GMP_H_HAVE_FILE 1 167 #endif 168 169 void mp_set_memory_functions _PROTO ((void *(*) (size_t), 170 void *(*) (void *, size_t, size_t), 171 void (*) (void *, size_t))); 172 extern const int mp_bits_per_limb; 173 174 /**************** Integer (i.e. Z) routines. ****************/ 175 176 #if defined (__cplusplus) 177 extern "C" { 178 #endif 179 void *_mpz_realloc _PROTO ((mpz_ptr, mp_size_t)); 180 181 void mpz_abs _PROTO ((mpz_ptr, mpz_srcptr)); 182 void mpz_add _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 183 void mpz_add_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 184 void mpz_and _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 185 void mpz_array_init _PROTO ((mpz_ptr, mp_size_t, mp_size_t)); 186 void mpz_cdiv_q _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 187 unsigned long int mpz_cdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 188 void mpz_cdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr)); 189 unsigned long int mpz_cdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int)); 190 void mpz_cdiv_r _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 191 unsigned long int mpz_cdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 192 unsigned long int mpz_cdiv_ui _PROTO ((mpz_srcptr, unsigned long int)); 193 void mpz_clear _PROTO ((mpz_ptr)); 194 void mpz_clrbit _PROTO ((mpz_ptr, unsigned long int)); 195 int mpz_cmp _PROTO ((mpz_srcptr, mpz_srcptr)); 196 int mpz_cmp_si _PROTO ((mpz_srcptr, signed long int)); 197 int mpz_cmp_ui _PROTO ((mpz_srcptr, unsigned long int)); 198 void mpz_com _PROTO ((mpz_ptr, mpz_srcptr)); 199 void mpz_divexact _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 200 void mpz_fac_ui _PROTO ((mpz_ptr, unsigned long int)); 201 void mpz_fdiv_q _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 202 void mpz_fdiv_q_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 203 unsigned long int mpz_fdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 204 void mpz_fdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr)); 205 unsigned long int mpz_fdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int)); 206 void mpz_fdiv_r _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 207 void mpz_fdiv_r_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 208 unsigned long int mpz_fdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 209 unsigned long int mpz_fdiv_ui _PROTO ((mpz_srcptr, unsigned long int)); 210 void mpz_gcd _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 211 unsigned long int mpz_gcd_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 212 void mpz_gcdext _PROTO ((mpz_ptr, mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr)); 213 /* signed */ long int mpz_get_si _PROTO ((mpz_srcptr)); 214 char *mpz_get_str _PROTO ((char *, int, mpz_srcptr)); 215 unsigned long int mpz_get_ui _PROTO ((mpz_srcptr)); 216 mp_limb_t mpz_getlimbn _PROTO ((mpz_srcptr, mp_size_t)); 217 unsigned long int mpz_hamdist _PROTO ((mpz_srcptr, mpz_srcptr)); 218 void mpz_init _PROTO ((mpz_ptr)); 219 #ifdef _GMP_H_HAVE_FILE 220 size_t mpz_inp_binary _PROTO ((mpz_ptr, FILE *)); 221 size_t mpz_inp_raw _PROTO ((mpz_ptr, FILE *)); 222 size_t mpz_inp_str _PROTO ((mpz_ptr, FILE *, int)); 223 #endif 224 void mpz_init_set _PROTO ((mpz_ptr, mpz_srcptr)); 225 void mpz_init_set_d _PROTO ((mpz_ptr, double)); 226 void mpz_init_set_si _PROTO ((mpz_ptr, signed long int)); 227 int mpz_init_set_str _PROTO ((mpz_ptr, const char *, int)); 228 void mpz_init_set_ui _PROTO ((mpz_ptr, unsigned long int)); 229 int mpz_invert _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 230 void mpz_ior _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 231 int mpz_jacobi _PROTO ((mpz_srcptr, mpz_srcptr)); 232 int mpz_legendre _PROTO ((mpz_srcptr, mpz_srcptr)); 233 void mpz_mod _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 234 void mpz_mul _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 235 void mpz_mul_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 236 void mpz_mul_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 237 void mpz_neg _PROTO ((mpz_ptr, mpz_srcptr)); 238 #ifdef _GMP_H_HAVE_FILE 239 size_t mpz_out_binary _PROTO ((FILE *, mpz_srcptr)); 240 size_t mpz_out_raw _PROTO ((FILE *, mpz_srcptr)); 241 size_t mpz_out_str _PROTO ((FILE *, int, mpz_srcptr)); 242 #endif 243 int mpz_perfect_square_p _PROTO ((mpz_srcptr)); 244 unsigned long int mpz_popcount _PROTO ((mpz_srcptr)); 245 void mpz_pow_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 246 void mpz_powm _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr, mpz_srcptr)); 247 void mpz_powm_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int, mpz_srcptr)); 248 int mpz_probab_prime_p _PROTO ((mpz_srcptr, int)); 249 void mpz_random _PROTO ((mpz_ptr, mp_size_t)); 250 void mpz_random2 _PROTO ((mpz_ptr, mp_size_t)); 251 unsigned long int mpz_scan0 _PROTO ((mpz_srcptr, unsigned long int)); 252 unsigned long int mpz_scan1 _PROTO ((mpz_srcptr, unsigned long int)); 253 void mpz_set _PROTO ((mpz_ptr, mpz_srcptr)); 254 void mpz_set_d _PROTO ((mpz_ptr, double)); 255 void mpz_set_si _PROTO ((mpz_ptr, signed long int)); 256 int mpz_set_str _PROTO ((mpz_ptr, const char *, int)); 257 void mpz_set_ui _PROTO ((mpz_ptr, unsigned long int)); 258 void mpz_setbit _PROTO ((mpz_ptr, unsigned long int)); 259 size_t mpz_size _PROTO ((mpz_srcptr)); 260 size_t mpz_sizeinbase _PROTO ((mpz_srcptr, int)); 261 void mpz_sqrt _PROTO ((mpz_ptr, mpz_srcptr)); 262 void mpz_sqrtrem _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr)); 263 void mpz_sub _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 264 void mpz_sub_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 265 void mpz_tdiv_q _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 266 void mpz_tdiv_q_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 267 void mpz_tdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 268 void mpz_tdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr)); 269 void mpz_tdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int)); 270 void mpz_tdiv_r _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr)); 271 void mpz_tdiv_r_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 272 void mpz_tdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int)); 273 void mpz_ui_pow_ui _PROTO ((mpz_ptr, unsigned long int, unsigned long int)); 274 275 /**************** Rational (i.e. Q) routines. ****************/ 276 277 void mpq_init _PROTO ((mpq_ptr)); 278 void mpq_clear _PROTO ((mpq_ptr)); 279 void mpq_set _PROTO ((mpq_ptr, mpq_srcptr)); 280 void mpq_set_ui _PROTO ((mpq_ptr, unsigned long int, unsigned long int)); 281 void mpq_set_si _PROTO ((mpq_ptr, signed long int, unsigned long int)); 282 void mpq_add _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr)); 283 void mpq_sub _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr)); 284 void mpq_mul _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr)); 285 void mpq_div _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr)); 286 void mpq_neg _PROTO ((mpq_ptr, mpq_srcptr)); 287 int mpq_cmp _PROTO ((mpq_srcptr, mpq_srcptr)); 288 int mpq_cmp_ui _PROTO ((mpq_srcptr, unsigned long int, unsigned long int)); 289 void mpq_inv _PROTO ((mpq_ptr, mpq_srcptr)); 290 void mpq_set_num _PROTO ((mpq_ptr, mpz_srcptr)); 291 void mpq_set_den _PROTO ((mpq_ptr, mpz_srcptr)); 292 void mpq_get_num _PROTO ((mpz_ptr, mpq_srcptr)); 293 void mpq_get_den _PROTO ((mpz_ptr, mpq_srcptr)); 294 double mpq_get_d _PROTO ((mpq_srcptr)); 295 void mpq_canonicalize _PROTO ((mpq_ptr)); 296 297 /**************** Float (i.e. F) routines. ****************/ 298 299 void mpf_abs _PROTO ((mpf_ptr, mpf_srcptr)); 300 void mpf_add _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr)); 301 void mpf_add_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int)); 302 void mpf_clear _PROTO ((mpf_ptr)); 303 int mpf_cmp _PROTO ((mpf_srcptr, mpf_srcptr)); 304 int mpf_cmp_si _PROTO ((mpf_srcptr, signed long int)); 305 int mpf_cmp_ui _PROTO ((mpf_srcptr, unsigned long int)); 306 void mpf_div _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr)); 307 void mpf_div_2exp _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int)); 308 void mpf_div_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int)); 309 void mpf_dump _PROTO ((mpf_srcptr)); 310 int mpf_eq _PROTO ((mpf_srcptr, mpf_srcptr, unsigned long int)); 311 unsigned long int mpf_get_prec _PROTO ((mpf_srcptr)); 312 char *mpf_get_str _PROTO ((char *, mp_exp_t *, int, size_t, mpf_srcptr)); 313 void mpf_init _PROTO ((mpf_ptr)); 314 void mpf_init2 _PROTO ((mpf_ptr, unsigned long int)); 315 #ifdef _GMP_H_HAVE_FILE 316 size_t mpf_inp_str _PROTO ((mpf_ptr, FILE *, int)); 317 #endif 318 void mpf_init_set _PROTO ((mpf_ptr, mpf_srcptr)); 319 void mpf_init_set_d _PROTO ((mpf_ptr, double)); 320 void mpf_init_set_si _PROTO ((mpf_ptr, signed long int)); 321 int mpf_init_set_str _PROTO ((mpf_ptr, char *, int)); 322 void mpf_init_set_ui _PROTO ((mpf_ptr, unsigned long int)); 323 void mpf_mul _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr)); 324 void mpf_mul_2exp _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int)); 325 void mpf_mul_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int)); 326 void mpf_neg _PROTO ((mpf_ptr, mpf_srcptr)); 327 #ifdef _GMP_H_HAVE_FILE 328 size_t mpf_out_str _PROTO ((FILE *, int, size_t, mpf_srcptr)); 329 #endif 330 void mpf_random2 _PROTO ((mpf_ptr, mp_size_t, mp_exp_t)); 331 void mpf_reldiff _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr)); 332 void mpf_set _PROTO ((mpf_ptr, mpf_srcptr)); 333 void mpf_set_d _PROTO ((mpf_ptr, double)); 334 void mpf_set_default_prec _PROTO ((unsigned long int)); 335 void mpf_set_prec _PROTO ((mpf_ptr, unsigned long int)); 336 void mpf_set_prec_raw _PROTO ((mpf_ptr, unsigned long int)); 337 void mpf_set_si _PROTO ((mpf_ptr, signed long int)); 338 int mpf_set_str _PROTO ((mpf_ptr, const char *, int)); 339 void mpf_set_ui _PROTO ((mpf_ptr, unsigned long int)); 340 size_t mpf_size _PROTO ((mpf_srcptr)); 341 void mpf_sqrt _PROTO ((mpf_ptr, mpf_srcptr)); 342 void mpf_sqrt_ui _PROTO ((mpf_ptr, unsigned long int)); 343 void mpf_sub _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr)); 344 void mpf_sub_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int)); 345 void mpf_ui_div _PROTO ((mpf_ptr, unsigned long int, mpf_srcptr)); 346 void mpf_ui_sub _PROTO ((mpf_ptr, unsigned long int, mpf_srcptr)); 347 #if defined (__cplusplus) 348 } 349 #endif 350 /************ Low level positive-integer (i.e. N) routines. ************/ 351 352 /* This is ugly, but we need to make usr calls reach the prefixed function. */ 353 #define mpn_add __MPN(add) 354 #define mpn_add_1 __MPN(add_1) 355 #define mpn_add_n __MPN(add_n) 356 #define mpn_addmul_1 __MPN(addmul_1) 357 #define mpn_bdivmod __MPN(bdivmod) 358 #define mpn_cmp __MPN(cmp) 359 #define mpn_divmod_1 __MPN(divmod_1) 360 #define mpn_divrem __MPN(divrem) 361 #define mpn_divrem_1 __MPN(divrem_1) 362 #define mpn_dump __MPN(dump) 363 #define mpn_gcd __MPN(gcd) 364 #define mpn_gcd_1 __MPN(gcd_1) 365 #define mpn_gcdext __MPN(gcdext) 366 #define mpn_get_str __MPN(get_str) 367 #define mpn_hamdist __MPN(hamdist) 368 #define mpn_lshift __MPN(lshift) 369 #define mpn_mod_1 __MPN(mod_1) 370 #define mpn_mul __MPN(mul) 371 #define mpn_mul_1 __MPN(mul_1) 372 #define mpn_mul_n __MPN(mul_n) 373 #define mpn_perfect_square_p __MPN(perfect_square_p) 374 #define mpn_popcount __MPN(popcount) 375 #define mpn_preinv_mod_1 __MPN(preinv_mod_1) 376 #define mpn_random2 __MPN(random2) 377 #define mpn_rshift __MPN(rshift) 378 #define mpn_scan0 __MPN(scan0) 379 #define mpn_scan1 __MPN(scan1) 380 #define mpn_set_str __MPN(set_str) 381 #define mpn_sqrtrem __MPN(sqrtrem) 382 #define mpn_sub __MPN(sub) 383 #define mpn_sub_1 __MPN(sub_1) 384 #define mpn_sub_n __MPN(sub_n) 385 #define mpn_submul_1 __MPN(submul_1) 386 #define mpn_udiv_w_sdiv __MPN(udiv_w_sdiv) 387 388 #if defined (__cplusplus) 389 extern "C" { 390 #endif 391 mp_limb_t mpn_add _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr,mp_size_t)); 392 mp_limb_t mpn_add_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t)); 393 mp_limb_t mpn_add_n _PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t)); 394 mp_limb_t mpn_addmul_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t)); 395 mp_limb_t mpn_bdivmod _PROTO ((mp_ptr, mp_ptr, mp_size_t, mp_srcptr, mp_size_t, unsigned long int)); 396 int mpn_cmp _PROTO ((mp_srcptr, mp_srcptr, mp_size_t)); 397 mp_limb_t mpn_divmod_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t)); 398 mp_limb_t mpn_divrem _PROTO ((mp_ptr, mp_size_t, mp_ptr, mp_size_t, mp_srcptr, mp_size_t)); 399 mp_limb_t mpn_divrem_1 _PROTO ((mp_ptr, mp_size_t, mp_srcptr, mp_size_t, mp_limb_t)); 400 void mpn_dump _PROTO ((mp_srcptr, mp_size_t)); 401 mp_size_t mpn_gcd _PROTO ((mp_ptr, mp_ptr, mp_size_t, mp_ptr, mp_size_t)); 402 mp_limb_t mpn_gcd_1 _PROTO ((mp_srcptr, mp_size_t, mp_limb_t)); 403 mp_size_t mpn_gcdext _PROTO ((mp_ptr, mp_ptr, mp_ptr, mp_size_t, mp_ptr, mp_size_t)); 404 size_t mpn_get_str _PROTO ((unsigned char *, int, mp_ptr, mp_size_t)); 405 unsigned long int mpn_hamdist _PROTO ((mp_srcptr, mp_srcptr, mp_size_t)); 406 mp_limb_t mpn_lshift _PROTO ((mp_ptr, mp_srcptr, mp_size_t, unsigned int)); 407 mp_limb_t mpn_mod_1 _PROTO ((mp_srcptr, mp_size_t, mp_limb_t)); 408 mp_limb_t mpn_mul _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t)); 409 mp_limb_t mpn_mul_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t)); 410 void mpn_mul_n _PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t)); 411 int mpn_perfect_square_p _PROTO ((mp_srcptr, mp_size_t)); 412 unsigned long int mpn_popcount _PROTO ((mp_srcptr, mp_size_t)); 413 mp_limb_t mpn_preinv_mod_1 _PROTO ((mp_srcptr, mp_size_t, mp_limb_t, mp_limb_t)); 414 void mpn_random2 _PROTO ((mp_ptr, mp_size_t)); 415 mp_limb_t mpn_rshift _PROTO ((mp_ptr, mp_srcptr, mp_size_t, unsigned int)); 416 unsigned long int mpn_scan0 _PROTO ((mp_srcptr, unsigned long int)); 417 unsigned long int mpn_scan1 _PROTO ((mp_srcptr, unsigned long int)); 418 mp_size_t mpn_set_str _PROTO ((mp_ptr, const unsigned char *, size_t, int)); 419 mp_size_t mpn_sqrtrem _PROTO ((mp_ptr, mp_ptr, mp_srcptr, mp_size_t)); 420 mp_limb_t mpn_sub _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr,mp_size_t)); 421 mp_limb_t mpn_sub_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t)); 422 mp_limb_t mpn_sub_n _PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t)); 423 mp_limb_t mpn_submul_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t)); 424 #if defined (__cplusplus) 425 } 426 #endif 427 428 #if defined (__GNUC__) || defined (_FORCE_INLINES) 429 _EXTERN_INLINE mp_limb_t 430 #if defined (__STDC__) || defined (__cplusplus) 431 mpn_add_1 (register mp_ptr res_ptr, 432 register mp_srcptr s1_ptr, 433 register mp_size_t s1_size, 434 register mp_limb_t s2_limb) 435 #else 436 mpn_add_1 (res_ptr, s1_ptr, s1_size, s2_limb) 437 register mp_ptr res_ptr; 438 register mp_srcptr s1_ptr; 439 register mp_size_t s1_size; 440 register mp_limb_t s2_limb; 441 #endif 442 { 443 register mp_limb_t x; 444 445 x = *s1_ptr++; 446 s2_limb = x + s2_limb; 447 *res_ptr++ = s2_limb; 448 if (s2_limb < x) 449 { 450 while (--s1_size != 0) 451 { 452 x = *s1_ptr++ + 1; 453 *res_ptr++ = x; 454 if (x != 0) 455 goto fin; 456 } 457 458 return 1; 459 } 460 461 fin: 462 if (res_ptr != s1_ptr) 463 { 464 mp_size_t i; 465 for (i = 0; i < s1_size - 1; i++) 466 res_ptr[i] = s1_ptr[i]; 467 } 468 return 0; 469 } 470 471 _EXTERN_INLINE mp_limb_t 472 #if defined (__STDC__) || defined (__cplusplus) 473 mpn_add (register mp_ptr res_ptr, 474 register mp_srcptr s1_ptr, 475 register mp_size_t s1_size, 476 register mp_srcptr s2_ptr, 477 register mp_size_t s2_size) 478 #else 479 mpn_add (res_ptr, s1_ptr, s1_size, s2_ptr, s2_size) 480 register mp_ptr res_ptr; 481 register mp_srcptr s1_ptr; 482 register mp_size_t s1_size; 483 register mp_srcptr s2_ptr; 484 register mp_size_t s2_size; 485 #endif 486 { 487 mp_limb_t cy_limb = 0; 488 489 if (s2_size != 0) 490 cy_limb = mpn_add_n (res_ptr, s1_ptr, s2_ptr, s2_size); 491 492 if (s1_size - s2_size != 0) 493 cy_limb = mpn_add_1 (res_ptr + s2_size, 494 s1_ptr + s2_size, 495 s1_size - s2_size, 496 cy_limb); 497 return cy_limb; 498 } 499 500 _EXTERN_INLINE mp_limb_t 501 #if defined (__STDC__) || defined (__cplusplus) 502 mpn_sub_1 (register mp_ptr res_ptr, 503 register mp_srcptr s1_ptr, 504 register mp_size_t s1_size, 505 register mp_limb_t s2_limb) 506 #else 507 mpn_sub_1 (res_ptr, s1_ptr, s1_size, s2_limb) 508 register mp_ptr res_ptr; 509 register mp_srcptr s1_ptr; 510 register mp_size_t s1_size; 511 register mp_limb_t s2_limb; 512 #endif 513 { 514 register mp_limb_t x; 515 516 x = *s1_ptr++; 517 s2_limb = x - s2_limb; 518 *res_ptr++ = s2_limb; 519 if (s2_limb > x) 520 { 521 while (--s1_size != 0) 522 { 523 x = *s1_ptr++; 524 *res_ptr++ = x - 1; 525 if (x != 0) 526 goto fin; 527 } 528 529 return 1; 530 } 531 532 fin: 533 if (res_ptr != s1_ptr) 534 { 535 mp_size_t i; 536 for (i = 0; i < s1_size - 1; i++) 537 res_ptr[i] = s1_ptr[i]; 538 } 539 return 0; 540 } 541 542 _EXTERN_INLINE mp_limb_t 543 #if defined (__STDC__) || defined (__cplusplus) 544 mpn_sub (register mp_ptr res_ptr, 545 register mp_srcptr s1_ptr, 546 register mp_size_t s1_size, 547 register mp_srcptr s2_ptr, 548 register mp_size_t s2_size) 549 #else 550 mpn_sub (res_ptr, s1_ptr, s1_size, s2_ptr, s2_size) 551 register mp_ptr res_ptr; 552 register mp_srcptr s1_ptr; 553 register mp_size_t s1_size; 554 register mp_srcptr s2_ptr; 555 register mp_size_t s2_size; 556 #endif 557 { 558 mp_limb_t cy_limb = 0; 559 560 if (s2_size != 0) 561 cy_limb = mpn_sub_n (res_ptr, s1_ptr, s2_ptr, s2_size); 562 563 if (s1_size - s2_size != 0) 564 cy_limb = mpn_sub_1 (res_ptr + s2_size, 565 s1_ptr + s2_size, 566 s1_size - s2_size, 567 cy_limb); 568 return cy_limb; 569 } 570 #endif /* __GNUC__ */ 571 572 /* Allow faster testing for negative, zero, and positive. */ 573 #define mpz_sgn(Z) ((Z)->_mp_size < 0 ? -1 : (Z)->_mp_size > 0) 574 #define mpf_sgn(F) ((F)->_mp_size < 0 ? -1 : (F)->_mp_size > 0) 575 #define mpq_sgn(Q) ((Q)->_mp_num._mp_size < 0 ? -1 : (Q)->_mp_num._mp_size > 0) 576 577 /* Allow direct user access to numerator and denominator of a mpq_t object. */ 578 #define mpq_numref(Q) (&((Q)->_mp_num)) 579 #define mpq_denref(Q) (&((Q)->_mp_den)) 580 581 /* When using GCC, optimize certain common comparisons. */ 582 #if defined (__GNUC__) 583 #define mpz_cmp_ui(Z,UI) \ 584 (__builtin_constant_p (UI) && (UI) == 0 \ 585 ? mpz_sgn (Z) : mpz_cmp_ui (Z,UI)) 586 #define mpz_cmp_si(Z,UI) \ 587 (__builtin_constant_p (UI) && (UI) == 0 ? mpz_sgn (Z) \ 588 : __builtin_constant_p (UI) && (UI) > 0 ? mpz_cmp_ui (Z,UI) \ 589 : mpz_cmp_si (Z,UI)) 590 #define mpq_cmp_ui(Q,NUI,DUI) \ 591 (__builtin_constant_p (NUI) && (NUI) == 0 \ 592 ? mpq_sgn (Q) : mpq_cmp_ui (Q,NUI,DUI)) 593 #endif 594 595 #define mpn_divmod(qp,np,nsize,dp,dsize) mpn_divrem (qp,0,np,nsize,dp,dsize) 596 #if 0 597 #define mpn_divmod_1(qp,np,nsize,dlimb) mpn_divrem_1 (qp,0,np,nsize,dlimb) 598 #endif 599 600 /* Compatibility with GMP 1. */ 601 #define mpz_mdiv mpz_fdiv_q 602 #define mpz_mdivmod mpz_fdiv_qr 603 #define mpz_mmod mpz_fdiv_r 604 #define mpz_mdiv_ui mpz_fdiv_q_ui 605 #define mpz_mdivmod_ui(q,r,n,d) \ 606 ((r == 0) ? mpz_fdiv_q_ui (q,n,d) : mpz_fdiv_qr_ui (q,r,n,d)) 607 #define mpz_mmod_ui(r,n,d) \ 608 ((r == 0) ? mpz_fdiv_ui (n,d) : mpz_fdiv_r_ui (r,n,d)) 609 610 /* Useful synonyms, but not quite compatible with GMP 1. */ 611 #define mpz_div mpz_fdiv_q 612 #define mpz_divmod mpz_fdiv_qr 613 #define mpz_div_ui mpz_fdiv_q_ui 614 #define mpz_divmod_ui mpz_fdiv_qr_ui 615 #define mpz_mod_ui mpz_fdiv_r_ui 616 #define mpz_div_2exp mpz_fdiv_q_2exp 617 #define mpz_mod_2exp mpz_fdiv_r_2exp 618 619 #define __GNU_MP_VERSION 2 620 #define __GNU_MP_VERSION_MINOR 0 621 #define __GMP_H__ 622 #endif /* __GMP_H__ */ 623