1 /* Definitions for GNU multiple precision functions. -*- mode: c -*- 2 3 Copyright 1991, 1993-1997, 1999-2016, 2020 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 either: 9 10 * the GNU Lesser General Public License as published by the Free 11 Software Foundation; either version 3 of the License, or (at your 12 option) any later version. 13 14 or 15 16 * the GNU General Public License as published by the Free Software 17 Foundation; either version 2 of the License, or (at your option) any 18 later version. 19 20 or both in parallel, as here. 21 22 The GNU MP Library is distributed in the hope that it will be useful, but 23 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 24 or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 25 for more details. 26 27 You should have received copies of the GNU General Public License and the 28 GNU Lesser General Public License along with the GNU MP Library. If not, 29 see https://www.gnu.org/licenses/. */ 30 31 #ifndef __GMP_H__ 32 33 #if defined (__cplusplus) 34 #include <iosfwd> /* for std::istream, std::ostream, std::string */ 35 #include <cstdio> 36 #endif 37 38 39 /* Instantiated by configure. */ 40 #if ! defined (__GMP_WITHIN_CONFIGURE) 41 #define __GMP_HAVE_HOST_CPU_FAMILY_power 0 42 #define __GMP_HAVE_HOST_CPU_FAMILY_powerpc 1 43 #define GMP_LIMB_BITS 32 44 #define GMP_NAIL_BITS 0 45 #endif 46 #define GMP_NUMB_BITS (GMP_LIMB_BITS - GMP_NAIL_BITS) 47 #define GMP_NUMB_MASK ((~ __GMP_CAST (mp_limb_t, 0)) >> GMP_NAIL_BITS) 48 #define GMP_NUMB_MAX GMP_NUMB_MASK 49 #define GMP_NAIL_MASK (~ GMP_NUMB_MASK) 50 51 52 #ifndef __GNU_MP__ 53 #define __GNU_MP__ 6 54 55 #include <stddef.h> /* for size_t */ 56 #include <limits.h> 57 58 /* Instantiated by configure. */ 59 #if ! defined (__GMP_WITHIN_CONFIGURE) 60 /* #undef _LONG_LONG_LIMB */ 61 #define __GMP_LIBGMP_DLL 0 62 #endif 63 64 65 /* __GMP_DECLSPEC supports Windows DLL versions of libgmp, and is empty in 66 all other circumstances. 67 68 When compiling objects for libgmp, __GMP_DECLSPEC is an export directive, 69 or when compiling for an application it's an import directive. The two 70 cases are differentiated by __GMP_WITHIN_GMP defined by the GMP Makefiles 71 (and not defined from an application). 72 73 __GMP_DECLSPEC_XX is similarly used for libgmpxx. __GMP_WITHIN_GMPXX 74 indicates when building libgmpxx, and in that case libgmpxx functions are 75 exports, but libgmp functions which might get called are imports. 76 77 Libtool DLL_EXPORT define is not used. 78 79 There's no attempt to support GMP built both static and DLL. Doing so 80 would mean applications would have to tell us which of the two is going 81 to be used when linking, and that seems very tedious and error prone if 82 using GMP by hand, and equally tedious from a package since autoconf and 83 automake don't give much help. 84 85 __GMP_DECLSPEC is required on all documented global functions and 86 variables, the various internals in gmp-impl.h etc can be left unadorned. 87 But internals used by the test programs or speed measuring programs 88 should have __GMP_DECLSPEC, and certainly constants or variables must 89 have it or the wrong address will be resolved. 90 91 In gcc __declspec can go at either the start or end of a prototype. 92 93 In Microsoft C __declspec must go at the start, or after the type like 94 void __declspec(...) *foo()". There's no __dllexport or anything to 95 guard against someone foolish #defining dllexport. _export used to be 96 available, but no longer. 97 98 In Borland C _export still exists, but needs to go after the type, like 99 "void _export foo();". Would have to change the __GMP_DECLSPEC syntax to 100 make use of that. Probably more trouble than it's worth. */ 101 102 #if defined (__GNUC__) 103 #define __GMP_DECLSPEC_EXPORT __declspec(__dllexport__) 104 #define __GMP_DECLSPEC_IMPORT __declspec(__dllimport__) 105 #endif 106 #if defined (_MSC_VER) || defined (__BORLANDC__) 107 #define __GMP_DECLSPEC_EXPORT __declspec(dllexport) 108 #define __GMP_DECLSPEC_IMPORT __declspec(dllimport) 109 #endif 110 #ifdef __WATCOMC__ 111 #define __GMP_DECLSPEC_EXPORT __export 112 #define __GMP_DECLSPEC_IMPORT __import 113 #endif 114 #ifdef __IBMC__ 115 #define __GMP_DECLSPEC_EXPORT _Export 116 #define __GMP_DECLSPEC_IMPORT _Import 117 #endif 118 119 #if __GMP_LIBGMP_DLL 120 #ifdef __GMP_WITHIN_GMP 121 /* compiling to go into a DLL libgmp */ 122 #define __GMP_DECLSPEC __GMP_DECLSPEC_EXPORT 123 #else 124 /* compiling to go into an application which will link to a DLL libgmp */ 125 #define __GMP_DECLSPEC __GMP_DECLSPEC_IMPORT 126 #endif 127 #else 128 /* all other cases */ 129 #define __GMP_DECLSPEC 130 #endif 131 132 133 #ifdef __GMP_SHORT_LIMB 134 typedef unsigned int mp_limb_t; 135 typedef int mp_limb_signed_t; 136 #else 137 #ifdef _LONG_LONG_LIMB 138 typedef unsigned long long int mp_limb_t; 139 typedef long long int mp_limb_signed_t; 140 #else 141 typedef unsigned long int mp_limb_t; 142 typedef long int mp_limb_signed_t; 143 #endif 144 #endif 145 typedef unsigned long int mp_bitcnt_t; 146 147 /* For reference, note that the name __mpz_struct gets into C++ mangled 148 function names, which means although the "__" suggests an internal, we 149 must leave this name for binary compatibility. */ 150 typedef struct 151 { 152 int _mp_alloc; /* Number of *limbs* allocated and pointed 153 to by the _mp_d field. */ 154 int _mp_size; /* abs(_mp_size) is the number of limbs the 155 last field points to. If _mp_size is 156 negative this is a negative number. */ 157 mp_limb_t *_mp_d; /* Pointer to the limbs. */ 158 } __mpz_struct; 159 160 #endif /* __GNU_MP__ */ 161 162 163 typedef __mpz_struct MP_INT; /* gmp 1 source compatibility */ 164 typedef __mpz_struct mpz_t[1]; 165 166 typedef mp_limb_t * mp_ptr; 167 typedef const mp_limb_t * mp_srcptr; 168 #if defined (_CRAY) && ! defined (_CRAYMPP) 169 /* plain `int' is much faster (48 bits) */ 170 #define __GMP_MP_SIZE_T_INT 1 171 typedef int mp_size_t; 172 typedef int mp_exp_t; 173 #else 174 #define __GMP_MP_SIZE_T_INT 0 175 typedef long int mp_size_t; 176 typedef long int mp_exp_t; 177 #endif 178 179 typedef struct 180 { 181 __mpz_struct _mp_num; 182 __mpz_struct _mp_den; 183 } __mpq_struct; 184 185 typedef __mpq_struct MP_RAT; /* gmp 1 source compatibility */ 186 typedef __mpq_struct mpq_t[1]; 187 188 typedef struct 189 { 190 int _mp_prec; /* Max precision, in number of `mp_limb_t's. 191 Set by mpf_init and modified by 192 mpf_set_prec. The area pointed to by the 193 _mp_d field contains `prec' + 1 limbs. */ 194 int _mp_size; /* abs(_mp_size) is the number of limbs the 195 last field points to. If _mp_size is 196 negative this is a negative number. */ 197 mp_exp_t _mp_exp; /* Exponent, in the base of `mp_limb_t'. */ 198 mp_limb_t *_mp_d; /* Pointer to the limbs. */ 199 } __mpf_struct; 200 201 /* typedef __mpf_struct MP_FLOAT; */ 202 typedef __mpf_struct mpf_t[1]; 203 204 /* Available random number generation algorithms. */ 205 typedef enum 206 { 207 GMP_RAND_ALG_DEFAULT = 0, 208 GMP_RAND_ALG_LC = GMP_RAND_ALG_DEFAULT /* Linear congruential. */ 209 } gmp_randalg_t; 210 211 /* Random state struct. */ 212 typedef struct 213 { 214 mpz_t _mp_seed; /* _mp_d member points to state of the generator. */ 215 gmp_randalg_t _mp_alg; /* Currently unused. */ 216 union { 217 void *_mp_lc; /* Pointer to function pointers structure. */ 218 } _mp_algdata; 219 } __gmp_randstate_struct; 220 typedef __gmp_randstate_struct gmp_randstate_t[1]; 221 222 /* Types for function declarations in gmp files. */ 223 /* ??? Should not pollute user name space with these ??? */ 224 typedef const __mpz_struct *mpz_srcptr; 225 typedef __mpz_struct *mpz_ptr; 226 typedef const __mpf_struct *mpf_srcptr; 227 typedef __mpf_struct *mpf_ptr; 228 typedef const __mpq_struct *mpq_srcptr; 229 typedef __mpq_struct *mpq_ptr; 230 231 232 #if __GMP_LIBGMP_DLL 233 #ifdef __GMP_WITHIN_GMPXX 234 /* compiling to go into a DLL libgmpxx */ 235 #define __GMP_DECLSPEC_XX __GMP_DECLSPEC_EXPORT 236 #else 237 /* compiling to go into a application which will link to a DLL libgmpxx */ 238 #define __GMP_DECLSPEC_XX __GMP_DECLSPEC_IMPORT 239 #endif 240 #else 241 /* all other cases */ 242 #define __GMP_DECLSPEC_XX 243 #endif 244 245 246 #ifndef __MPN 247 #define __MPN(x) __gmpn_##x 248 #endif 249 250 /* For reference, "defined(EOF)" cannot be used here. In g++ 2.95.4, 251 <iostream> defines EOF but not FILE. */ 252 #if defined (FILE) \ 253 || defined (H_STDIO) \ 254 || defined (_H_STDIO) /* AIX */ \ 255 || defined (_STDIO_H) /* glibc, Sun, SCO */ \ 256 || defined (_STDIO_H_) /* BSD, OSF */ \ 257 || defined (__STDIO_H) /* Borland */ \ 258 || defined (__STDIO_H__) /* IRIX */ \ 259 || defined (_STDIO_INCLUDED) /* HPUX */ \ 260 || defined (__dj_include_stdio_h_) /* DJGPP */ \ 261 || defined (_FILE_DEFINED) /* Microsoft */ \ 262 || defined (__STDIO__) /* Apple MPW MrC */ \ 263 || defined (_MSL_STDIO_H) /* Metrowerks */ \ 264 || defined (_STDIO_H_INCLUDED) /* QNX4 */ \ 265 || defined (_ISO_STDIO_ISO_H) /* Sun C++ */ \ 266 || defined (__STDIO_LOADED) /* VMS */ \ 267 || defined (__DEFINED_FILE) /* musl */ 268 #define _GMP_H_HAVE_FILE 1 269 #endif 270 271 /* In ISO C, if a prototype involving "struct obstack *" is given without 272 that structure defined, then the struct is scoped down to just the 273 prototype, causing a conflict if it's subsequently defined for real. So 274 only give prototypes if we've got obstack.h. */ 275 #if defined (_OBSTACK_H) /* glibc <obstack.h> */ 276 #define _GMP_H_HAVE_OBSTACK 1 277 #endif 278 279 /* The prototypes for gmp_vprintf etc are provided only if va_list is defined, 280 via an application having included <stdarg.h>. Usually va_list is a typedef 281 so can't be tested directly, but C99 specifies that va_start is a macro. 282 283 <stdio.h> will define some sort of va_list for vprintf and vfprintf, but 284 let's not bother trying to use that since it's not standard and since 285 application uses for gmp_vprintf etc will almost certainly require the 286 whole <stdarg.h> anyway. */ 287 288 #ifdef va_start 289 #define _GMP_H_HAVE_VA_LIST 1 290 #endif 291 292 /* Test for gcc >= maj.min, as per __GNUC_PREREQ in glibc */ 293 #if defined (__GNUC__) && defined (__GNUC_MINOR__) 294 #define __GMP_GNUC_PREREQ(maj, min) \ 295 ((__GNUC__ << 16) + __GNUC_MINOR__ >= ((maj) << 16) + (min)) 296 #else 297 #define __GMP_GNUC_PREREQ(maj, min) 0 298 #endif 299 300 /* "pure" is in gcc 2.96 and up, see "(gcc)Function Attributes". Basically 301 it means a function does nothing but examine its arguments and memory 302 (global or via arguments) to generate a return value, but changes nothing 303 and has no side-effects. __GMP_NO_ATTRIBUTE_CONST_PURE lets 304 tune/common.c etc turn this off when trying to write timing loops. */ 305 #if __GMP_GNUC_PREREQ (2,96) && ! defined (__GMP_NO_ATTRIBUTE_CONST_PURE) 306 #define __GMP_ATTRIBUTE_PURE __attribute__ ((__pure__)) 307 #else 308 #define __GMP_ATTRIBUTE_PURE 309 #endif 310 311 312 /* __GMP_CAST allows us to use static_cast in C++, so our macros are clean 313 to "g++ -Wold-style-cast". 314 315 Casts in "extern inline" code within an extern "C" block don't induce 316 these warnings, so __GMP_CAST only needs to be used on documented 317 macros. */ 318 319 #ifdef __cplusplus 320 #define __GMP_CAST(type, expr) (static_cast<type> (expr)) 321 #else 322 #define __GMP_CAST(type, expr) ((type) (expr)) 323 #endif 324 325 326 /* An empty "throw ()" means the function doesn't throw any C++ exceptions, 327 this can save some stack frame info in applications. 328 329 Currently it's given only on functions which never divide-by-zero etc, 330 don't allocate memory, and are expected to never need to allocate memory. 331 This leaves open the possibility of a C++ throw from a future GMP 332 exceptions scheme. 333 334 mpz_set_ui etc are omitted to leave open the lazy allocation scheme 335 described in doc/tasks.html. mpz_get_d etc are omitted to leave open 336 exceptions for float overflows. 337 338 Note that __GMP_NOTHROW must be given on any inlines the same as on their 339 prototypes (for g++ at least, where they're used together). Note also 340 that g++ 3.0 demands that __GMP_NOTHROW is before other attributes like 341 __GMP_ATTRIBUTE_PURE. */ 342 343 #if defined (__cplusplus) 344 #if __cplusplus >= 201103L 345 #define __GMP_NOTHROW noexcept 346 #else 347 #define __GMP_NOTHROW throw () 348 #endif 349 #else 350 #define __GMP_NOTHROW 351 #endif 352 353 354 /* PORTME: What other compilers have a useful "extern inline"? "static 355 inline" would be an acceptable substitute if the compiler (or linker) 356 discards unused statics. */ 357 358 /* gcc has __inline__ in all modes, including strict ansi. Give a prototype 359 for an inline too, so as to correctly specify "dllimport" on windows, in 360 case the function is called rather than inlined. 361 GCC 4.3 and above with -std=c99 or -std=gnu99 implements ISO C99 362 inline semantics, unless -fgnu89-inline is used. */ 363 #ifdef __GNUC__ 364 #if (defined __GNUC_STDC_INLINE__) || (__GNUC__ == 4 && __GNUC_MINOR__ == 2) \ 365 || (defined __GNUC_GNU_INLINE__ && defined __cplusplus) 366 #define __GMP_EXTERN_INLINE extern __inline__ __attribute__ ((__gnu_inline__)) 367 #else 368 #define __GMP_EXTERN_INLINE extern __inline__ 369 #endif 370 #define __GMP_INLINE_PROTOTYPES 1 371 #endif 372 373 /* DEC C (eg. version 5.9) supports "static __inline foo()", even in -std1 374 strict ANSI mode. Inlining is done even when not optimizing (ie. -O0 375 mode, which is the default), but an unnecessary local copy of foo is 376 emitted unless -O is used. "extern __inline" is accepted, but the 377 "extern" appears to be ignored, ie. it becomes a plain global function 378 but which is inlined within its file. Don't know if all old versions of 379 DEC C supported __inline, but as a start let's do the right thing for 380 current versions. */ 381 #ifdef __DECC 382 #define __GMP_EXTERN_INLINE static __inline 383 #endif 384 385 /* SCO OpenUNIX 8 cc supports "static inline foo()" but not in -Xc strict 386 ANSI mode (__STDC__ is 1 in that mode). Inlining only actually takes 387 place under -O. Without -O "foo" seems to be emitted whether it's used 388 or not, which is wasteful. "extern inline foo()" isn't useful, the 389 "extern" is apparently ignored, so foo is inlined if possible but also 390 emitted as a global, which causes multiple definition errors when 391 building a shared libgmp. */ 392 #ifdef __SCO_VERSION__ 393 #if __SCO_VERSION__ > 400000000 && __STDC__ != 1 \ 394 && ! defined (__GMP_EXTERN_INLINE) 395 #define __GMP_EXTERN_INLINE static inline 396 #endif 397 #endif 398 399 /* Microsoft's C compiler accepts __inline */ 400 #ifdef _MSC_VER 401 #define __GMP_EXTERN_INLINE __inline 402 #endif 403 404 /* Recent enough Sun C compilers want "inline" */ 405 #if defined (__SUNPRO_C) && __SUNPRO_C >= 0x560 \ 406 && ! defined (__GMP_EXTERN_INLINE) 407 #define __GMP_EXTERN_INLINE inline 408 #endif 409 410 /* Somewhat older Sun C compilers want "static inline" */ 411 #if defined (__SUNPRO_C) && __SUNPRO_C >= 0x540 \ 412 && ! defined (__GMP_EXTERN_INLINE) 413 #define __GMP_EXTERN_INLINE static inline 414 #endif 415 416 417 /* C++ always has "inline" and since it's a normal feature the linker should 418 discard duplicate non-inlined copies, or if it doesn't then that's a 419 problem for everyone, not just GMP. */ 420 #if defined (__cplusplus) && ! defined (__GMP_EXTERN_INLINE) 421 #define __GMP_EXTERN_INLINE inline 422 #endif 423 424 /* Don't do any inlining within a configure run, since if the compiler ends 425 up emitting copies of the code into the object file it can end up 426 demanding the various support routines (like mpn_popcount) for linking, 427 making the "alloca" test and perhaps others fail. And on hppa ia64 a 428 pre-release gcc 3.2 was seen not respecting the "extern" in "extern 429 __inline__", triggering this problem too. */ 430 #if defined (__GMP_WITHIN_CONFIGURE) && ! __GMP_WITHIN_CONFIGURE_INLINE 431 #undef __GMP_EXTERN_INLINE 432 #endif 433 434 /* By default, don't give a prototype when there's going to be an inline 435 version. Note in particular that Cray C++ objects to the combination of 436 prototype and inline. */ 437 #ifdef __GMP_EXTERN_INLINE 438 #ifndef __GMP_INLINE_PROTOTYPES 439 #define __GMP_INLINE_PROTOTYPES 0 440 #endif 441 #else 442 #define __GMP_INLINE_PROTOTYPES 1 443 #endif 444 445 446 #define __GMP_ABS(x) ((x) >= 0 ? (x) : -(x)) 447 #define __GMP_MAX(h,i) ((h) > (i) ? (h) : (i)) 448 449 450 /* __builtin_expect is in gcc 3.0, and not in 2.95. */ 451 #if __GMP_GNUC_PREREQ (3,0) 452 #define __GMP_LIKELY(cond) __builtin_expect ((cond) != 0, 1) 453 #define __GMP_UNLIKELY(cond) __builtin_expect ((cond) != 0, 0) 454 #else 455 #define __GMP_LIKELY(cond) (cond) 456 #define __GMP_UNLIKELY(cond) (cond) 457 #endif 458 459 #ifdef _CRAY 460 #define __GMP_CRAY_Pragma(str) _Pragma (str) 461 #else 462 #define __GMP_CRAY_Pragma(str) 463 #endif 464 465 466 /* Allow direct user access to numerator and denominator of an mpq_t object. */ 467 #define mpq_numref(Q) (&((Q)->_mp_num)) 468 #define mpq_denref(Q) (&((Q)->_mp_den)) 469 470 471 #if defined (__cplusplus) 472 extern "C" { 473 using std::FILE; 474 #endif 475 476 #define mp_set_memory_functions __gmp_set_memory_functions 477 __GMP_DECLSPEC void mp_set_memory_functions (void *(*) (size_t), 478 void *(*) (void *, size_t, size_t), 479 void (*) (void *, size_t)) __GMP_NOTHROW; 480 481 #define mp_get_memory_functions __gmp_get_memory_functions 482 __GMP_DECLSPEC void mp_get_memory_functions (void *(**) (size_t), 483 void *(**) (void *, size_t, size_t), 484 void (**) (void *, size_t)) __GMP_NOTHROW; 485 486 #define mp_bits_per_limb __gmp_bits_per_limb 487 __GMP_DECLSPEC extern const int mp_bits_per_limb; 488 489 #define gmp_errno __gmp_errno 490 __GMP_DECLSPEC extern int gmp_errno; 491 492 #define gmp_version __gmp_version 493 __GMP_DECLSPEC extern const char * const gmp_version; 494 495 496 /**************** Random number routines. ****************/ 497 498 /* obsolete */ 499 #define gmp_randinit __gmp_randinit 500 __GMP_DECLSPEC void gmp_randinit (gmp_randstate_t, gmp_randalg_t, ...); 501 502 #define gmp_randinit_default __gmp_randinit_default 503 __GMP_DECLSPEC void gmp_randinit_default (gmp_randstate_t); 504 505 #define gmp_randinit_lc_2exp __gmp_randinit_lc_2exp 506 __GMP_DECLSPEC void gmp_randinit_lc_2exp (gmp_randstate_t, mpz_srcptr, unsigned long int, mp_bitcnt_t); 507 508 #define gmp_randinit_lc_2exp_size __gmp_randinit_lc_2exp_size 509 __GMP_DECLSPEC int gmp_randinit_lc_2exp_size (gmp_randstate_t, mp_bitcnt_t); 510 511 #define gmp_randinit_mt __gmp_randinit_mt 512 __GMP_DECLSPEC void gmp_randinit_mt (gmp_randstate_t); 513 514 #define gmp_randinit_set __gmp_randinit_set 515 __GMP_DECLSPEC void gmp_randinit_set (gmp_randstate_t, const __gmp_randstate_struct *); 516 517 #define gmp_randseed __gmp_randseed 518 __GMP_DECLSPEC void gmp_randseed (gmp_randstate_t, mpz_srcptr); 519 520 #define gmp_randseed_ui __gmp_randseed_ui 521 __GMP_DECLSPEC void gmp_randseed_ui (gmp_randstate_t, unsigned long int); 522 523 #define gmp_randclear __gmp_randclear 524 __GMP_DECLSPEC void gmp_randclear (gmp_randstate_t); 525 526 #define gmp_urandomb_ui __gmp_urandomb_ui 527 __GMP_DECLSPEC unsigned long gmp_urandomb_ui (gmp_randstate_t, unsigned long); 528 529 #define gmp_urandomm_ui __gmp_urandomm_ui 530 __GMP_DECLSPEC unsigned long gmp_urandomm_ui (gmp_randstate_t, unsigned long); 531 532 533 /**************** Formatted output routines. ****************/ 534 535 #define gmp_asprintf __gmp_asprintf 536 __GMP_DECLSPEC int gmp_asprintf (char **, const char *, ...); 537 538 #define gmp_fprintf __gmp_fprintf 539 #ifdef _GMP_H_HAVE_FILE 540 __GMP_DECLSPEC int gmp_fprintf (FILE *, const char *, ...); 541 #endif 542 543 #define gmp_obstack_printf __gmp_obstack_printf 544 #if defined (_GMP_H_HAVE_OBSTACK) 545 __GMP_DECLSPEC int gmp_obstack_printf (struct obstack *, const char *, ...); 546 #endif 547 548 #define gmp_obstack_vprintf __gmp_obstack_vprintf 549 #if defined (_GMP_H_HAVE_OBSTACK) && defined (_GMP_H_HAVE_VA_LIST) 550 __GMP_DECLSPEC int gmp_obstack_vprintf (struct obstack *, const char *, va_list); 551 #endif 552 553 #define gmp_printf __gmp_printf 554 __GMP_DECLSPEC int gmp_printf (const char *, ...); 555 556 #define gmp_snprintf __gmp_snprintf 557 __GMP_DECLSPEC int gmp_snprintf (char *, size_t, const char *, ...); 558 559 #define gmp_sprintf __gmp_sprintf 560 __GMP_DECLSPEC int gmp_sprintf (char *, const char *, ...); 561 562 #define gmp_vasprintf __gmp_vasprintf 563 #if defined (_GMP_H_HAVE_VA_LIST) 564 __GMP_DECLSPEC int gmp_vasprintf (char **, const char *, va_list); 565 #endif 566 567 #define gmp_vfprintf __gmp_vfprintf 568 #if defined (_GMP_H_HAVE_FILE) && defined (_GMP_H_HAVE_VA_LIST) 569 __GMP_DECLSPEC int gmp_vfprintf (FILE *, const char *, va_list); 570 #endif 571 572 #define gmp_vprintf __gmp_vprintf 573 #if defined (_GMP_H_HAVE_VA_LIST) 574 __GMP_DECLSPEC int gmp_vprintf (const char *, va_list); 575 #endif 576 577 #define gmp_vsnprintf __gmp_vsnprintf 578 #if defined (_GMP_H_HAVE_VA_LIST) 579 __GMP_DECLSPEC int gmp_vsnprintf (char *, size_t, const char *, va_list); 580 #endif 581 582 #define gmp_vsprintf __gmp_vsprintf 583 #if defined (_GMP_H_HAVE_VA_LIST) 584 __GMP_DECLSPEC int gmp_vsprintf (char *, const char *, va_list); 585 #endif 586 587 588 /**************** Formatted input routines. ****************/ 589 590 #define gmp_fscanf __gmp_fscanf 591 #ifdef _GMP_H_HAVE_FILE 592 __GMP_DECLSPEC int gmp_fscanf (FILE *, const char *, ...); 593 #endif 594 595 #define gmp_scanf __gmp_scanf 596 __GMP_DECLSPEC int gmp_scanf (const char *, ...); 597 598 #define gmp_sscanf __gmp_sscanf 599 __GMP_DECLSPEC int gmp_sscanf (const char *, const char *, ...); 600 601 #define gmp_vfscanf __gmp_vfscanf 602 #if defined (_GMP_H_HAVE_FILE) && defined (_GMP_H_HAVE_VA_LIST) 603 __GMP_DECLSPEC int gmp_vfscanf (FILE *, const char *, va_list); 604 #endif 605 606 #define gmp_vscanf __gmp_vscanf 607 #if defined (_GMP_H_HAVE_VA_LIST) 608 __GMP_DECLSPEC int gmp_vscanf (const char *, va_list); 609 #endif 610 611 #define gmp_vsscanf __gmp_vsscanf 612 #if defined (_GMP_H_HAVE_VA_LIST) 613 __GMP_DECLSPEC int gmp_vsscanf (const char *, const char *, va_list); 614 #endif 615 616 617 /**************** Integer (i.e. Z) routines. ****************/ 618 619 #define _mpz_realloc __gmpz_realloc 620 #define mpz_realloc __gmpz_realloc 621 __GMP_DECLSPEC void *_mpz_realloc (mpz_ptr, mp_size_t); 622 623 #define mpz_abs __gmpz_abs 624 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_abs) 625 __GMP_DECLSPEC void mpz_abs (mpz_ptr, mpz_srcptr); 626 #endif 627 628 #define mpz_add __gmpz_add 629 __GMP_DECLSPEC void mpz_add (mpz_ptr, mpz_srcptr, mpz_srcptr); 630 631 #define mpz_add_ui __gmpz_add_ui 632 __GMP_DECLSPEC void mpz_add_ui (mpz_ptr, mpz_srcptr, unsigned long int); 633 634 #define mpz_addmul __gmpz_addmul 635 __GMP_DECLSPEC void mpz_addmul (mpz_ptr, mpz_srcptr, mpz_srcptr); 636 637 #define mpz_addmul_ui __gmpz_addmul_ui 638 __GMP_DECLSPEC void mpz_addmul_ui (mpz_ptr, mpz_srcptr, unsigned long int); 639 640 #define mpz_and __gmpz_and 641 __GMP_DECLSPEC void mpz_and (mpz_ptr, mpz_srcptr, mpz_srcptr); 642 643 #define mpz_array_init __gmpz_array_init 644 __GMP_DECLSPEC void mpz_array_init (mpz_ptr, mp_size_t, mp_size_t); 645 646 #define mpz_bin_ui __gmpz_bin_ui 647 __GMP_DECLSPEC void mpz_bin_ui (mpz_ptr, mpz_srcptr, unsigned long int); 648 649 #define mpz_bin_uiui __gmpz_bin_uiui 650 __GMP_DECLSPEC void mpz_bin_uiui (mpz_ptr, unsigned long int, unsigned long int); 651 652 #define mpz_cdiv_q __gmpz_cdiv_q 653 __GMP_DECLSPEC void mpz_cdiv_q (mpz_ptr, mpz_srcptr, mpz_srcptr); 654 655 #define mpz_cdiv_q_2exp __gmpz_cdiv_q_2exp 656 __GMP_DECLSPEC void mpz_cdiv_q_2exp (mpz_ptr, mpz_srcptr, mp_bitcnt_t); 657 658 #define mpz_cdiv_q_ui __gmpz_cdiv_q_ui 659 __GMP_DECLSPEC unsigned long int mpz_cdiv_q_ui (mpz_ptr, mpz_srcptr, unsigned long int); 660 661 #define mpz_cdiv_qr __gmpz_cdiv_qr 662 __GMP_DECLSPEC void mpz_cdiv_qr (mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr); 663 664 #define mpz_cdiv_qr_ui __gmpz_cdiv_qr_ui 665 __GMP_DECLSPEC unsigned long int mpz_cdiv_qr_ui (mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int); 666 667 #define mpz_cdiv_r __gmpz_cdiv_r 668 __GMP_DECLSPEC void mpz_cdiv_r (mpz_ptr, mpz_srcptr, mpz_srcptr); 669 670 #define mpz_cdiv_r_2exp __gmpz_cdiv_r_2exp 671 __GMP_DECLSPEC void mpz_cdiv_r_2exp (mpz_ptr, mpz_srcptr, mp_bitcnt_t); 672 673 #define mpz_cdiv_r_ui __gmpz_cdiv_r_ui 674 __GMP_DECLSPEC unsigned long int mpz_cdiv_r_ui (mpz_ptr, mpz_srcptr, unsigned long int); 675 676 #define mpz_cdiv_ui __gmpz_cdiv_ui 677 __GMP_DECLSPEC unsigned long int mpz_cdiv_ui (mpz_srcptr, unsigned long int) __GMP_ATTRIBUTE_PURE; 678 679 #define mpz_clear __gmpz_clear 680 __GMP_DECLSPEC void mpz_clear (mpz_ptr); 681 682 #define mpz_clears __gmpz_clears 683 __GMP_DECLSPEC void mpz_clears (mpz_ptr, ...); 684 685 #define mpz_clrbit __gmpz_clrbit 686 __GMP_DECLSPEC void mpz_clrbit (mpz_ptr, mp_bitcnt_t); 687 688 #define mpz_cmp __gmpz_cmp 689 __GMP_DECLSPEC int mpz_cmp (mpz_srcptr, mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 690 691 #define mpz_cmp_d __gmpz_cmp_d 692 __GMP_DECLSPEC int mpz_cmp_d (mpz_srcptr, double) __GMP_ATTRIBUTE_PURE; 693 694 #define _mpz_cmp_si __gmpz_cmp_si 695 __GMP_DECLSPEC int _mpz_cmp_si (mpz_srcptr, signed long int) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 696 697 #define _mpz_cmp_ui __gmpz_cmp_ui 698 __GMP_DECLSPEC int _mpz_cmp_ui (mpz_srcptr, unsigned long int) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 699 700 #define mpz_cmpabs __gmpz_cmpabs 701 __GMP_DECLSPEC int mpz_cmpabs (mpz_srcptr, mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 702 703 #define mpz_cmpabs_d __gmpz_cmpabs_d 704 __GMP_DECLSPEC int mpz_cmpabs_d (mpz_srcptr, double) __GMP_ATTRIBUTE_PURE; 705 706 #define mpz_cmpabs_ui __gmpz_cmpabs_ui 707 __GMP_DECLSPEC int mpz_cmpabs_ui (mpz_srcptr, unsigned long int) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 708 709 #define mpz_com __gmpz_com 710 __GMP_DECLSPEC void mpz_com (mpz_ptr, mpz_srcptr); 711 712 #define mpz_combit __gmpz_combit 713 __GMP_DECLSPEC void mpz_combit (mpz_ptr, mp_bitcnt_t); 714 715 #define mpz_congruent_p __gmpz_congruent_p 716 __GMP_DECLSPEC int mpz_congruent_p (mpz_srcptr, mpz_srcptr, mpz_srcptr) __GMP_ATTRIBUTE_PURE; 717 718 #define mpz_congruent_2exp_p __gmpz_congruent_2exp_p 719 __GMP_DECLSPEC int mpz_congruent_2exp_p (mpz_srcptr, mpz_srcptr, mp_bitcnt_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 720 721 #define mpz_congruent_ui_p __gmpz_congruent_ui_p 722 __GMP_DECLSPEC int mpz_congruent_ui_p (mpz_srcptr, unsigned long, unsigned long) __GMP_ATTRIBUTE_PURE; 723 724 #define mpz_divexact __gmpz_divexact 725 __GMP_DECLSPEC void mpz_divexact (mpz_ptr, mpz_srcptr, mpz_srcptr); 726 727 #define mpz_divexact_ui __gmpz_divexact_ui 728 __GMP_DECLSPEC void mpz_divexact_ui (mpz_ptr, mpz_srcptr, unsigned long); 729 730 #define mpz_divisible_p __gmpz_divisible_p 731 __GMP_DECLSPEC int mpz_divisible_p (mpz_srcptr, mpz_srcptr) __GMP_ATTRIBUTE_PURE; 732 733 #define mpz_divisible_ui_p __gmpz_divisible_ui_p 734 __GMP_DECLSPEC int mpz_divisible_ui_p (mpz_srcptr, unsigned long) __GMP_ATTRIBUTE_PURE; 735 736 #define mpz_divisible_2exp_p __gmpz_divisible_2exp_p 737 __GMP_DECLSPEC int mpz_divisible_2exp_p (mpz_srcptr, mp_bitcnt_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 738 739 #define mpz_dump __gmpz_dump 740 __GMP_DECLSPEC void mpz_dump (mpz_srcptr); 741 742 #define mpz_export __gmpz_export 743 __GMP_DECLSPEC void *mpz_export (void *, size_t *, int, size_t, int, size_t, mpz_srcptr); 744 745 #define mpz_fac_ui __gmpz_fac_ui 746 __GMP_DECLSPEC void mpz_fac_ui (mpz_ptr, unsigned long int); 747 748 #define mpz_2fac_ui __gmpz_2fac_ui 749 __GMP_DECLSPEC void mpz_2fac_ui (mpz_ptr, unsigned long int); 750 751 #define mpz_mfac_uiui __gmpz_mfac_uiui 752 __GMP_DECLSPEC void mpz_mfac_uiui (mpz_ptr, unsigned long int, unsigned long int); 753 754 #define mpz_primorial_ui __gmpz_primorial_ui 755 __GMP_DECLSPEC void mpz_primorial_ui (mpz_ptr, unsigned long int); 756 757 #define mpz_fdiv_q __gmpz_fdiv_q 758 __GMP_DECLSPEC void mpz_fdiv_q (mpz_ptr, mpz_srcptr, mpz_srcptr); 759 760 #define mpz_fdiv_q_2exp __gmpz_fdiv_q_2exp 761 __GMP_DECLSPEC void mpz_fdiv_q_2exp (mpz_ptr, mpz_srcptr, mp_bitcnt_t); 762 763 #define mpz_fdiv_q_ui __gmpz_fdiv_q_ui 764 __GMP_DECLSPEC unsigned long int mpz_fdiv_q_ui (mpz_ptr, mpz_srcptr, unsigned long int); 765 766 #define mpz_fdiv_qr __gmpz_fdiv_qr 767 __GMP_DECLSPEC void mpz_fdiv_qr (mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr); 768 769 #define mpz_fdiv_qr_ui __gmpz_fdiv_qr_ui 770 __GMP_DECLSPEC unsigned long int mpz_fdiv_qr_ui (mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int); 771 772 #define mpz_fdiv_r __gmpz_fdiv_r 773 __GMP_DECLSPEC void mpz_fdiv_r (mpz_ptr, mpz_srcptr, mpz_srcptr); 774 775 #define mpz_fdiv_r_2exp __gmpz_fdiv_r_2exp 776 __GMP_DECLSPEC void mpz_fdiv_r_2exp (mpz_ptr, mpz_srcptr, mp_bitcnt_t); 777 778 #define mpz_fdiv_r_ui __gmpz_fdiv_r_ui 779 __GMP_DECLSPEC unsigned long int mpz_fdiv_r_ui (mpz_ptr, mpz_srcptr, unsigned long int); 780 781 #define mpz_fdiv_ui __gmpz_fdiv_ui 782 __GMP_DECLSPEC unsigned long int mpz_fdiv_ui (mpz_srcptr, unsigned long int) __GMP_ATTRIBUTE_PURE; 783 784 #define mpz_fib_ui __gmpz_fib_ui 785 __GMP_DECLSPEC void mpz_fib_ui (mpz_ptr, unsigned long int); 786 787 #define mpz_fib2_ui __gmpz_fib2_ui 788 __GMP_DECLSPEC void mpz_fib2_ui (mpz_ptr, mpz_ptr, unsigned long int); 789 790 #define mpz_fits_sint_p __gmpz_fits_sint_p 791 __GMP_DECLSPEC int mpz_fits_sint_p (mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 792 793 #define mpz_fits_slong_p __gmpz_fits_slong_p 794 __GMP_DECLSPEC int mpz_fits_slong_p (mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 795 796 #define mpz_fits_sshort_p __gmpz_fits_sshort_p 797 __GMP_DECLSPEC int mpz_fits_sshort_p (mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 798 799 #define mpz_fits_uint_p __gmpz_fits_uint_p 800 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_fits_uint_p) 801 __GMP_DECLSPEC int mpz_fits_uint_p (mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 802 #endif 803 804 #define mpz_fits_ulong_p __gmpz_fits_ulong_p 805 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_fits_ulong_p) 806 __GMP_DECLSPEC int mpz_fits_ulong_p (mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 807 #endif 808 809 #define mpz_fits_ushort_p __gmpz_fits_ushort_p 810 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_fits_ushort_p) 811 __GMP_DECLSPEC int mpz_fits_ushort_p (mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 812 #endif 813 814 #define mpz_gcd __gmpz_gcd 815 __GMP_DECLSPEC void mpz_gcd (mpz_ptr, mpz_srcptr, mpz_srcptr); 816 817 #define mpz_gcd_ui __gmpz_gcd_ui 818 __GMP_DECLSPEC unsigned long int mpz_gcd_ui (mpz_ptr, mpz_srcptr, unsigned long int); 819 820 #define mpz_gcdext __gmpz_gcdext 821 __GMP_DECLSPEC void mpz_gcdext (mpz_ptr, mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr); 822 823 #define mpz_get_d __gmpz_get_d 824 __GMP_DECLSPEC double mpz_get_d (mpz_srcptr) __GMP_ATTRIBUTE_PURE; 825 826 #define mpz_get_d_2exp __gmpz_get_d_2exp 827 __GMP_DECLSPEC double mpz_get_d_2exp (signed long int *, mpz_srcptr); 828 829 #define mpz_get_si __gmpz_get_si 830 __GMP_DECLSPEC /* signed */ long int mpz_get_si (mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 831 832 #define mpz_get_str __gmpz_get_str 833 __GMP_DECLSPEC char *mpz_get_str (char *, int, mpz_srcptr); 834 835 #define mpz_get_ui __gmpz_get_ui 836 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_get_ui) 837 __GMP_DECLSPEC unsigned long int mpz_get_ui (mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 838 #endif 839 840 #define mpz_getlimbn __gmpz_getlimbn 841 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_getlimbn) 842 __GMP_DECLSPEC mp_limb_t mpz_getlimbn (mpz_srcptr, mp_size_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 843 #endif 844 845 #define mpz_hamdist __gmpz_hamdist 846 __GMP_DECLSPEC mp_bitcnt_t mpz_hamdist (mpz_srcptr, mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 847 848 #define mpz_import __gmpz_import 849 __GMP_DECLSPEC void mpz_import (mpz_ptr, size_t, int, size_t, int, size_t, const void *); 850 851 #define mpz_init __gmpz_init 852 __GMP_DECLSPEC void mpz_init (mpz_ptr) __GMP_NOTHROW; 853 854 #define mpz_init2 __gmpz_init2 855 __GMP_DECLSPEC void mpz_init2 (mpz_ptr, mp_bitcnt_t); 856 857 #define mpz_inits __gmpz_inits 858 __GMP_DECLSPEC void mpz_inits (mpz_ptr, ...) __GMP_NOTHROW; 859 860 #define mpz_init_set __gmpz_init_set 861 __GMP_DECLSPEC void mpz_init_set (mpz_ptr, mpz_srcptr); 862 863 #define mpz_init_set_d __gmpz_init_set_d 864 __GMP_DECLSPEC void mpz_init_set_d (mpz_ptr, double); 865 866 #define mpz_init_set_si __gmpz_init_set_si 867 __GMP_DECLSPEC void mpz_init_set_si (mpz_ptr, signed long int); 868 869 #define mpz_init_set_str __gmpz_init_set_str 870 __GMP_DECLSPEC int mpz_init_set_str (mpz_ptr, const char *, int); 871 872 #define mpz_init_set_ui __gmpz_init_set_ui 873 __GMP_DECLSPEC void mpz_init_set_ui (mpz_ptr, unsigned long int); 874 875 #define mpz_inp_raw __gmpz_inp_raw 876 #ifdef _GMP_H_HAVE_FILE 877 __GMP_DECLSPEC size_t mpz_inp_raw (mpz_ptr, FILE *); 878 #endif 879 880 #define mpz_inp_str __gmpz_inp_str 881 #ifdef _GMP_H_HAVE_FILE 882 __GMP_DECLSPEC size_t mpz_inp_str (mpz_ptr, FILE *, int); 883 #endif 884 885 #define mpz_invert __gmpz_invert 886 __GMP_DECLSPEC int mpz_invert (mpz_ptr, mpz_srcptr, mpz_srcptr); 887 888 #define mpz_ior __gmpz_ior 889 __GMP_DECLSPEC void mpz_ior (mpz_ptr, mpz_srcptr, mpz_srcptr); 890 891 #define mpz_jacobi __gmpz_jacobi 892 __GMP_DECLSPEC int mpz_jacobi (mpz_srcptr, mpz_srcptr) __GMP_ATTRIBUTE_PURE; 893 894 #define mpz_kronecker mpz_jacobi /* alias */ 895 896 #define mpz_kronecker_si __gmpz_kronecker_si 897 __GMP_DECLSPEC int mpz_kronecker_si (mpz_srcptr, long) __GMP_ATTRIBUTE_PURE; 898 899 #define mpz_kronecker_ui __gmpz_kronecker_ui 900 __GMP_DECLSPEC int mpz_kronecker_ui (mpz_srcptr, unsigned long) __GMP_ATTRIBUTE_PURE; 901 902 #define mpz_si_kronecker __gmpz_si_kronecker 903 __GMP_DECLSPEC int mpz_si_kronecker (long, mpz_srcptr) __GMP_ATTRIBUTE_PURE; 904 905 #define mpz_ui_kronecker __gmpz_ui_kronecker 906 __GMP_DECLSPEC int mpz_ui_kronecker (unsigned long, mpz_srcptr) __GMP_ATTRIBUTE_PURE; 907 908 #define mpz_lcm __gmpz_lcm 909 __GMP_DECLSPEC void mpz_lcm (mpz_ptr, mpz_srcptr, mpz_srcptr); 910 911 #define mpz_lcm_ui __gmpz_lcm_ui 912 __GMP_DECLSPEC void mpz_lcm_ui (mpz_ptr, mpz_srcptr, unsigned long); 913 914 #define mpz_legendre mpz_jacobi /* alias */ 915 916 #define mpz_lucnum_ui __gmpz_lucnum_ui 917 __GMP_DECLSPEC void mpz_lucnum_ui (mpz_ptr, unsigned long int); 918 919 #define mpz_lucnum2_ui __gmpz_lucnum2_ui 920 __GMP_DECLSPEC void mpz_lucnum2_ui (mpz_ptr, mpz_ptr, unsigned long int); 921 922 #define mpz_millerrabin __gmpz_millerrabin 923 __GMP_DECLSPEC int mpz_millerrabin (mpz_srcptr, int) __GMP_ATTRIBUTE_PURE; 924 925 #define mpz_mod __gmpz_mod 926 __GMP_DECLSPEC void mpz_mod (mpz_ptr, mpz_srcptr, mpz_srcptr); 927 928 #define mpz_mod_ui mpz_fdiv_r_ui /* same as fdiv_r because divisor unsigned */ 929 930 #define mpz_mul __gmpz_mul 931 __GMP_DECLSPEC void mpz_mul (mpz_ptr, mpz_srcptr, mpz_srcptr); 932 933 #define mpz_mul_2exp __gmpz_mul_2exp 934 __GMP_DECLSPEC void mpz_mul_2exp (mpz_ptr, mpz_srcptr, mp_bitcnt_t); 935 936 #define mpz_mul_si __gmpz_mul_si 937 __GMP_DECLSPEC void mpz_mul_si (mpz_ptr, mpz_srcptr, long int); 938 939 #define mpz_mul_ui __gmpz_mul_ui 940 __GMP_DECLSPEC void mpz_mul_ui (mpz_ptr, mpz_srcptr, unsigned long int); 941 942 #define mpz_neg __gmpz_neg 943 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_neg) 944 __GMP_DECLSPEC void mpz_neg (mpz_ptr, mpz_srcptr); 945 #endif 946 947 #define mpz_nextprime __gmpz_nextprime 948 __GMP_DECLSPEC void mpz_nextprime (mpz_ptr, mpz_srcptr); 949 950 #define mpz_out_raw __gmpz_out_raw 951 #ifdef _GMP_H_HAVE_FILE 952 __GMP_DECLSPEC size_t mpz_out_raw (FILE *, mpz_srcptr); 953 #endif 954 955 #define mpz_out_str __gmpz_out_str 956 #ifdef _GMP_H_HAVE_FILE 957 __GMP_DECLSPEC size_t mpz_out_str (FILE *, int, mpz_srcptr); 958 #endif 959 960 #define mpz_perfect_power_p __gmpz_perfect_power_p 961 __GMP_DECLSPEC int mpz_perfect_power_p (mpz_srcptr) __GMP_ATTRIBUTE_PURE; 962 963 #define mpz_perfect_square_p __gmpz_perfect_square_p 964 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_perfect_square_p) 965 __GMP_DECLSPEC int mpz_perfect_square_p (mpz_srcptr) __GMP_ATTRIBUTE_PURE; 966 #endif 967 968 #define mpz_popcount __gmpz_popcount 969 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_popcount) 970 __GMP_DECLSPEC mp_bitcnt_t mpz_popcount (mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 971 #endif 972 973 #define mpz_pow_ui __gmpz_pow_ui 974 __GMP_DECLSPEC void mpz_pow_ui (mpz_ptr, mpz_srcptr, unsigned long int); 975 976 #define mpz_powm __gmpz_powm 977 __GMP_DECLSPEC void mpz_powm (mpz_ptr, mpz_srcptr, mpz_srcptr, mpz_srcptr); 978 979 #define mpz_powm_sec __gmpz_powm_sec 980 __GMP_DECLSPEC void mpz_powm_sec (mpz_ptr, mpz_srcptr, mpz_srcptr, mpz_srcptr); 981 982 #define mpz_powm_ui __gmpz_powm_ui 983 __GMP_DECLSPEC void mpz_powm_ui (mpz_ptr, mpz_srcptr, unsigned long int, mpz_srcptr); 984 985 #define mpz_probab_prime_p __gmpz_probab_prime_p 986 __GMP_DECLSPEC int mpz_probab_prime_p (mpz_srcptr, int) __GMP_ATTRIBUTE_PURE; 987 988 #define mpz_random __gmpz_random 989 __GMP_DECLSPEC void mpz_random (mpz_ptr, mp_size_t); 990 991 #define mpz_random2 __gmpz_random2 992 __GMP_DECLSPEC void mpz_random2 (mpz_ptr, mp_size_t); 993 994 #define mpz_realloc2 __gmpz_realloc2 995 __GMP_DECLSPEC void mpz_realloc2 (mpz_ptr, mp_bitcnt_t); 996 997 #define mpz_remove __gmpz_remove 998 __GMP_DECLSPEC mp_bitcnt_t mpz_remove (mpz_ptr, mpz_srcptr, mpz_srcptr); 999 1000 #define mpz_root __gmpz_root 1001 __GMP_DECLSPEC int mpz_root (mpz_ptr, mpz_srcptr, unsigned long int); 1002 1003 #define mpz_rootrem __gmpz_rootrem 1004 __GMP_DECLSPEC void mpz_rootrem (mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int); 1005 1006 #define mpz_rrandomb __gmpz_rrandomb 1007 __GMP_DECLSPEC void mpz_rrandomb (mpz_ptr, gmp_randstate_t, mp_bitcnt_t); 1008 1009 #define mpz_scan0 __gmpz_scan0 1010 __GMP_DECLSPEC mp_bitcnt_t mpz_scan0 (mpz_srcptr, mp_bitcnt_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1011 1012 #define mpz_scan1 __gmpz_scan1 1013 __GMP_DECLSPEC mp_bitcnt_t mpz_scan1 (mpz_srcptr, mp_bitcnt_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1014 1015 #define mpz_set __gmpz_set 1016 __GMP_DECLSPEC void mpz_set (mpz_ptr, mpz_srcptr); 1017 1018 #define mpz_set_d __gmpz_set_d 1019 __GMP_DECLSPEC void mpz_set_d (mpz_ptr, double); 1020 1021 #define mpz_set_f __gmpz_set_f 1022 __GMP_DECLSPEC void mpz_set_f (mpz_ptr, mpf_srcptr); 1023 1024 #define mpz_set_q __gmpz_set_q 1025 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_set_q) 1026 __GMP_DECLSPEC void mpz_set_q (mpz_ptr, mpq_srcptr); 1027 #endif 1028 1029 #define mpz_set_si __gmpz_set_si 1030 __GMP_DECLSPEC void mpz_set_si (mpz_ptr, signed long int); 1031 1032 #define mpz_set_str __gmpz_set_str 1033 __GMP_DECLSPEC int mpz_set_str (mpz_ptr, const char *, int); 1034 1035 #define mpz_set_ui __gmpz_set_ui 1036 __GMP_DECLSPEC void mpz_set_ui (mpz_ptr, unsigned long int); 1037 1038 #define mpz_setbit __gmpz_setbit 1039 __GMP_DECLSPEC void mpz_setbit (mpz_ptr, mp_bitcnt_t); 1040 1041 #define mpz_size __gmpz_size 1042 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_size) 1043 __GMP_DECLSPEC size_t mpz_size (mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1044 #endif 1045 1046 #define mpz_sizeinbase __gmpz_sizeinbase 1047 __GMP_DECLSPEC size_t mpz_sizeinbase (mpz_srcptr, int) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1048 1049 #define mpz_sqrt __gmpz_sqrt 1050 __GMP_DECLSPEC void mpz_sqrt (mpz_ptr, mpz_srcptr); 1051 1052 #define mpz_sqrtrem __gmpz_sqrtrem 1053 __GMP_DECLSPEC void mpz_sqrtrem (mpz_ptr, mpz_ptr, mpz_srcptr); 1054 1055 #define mpz_sub __gmpz_sub 1056 __GMP_DECLSPEC void mpz_sub (mpz_ptr, mpz_srcptr, mpz_srcptr); 1057 1058 #define mpz_sub_ui __gmpz_sub_ui 1059 __GMP_DECLSPEC void mpz_sub_ui (mpz_ptr, mpz_srcptr, unsigned long int); 1060 1061 #define mpz_ui_sub __gmpz_ui_sub 1062 __GMP_DECLSPEC void mpz_ui_sub (mpz_ptr, unsigned long int, mpz_srcptr); 1063 1064 #define mpz_submul __gmpz_submul 1065 __GMP_DECLSPEC void mpz_submul (mpz_ptr, mpz_srcptr, mpz_srcptr); 1066 1067 #define mpz_submul_ui __gmpz_submul_ui 1068 __GMP_DECLSPEC void mpz_submul_ui (mpz_ptr, mpz_srcptr, unsigned long int); 1069 1070 #define mpz_swap __gmpz_swap 1071 __GMP_DECLSPEC void mpz_swap (mpz_ptr, mpz_ptr) __GMP_NOTHROW; 1072 1073 #define mpz_tdiv_ui __gmpz_tdiv_ui 1074 __GMP_DECLSPEC unsigned long int mpz_tdiv_ui (mpz_srcptr, unsigned long int) __GMP_ATTRIBUTE_PURE; 1075 1076 #define mpz_tdiv_q __gmpz_tdiv_q 1077 __GMP_DECLSPEC void mpz_tdiv_q (mpz_ptr, mpz_srcptr, mpz_srcptr); 1078 1079 #define mpz_tdiv_q_2exp __gmpz_tdiv_q_2exp 1080 __GMP_DECLSPEC void mpz_tdiv_q_2exp (mpz_ptr, mpz_srcptr, mp_bitcnt_t); 1081 1082 #define mpz_tdiv_q_ui __gmpz_tdiv_q_ui 1083 __GMP_DECLSPEC unsigned long int mpz_tdiv_q_ui (mpz_ptr, mpz_srcptr, unsigned long int); 1084 1085 #define mpz_tdiv_qr __gmpz_tdiv_qr 1086 __GMP_DECLSPEC void mpz_tdiv_qr (mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr); 1087 1088 #define mpz_tdiv_qr_ui __gmpz_tdiv_qr_ui 1089 __GMP_DECLSPEC unsigned long int mpz_tdiv_qr_ui (mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int); 1090 1091 #define mpz_tdiv_r __gmpz_tdiv_r 1092 __GMP_DECLSPEC void mpz_tdiv_r (mpz_ptr, mpz_srcptr, mpz_srcptr); 1093 1094 #define mpz_tdiv_r_2exp __gmpz_tdiv_r_2exp 1095 __GMP_DECLSPEC void mpz_tdiv_r_2exp (mpz_ptr, mpz_srcptr, mp_bitcnt_t); 1096 1097 #define mpz_tdiv_r_ui __gmpz_tdiv_r_ui 1098 __GMP_DECLSPEC unsigned long int mpz_tdiv_r_ui (mpz_ptr, mpz_srcptr, unsigned long int); 1099 1100 #define mpz_tstbit __gmpz_tstbit 1101 __GMP_DECLSPEC int mpz_tstbit (mpz_srcptr, mp_bitcnt_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1102 1103 #define mpz_ui_pow_ui __gmpz_ui_pow_ui 1104 __GMP_DECLSPEC void mpz_ui_pow_ui (mpz_ptr, unsigned long int, unsigned long int); 1105 1106 #define mpz_urandomb __gmpz_urandomb 1107 __GMP_DECLSPEC void mpz_urandomb (mpz_ptr, gmp_randstate_t, mp_bitcnt_t); 1108 1109 #define mpz_urandomm __gmpz_urandomm 1110 __GMP_DECLSPEC void mpz_urandomm (mpz_ptr, gmp_randstate_t, mpz_srcptr); 1111 1112 #define mpz_xor __gmpz_xor 1113 #define mpz_eor __gmpz_xor 1114 __GMP_DECLSPEC void mpz_xor (mpz_ptr, mpz_srcptr, mpz_srcptr); 1115 1116 #define mpz_limbs_read __gmpz_limbs_read 1117 __GMP_DECLSPEC mp_srcptr mpz_limbs_read (mpz_srcptr); 1118 1119 #define mpz_limbs_write __gmpz_limbs_write 1120 __GMP_DECLSPEC mp_ptr mpz_limbs_write (mpz_ptr, mp_size_t); 1121 1122 #define mpz_limbs_modify __gmpz_limbs_modify 1123 __GMP_DECLSPEC mp_ptr mpz_limbs_modify (mpz_ptr, mp_size_t); 1124 1125 #define mpz_limbs_finish __gmpz_limbs_finish 1126 __GMP_DECLSPEC void mpz_limbs_finish (mpz_ptr, mp_size_t); 1127 1128 #define mpz_roinit_n __gmpz_roinit_n 1129 __GMP_DECLSPEC mpz_srcptr mpz_roinit_n (mpz_ptr, mp_srcptr, mp_size_t); 1130 1131 #define MPZ_ROINIT_N(xp, xs) {{0, (xs),(xp) }} 1132 1133 /**************** Rational (i.e. Q) routines. ****************/ 1134 1135 #define mpq_abs __gmpq_abs 1136 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpq_abs) 1137 __GMP_DECLSPEC void mpq_abs (mpq_ptr, mpq_srcptr); 1138 #endif 1139 1140 #define mpq_add __gmpq_add 1141 __GMP_DECLSPEC void mpq_add (mpq_ptr, mpq_srcptr, mpq_srcptr); 1142 1143 #define mpq_canonicalize __gmpq_canonicalize 1144 __GMP_DECLSPEC void mpq_canonicalize (mpq_ptr); 1145 1146 #define mpq_clear __gmpq_clear 1147 __GMP_DECLSPEC void mpq_clear (mpq_ptr); 1148 1149 #define mpq_clears __gmpq_clears 1150 __GMP_DECLSPEC void mpq_clears (mpq_ptr, ...); 1151 1152 #define mpq_cmp __gmpq_cmp 1153 __GMP_DECLSPEC int mpq_cmp (mpq_srcptr, mpq_srcptr) __GMP_ATTRIBUTE_PURE; 1154 1155 #define _mpq_cmp_si __gmpq_cmp_si 1156 __GMP_DECLSPEC int _mpq_cmp_si (mpq_srcptr, long, unsigned long) __GMP_ATTRIBUTE_PURE; 1157 1158 #define _mpq_cmp_ui __gmpq_cmp_ui 1159 __GMP_DECLSPEC int _mpq_cmp_ui (mpq_srcptr, unsigned long int, unsigned long int) __GMP_ATTRIBUTE_PURE; 1160 1161 #define mpq_cmp_z __gmpq_cmp_z 1162 __GMP_DECLSPEC int mpq_cmp_z (mpq_srcptr, mpz_srcptr) __GMP_ATTRIBUTE_PURE; 1163 1164 #define mpq_div __gmpq_div 1165 __GMP_DECLSPEC void mpq_div (mpq_ptr, mpq_srcptr, mpq_srcptr); 1166 1167 #define mpq_div_2exp __gmpq_div_2exp 1168 __GMP_DECLSPEC void mpq_div_2exp (mpq_ptr, mpq_srcptr, mp_bitcnt_t); 1169 1170 #define mpq_equal __gmpq_equal 1171 __GMP_DECLSPEC int mpq_equal (mpq_srcptr, mpq_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1172 1173 #define mpq_get_num __gmpq_get_num 1174 __GMP_DECLSPEC void mpq_get_num (mpz_ptr, mpq_srcptr); 1175 1176 #define mpq_get_den __gmpq_get_den 1177 __GMP_DECLSPEC void mpq_get_den (mpz_ptr, mpq_srcptr); 1178 1179 #define mpq_get_d __gmpq_get_d 1180 __GMP_DECLSPEC double mpq_get_d (mpq_srcptr) __GMP_ATTRIBUTE_PURE; 1181 1182 #define mpq_get_str __gmpq_get_str 1183 __GMP_DECLSPEC char *mpq_get_str (char *, int, mpq_srcptr); 1184 1185 #define mpq_init __gmpq_init 1186 __GMP_DECLSPEC void mpq_init (mpq_ptr); 1187 1188 #define mpq_inits __gmpq_inits 1189 __GMP_DECLSPEC void mpq_inits (mpq_ptr, ...); 1190 1191 #define mpq_inp_str __gmpq_inp_str 1192 #ifdef _GMP_H_HAVE_FILE 1193 __GMP_DECLSPEC size_t mpq_inp_str (mpq_ptr, FILE *, int); 1194 #endif 1195 1196 #define mpq_inv __gmpq_inv 1197 __GMP_DECLSPEC void mpq_inv (mpq_ptr, mpq_srcptr); 1198 1199 #define mpq_mul __gmpq_mul 1200 __GMP_DECLSPEC void mpq_mul (mpq_ptr, mpq_srcptr, mpq_srcptr); 1201 1202 #define mpq_mul_2exp __gmpq_mul_2exp 1203 __GMP_DECLSPEC void mpq_mul_2exp (mpq_ptr, mpq_srcptr, mp_bitcnt_t); 1204 1205 #define mpq_neg __gmpq_neg 1206 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpq_neg) 1207 __GMP_DECLSPEC void mpq_neg (mpq_ptr, mpq_srcptr); 1208 #endif 1209 1210 #define mpq_out_str __gmpq_out_str 1211 #ifdef _GMP_H_HAVE_FILE 1212 __GMP_DECLSPEC size_t mpq_out_str (FILE *, int, mpq_srcptr); 1213 #endif 1214 1215 #define mpq_set __gmpq_set 1216 __GMP_DECLSPEC void mpq_set (mpq_ptr, mpq_srcptr); 1217 1218 #define mpq_set_d __gmpq_set_d 1219 __GMP_DECLSPEC void mpq_set_d (mpq_ptr, double); 1220 1221 #define mpq_set_den __gmpq_set_den 1222 __GMP_DECLSPEC void mpq_set_den (mpq_ptr, mpz_srcptr); 1223 1224 #define mpq_set_f __gmpq_set_f 1225 __GMP_DECLSPEC void mpq_set_f (mpq_ptr, mpf_srcptr); 1226 1227 #define mpq_set_num __gmpq_set_num 1228 __GMP_DECLSPEC void mpq_set_num (mpq_ptr, mpz_srcptr); 1229 1230 #define mpq_set_si __gmpq_set_si 1231 __GMP_DECLSPEC void mpq_set_si (mpq_ptr, signed long int, unsigned long int); 1232 1233 #define mpq_set_str __gmpq_set_str 1234 __GMP_DECLSPEC int mpq_set_str (mpq_ptr, const char *, int); 1235 1236 #define mpq_set_ui __gmpq_set_ui 1237 __GMP_DECLSPEC void mpq_set_ui (mpq_ptr, unsigned long int, unsigned long int); 1238 1239 #define mpq_set_z __gmpq_set_z 1240 __GMP_DECLSPEC void mpq_set_z (mpq_ptr, mpz_srcptr); 1241 1242 #define mpq_sub __gmpq_sub 1243 __GMP_DECLSPEC void mpq_sub (mpq_ptr, mpq_srcptr, mpq_srcptr); 1244 1245 #define mpq_swap __gmpq_swap 1246 __GMP_DECLSPEC void mpq_swap (mpq_ptr, mpq_ptr) __GMP_NOTHROW; 1247 1248 1249 /**************** Float (i.e. F) routines. ****************/ 1250 1251 #define mpf_abs __gmpf_abs 1252 __GMP_DECLSPEC void mpf_abs (mpf_ptr, mpf_srcptr); 1253 1254 #define mpf_add __gmpf_add 1255 __GMP_DECLSPEC void mpf_add (mpf_ptr, mpf_srcptr, mpf_srcptr); 1256 1257 #define mpf_add_ui __gmpf_add_ui 1258 __GMP_DECLSPEC void mpf_add_ui (mpf_ptr, mpf_srcptr, unsigned long int); 1259 #define mpf_ceil __gmpf_ceil 1260 __GMP_DECLSPEC void mpf_ceil (mpf_ptr, mpf_srcptr); 1261 1262 #define mpf_clear __gmpf_clear 1263 __GMP_DECLSPEC void mpf_clear (mpf_ptr); 1264 1265 #define mpf_clears __gmpf_clears 1266 __GMP_DECLSPEC void mpf_clears (mpf_ptr, ...); 1267 1268 #define mpf_cmp __gmpf_cmp 1269 __GMP_DECLSPEC int mpf_cmp (mpf_srcptr, mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1270 1271 #define mpf_cmp_z __gmpf_cmp_z 1272 __GMP_DECLSPEC int mpf_cmp_z (mpf_srcptr, mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1273 1274 #define mpf_cmp_d __gmpf_cmp_d 1275 __GMP_DECLSPEC int mpf_cmp_d (mpf_srcptr, double) __GMP_ATTRIBUTE_PURE; 1276 1277 #define mpf_cmp_si __gmpf_cmp_si 1278 __GMP_DECLSPEC int mpf_cmp_si (mpf_srcptr, signed long int) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1279 1280 #define mpf_cmp_ui __gmpf_cmp_ui 1281 __GMP_DECLSPEC int mpf_cmp_ui (mpf_srcptr, unsigned long int) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1282 1283 #define mpf_div __gmpf_div 1284 __GMP_DECLSPEC void mpf_div (mpf_ptr, mpf_srcptr, mpf_srcptr); 1285 1286 #define mpf_div_2exp __gmpf_div_2exp 1287 __GMP_DECLSPEC void mpf_div_2exp (mpf_ptr, mpf_srcptr, mp_bitcnt_t); 1288 1289 #define mpf_div_ui __gmpf_div_ui 1290 __GMP_DECLSPEC void mpf_div_ui (mpf_ptr, mpf_srcptr, unsigned long int); 1291 1292 #define mpf_dump __gmpf_dump 1293 __GMP_DECLSPEC void mpf_dump (mpf_srcptr); 1294 1295 #define mpf_eq __gmpf_eq 1296 __GMP_DECLSPEC int mpf_eq (mpf_srcptr, mpf_srcptr, mp_bitcnt_t) __GMP_ATTRIBUTE_PURE; 1297 1298 #define mpf_fits_sint_p __gmpf_fits_sint_p 1299 __GMP_DECLSPEC int mpf_fits_sint_p (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1300 1301 #define mpf_fits_slong_p __gmpf_fits_slong_p 1302 __GMP_DECLSPEC int mpf_fits_slong_p (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1303 1304 #define mpf_fits_sshort_p __gmpf_fits_sshort_p 1305 __GMP_DECLSPEC int mpf_fits_sshort_p (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1306 1307 #define mpf_fits_uint_p __gmpf_fits_uint_p 1308 __GMP_DECLSPEC int mpf_fits_uint_p (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1309 1310 #define mpf_fits_ulong_p __gmpf_fits_ulong_p 1311 __GMP_DECLSPEC int mpf_fits_ulong_p (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1312 1313 #define mpf_fits_ushort_p __gmpf_fits_ushort_p 1314 __GMP_DECLSPEC int mpf_fits_ushort_p (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1315 1316 #define mpf_floor __gmpf_floor 1317 __GMP_DECLSPEC void mpf_floor (mpf_ptr, mpf_srcptr); 1318 1319 #define mpf_get_d __gmpf_get_d 1320 __GMP_DECLSPEC double mpf_get_d (mpf_srcptr) __GMP_ATTRIBUTE_PURE; 1321 1322 #define mpf_get_d_2exp __gmpf_get_d_2exp 1323 __GMP_DECLSPEC double mpf_get_d_2exp (signed long int *, mpf_srcptr); 1324 1325 #define mpf_get_default_prec __gmpf_get_default_prec 1326 __GMP_DECLSPEC mp_bitcnt_t mpf_get_default_prec (void) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1327 1328 #define mpf_get_prec __gmpf_get_prec 1329 __GMP_DECLSPEC mp_bitcnt_t mpf_get_prec (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1330 1331 #define mpf_get_si __gmpf_get_si 1332 __GMP_DECLSPEC long mpf_get_si (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1333 1334 #define mpf_get_str __gmpf_get_str 1335 __GMP_DECLSPEC char *mpf_get_str (char *, mp_exp_t *, int, size_t, mpf_srcptr); 1336 1337 #define mpf_get_ui __gmpf_get_ui 1338 __GMP_DECLSPEC unsigned long mpf_get_ui (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1339 1340 #define mpf_init __gmpf_init 1341 __GMP_DECLSPEC void mpf_init (mpf_ptr); 1342 1343 #define mpf_init2 __gmpf_init2 1344 __GMP_DECLSPEC void mpf_init2 (mpf_ptr, mp_bitcnt_t); 1345 1346 #define mpf_inits __gmpf_inits 1347 __GMP_DECLSPEC void mpf_inits (mpf_ptr, ...); 1348 1349 #define mpf_init_set __gmpf_init_set 1350 __GMP_DECLSPEC void mpf_init_set (mpf_ptr, mpf_srcptr); 1351 1352 #define mpf_init_set_d __gmpf_init_set_d 1353 __GMP_DECLSPEC void mpf_init_set_d (mpf_ptr, double); 1354 1355 #define mpf_init_set_si __gmpf_init_set_si 1356 __GMP_DECLSPEC void mpf_init_set_si (mpf_ptr, signed long int); 1357 1358 #define mpf_init_set_str __gmpf_init_set_str 1359 __GMP_DECLSPEC int mpf_init_set_str (mpf_ptr, const char *, int); 1360 1361 #define mpf_init_set_ui __gmpf_init_set_ui 1362 __GMP_DECLSPEC void mpf_init_set_ui (mpf_ptr, unsigned long int); 1363 1364 #define mpf_inp_str __gmpf_inp_str 1365 #ifdef _GMP_H_HAVE_FILE 1366 __GMP_DECLSPEC size_t mpf_inp_str (mpf_ptr, FILE *, int); 1367 #endif 1368 1369 #define mpf_integer_p __gmpf_integer_p 1370 __GMP_DECLSPEC int mpf_integer_p (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1371 1372 #define mpf_mul __gmpf_mul 1373 __GMP_DECLSPEC void mpf_mul (mpf_ptr, mpf_srcptr, mpf_srcptr); 1374 1375 #define mpf_mul_2exp __gmpf_mul_2exp 1376 __GMP_DECLSPEC void mpf_mul_2exp (mpf_ptr, mpf_srcptr, mp_bitcnt_t); 1377 1378 #define mpf_mul_ui __gmpf_mul_ui 1379 __GMP_DECLSPEC void mpf_mul_ui (mpf_ptr, mpf_srcptr, unsigned long int); 1380 1381 #define mpf_neg __gmpf_neg 1382 __GMP_DECLSPEC void mpf_neg (mpf_ptr, mpf_srcptr); 1383 1384 #define mpf_out_str __gmpf_out_str 1385 #ifdef _GMP_H_HAVE_FILE 1386 __GMP_DECLSPEC size_t mpf_out_str (FILE *, int, size_t, mpf_srcptr); 1387 #endif 1388 1389 #define mpf_pow_ui __gmpf_pow_ui 1390 __GMP_DECLSPEC void mpf_pow_ui (mpf_ptr, mpf_srcptr, unsigned long int); 1391 1392 #define mpf_random2 __gmpf_random2 1393 __GMP_DECLSPEC void mpf_random2 (mpf_ptr, mp_size_t, mp_exp_t); 1394 1395 #define mpf_reldiff __gmpf_reldiff 1396 __GMP_DECLSPEC void mpf_reldiff (mpf_ptr, mpf_srcptr, mpf_srcptr); 1397 1398 #define mpf_set __gmpf_set 1399 __GMP_DECLSPEC void mpf_set (mpf_ptr, mpf_srcptr); 1400 1401 #define mpf_set_d __gmpf_set_d 1402 __GMP_DECLSPEC void mpf_set_d (mpf_ptr, double); 1403 1404 #define mpf_set_default_prec __gmpf_set_default_prec 1405 __GMP_DECLSPEC void mpf_set_default_prec (mp_bitcnt_t) __GMP_NOTHROW; 1406 1407 #define mpf_set_prec __gmpf_set_prec 1408 __GMP_DECLSPEC void mpf_set_prec (mpf_ptr, mp_bitcnt_t); 1409 1410 #define mpf_set_prec_raw __gmpf_set_prec_raw 1411 __GMP_DECLSPEC void mpf_set_prec_raw (mpf_ptr, mp_bitcnt_t) __GMP_NOTHROW; 1412 1413 #define mpf_set_q __gmpf_set_q 1414 __GMP_DECLSPEC void mpf_set_q (mpf_ptr, mpq_srcptr); 1415 1416 #define mpf_set_si __gmpf_set_si 1417 __GMP_DECLSPEC void mpf_set_si (mpf_ptr, signed long int); 1418 1419 #define mpf_set_str __gmpf_set_str 1420 __GMP_DECLSPEC int mpf_set_str (mpf_ptr, const char *, int); 1421 1422 #define mpf_set_ui __gmpf_set_ui 1423 __GMP_DECLSPEC void mpf_set_ui (mpf_ptr, unsigned long int); 1424 1425 #define mpf_set_z __gmpf_set_z 1426 __GMP_DECLSPEC void mpf_set_z (mpf_ptr, mpz_srcptr); 1427 1428 #define mpf_size __gmpf_size 1429 __GMP_DECLSPEC size_t mpf_size (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1430 1431 #define mpf_sqrt __gmpf_sqrt 1432 __GMP_DECLSPEC void mpf_sqrt (mpf_ptr, mpf_srcptr); 1433 1434 #define mpf_sqrt_ui __gmpf_sqrt_ui 1435 __GMP_DECLSPEC void mpf_sqrt_ui (mpf_ptr, unsigned long int); 1436 1437 #define mpf_sub __gmpf_sub 1438 __GMP_DECLSPEC void mpf_sub (mpf_ptr, mpf_srcptr, mpf_srcptr); 1439 1440 #define mpf_sub_ui __gmpf_sub_ui 1441 __GMP_DECLSPEC void mpf_sub_ui (mpf_ptr, mpf_srcptr, unsigned long int); 1442 1443 #define mpf_swap __gmpf_swap 1444 __GMP_DECLSPEC void mpf_swap (mpf_ptr, mpf_ptr) __GMP_NOTHROW; 1445 1446 #define mpf_trunc __gmpf_trunc 1447 __GMP_DECLSPEC void mpf_trunc (mpf_ptr, mpf_srcptr); 1448 1449 #define mpf_ui_div __gmpf_ui_div 1450 __GMP_DECLSPEC void mpf_ui_div (mpf_ptr, unsigned long int, mpf_srcptr); 1451 1452 #define mpf_ui_sub __gmpf_ui_sub 1453 __GMP_DECLSPEC void mpf_ui_sub (mpf_ptr, unsigned long int, mpf_srcptr); 1454 1455 #define mpf_urandomb __gmpf_urandomb 1456 __GMP_DECLSPEC void mpf_urandomb (mpf_t, gmp_randstate_t, mp_bitcnt_t); 1457 1458 1459 /************ Low level positive-integer (i.e. N) routines. ************/ 1460 1461 /* This is ugly, but we need to make user calls reach the prefixed function. */ 1462 1463 #define mpn_add __MPN(add) 1464 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpn_add) 1465 __GMP_DECLSPEC mp_limb_t mpn_add (mp_ptr, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t); 1466 #endif 1467 1468 #define mpn_add_1 __MPN(add_1) 1469 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpn_add_1) 1470 __GMP_DECLSPEC mp_limb_t mpn_add_1 (mp_ptr, mp_srcptr, mp_size_t, mp_limb_t) __GMP_NOTHROW; 1471 #endif 1472 1473 #define mpn_add_n __MPN(add_n) 1474 __GMP_DECLSPEC mp_limb_t mpn_add_n (mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); 1475 1476 #define mpn_addmul_1 __MPN(addmul_1) 1477 __GMP_DECLSPEC mp_limb_t mpn_addmul_1 (mp_ptr, mp_srcptr, mp_size_t, mp_limb_t); 1478 1479 #define mpn_cmp __MPN(cmp) 1480 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpn_cmp) 1481 __GMP_DECLSPEC int mpn_cmp (mp_srcptr, mp_srcptr, mp_size_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1482 #endif 1483 1484 #define mpn_zero_p __MPN(zero_p) 1485 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpn_zero_p) 1486 __GMP_DECLSPEC int mpn_zero_p (mp_srcptr, mp_size_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1487 #endif 1488 1489 #define mpn_divexact_1 __MPN(divexact_1) 1490 __GMP_DECLSPEC void mpn_divexact_1 (mp_ptr, mp_srcptr, mp_size_t, mp_limb_t); 1491 1492 #define mpn_divexact_by3(dst,src,size) \ 1493 mpn_divexact_by3c (dst, src, size, __GMP_CAST (mp_limb_t, 0)) 1494 1495 #define mpn_divexact_by3c __MPN(divexact_by3c) 1496 __GMP_DECLSPEC mp_limb_t mpn_divexact_by3c (mp_ptr, mp_srcptr, mp_size_t, mp_limb_t); 1497 1498 #define mpn_divmod_1(qp,np,nsize,dlimb) \ 1499 mpn_divrem_1 (qp, __GMP_CAST (mp_size_t, 0), np, nsize, dlimb) 1500 1501 #define mpn_divrem __MPN(divrem) 1502 __GMP_DECLSPEC mp_limb_t mpn_divrem (mp_ptr, mp_size_t, mp_ptr, mp_size_t, mp_srcptr, mp_size_t); 1503 1504 #define mpn_divrem_1 __MPN(divrem_1) 1505 __GMP_DECLSPEC mp_limb_t mpn_divrem_1 (mp_ptr, mp_size_t, mp_srcptr, mp_size_t, mp_limb_t); 1506 1507 #define mpn_divrem_2 __MPN(divrem_2) 1508 __GMP_DECLSPEC mp_limb_t mpn_divrem_2 (mp_ptr, mp_size_t, mp_ptr, mp_size_t, mp_srcptr); 1509 1510 #define mpn_div_qr_1 __MPN(div_qr_1) 1511 __GMP_DECLSPEC mp_limb_t mpn_div_qr_1 (mp_ptr, mp_limb_t *, mp_srcptr, mp_size_t, mp_limb_t); 1512 1513 #define mpn_div_qr_2 __MPN(div_qr_2) 1514 __GMP_DECLSPEC mp_limb_t mpn_div_qr_2 (mp_ptr, mp_ptr, mp_srcptr, mp_size_t, mp_srcptr); 1515 1516 #define mpn_gcd __MPN(gcd) 1517 __GMP_DECLSPEC mp_size_t mpn_gcd (mp_ptr, mp_ptr, mp_size_t, mp_ptr, mp_size_t); 1518 1519 #define mpn_gcd_11 __MPN(gcd_11) 1520 __GMP_DECLSPEC mp_limb_t mpn_gcd_11 (mp_limb_t, mp_limb_t) __GMP_ATTRIBUTE_PURE; 1521 1522 #define mpn_gcd_1 __MPN(gcd_1) 1523 __GMP_DECLSPEC mp_limb_t mpn_gcd_1 (mp_srcptr, mp_size_t, mp_limb_t) __GMP_ATTRIBUTE_PURE; 1524 1525 #define mpn_gcdext_1 __MPN(gcdext_1) 1526 __GMP_DECLSPEC mp_limb_t mpn_gcdext_1 (mp_limb_signed_t *, mp_limb_signed_t *, mp_limb_t, mp_limb_t); 1527 1528 #define mpn_gcdext __MPN(gcdext) 1529 __GMP_DECLSPEC mp_size_t mpn_gcdext (mp_ptr, mp_ptr, mp_size_t *, mp_ptr, mp_size_t, mp_ptr, mp_size_t); 1530 1531 #define mpn_get_str __MPN(get_str) 1532 __GMP_DECLSPEC size_t mpn_get_str (unsigned char *, int, mp_ptr, mp_size_t); 1533 1534 #define mpn_hamdist __MPN(hamdist) 1535 __GMP_DECLSPEC mp_bitcnt_t mpn_hamdist (mp_srcptr, mp_srcptr, mp_size_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1536 1537 #define mpn_lshift __MPN(lshift) 1538 __GMP_DECLSPEC mp_limb_t mpn_lshift (mp_ptr, mp_srcptr, mp_size_t, unsigned int); 1539 1540 #define mpn_mod_1 __MPN(mod_1) 1541 __GMP_DECLSPEC mp_limb_t mpn_mod_1 (mp_srcptr, mp_size_t, mp_limb_t) __GMP_ATTRIBUTE_PURE; 1542 1543 #define mpn_mul __MPN(mul) 1544 __GMP_DECLSPEC mp_limb_t mpn_mul (mp_ptr, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t); 1545 1546 #define mpn_mul_1 __MPN(mul_1) 1547 __GMP_DECLSPEC mp_limb_t mpn_mul_1 (mp_ptr, mp_srcptr, mp_size_t, mp_limb_t); 1548 1549 #define mpn_mul_n __MPN(mul_n) 1550 __GMP_DECLSPEC void mpn_mul_n (mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); 1551 1552 #define mpn_sqr __MPN(sqr) 1553 __GMP_DECLSPEC void mpn_sqr (mp_ptr, mp_srcptr, mp_size_t); 1554 1555 #define mpn_neg __MPN(neg) 1556 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpn_neg) 1557 __GMP_DECLSPEC mp_limb_t mpn_neg (mp_ptr, mp_srcptr, mp_size_t); 1558 #endif 1559 1560 #define mpn_com __MPN(com) 1561 __GMP_DECLSPEC void mpn_com (mp_ptr, mp_srcptr, mp_size_t); 1562 1563 #define mpn_perfect_square_p __MPN(perfect_square_p) 1564 __GMP_DECLSPEC int mpn_perfect_square_p (mp_srcptr, mp_size_t) __GMP_ATTRIBUTE_PURE; 1565 1566 #define mpn_perfect_power_p __MPN(perfect_power_p) 1567 __GMP_DECLSPEC int mpn_perfect_power_p (mp_srcptr, mp_size_t) __GMP_ATTRIBUTE_PURE; 1568 1569 #define mpn_popcount __MPN(popcount) 1570 __GMP_DECLSPEC mp_bitcnt_t mpn_popcount (mp_srcptr, mp_size_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE; 1571 1572 #define mpn_pow_1 __MPN(pow_1) 1573 __GMP_DECLSPEC mp_size_t mpn_pow_1 (mp_ptr, mp_srcptr, mp_size_t, mp_limb_t, mp_ptr); 1574 1575 /* undocumented now, but retained here for upward compatibility */ 1576 #define mpn_preinv_mod_1 __MPN(preinv_mod_1) 1577 __GMP_DECLSPEC mp_limb_t mpn_preinv_mod_1 (mp_srcptr, mp_size_t, mp_limb_t, mp_limb_t) __GMP_ATTRIBUTE_PURE; 1578 1579 #define mpn_random __MPN(random) 1580 __GMP_DECLSPEC void mpn_random (mp_ptr, mp_size_t); 1581 1582 #define mpn_random2 __MPN(random2) 1583 __GMP_DECLSPEC void mpn_random2 (mp_ptr, mp_size_t); 1584 1585 #define mpn_rshift __MPN(rshift) 1586 __GMP_DECLSPEC mp_limb_t mpn_rshift (mp_ptr, mp_srcptr, mp_size_t, unsigned int); 1587 1588 #define mpn_scan0 __MPN(scan0) 1589 __GMP_DECLSPEC mp_bitcnt_t mpn_scan0 (mp_srcptr, mp_bitcnt_t) __GMP_ATTRIBUTE_PURE; 1590 1591 #define mpn_scan1 __MPN(scan1) 1592 __GMP_DECLSPEC mp_bitcnt_t mpn_scan1 (mp_srcptr, mp_bitcnt_t) __GMP_ATTRIBUTE_PURE; 1593 1594 #define mpn_set_str __MPN(set_str) 1595 __GMP_DECLSPEC mp_size_t mpn_set_str (mp_ptr, const unsigned char *, size_t, int); 1596 1597 #define mpn_sizeinbase __MPN(sizeinbase) 1598 __GMP_DECLSPEC size_t mpn_sizeinbase (mp_srcptr, mp_size_t, int); 1599 1600 #define mpn_sqrtrem __MPN(sqrtrem) 1601 __GMP_DECLSPEC mp_size_t mpn_sqrtrem (mp_ptr, mp_ptr, mp_srcptr, mp_size_t); 1602 1603 #define mpn_sub __MPN(sub) 1604 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpn_sub) 1605 __GMP_DECLSPEC mp_limb_t mpn_sub (mp_ptr, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t); 1606 #endif 1607 1608 #define mpn_sub_1 __MPN(sub_1) 1609 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpn_sub_1) 1610 __GMP_DECLSPEC mp_limb_t mpn_sub_1 (mp_ptr, mp_srcptr, mp_size_t, mp_limb_t) __GMP_NOTHROW; 1611 #endif 1612 1613 #define mpn_sub_n __MPN(sub_n) 1614 __GMP_DECLSPEC mp_limb_t mpn_sub_n (mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); 1615 1616 #define mpn_submul_1 __MPN(submul_1) 1617 __GMP_DECLSPEC mp_limb_t mpn_submul_1 (mp_ptr, mp_srcptr, mp_size_t, mp_limb_t); 1618 1619 #define mpn_tdiv_qr __MPN(tdiv_qr) 1620 __GMP_DECLSPEC void mpn_tdiv_qr (mp_ptr, mp_ptr, mp_size_t, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t); 1621 1622 #define mpn_and_n __MPN(and_n) 1623 __GMP_DECLSPEC void mpn_and_n (mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); 1624 #define mpn_andn_n __MPN(andn_n) 1625 __GMP_DECLSPEC void mpn_andn_n (mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); 1626 #define mpn_nand_n __MPN(nand_n) 1627 __GMP_DECLSPEC void mpn_nand_n (mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); 1628 #define mpn_ior_n __MPN(ior_n) 1629 __GMP_DECLSPEC void mpn_ior_n (mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); 1630 #define mpn_iorn_n __MPN(iorn_n) 1631 __GMP_DECLSPEC void mpn_iorn_n (mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); 1632 #define mpn_nior_n __MPN(nior_n) 1633 __GMP_DECLSPEC void mpn_nior_n (mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); 1634 #define mpn_xor_n __MPN(xor_n) 1635 __GMP_DECLSPEC void mpn_xor_n (mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); 1636 #define mpn_xnor_n __MPN(xnor_n) 1637 __GMP_DECLSPEC void mpn_xnor_n (mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); 1638 1639 #define mpn_copyi __MPN(copyi) 1640 __GMP_DECLSPEC void mpn_copyi (mp_ptr, mp_srcptr, mp_size_t); 1641 #define mpn_copyd __MPN(copyd) 1642 __GMP_DECLSPEC void mpn_copyd (mp_ptr, mp_srcptr, mp_size_t); 1643 #define mpn_zero __MPN(zero) 1644 __GMP_DECLSPEC void mpn_zero (mp_ptr, mp_size_t); 1645 1646 #define mpn_cnd_add_n __MPN(cnd_add_n) 1647 __GMP_DECLSPEC mp_limb_t mpn_cnd_add_n (mp_limb_t, mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); 1648 #define mpn_cnd_sub_n __MPN(cnd_sub_n) 1649 __GMP_DECLSPEC mp_limb_t mpn_cnd_sub_n (mp_limb_t, mp_ptr, mp_srcptr, mp_srcptr, mp_size_t); 1650 1651 #define mpn_sec_add_1 __MPN(sec_add_1) 1652 __GMP_DECLSPEC mp_limb_t mpn_sec_add_1 (mp_ptr, mp_srcptr, mp_size_t, mp_limb_t, mp_ptr); 1653 #define mpn_sec_add_1_itch __MPN(sec_add_1_itch) 1654 __GMP_DECLSPEC mp_size_t mpn_sec_add_1_itch (mp_size_t) __GMP_ATTRIBUTE_PURE; 1655 1656 #define mpn_sec_sub_1 __MPN(sec_sub_1) 1657 __GMP_DECLSPEC mp_limb_t mpn_sec_sub_1 (mp_ptr, mp_srcptr, mp_size_t, mp_limb_t, mp_ptr); 1658 #define mpn_sec_sub_1_itch __MPN(sec_sub_1_itch) 1659 __GMP_DECLSPEC mp_size_t mpn_sec_sub_1_itch (mp_size_t) __GMP_ATTRIBUTE_PURE; 1660 1661 #define mpn_cnd_swap __MPN(cnd_swap) 1662 __GMP_DECLSPEC void mpn_cnd_swap (mp_limb_t, volatile mp_limb_t *, volatile mp_limb_t *, mp_size_t); 1663 1664 #define mpn_sec_mul __MPN(sec_mul) 1665 __GMP_DECLSPEC void mpn_sec_mul (mp_ptr, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t, mp_ptr); 1666 #define mpn_sec_mul_itch __MPN(sec_mul_itch) 1667 __GMP_DECLSPEC mp_size_t mpn_sec_mul_itch (mp_size_t, mp_size_t) __GMP_ATTRIBUTE_PURE; 1668 1669 #define mpn_sec_sqr __MPN(sec_sqr) 1670 __GMP_DECLSPEC void mpn_sec_sqr (mp_ptr, mp_srcptr, mp_size_t, mp_ptr); 1671 #define mpn_sec_sqr_itch __MPN(sec_sqr_itch) 1672 __GMP_DECLSPEC mp_size_t mpn_sec_sqr_itch (mp_size_t) __GMP_ATTRIBUTE_PURE; 1673 1674 #define mpn_sec_powm __MPN(sec_powm) 1675 __GMP_DECLSPEC void mpn_sec_powm (mp_ptr, mp_srcptr, mp_size_t, mp_srcptr, mp_bitcnt_t, mp_srcptr, mp_size_t, mp_ptr); 1676 #define mpn_sec_powm_itch __MPN(sec_powm_itch) 1677 __GMP_DECLSPEC mp_size_t mpn_sec_powm_itch (mp_size_t, mp_bitcnt_t, mp_size_t) __GMP_ATTRIBUTE_PURE; 1678 1679 #define mpn_sec_tabselect __MPN(sec_tabselect) 1680 __GMP_DECLSPEC void mpn_sec_tabselect (volatile mp_limb_t *, volatile const mp_limb_t *, mp_size_t, mp_size_t, mp_size_t); 1681 1682 #define mpn_sec_div_qr __MPN(sec_div_qr) 1683 __GMP_DECLSPEC mp_limb_t mpn_sec_div_qr (mp_ptr, mp_ptr, mp_size_t, mp_srcptr, mp_size_t, mp_ptr); 1684 #define mpn_sec_div_qr_itch __MPN(sec_div_qr_itch) 1685 __GMP_DECLSPEC mp_size_t mpn_sec_div_qr_itch (mp_size_t, mp_size_t) __GMP_ATTRIBUTE_PURE; 1686 #define mpn_sec_div_r __MPN(sec_div_r) 1687 __GMP_DECLSPEC void mpn_sec_div_r (mp_ptr, mp_size_t, mp_srcptr, mp_size_t, mp_ptr); 1688 #define mpn_sec_div_r_itch __MPN(sec_div_r_itch) 1689 __GMP_DECLSPEC mp_size_t mpn_sec_div_r_itch (mp_size_t, mp_size_t) __GMP_ATTRIBUTE_PURE; 1690 1691 #define mpn_sec_invert __MPN(sec_invert) 1692 __GMP_DECLSPEC int mpn_sec_invert (mp_ptr, mp_ptr, mp_srcptr, mp_size_t, mp_bitcnt_t, mp_ptr); 1693 #define mpn_sec_invert_itch __MPN(sec_invert_itch) 1694 __GMP_DECLSPEC mp_size_t mpn_sec_invert_itch (mp_size_t) __GMP_ATTRIBUTE_PURE; 1695 1696 1697 /**************** mpz inlines ****************/ 1698 1699 /* The following are provided as inlines where possible, but always exist as 1700 library functions too, for binary compatibility. 1701 1702 Within gmp itself this inlining generally isn't relied on, since it 1703 doesn't get done for all compilers, whereas if something is worth 1704 inlining then it's worth arranging always. 1705 1706 There are two styles of inlining here. When the same bit of code is 1707 wanted for the inline as for the library version, then __GMP_FORCE_foo 1708 arranges for that code to be emitted and the __GMP_EXTERN_INLINE 1709 directive suppressed, eg. mpz_fits_uint_p. When a different bit of code 1710 is wanted for the inline than for the library version, then 1711 __GMP_FORCE_foo arranges the inline to be suppressed, eg. mpz_abs. */ 1712 1713 #if defined (__GMP_EXTERN_INLINE) && ! defined (__GMP_FORCE_mpz_abs) 1714 __GMP_EXTERN_INLINE void 1715 mpz_abs (mpz_ptr __gmp_w, mpz_srcptr __gmp_u) 1716 { 1717 if (__gmp_w != __gmp_u) 1718 mpz_set (__gmp_w, __gmp_u); 1719 __gmp_w->_mp_size = __GMP_ABS (__gmp_w->_mp_size); 1720 } 1721 #endif 1722 1723 #if GMP_NAIL_BITS == 0 1724 #define __GMPZ_FITS_UTYPE_P(z,maxval) \ 1725 mp_size_t __gmp_n = z->_mp_size; \ 1726 mp_ptr __gmp_p = z->_mp_d; \ 1727 return (__gmp_n == 0 || (__gmp_n == 1 && __gmp_p[0] <= maxval)); 1728 #else 1729 #define __GMPZ_FITS_UTYPE_P(z,maxval) \ 1730 mp_size_t __gmp_n = z->_mp_size; \ 1731 mp_ptr __gmp_p = z->_mp_d; \ 1732 return (__gmp_n == 0 || (__gmp_n == 1 && __gmp_p[0] <= maxval) \ 1733 || (__gmp_n == 2 && __gmp_p[1] <= ((mp_limb_t) maxval >> GMP_NUMB_BITS))); 1734 #endif 1735 1736 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_fits_uint_p) 1737 #if ! defined (__GMP_FORCE_mpz_fits_uint_p) 1738 __GMP_EXTERN_INLINE 1739 #endif 1740 int 1741 mpz_fits_uint_p (mpz_srcptr __gmp_z) __GMP_NOTHROW 1742 { 1743 __GMPZ_FITS_UTYPE_P (__gmp_z, UINT_MAX); 1744 } 1745 #endif 1746 1747 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_fits_ulong_p) 1748 #if ! defined (__GMP_FORCE_mpz_fits_ulong_p) 1749 __GMP_EXTERN_INLINE 1750 #endif 1751 int 1752 mpz_fits_ulong_p (mpz_srcptr __gmp_z) __GMP_NOTHROW 1753 { 1754 __GMPZ_FITS_UTYPE_P (__gmp_z, ULONG_MAX); 1755 } 1756 #endif 1757 1758 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_fits_ushort_p) 1759 #if ! defined (__GMP_FORCE_mpz_fits_ushort_p) 1760 __GMP_EXTERN_INLINE 1761 #endif 1762 int 1763 mpz_fits_ushort_p (mpz_srcptr __gmp_z) __GMP_NOTHROW 1764 { 1765 __GMPZ_FITS_UTYPE_P (__gmp_z, USHRT_MAX); 1766 } 1767 #endif 1768 1769 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_get_ui) 1770 #if ! defined (__GMP_FORCE_mpz_get_ui) 1771 __GMP_EXTERN_INLINE 1772 #endif 1773 unsigned long 1774 mpz_get_ui (mpz_srcptr __gmp_z) __GMP_NOTHROW 1775 { 1776 mp_ptr __gmp_p = __gmp_z->_mp_d; 1777 mp_size_t __gmp_n = __gmp_z->_mp_size; 1778 mp_limb_t __gmp_l = __gmp_p[0]; 1779 /* This is a "#if" rather than a plain "if" so as to avoid gcc warnings 1780 about "<< GMP_NUMB_BITS" exceeding the type size, and to avoid Borland 1781 C++ 6.0 warnings about condition always true for something like 1782 "ULONG_MAX < GMP_NUMB_MASK". */ 1783 #if GMP_NAIL_BITS == 0 || defined (_LONG_LONG_LIMB) 1784 /* limb==long and no nails, or limb==longlong, one limb is enough */ 1785 return (__gmp_n != 0 ? __gmp_l : 0); 1786 #else 1787 /* limb==long and nails, need two limbs when available */ 1788 __gmp_n = __GMP_ABS (__gmp_n); 1789 if (__gmp_n <= 1) 1790 return (__gmp_n != 0 ? __gmp_l : 0); 1791 else 1792 return __gmp_l + (__gmp_p[1] << GMP_NUMB_BITS); 1793 #endif 1794 } 1795 #endif 1796 1797 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_getlimbn) 1798 #if ! defined (__GMP_FORCE_mpz_getlimbn) 1799 __GMP_EXTERN_INLINE 1800 #endif 1801 mp_limb_t 1802 mpz_getlimbn (mpz_srcptr __gmp_z, mp_size_t __gmp_n) __GMP_NOTHROW 1803 { 1804 mp_limb_t __gmp_result = 0; 1805 if (__GMP_LIKELY (__gmp_n >= 0 && __gmp_n < __GMP_ABS (__gmp_z->_mp_size))) 1806 __gmp_result = __gmp_z->_mp_d[__gmp_n]; 1807 return __gmp_result; 1808 } 1809 #endif 1810 1811 #if defined (__GMP_EXTERN_INLINE) && ! defined (__GMP_FORCE_mpz_neg) 1812 __GMP_EXTERN_INLINE void 1813 mpz_neg (mpz_ptr __gmp_w, mpz_srcptr __gmp_u) 1814 { 1815 if (__gmp_w != __gmp_u) 1816 mpz_set (__gmp_w, __gmp_u); 1817 __gmp_w->_mp_size = - __gmp_w->_mp_size; 1818 } 1819 #endif 1820 1821 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_perfect_square_p) 1822 #if ! defined (__GMP_FORCE_mpz_perfect_square_p) 1823 __GMP_EXTERN_INLINE 1824 #endif 1825 int 1826 mpz_perfect_square_p (mpz_srcptr __gmp_a) 1827 { 1828 mp_size_t __gmp_asize; 1829 int __gmp_result; 1830 1831 __gmp_asize = __gmp_a->_mp_size; 1832 __gmp_result = (__gmp_asize >= 0); /* zero is a square, negatives are not */ 1833 if (__GMP_LIKELY (__gmp_asize > 0)) 1834 __gmp_result = mpn_perfect_square_p (__gmp_a->_mp_d, __gmp_asize); 1835 return __gmp_result; 1836 } 1837 #endif 1838 1839 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_popcount) 1840 #if ! defined (__GMP_FORCE_mpz_popcount) 1841 __GMP_EXTERN_INLINE 1842 #endif 1843 mp_bitcnt_t 1844 mpz_popcount (mpz_srcptr __gmp_u) __GMP_NOTHROW 1845 { 1846 mp_size_t __gmp_usize; 1847 mp_bitcnt_t __gmp_result; 1848 1849 __gmp_usize = __gmp_u->_mp_size; 1850 __gmp_result = (__gmp_usize < 0 ? ~ __GMP_CAST (mp_bitcnt_t, 0) : __GMP_CAST (mp_bitcnt_t, 0)); 1851 if (__GMP_LIKELY (__gmp_usize > 0)) 1852 __gmp_result = mpn_popcount (__gmp_u->_mp_d, __gmp_usize); 1853 return __gmp_result; 1854 } 1855 #endif 1856 1857 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_set_q) 1858 #if ! defined (__GMP_FORCE_mpz_set_q) 1859 __GMP_EXTERN_INLINE 1860 #endif 1861 void 1862 mpz_set_q (mpz_ptr __gmp_w, mpq_srcptr __gmp_u) 1863 { 1864 mpz_tdiv_q (__gmp_w, mpq_numref (__gmp_u), mpq_denref (__gmp_u)); 1865 } 1866 #endif 1867 1868 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_size) 1869 #if ! defined (__GMP_FORCE_mpz_size) 1870 __GMP_EXTERN_INLINE 1871 #endif 1872 size_t 1873 mpz_size (mpz_srcptr __gmp_z) __GMP_NOTHROW 1874 { 1875 return __GMP_ABS (__gmp_z->_mp_size); 1876 } 1877 #endif 1878 1879 1880 /**************** mpq inlines ****************/ 1881 1882 #if defined (__GMP_EXTERN_INLINE) && ! defined (__GMP_FORCE_mpq_abs) 1883 __GMP_EXTERN_INLINE void 1884 mpq_abs (mpq_ptr __gmp_w, mpq_srcptr __gmp_u) 1885 { 1886 if (__gmp_w != __gmp_u) 1887 mpq_set (__gmp_w, __gmp_u); 1888 __gmp_w->_mp_num._mp_size = __GMP_ABS (__gmp_w->_mp_num._mp_size); 1889 } 1890 #endif 1891 1892 #if defined (__GMP_EXTERN_INLINE) && ! defined (__GMP_FORCE_mpq_neg) 1893 __GMP_EXTERN_INLINE void 1894 mpq_neg (mpq_ptr __gmp_w, mpq_srcptr __gmp_u) 1895 { 1896 if (__gmp_w != __gmp_u) 1897 mpq_set (__gmp_w, __gmp_u); 1898 __gmp_w->_mp_num._mp_size = - __gmp_w->_mp_num._mp_size; 1899 } 1900 #endif 1901 1902 1903 /**************** mpn inlines ****************/ 1904 1905 /* The comments with __GMPN_ADD_1 below apply here too. 1906 1907 The test for FUNCTION returning 0 should predict well. If it's assumed 1908 {yp,ysize} will usually have a random number of bits then the high limb 1909 won't be full and a carry out will occur a good deal less than 50% of the 1910 time. 1911 1912 ysize==0 isn't a documented feature, but is used internally in a few 1913 places. 1914 1915 Producing cout last stops it using up a register during the main part of 1916 the calculation, though gcc (as of 3.0) on an "if (mpn_add (...))" 1917 doesn't seem able to move the true and false legs of the conditional up 1918 to the two places cout is generated. */ 1919 1920 #define __GMPN_AORS(cout, wp, xp, xsize, yp, ysize, FUNCTION, TEST) \ 1921 do { \ 1922 mp_size_t __gmp_i; \ 1923 mp_limb_t __gmp_x; \ 1924 \ 1925 /* ASSERT ((ysize) >= 0); */ \ 1926 /* ASSERT ((xsize) >= (ysize)); */ \ 1927 /* ASSERT (MPN_SAME_OR_SEPARATE2_P (wp, xsize, xp, xsize)); */ \ 1928 /* ASSERT (MPN_SAME_OR_SEPARATE2_P (wp, xsize, yp, ysize)); */ \ 1929 \ 1930 __gmp_i = (ysize); \ 1931 if (__gmp_i != 0) \ 1932 { \ 1933 if (FUNCTION (wp, xp, yp, __gmp_i)) \ 1934 { \ 1935 do \ 1936 { \ 1937 if (__gmp_i >= (xsize)) \ 1938 { \ 1939 (cout) = 1; \ 1940 goto __gmp_done; \ 1941 } \ 1942 __gmp_x = (xp)[__gmp_i]; \ 1943 } \ 1944 while (TEST); \ 1945 } \ 1946 } \ 1947 if ((wp) != (xp)) \ 1948 __GMPN_COPY_REST (wp, xp, xsize, __gmp_i); \ 1949 (cout) = 0; \ 1950 __gmp_done: \ 1951 ; \ 1952 } while (0) 1953 1954 #define __GMPN_ADD(cout, wp, xp, xsize, yp, ysize) \ 1955 __GMPN_AORS (cout, wp, xp, xsize, yp, ysize, mpn_add_n, \ 1956 (((wp)[__gmp_i++] = (__gmp_x + 1) & GMP_NUMB_MASK) == 0)) 1957 #define __GMPN_SUB(cout, wp, xp, xsize, yp, ysize) \ 1958 __GMPN_AORS (cout, wp, xp, xsize, yp, ysize, mpn_sub_n, \ 1959 (((wp)[__gmp_i++] = (__gmp_x - 1) & GMP_NUMB_MASK), __gmp_x == 0)) 1960 1961 1962 /* The use of __gmp_i indexing is designed to ensure a compile time src==dst 1963 remains nice and clear to the compiler, so that __GMPN_COPY_REST can 1964 disappear, and the load/add/store gets a chance to become a 1965 read-modify-write on CISC CPUs. 1966 1967 Alternatives: 1968 1969 Using a pair of pointers instead of indexing would be possible, but gcc 1970 isn't able to recognise compile-time src==dst in that case, even when the 1971 pointers are incremented more or less together. Other compilers would 1972 very likely have similar difficulty. 1973 1974 gcc could use "if (__builtin_constant_p(src==dst) && src==dst)" or 1975 similar to detect a compile-time src==dst. This works nicely on gcc 1976 2.95.x, it's not good on gcc 3.0 where __builtin_constant_p(p==p) seems 1977 to be always false, for a pointer p. But the current code form seems 1978 good enough for src==dst anyway. 1979 1980 gcc on x86 as usual doesn't give particularly good flags handling for the 1981 carry/borrow detection. It's tempting to want some multi instruction asm 1982 blocks to help it, and this was tried, but in truth there's only a few 1983 instructions to save and any gain is all too easily lost by register 1984 juggling setting up for the asm. */ 1985 1986 #if GMP_NAIL_BITS == 0 1987 #define __GMPN_AORS_1(cout, dst, src, n, v, OP, CB) \ 1988 do { \ 1989 mp_size_t __gmp_i; \ 1990 mp_limb_t __gmp_x, __gmp_r; \ 1991 \ 1992 /* ASSERT ((n) >= 1); */ \ 1993 /* ASSERT (MPN_SAME_OR_SEPARATE_P (dst, src, n)); */ \ 1994 \ 1995 __gmp_x = (src)[0]; \ 1996 __gmp_r = __gmp_x OP (v); \ 1997 (dst)[0] = __gmp_r; \ 1998 if (CB (__gmp_r, __gmp_x, (v))) \ 1999 { \ 2000 (cout) = 1; \ 2001 for (__gmp_i = 1; __gmp_i < (n);) \ 2002 { \ 2003 __gmp_x = (src)[__gmp_i]; \ 2004 __gmp_r = __gmp_x OP 1; \ 2005 (dst)[__gmp_i] = __gmp_r; \ 2006 ++__gmp_i; \ 2007 if (!CB (__gmp_r, __gmp_x, 1)) \ 2008 { \ 2009 if ((src) != (dst)) \ 2010 __GMPN_COPY_REST (dst, src, n, __gmp_i); \ 2011 (cout) = 0; \ 2012 break; \ 2013 } \ 2014 } \ 2015 } \ 2016 else \ 2017 { \ 2018 if ((src) != (dst)) \ 2019 __GMPN_COPY_REST (dst, src, n, 1); \ 2020 (cout) = 0; \ 2021 } \ 2022 } while (0) 2023 #endif 2024 2025 #if GMP_NAIL_BITS >= 1 2026 #define __GMPN_AORS_1(cout, dst, src, n, v, OP, CB) \ 2027 do { \ 2028 mp_size_t __gmp_i; \ 2029 mp_limb_t __gmp_x, __gmp_r; \ 2030 \ 2031 /* ASSERT ((n) >= 1); */ \ 2032 /* ASSERT (MPN_SAME_OR_SEPARATE_P (dst, src, n)); */ \ 2033 \ 2034 __gmp_x = (src)[0]; \ 2035 __gmp_r = __gmp_x OP (v); \ 2036 (dst)[0] = __gmp_r & GMP_NUMB_MASK; \ 2037 if (__gmp_r >> GMP_NUMB_BITS != 0) \ 2038 { \ 2039 (cout) = 1; \ 2040 for (__gmp_i = 1; __gmp_i < (n);) \ 2041 { \ 2042 __gmp_x = (src)[__gmp_i]; \ 2043 __gmp_r = __gmp_x OP 1; \ 2044 (dst)[__gmp_i] = __gmp_r & GMP_NUMB_MASK; \ 2045 ++__gmp_i; \ 2046 if (__gmp_r >> GMP_NUMB_BITS == 0) \ 2047 { \ 2048 if ((src) != (dst)) \ 2049 __GMPN_COPY_REST (dst, src, n, __gmp_i); \ 2050 (cout) = 0; \ 2051 break; \ 2052 } \ 2053 } \ 2054 } \ 2055 else \ 2056 { \ 2057 if ((src) != (dst)) \ 2058 __GMPN_COPY_REST (dst, src, n, 1); \ 2059 (cout) = 0; \ 2060 } \ 2061 } while (0) 2062 #endif 2063 2064 #define __GMPN_ADDCB(r,x,y) ((r) < (y)) 2065 #define __GMPN_SUBCB(r,x,y) ((x) < (y)) 2066 2067 #define __GMPN_ADD_1(cout, dst, src, n, v) \ 2068 __GMPN_AORS_1(cout, dst, src, n, v, +, __GMPN_ADDCB) 2069 #define __GMPN_SUB_1(cout, dst, src, n, v) \ 2070 __GMPN_AORS_1(cout, dst, src, n, v, -, __GMPN_SUBCB) 2071 2072 2073 /* Compare {xp,size} and {yp,size}, setting "result" to positive, zero or 2074 negative. size==0 is allowed. On random data usually only one limb will 2075 need to be examined to get a result, so it's worth having it inline. */ 2076 #define __GMPN_CMP(result, xp, yp, size) \ 2077 do { \ 2078 mp_size_t __gmp_i; \ 2079 mp_limb_t __gmp_x, __gmp_y; \ 2080 \ 2081 /* ASSERT ((size) >= 0); */ \ 2082 \ 2083 (result) = 0; \ 2084 __gmp_i = (size); \ 2085 while (--__gmp_i >= 0) \ 2086 { \ 2087 __gmp_x = (xp)[__gmp_i]; \ 2088 __gmp_y = (yp)[__gmp_i]; \ 2089 if (__gmp_x != __gmp_y) \ 2090 { \ 2091 /* Cannot use __gmp_x - __gmp_y, may overflow an "int" */ \ 2092 (result) = (__gmp_x > __gmp_y ? 1 : -1); \ 2093 break; \ 2094 } \ 2095 } \ 2096 } while (0) 2097 2098 2099 #if defined (__GMPN_COPY) && ! defined (__GMPN_COPY_REST) 2100 #define __GMPN_COPY_REST(dst, src, size, start) \ 2101 do { \ 2102 /* ASSERT ((start) >= 0); */ \ 2103 /* ASSERT ((start) <= (size)); */ \ 2104 __GMPN_COPY ((dst)+(start), (src)+(start), (size)-(start)); \ 2105 } while (0) 2106 #endif 2107 2108 /* Copy {src,size} to {dst,size}, starting at "start". This is designed to 2109 keep the indexing dst[j] and src[j] nice and simple for __GMPN_ADD_1, 2110 __GMPN_ADD, etc. */ 2111 #if ! defined (__GMPN_COPY_REST) 2112 #define __GMPN_COPY_REST(dst, src, size, start) \ 2113 do { \ 2114 mp_size_t __gmp_j; \ 2115 /* ASSERT ((size) >= 0); */ \ 2116 /* ASSERT ((start) >= 0); */ \ 2117 /* ASSERT ((start) <= (size)); */ \ 2118 /* ASSERT (MPN_SAME_OR_SEPARATE_P (dst, src, size)); */ \ 2119 __GMP_CRAY_Pragma ("_CRI ivdep"); \ 2120 for (__gmp_j = (start); __gmp_j < (size); __gmp_j++) \ 2121 (dst)[__gmp_j] = (src)[__gmp_j]; \ 2122 } while (0) 2123 #endif 2124 2125 /* Enhancement: Use some of the smarter code from gmp-impl.h. Maybe use 2126 mpn_copyi if there's a native version, and if we don't mind demanding 2127 binary compatibility for it (on targets which use it). */ 2128 2129 #if ! defined (__GMPN_COPY) 2130 #define __GMPN_COPY(dst, src, size) __GMPN_COPY_REST (dst, src, size, 0) 2131 #endif 2132 2133 2134 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_add) 2135 #if ! defined (__GMP_FORCE_mpn_add) 2136 __GMP_EXTERN_INLINE 2137 #endif 2138 mp_limb_t 2139 mpn_add (mp_ptr __gmp_wp, mp_srcptr __gmp_xp, mp_size_t __gmp_xsize, mp_srcptr __gmp_yp, mp_size_t __gmp_ysize) 2140 { 2141 mp_limb_t __gmp_c; 2142 __GMPN_ADD (__gmp_c, __gmp_wp, __gmp_xp, __gmp_xsize, __gmp_yp, __gmp_ysize); 2143 return __gmp_c; 2144 } 2145 #endif 2146 2147 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_add_1) 2148 #if ! defined (__GMP_FORCE_mpn_add_1) 2149 __GMP_EXTERN_INLINE 2150 #endif 2151 mp_limb_t 2152 mpn_add_1 (mp_ptr __gmp_dst, mp_srcptr __gmp_src, mp_size_t __gmp_size, mp_limb_t __gmp_n) __GMP_NOTHROW 2153 { 2154 mp_limb_t __gmp_c; 2155 __GMPN_ADD_1 (__gmp_c, __gmp_dst, __gmp_src, __gmp_size, __gmp_n); 2156 return __gmp_c; 2157 } 2158 #endif 2159 2160 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_cmp) 2161 #if ! defined (__GMP_FORCE_mpn_cmp) 2162 __GMP_EXTERN_INLINE 2163 #endif 2164 int 2165 mpn_cmp (mp_srcptr __gmp_xp, mp_srcptr __gmp_yp, mp_size_t __gmp_size) __GMP_NOTHROW 2166 { 2167 int __gmp_result; 2168 __GMPN_CMP (__gmp_result, __gmp_xp, __gmp_yp, __gmp_size); 2169 return __gmp_result; 2170 } 2171 #endif 2172 2173 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_zero_p) 2174 #if ! defined (__GMP_FORCE_mpn_zero_p) 2175 __GMP_EXTERN_INLINE 2176 #endif 2177 int 2178 mpn_zero_p (mp_srcptr __gmp_p, mp_size_t __gmp_n) __GMP_NOTHROW 2179 { 2180 /* if (__GMP_LIKELY (__gmp_n > 0)) */ 2181 do { 2182 if (__gmp_p[--__gmp_n] != 0) 2183 return 0; 2184 } while (__gmp_n != 0); 2185 return 1; 2186 } 2187 #endif 2188 2189 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_sub) 2190 #if ! defined (__GMP_FORCE_mpn_sub) 2191 __GMP_EXTERN_INLINE 2192 #endif 2193 mp_limb_t 2194 mpn_sub (mp_ptr __gmp_wp, mp_srcptr __gmp_xp, mp_size_t __gmp_xsize, mp_srcptr __gmp_yp, mp_size_t __gmp_ysize) 2195 { 2196 mp_limb_t __gmp_c; 2197 __GMPN_SUB (__gmp_c, __gmp_wp, __gmp_xp, __gmp_xsize, __gmp_yp, __gmp_ysize); 2198 return __gmp_c; 2199 } 2200 #endif 2201 2202 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_sub_1) 2203 #if ! defined (__GMP_FORCE_mpn_sub_1) 2204 __GMP_EXTERN_INLINE 2205 #endif 2206 mp_limb_t 2207 mpn_sub_1 (mp_ptr __gmp_dst, mp_srcptr __gmp_src, mp_size_t __gmp_size, mp_limb_t __gmp_n) __GMP_NOTHROW 2208 { 2209 mp_limb_t __gmp_c; 2210 __GMPN_SUB_1 (__gmp_c, __gmp_dst, __gmp_src, __gmp_size, __gmp_n); 2211 return __gmp_c; 2212 } 2213 #endif 2214 2215 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_neg) 2216 #if ! defined (__GMP_FORCE_mpn_neg) 2217 __GMP_EXTERN_INLINE 2218 #endif 2219 mp_limb_t 2220 mpn_neg (mp_ptr __gmp_rp, mp_srcptr __gmp_up, mp_size_t __gmp_n) 2221 { 2222 while (*__gmp_up == 0) /* Low zero limbs are unchanged by negation. */ 2223 { 2224 *__gmp_rp = 0; 2225 if (!--__gmp_n) /* All zero */ 2226 return 0; 2227 ++__gmp_up; ++__gmp_rp; 2228 } 2229 2230 *__gmp_rp = (- *__gmp_up) & GMP_NUMB_MASK; 2231 2232 if (--__gmp_n) /* Higher limbs get complemented. */ 2233 mpn_com (++__gmp_rp, ++__gmp_up, __gmp_n); 2234 2235 return 1; 2236 } 2237 #endif 2238 2239 #if defined (__cplusplus) 2240 } 2241 #endif 2242 2243 2244 /* Allow faster testing for negative, zero, and positive. */ 2245 #define mpz_sgn(Z) ((Z)->_mp_size < 0 ? -1 : (Z)->_mp_size > 0) 2246 #define mpf_sgn(F) ((F)->_mp_size < 0 ? -1 : (F)->_mp_size > 0) 2247 #define mpq_sgn(Q) ((Q)->_mp_num._mp_size < 0 ? -1 : (Q)->_mp_num._mp_size > 0) 2248 2249 /* When using GCC, optimize certain common comparisons. */ 2250 #if defined (__GNUC__) && __GNUC__ >= 2 2251 #define mpz_cmp_ui(Z,UI) \ 2252 (__builtin_constant_p (UI) && (UI) == 0 \ 2253 ? mpz_sgn (Z) : _mpz_cmp_ui (Z,UI)) 2254 #define mpz_cmp_si(Z,SI) \ 2255 (__builtin_constant_p ((SI) >= 0) && (SI) >= 0 \ 2256 ? mpz_cmp_ui (Z, __GMP_CAST (unsigned long, SI)) \ 2257 : _mpz_cmp_si (Z,SI)) 2258 #define mpq_cmp_ui(Q,NUI,DUI) \ 2259 (__builtin_constant_p (NUI) && (NUI) == 0 ? mpq_sgn (Q) \ 2260 : __builtin_constant_p ((NUI) == (DUI)) && (NUI) == (DUI) \ 2261 ? mpz_cmp (mpq_numref (Q), mpq_denref (Q)) \ 2262 : _mpq_cmp_ui (Q,NUI,DUI)) 2263 #define mpq_cmp_si(q,n,d) \ 2264 (__builtin_constant_p ((n) >= 0) && (n) >= 0 \ 2265 ? mpq_cmp_ui (q, __GMP_CAST (unsigned long, n), d) \ 2266 : _mpq_cmp_si (q, n, d)) 2267 #else 2268 #define mpz_cmp_ui(Z,UI) _mpz_cmp_ui (Z,UI) 2269 #define mpz_cmp_si(Z,UI) _mpz_cmp_si (Z,UI) 2270 #define mpq_cmp_ui(Q,NUI,DUI) _mpq_cmp_ui (Q,NUI,DUI) 2271 #define mpq_cmp_si(q,n,d) _mpq_cmp_si(q,n,d) 2272 #endif 2273 2274 2275 /* Using "&" rather than "&&" means these can come out branch-free. Every 2276 mpz_t has at least one limb allocated, so fetching the low limb is always 2277 allowed. */ 2278 #define mpz_odd_p(z) (((z)->_mp_size != 0) & __GMP_CAST (int, (z)->_mp_d[0])) 2279 #define mpz_even_p(z) (! mpz_odd_p (z)) 2280 2281 2282 /**************** C++ routines ****************/ 2283 2284 #ifdef __cplusplus 2285 __GMP_DECLSPEC_XX std::ostream& operator<< (std::ostream &, mpz_srcptr); 2286 __GMP_DECLSPEC_XX std::ostream& operator<< (std::ostream &, mpq_srcptr); 2287 __GMP_DECLSPEC_XX std::ostream& operator<< (std::ostream &, mpf_srcptr); 2288 __GMP_DECLSPEC_XX std::istream& operator>> (std::istream &, mpz_ptr); 2289 __GMP_DECLSPEC_XX std::istream& operator>> (std::istream &, mpq_ptr); 2290 __GMP_DECLSPEC_XX std::istream& operator>> (std::istream &, mpf_ptr); 2291 #endif 2292 2293 2294 /* Source-level compatibility with GMP 2 and earlier. */ 2295 #define mpn_divmod(qp,np,nsize,dp,dsize) \ 2296 mpn_divrem (qp, __GMP_CAST (mp_size_t, 0), np, nsize, dp, dsize) 2297 2298 /* Source-level compatibility with GMP 1. */ 2299 #define mpz_mdiv mpz_fdiv_q 2300 #define mpz_mdivmod mpz_fdiv_qr 2301 #define mpz_mmod mpz_fdiv_r 2302 #define mpz_mdiv_ui mpz_fdiv_q_ui 2303 #define mpz_mdivmod_ui(q,r,n,d) \ 2304 (((r) == 0) ? mpz_fdiv_q_ui (q,n,d) : mpz_fdiv_qr_ui (q,r,n,d)) 2305 #define mpz_mmod_ui(r,n,d) \ 2306 (((r) == 0) ? mpz_fdiv_ui (n,d) : mpz_fdiv_r_ui (r,n,d)) 2307 2308 /* Useful synonyms, but not quite compatible with GMP 1. */ 2309 #define mpz_div mpz_fdiv_q 2310 #define mpz_divmod mpz_fdiv_qr 2311 #define mpz_div_ui mpz_fdiv_q_ui 2312 #define mpz_divmod_ui mpz_fdiv_qr_ui 2313 #define mpz_div_2exp mpz_fdiv_q_2exp 2314 #define mpz_mod_2exp mpz_fdiv_r_2exp 2315 2316 enum 2317 { 2318 GMP_ERROR_NONE = 0, 2319 GMP_ERROR_UNSUPPORTED_ARGUMENT = 1, 2320 GMP_ERROR_DIVISION_BY_ZERO = 2, 2321 GMP_ERROR_SQRT_OF_NEGATIVE = 4, 2322 GMP_ERROR_INVALID_ARGUMENT = 8 2323 }; 2324 2325 /* Define CC and CFLAGS which were used to build this version of GMP */ 2326 #define __GMP_CC "gcc" 2327 #define __GMP_CFLAGS "-m32 -O2 -pedantic" 2328 2329 /* Major version number is the value of __GNU_MP__ too, above. */ 2330 #define __GNU_MP_VERSION 6 2331 #define __GNU_MP_VERSION_MINOR 2 2332 #define __GNU_MP_VERSION_PATCHLEVEL 0 2333 #define __GNU_MP_RELEASE (__GNU_MP_VERSION * 10000 + __GNU_MP_VERSION_MINOR * 100 + __GNU_MP_VERSION_PATCHLEVEL) 2334 2335 #define __GMP_H__ 2336 #endif /* __GMP_H__ */ 2337