1 /* Definitions for GNU multiple precision functions. -*- mode: c -*-
3 Copyright 1991, 1993-1997, 1999-2016, 2020 Free Software Foundation, Inc.
5 This file is part of the GNU MP Library.
7 The GNU MP Library is free software; you can redistribute it and/or modify
8 it under the terms of either:
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.
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
20 or both in parallel, as here.
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
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/. */
33 #if defined (__cplusplus)
34 #include <iosfwd> /* for std::istream, std::ostream, std::string */
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 0
43 #define GMP_LIMB_BITS 64
44 #define GMP_NAIL_BITS 0
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)
55 #include <stddef.h> /* for size_t */
58 /* Instantiated by configure. */
59 #if ! defined (__GMP_WITHIN_CONFIGURE)
60 /* #undef _LONG_LONG_LIMB */
61 #define __GMP_LIBGMP_DLL 0
65 /* __GMP_DECLSPEC supports Windows DLL versions of libgmp, and is empty in
66 all other circumstances.
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).
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.
77 Libtool DLL_EXPORT define is not used.
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.
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.
91 In gcc __declspec can go at either the start or end of a prototype.
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.
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. */
102 #if defined (__GNUC__)
103 #define __GMP_DECLSPEC_EXPORT __declspec(__dllexport__)
104 #define __GMP_DECLSPEC_IMPORT __declspec(__dllimport__)
106 #if defined (_MSC_VER) || defined (__BORLANDC__)
107 #define __GMP_DECLSPEC_EXPORT __declspec(dllexport)
108 #define __GMP_DECLSPEC_IMPORT __declspec(dllimport)
111 #define __GMP_DECLSPEC_EXPORT __export
112 #define __GMP_DECLSPEC_IMPORT __import
115 #define __GMP_DECLSPEC_EXPORT _Export
116 #define __GMP_DECLSPEC_IMPORT _Import
120 #ifdef __GMP_WITHIN_GMP
121 /* compiling to go into a DLL libgmp */
122 #define __GMP_DECLSPEC __GMP_DECLSPEC_EXPORT
124 /* compiling to go into an application which will link to a DLL libgmp */
125 #define __GMP_DECLSPEC __GMP_DECLSPEC_IMPORT
128 /* all other cases */
129 #define __GMP_DECLSPEC
133 #ifdef __GMP_SHORT_LIMB
134 typedef unsigned int mp_limb_t;
135 typedef int mp_limb_signed_t;
137 #ifdef _LONG_LONG_LIMB
138 typedef unsigned long long int mp_limb_t;
139 typedef long long int mp_limb_signed_t;
141 typedef unsigned long int mp_limb_t;
142 typedef long int mp_limb_signed_t;
145 typedef unsigned long int mp_bitcnt_t;
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. */
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. */
160 #endif /* __GNU_MP__ */
163 typedef __mpz_struct MP_INT; /* gmp 1 source compatibility */
164 typedef __mpz_struct mpz_t[1];
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;
174 #define __GMP_MP_SIZE_T_INT 0
175 typedef long int mp_size_t;
176 typedef long int mp_exp_t;
181 __mpz_struct _mp_num;
182 __mpz_struct _mp_den;
185 typedef __mpq_struct MP_RAT; /* gmp 1 source compatibility */
186 typedef __mpq_struct mpq_t[1];
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. */
201 /* typedef __mpf_struct MP_FLOAT; */
202 typedef __mpf_struct mpf_t[1];
204 /* Available random number generation algorithms. */
207 GMP_RAND_ALG_DEFAULT = 0,
208 GMP_RAND_ALG_LC = GMP_RAND_ALG_DEFAULT /* Linear congruential. */
211 /* Random state struct. */
214 mpz_t _mp_seed; /* _mp_d member points to state of the generator. */
215 gmp_randalg_t _mp_alg; /* Currently unused. */
217 void *_mp_lc; /* Pointer to function pointers structure. */
219 } __gmp_randstate_struct;
220 typedef __gmp_randstate_struct gmp_randstate_t[1];
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;
233 #ifdef __GMP_WITHIN_GMPXX
234 /* compiling to go into a DLL libgmpxx */
235 #define __GMP_DECLSPEC_XX __GMP_DECLSPEC_EXPORT
237 /* compiling to go into a application which will link to a DLL libgmpxx */
238 #define __GMP_DECLSPEC_XX __GMP_DECLSPEC_IMPORT
241 /* all other cases */
242 #define __GMP_DECLSPEC_XX
247 #define __MPN(x) __gmpn_##x
250 /* For reference, "defined(EOF)" cannot be used here. In g++ 2.95.4,
251 <iostream> defines EOF but not 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
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
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.
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. */
289 #define _GMP_H_HAVE_VA_LIST 1
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))
297 #define __GMP_GNUC_PREREQ(maj, min) 0
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__))
308 #define __GMP_ATTRIBUTE_PURE
312 /* __GMP_CAST allows us to use static_cast in C++, so our macros are clean
313 to "g++ -Wold-style-cast".
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
320 #define __GMP_CAST(type, expr) (static_cast<type> (expr))
322 #define __GMP_CAST(type, expr) ((type) (expr))
326 /* An empty "throw ()" means the function doesn't throw any C++ exceptions,
327 this can save some stack frame info in applications.
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
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.
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. */
343 #if defined (__cplusplus)
344 #if __cplusplus >= 201103L
345 #define __GMP_NOTHROW noexcept
347 #define __GMP_NOTHROW throw ()
350 #define __GMP_NOTHROW
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. */
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. */
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__))
368 #define __GMP_EXTERN_INLINE extern __inline__
370 #define __GMP_INLINE_PROTOTYPES 1
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
382 #define __GMP_EXTERN_INLINE static __inline
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
399 /* Microsoft's C compiler accepts __inline */
401 #define __GMP_EXTERN_INLINE __inline
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
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
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
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
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
442 #define __GMP_INLINE_PROTOTYPES 1
446 #define __GMP_ABS(x) ((x) >= 0 ? (x) : -(x))
447 #define __GMP_MAX(h,i) ((h) > (i) ? (h) : (i))
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)
455 #define __GMP_LIKELY(cond) (cond)
456 #define __GMP_UNLIKELY(cond) (cond)
460 #define __GMP_CRAY_Pragma(str) _Pragma (str)
462 #define __GMP_CRAY_Pragma(str)
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))
471 #if defined (__cplusplus)
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;
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;
486 #define mp_bits_per_limb __gmp_bits_per_limb
487 __GMP_DECLSPEC extern const int mp_bits_per_limb;
489 #define gmp_errno __gmp_errno
490 __GMP_DECLSPEC extern int gmp_errno;
492 #define gmp_version __gmp_version
493 __GMP_DECLSPEC extern const char * const gmp_version;
496 /**************** Random number routines. ****************/
499 #define gmp_randinit __gmp_randinit
500 __GMP_DECLSPEC void gmp_randinit (gmp_randstate_t, gmp_randalg_t, ...);
502 #define gmp_randinit_default __gmp_randinit_default
503 __GMP_DECLSPEC void gmp_randinit_default (gmp_randstate_t);
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);
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);
511 #define gmp_randinit_mt __gmp_randinit_mt
512 __GMP_DECLSPEC void gmp_randinit_mt (gmp_randstate_t);
514 #define gmp_randinit_set __gmp_randinit_set
515 __GMP_DECLSPEC void gmp_randinit_set (gmp_randstate_t, const __gmp_randstate_struct *);
517 #define gmp_randseed __gmp_randseed
518 __GMP_DECLSPEC void gmp_randseed (gmp_randstate_t, mpz_srcptr);
520 #define gmp_randseed_ui __gmp_randseed_ui
521 __GMP_DECLSPEC void gmp_randseed_ui (gmp_randstate_t, unsigned long int);
523 #define gmp_randclear __gmp_randclear
524 __GMP_DECLSPEC void gmp_randclear (gmp_randstate_t);
526 #define gmp_urandomb_ui __gmp_urandomb_ui
527 __GMP_DECLSPEC unsigned long gmp_urandomb_ui (gmp_randstate_t, unsigned long);
529 #define gmp_urandomm_ui __gmp_urandomm_ui
530 __GMP_DECLSPEC unsigned long gmp_urandomm_ui (gmp_randstate_t, unsigned long);
533 /**************** Formatted output routines. ****************/
535 #define gmp_asprintf __gmp_asprintf
536 __GMP_DECLSPEC int gmp_asprintf (char **, const char *, ...);
538 #define gmp_fprintf __gmp_fprintf
539 #ifdef _GMP_H_HAVE_FILE
540 __GMP_DECLSPEC int gmp_fprintf (FILE *, const char *, ...);
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 *, ...);
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);
553 #define gmp_printf __gmp_printf
554 __GMP_DECLSPEC int gmp_printf (const char *, ...);
556 #define gmp_snprintf __gmp_snprintf
557 __GMP_DECLSPEC int gmp_snprintf (char *, size_t, const char *, ...);
559 #define gmp_sprintf __gmp_sprintf
560 __GMP_DECLSPEC int gmp_sprintf (char *, const char *, ...);
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);
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);
572 #define gmp_vprintf __gmp_vprintf
573 #if defined (_GMP_H_HAVE_VA_LIST)
574 __GMP_DECLSPEC int gmp_vprintf (const char *, va_list);
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);
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);
588 /**************** Formatted input routines. ****************/
590 #define gmp_fscanf __gmp_fscanf
591 #ifdef _GMP_H_HAVE_FILE
592 __GMP_DECLSPEC int gmp_fscanf (FILE *, const char *, ...);
595 #define gmp_scanf __gmp_scanf
596 __GMP_DECLSPEC int gmp_scanf (const char *, ...);
598 #define gmp_sscanf __gmp_sscanf
599 __GMP_DECLSPEC int gmp_sscanf (const char *, const char *, ...);
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);
606 #define gmp_vscanf __gmp_vscanf
607 #if defined (_GMP_H_HAVE_VA_LIST)
608 __GMP_DECLSPEC int gmp_vscanf (const char *, va_list);
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);
617 /**************** Integer (i.e. Z) routines. ****************/
619 #define _mpz_realloc __gmpz_realloc
620 #define mpz_realloc __gmpz_realloc
621 __GMP_DECLSPEC void *_mpz_realloc (mpz_ptr, mp_size_t);
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);
628 #define mpz_add __gmpz_add
629 __GMP_DECLSPEC void mpz_add (mpz_ptr, mpz_srcptr, mpz_srcptr);
631 #define mpz_add_ui __gmpz_add_ui
632 __GMP_DECLSPEC void mpz_add_ui (mpz_ptr, mpz_srcptr, unsigned long int);
634 #define mpz_addmul __gmpz_addmul
635 __GMP_DECLSPEC void mpz_addmul (mpz_ptr, mpz_srcptr, mpz_srcptr);
637 #define mpz_addmul_ui __gmpz_addmul_ui
638 __GMP_DECLSPEC void mpz_addmul_ui (mpz_ptr, mpz_srcptr, unsigned long int);
640 #define mpz_and __gmpz_and
641 __GMP_DECLSPEC void mpz_and (mpz_ptr, mpz_srcptr, mpz_srcptr);
643 #define mpz_array_init __gmpz_array_init
644 __GMP_DECLSPEC void mpz_array_init (mpz_ptr, mp_size_t, mp_size_t);
646 #define mpz_bin_ui __gmpz_bin_ui
647 __GMP_DECLSPEC void mpz_bin_ui (mpz_ptr, mpz_srcptr, unsigned long int);
649 #define mpz_bin_uiui __gmpz_bin_uiui
650 __GMP_DECLSPEC void mpz_bin_uiui (mpz_ptr, unsigned long int, unsigned long int);
652 #define mpz_cdiv_q __gmpz_cdiv_q
653 __GMP_DECLSPEC void mpz_cdiv_q (mpz_ptr, mpz_srcptr, mpz_srcptr);
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);
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);
661 #define mpz_cdiv_qr __gmpz_cdiv_qr
662 __GMP_DECLSPEC void mpz_cdiv_qr (mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr);
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);
667 #define mpz_cdiv_r __gmpz_cdiv_r
668 __GMP_DECLSPEC void mpz_cdiv_r (mpz_ptr, mpz_srcptr, mpz_srcptr);
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);
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);
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;
679 #define mpz_clear __gmpz_clear
680 __GMP_DECLSPEC void mpz_clear (mpz_ptr);
682 #define mpz_clears __gmpz_clears
683 __GMP_DECLSPEC void mpz_clears (mpz_ptr, ...);
685 #define mpz_clrbit __gmpz_clrbit
686 __GMP_DECLSPEC void mpz_clrbit (mpz_ptr, mp_bitcnt_t);
688 #define mpz_cmp __gmpz_cmp
689 __GMP_DECLSPEC int mpz_cmp (mpz_srcptr, mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
691 #define mpz_cmp_d __gmpz_cmp_d
692 __GMP_DECLSPEC int mpz_cmp_d (mpz_srcptr, double) __GMP_ATTRIBUTE_PURE;
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;
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;
700 #define mpz_cmpabs __gmpz_cmpabs
701 __GMP_DECLSPEC int mpz_cmpabs (mpz_srcptr, mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
703 #define mpz_cmpabs_d __gmpz_cmpabs_d
704 __GMP_DECLSPEC int mpz_cmpabs_d (mpz_srcptr, double) __GMP_ATTRIBUTE_PURE;
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;
709 #define mpz_com __gmpz_com
710 __GMP_DECLSPEC void mpz_com (mpz_ptr, mpz_srcptr);
712 #define mpz_combit __gmpz_combit
713 __GMP_DECLSPEC void mpz_combit (mpz_ptr, mp_bitcnt_t);
715 #define mpz_congruent_p __gmpz_congruent_p
716 __GMP_DECLSPEC int mpz_congruent_p (mpz_srcptr, mpz_srcptr, mpz_srcptr) __GMP_ATTRIBUTE_PURE;
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;
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;
724 #define mpz_divexact __gmpz_divexact
725 __GMP_DECLSPEC void mpz_divexact (mpz_ptr, mpz_srcptr, mpz_srcptr);
727 #define mpz_divexact_ui __gmpz_divexact_ui
728 __GMP_DECLSPEC void mpz_divexact_ui (mpz_ptr, mpz_srcptr, unsigned long);
730 #define mpz_divisible_p __gmpz_divisible_p
731 __GMP_DECLSPEC int mpz_divisible_p (mpz_srcptr, mpz_srcptr) __GMP_ATTRIBUTE_PURE;
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;
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;
739 #define mpz_dump __gmpz_dump
740 __GMP_DECLSPEC void mpz_dump (mpz_srcptr);
742 #define mpz_export __gmpz_export
743 __GMP_DECLSPEC void *mpz_export (void *, size_t *, int, size_t, int, size_t, mpz_srcptr);
745 #define mpz_fac_ui __gmpz_fac_ui
746 __GMP_DECLSPEC void mpz_fac_ui (mpz_ptr, unsigned long int);
748 #define mpz_2fac_ui __gmpz_2fac_ui
749 __GMP_DECLSPEC void mpz_2fac_ui (mpz_ptr, unsigned long int);
751 #define mpz_mfac_uiui __gmpz_mfac_uiui
752 __GMP_DECLSPEC void mpz_mfac_uiui (mpz_ptr, unsigned long int, unsigned long int);
754 #define mpz_primorial_ui __gmpz_primorial_ui
755 __GMP_DECLSPEC void mpz_primorial_ui (mpz_ptr, unsigned long int);
757 #define mpz_fdiv_q __gmpz_fdiv_q
758 __GMP_DECLSPEC void mpz_fdiv_q (mpz_ptr, mpz_srcptr, mpz_srcptr);
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);
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);
766 #define mpz_fdiv_qr __gmpz_fdiv_qr
767 __GMP_DECLSPEC void mpz_fdiv_qr (mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr);
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);
772 #define mpz_fdiv_r __gmpz_fdiv_r
773 __GMP_DECLSPEC void mpz_fdiv_r (mpz_ptr, mpz_srcptr, mpz_srcptr);
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);
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);
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;
784 #define mpz_fib_ui __gmpz_fib_ui
785 __GMP_DECLSPEC void mpz_fib_ui (mpz_ptr, unsigned long int);
787 #define mpz_fib2_ui __gmpz_fib2_ui
788 __GMP_DECLSPEC void mpz_fib2_ui (mpz_ptr, mpz_ptr, unsigned long int);
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;
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;
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;
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;
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;
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;
814 #define mpz_gcd __gmpz_gcd
815 __GMP_DECLSPEC void mpz_gcd (mpz_ptr, mpz_srcptr, mpz_srcptr);
817 #define mpz_gcd_ui __gmpz_gcd_ui
818 __GMP_DECLSPEC unsigned long int mpz_gcd_ui (mpz_ptr, mpz_srcptr, unsigned long int);
820 #define mpz_gcdext __gmpz_gcdext
821 __GMP_DECLSPEC void mpz_gcdext (mpz_ptr, mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr);
823 #define mpz_get_d __gmpz_get_d
824 __GMP_DECLSPEC double mpz_get_d (mpz_srcptr) __GMP_ATTRIBUTE_PURE;
826 #define mpz_get_d_2exp __gmpz_get_d_2exp
827 __GMP_DECLSPEC double mpz_get_d_2exp (signed long int *, mpz_srcptr);
829 #define mpz_get_si __gmpz_get_si
830 __GMP_DECLSPEC /* signed */ long int mpz_get_si (mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
832 #define mpz_get_str __gmpz_get_str
833 __GMP_DECLSPEC char *mpz_get_str (char *, int, mpz_srcptr);
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;
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;
845 #define mpz_hamdist __gmpz_hamdist
846 __GMP_DECLSPEC mp_bitcnt_t mpz_hamdist (mpz_srcptr, mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
848 #define mpz_import __gmpz_import
849 __GMP_DECLSPEC void mpz_import (mpz_ptr, size_t, int, size_t, int, size_t, const void *);
851 #define mpz_init __gmpz_init
852 __GMP_DECLSPEC void mpz_init (mpz_ptr) __GMP_NOTHROW;
854 #define mpz_init2 __gmpz_init2
855 __GMP_DECLSPEC void mpz_init2 (mpz_ptr, mp_bitcnt_t);
857 #define mpz_inits __gmpz_inits
858 __GMP_DECLSPEC void mpz_inits (mpz_ptr, ...) __GMP_NOTHROW;
860 #define mpz_init_set __gmpz_init_set
861 __GMP_DECLSPEC void mpz_init_set (mpz_ptr, mpz_srcptr);
863 #define mpz_init_set_d __gmpz_init_set_d
864 __GMP_DECLSPEC void mpz_init_set_d (mpz_ptr, double);
866 #define mpz_init_set_si __gmpz_init_set_si
867 __GMP_DECLSPEC void mpz_init_set_si (mpz_ptr, signed long int);
869 #define mpz_init_set_str __gmpz_init_set_str
870 __GMP_DECLSPEC int mpz_init_set_str (mpz_ptr, const char *, int);
872 #define mpz_init_set_ui __gmpz_init_set_ui
873 __GMP_DECLSPEC void mpz_init_set_ui (mpz_ptr, unsigned long int);
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 *);
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);
885 #define mpz_invert __gmpz_invert
886 __GMP_DECLSPEC int mpz_invert (mpz_ptr, mpz_srcptr, mpz_srcptr);
888 #define mpz_ior __gmpz_ior
889 __GMP_DECLSPEC void mpz_ior (mpz_ptr, mpz_srcptr, mpz_srcptr);
891 #define mpz_jacobi __gmpz_jacobi
892 __GMP_DECLSPEC int mpz_jacobi (mpz_srcptr, mpz_srcptr) __GMP_ATTRIBUTE_PURE;
894 #define mpz_kronecker mpz_jacobi /* alias */
896 #define mpz_kronecker_si __gmpz_kronecker_si
897 __GMP_DECLSPEC int mpz_kronecker_si (mpz_srcptr, long) __GMP_ATTRIBUTE_PURE;
899 #define mpz_kronecker_ui __gmpz_kronecker_ui
900 __GMP_DECLSPEC int mpz_kronecker_ui (mpz_srcptr, unsigned long) __GMP_ATTRIBUTE_PURE;
902 #define mpz_si_kronecker __gmpz_si_kronecker
903 __GMP_DECLSPEC int mpz_si_kronecker (long, mpz_srcptr) __GMP_ATTRIBUTE_PURE;
905 #define mpz_ui_kronecker __gmpz_ui_kronecker
906 __GMP_DECLSPEC int mpz_ui_kronecker (unsigned long, mpz_srcptr) __GMP_ATTRIBUTE_PURE;
908 #define mpz_lcm __gmpz_lcm
909 __GMP_DECLSPEC void mpz_lcm (mpz_ptr, mpz_srcptr, mpz_srcptr);
911 #define mpz_lcm_ui __gmpz_lcm_ui
912 __GMP_DECLSPEC void mpz_lcm_ui (mpz_ptr, mpz_srcptr, unsigned long);
914 #define mpz_legendre mpz_jacobi /* alias */
916 #define mpz_lucnum_ui __gmpz_lucnum_ui
917 __GMP_DECLSPEC void mpz_lucnum_ui (mpz_ptr, unsigned long int);
919 #define mpz_lucnum2_ui __gmpz_lucnum2_ui
920 __GMP_DECLSPEC void mpz_lucnum2_ui (mpz_ptr, mpz_ptr, unsigned long int);
922 #define mpz_millerrabin __gmpz_millerrabin
923 __GMP_DECLSPEC int mpz_millerrabin (mpz_srcptr, int) __GMP_ATTRIBUTE_PURE;
925 #define mpz_mod __gmpz_mod
926 __GMP_DECLSPEC void mpz_mod (mpz_ptr, mpz_srcptr, mpz_srcptr);
928 #define mpz_mod_ui mpz_fdiv_r_ui /* same as fdiv_r because divisor unsigned */
930 #define mpz_mul __gmpz_mul
931 __GMP_DECLSPEC void mpz_mul (mpz_ptr, mpz_srcptr, mpz_srcptr);
933 #define mpz_mul_2exp __gmpz_mul_2exp
934 __GMP_DECLSPEC void mpz_mul_2exp (mpz_ptr, mpz_srcptr, mp_bitcnt_t);
936 #define mpz_mul_si __gmpz_mul_si
937 __GMP_DECLSPEC void mpz_mul_si (mpz_ptr, mpz_srcptr, long int);
939 #define mpz_mul_ui __gmpz_mul_ui
940 __GMP_DECLSPEC void mpz_mul_ui (mpz_ptr, mpz_srcptr, unsigned long int);
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);
947 #define mpz_nextprime __gmpz_nextprime
948 __GMP_DECLSPEC void mpz_nextprime (mpz_ptr, mpz_srcptr);
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);
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);
960 #define mpz_perfect_power_p __gmpz_perfect_power_p
961 __GMP_DECLSPEC int mpz_perfect_power_p (mpz_srcptr) __GMP_ATTRIBUTE_PURE;
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;
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;
973 #define mpz_pow_ui __gmpz_pow_ui
974 __GMP_DECLSPEC void mpz_pow_ui (mpz_ptr, mpz_srcptr, unsigned long int);
976 #define mpz_powm __gmpz_powm
977 __GMP_DECLSPEC void mpz_powm (mpz_ptr, mpz_srcptr, mpz_srcptr, mpz_srcptr);
979 #define mpz_powm_sec __gmpz_powm_sec
980 __GMP_DECLSPEC void mpz_powm_sec (mpz_ptr, mpz_srcptr, mpz_srcptr, mpz_srcptr);
982 #define mpz_powm_ui __gmpz_powm_ui
983 __GMP_DECLSPEC void mpz_powm_ui (mpz_ptr, mpz_srcptr, unsigned long int, mpz_srcptr);
985 #define mpz_probab_prime_p __gmpz_probab_prime_p
986 __GMP_DECLSPEC int mpz_probab_prime_p (mpz_srcptr, int) __GMP_ATTRIBUTE_PURE;
988 #define mpz_random __gmpz_random
989 __GMP_DECLSPEC void mpz_random (mpz_ptr, mp_size_t);
991 #define mpz_random2 __gmpz_random2
992 __GMP_DECLSPEC void mpz_random2 (mpz_ptr, mp_size_t);
994 #define mpz_realloc2 __gmpz_realloc2
995 __GMP_DECLSPEC void mpz_realloc2 (mpz_ptr, mp_bitcnt_t);
997 #define mpz_remove __gmpz_remove
998 __GMP_DECLSPEC mp_bitcnt_t mpz_remove (mpz_ptr, mpz_srcptr, mpz_srcptr);
1000 #define mpz_root __gmpz_root
1001 __GMP_DECLSPEC int mpz_root (mpz_ptr, mpz_srcptr, unsigned long int);
1003 #define mpz_rootrem __gmpz_rootrem
1004 __GMP_DECLSPEC void mpz_rootrem (mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int);
1006 #define mpz_rrandomb __gmpz_rrandomb
1007 __GMP_DECLSPEC void mpz_rrandomb (mpz_ptr, gmp_randstate_t, mp_bitcnt_t);
1009 #define mpz_scan0 __gmpz_scan0
1010 __GMP_DECLSPEC mp_bitcnt_t mpz_scan0 (mpz_srcptr, mp_bitcnt_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1012 #define mpz_scan1 __gmpz_scan1
1013 __GMP_DECLSPEC mp_bitcnt_t mpz_scan1 (mpz_srcptr, mp_bitcnt_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1015 #define mpz_set __gmpz_set
1016 __GMP_DECLSPEC void mpz_set (mpz_ptr, mpz_srcptr);
1018 #define mpz_set_d __gmpz_set_d
1019 __GMP_DECLSPEC void mpz_set_d (mpz_ptr, double);
1021 #define mpz_set_f __gmpz_set_f
1022 __GMP_DECLSPEC void mpz_set_f (mpz_ptr, mpf_srcptr);
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);
1029 #define mpz_set_si __gmpz_set_si
1030 __GMP_DECLSPEC void mpz_set_si (mpz_ptr, signed long int);
1032 #define mpz_set_str __gmpz_set_str
1033 __GMP_DECLSPEC int mpz_set_str (mpz_ptr, const char *, int);
1035 #define mpz_set_ui __gmpz_set_ui
1036 __GMP_DECLSPEC void mpz_set_ui (mpz_ptr, unsigned long int);
1038 #define mpz_setbit __gmpz_setbit
1039 __GMP_DECLSPEC void mpz_setbit (mpz_ptr, mp_bitcnt_t);
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;
1046 #define mpz_sizeinbase __gmpz_sizeinbase
1047 __GMP_DECLSPEC size_t mpz_sizeinbase (mpz_srcptr, int) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1049 #define mpz_sqrt __gmpz_sqrt
1050 __GMP_DECLSPEC void mpz_sqrt (mpz_ptr, mpz_srcptr);
1052 #define mpz_sqrtrem __gmpz_sqrtrem
1053 __GMP_DECLSPEC void mpz_sqrtrem (mpz_ptr, mpz_ptr, mpz_srcptr);
1055 #define mpz_sub __gmpz_sub
1056 __GMP_DECLSPEC void mpz_sub (mpz_ptr, mpz_srcptr, mpz_srcptr);
1058 #define mpz_sub_ui __gmpz_sub_ui
1059 __GMP_DECLSPEC void mpz_sub_ui (mpz_ptr, mpz_srcptr, unsigned long int);
1061 #define mpz_ui_sub __gmpz_ui_sub
1062 __GMP_DECLSPEC void mpz_ui_sub (mpz_ptr, unsigned long int, mpz_srcptr);
1064 #define mpz_submul __gmpz_submul
1065 __GMP_DECLSPEC void mpz_submul (mpz_ptr, mpz_srcptr, mpz_srcptr);
1067 #define mpz_submul_ui __gmpz_submul_ui
1068 __GMP_DECLSPEC void mpz_submul_ui (mpz_ptr, mpz_srcptr, unsigned long int);
1070 #define mpz_swap __gmpz_swap
1071 __GMP_DECLSPEC void mpz_swap (mpz_ptr, mpz_ptr) __GMP_NOTHROW;
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;
1076 #define mpz_tdiv_q __gmpz_tdiv_q
1077 __GMP_DECLSPEC void mpz_tdiv_q (mpz_ptr, mpz_srcptr, mpz_srcptr);
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);
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);
1085 #define mpz_tdiv_qr __gmpz_tdiv_qr
1086 __GMP_DECLSPEC void mpz_tdiv_qr (mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr);
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);
1091 #define mpz_tdiv_r __gmpz_tdiv_r
1092 __GMP_DECLSPEC void mpz_tdiv_r (mpz_ptr, mpz_srcptr, mpz_srcptr);
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);
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);
1100 #define mpz_tstbit __gmpz_tstbit
1101 __GMP_DECLSPEC int mpz_tstbit (mpz_srcptr, mp_bitcnt_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
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);
1106 #define mpz_urandomb __gmpz_urandomb
1107 __GMP_DECLSPEC void mpz_urandomb (mpz_ptr, gmp_randstate_t, mp_bitcnt_t);
1109 #define mpz_urandomm __gmpz_urandomm
1110 __GMP_DECLSPEC void mpz_urandomm (mpz_ptr, gmp_randstate_t, mpz_srcptr);
1112 #define mpz_xor __gmpz_xor
1113 #define mpz_eor __gmpz_xor
1114 __GMP_DECLSPEC void mpz_xor (mpz_ptr, mpz_srcptr, mpz_srcptr);
1116 #define mpz_limbs_read __gmpz_limbs_read
1117 __GMP_DECLSPEC mp_srcptr mpz_limbs_read (mpz_srcptr);
1119 #define mpz_limbs_write __gmpz_limbs_write
1120 __GMP_DECLSPEC mp_ptr mpz_limbs_write (mpz_ptr, mp_size_t);
1122 #define mpz_limbs_modify __gmpz_limbs_modify
1123 __GMP_DECLSPEC mp_ptr mpz_limbs_modify (mpz_ptr, mp_size_t);
1125 #define mpz_limbs_finish __gmpz_limbs_finish
1126 __GMP_DECLSPEC void mpz_limbs_finish (mpz_ptr, mp_size_t);
1128 #define mpz_roinit_n __gmpz_roinit_n
1129 __GMP_DECLSPEC mpz_srcptr mpz_roinit_n (mpz_ptr, mp_srcptr, mp_size_t);
1131 #define MPZ_ROINIT_N(xp, xs) {{0, (xs),(xp) }}
1133 /**************** Rational (i.e. Q) routines. ****************/
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);
1140 #define mpq_add __gmpq_add
1141 __GMP_DECLSPEC void mpq_add (mpq_ptr, mpq_srcptr, mpq_srcptr);
1143 #define mpq_canonicalize __gmpq_canonicalize
1144 __GMP_DECLSPEC void mpq_canonicalize (mpq_ptr);
1146 #define mpq_clear __gmpq_clear
1147 __GMP_DECLSPEC void mpq_clear (mpq_ptr);
1149 #define mpq_clears __gmpq_clears
1150 __GMP_DECLSPEC void mpq_clears (mpq_ptr, ...);
1152 #define mpq_cmp __gmpq_cmp
1153 __GMP_DECLSPEC int mpq_cmp (mpq_srcptr, mpq_srcptr) __GMP_ATTRIBUTE_PURE;
1155 #define _mpq_cmp_si __gmpq_cmp_si
1156 __GMP_DECLSPEC int _mpq_cmp_si (mpq_srcptr, long, unsigned long) __GMP_ATTRIBUTE_PURE;
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;
1161 #define mpq_cmp_z __gmpq_cmp_z
1162 __GMP_DECLSPEC int mpq_cmp_z (mpq_srcptr, mpz_srcptr) __GMP_ATTRIBUTE_PURE;
1164 #define mpq_div __gmpq_div
1165 __GMP_DECLSPEC void mpq_div (mpq_ptr, mpq_srcptr, mpq_srcptr);
1167 #define mpq_div_2exp __gmpq_div_2exp
1168 __GMP_DECLSPEC void mpq_div_2exp (mpq_ptr, mpq_srcptr, mp_bitcnt_t);
1170 #define mpq_equal __gmpq_equal
1171 __GMP_DECLSPEC int mpq_equal (mpq_srcptr, mpq_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1173 #define mpq_get_num __gmpq_get_num
1174 __GMP_DECLSPEC void mpq_get_num (mpz_ptr, mpq_srcptr);
1176 #define mpq_get_den __gmpq_get_den
1177 __GMP_DECLSPEC void mpq_get_den (mpz_ptr, mpq_srcptr);
1179 #define mpq_get_d __gmpq_get_d
1180 __GMP_DECLSPEC double mpq_get_d (mpq_srcptr) __GMP_ATTRIBUTE_PURE;
1182 #define mpq_get_str __gmpq_get_str
1183 __GMP_DECLSPEC char *mpq_get_str (char *, int, mpq_srcptr);
1185 #define mpq_init __gmpq_init
1186 __GMP_DECLSPEC void mpq_init (mpq_ptr);
1188 #define mpq_inits __gmpq_inits
1189 __GMP_DECLSPEC void mpq_inits (mpq_ptr, ...);
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);
1196 #define mpq_inv __gmpq_inv
1197 __GMP_DECLSPEC void mpq_inv (mpq_ptr, mpq_srcptr);
1199 #define mpq_mul __gmpq_mul
1200 __GMP_DECLSPEC void mpq_mul (mpq_ptr, mpq_srcptr, mpq_srcptr);
1202 #define mpq_mul_2exp __gmpq_mul_2exp
1203 __GMP_DECLSPEC void mpq_mul_2exp (mpq_ptr, mpq_srcptr, mp_bitcnt_t);
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);
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);
1215 #define mpq_set __gmpq_set
1216 __GMP_DECLSPEC void mpq_set (mpq_ptr, mpq_srcptr);
1218 #define mpq_set_d __gmpq_set_d
1219 __GMP_DECLSPEC void mpq_set_d (mpq_ptr, double);
1221 #define mpq_set_den __gmpq_set_den
1222 __GMP_DECLSPEC void mpq_set_den (mpq_ptr, mpz_srcptr);
1224 #define mpq_set_f __gmpq_set_f
1225 __GMP_DECLSPEC void mpq_set_f (mpq_ptr, mpf_srcptr);
1227 #define mpq_set_num __gmpq_set_num
1228 __GMP_DECLSPEC void mpq_set_num (mpq_ptr, mpz_srcptr);
1230 #define mpq_set_si __gmpq_set_si
1231 __GMP_DECLSPEC void mpq_set_si (mpq_ptr, signed long int, unsigned long int);
1233 #define mpq_set_str __gmpq_set_str
1234 __GMP_DECLSPEC int mpq_set_str (mpq_ptr, const char *, int);
1236 #define mpq_set_ui __gmpq_set_ui
1237 __GMP_DECLSPEC void mpq_set_ui (mpq_ptr, unsigned long int, unsigned long int);
1239 #define mpq_set_z __gmpq_set_z
1240 __GMP_DECLSPEC void mpq_set_z (mpq_ptr, mpz_srcptr);
1242 #define mpq_sub __gmpq_sub
1243 __GMP_DECLSPEC void mpq_sub (mpq_ptr, mpq_srcptr, mpq_srcptr);
1245 #define mpq_swap __gmpq_swap
1246 __GMP_DECLSPEC void mpq_swap (mpq_ptr, mpq_ptr) __GMP_NOTHROW;
1249 /**************** Float (i.e. F) routines. ****************/
1251 #define mpf_abs __gmpf_abs
1252 __GMP_DECLSPEC void mpf_abs (mpf_ptr, mpf_srcptr);
1254 #define mpf_add __gmpf_add
1255 __GMP_DECLSPEC void mpf_add (mpf_ptr, mpf_srcptr, mpf_srcptr);
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);
1262 #define mpf_clear __gmpf_clear
1263 __GMP_DECLSPEC void mpf_clear (mpf_ptr);
1265 #define mpf_clears __gmpf_clears
1266 __GMP_DECLSPEC void mpf_clears (mpf_ptr, ...);
1268 #define mpf_cmp __gmpf_cmp
1269 __GMP_DECLSPEC int mpf_cmp (mpf_srcptr, mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1271 #define mpf_cmp_z __gmpf_cmp_z
1272 __GMP_DECLSPEC int mpf_cmp_z (mpf_srcptr, mpz_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1274 #define mpf_cmp_d __gmpf_cmp_d
1275 __GMP_DECLSPEC int mpf_cmp_d (mpf_srcptr, double) __GMP_ATTRIBUTE_PURE;
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;
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;
1283 #define mpf_div __gmpf_div
1284 __GMP_DECLSPEC void mpf_div (mpf_ptr, mpf_srcptr, mpf_srcptr);
1286 #define mpf_div_2exp __gmpf_div_2exp
1287 __GMP_DECLSPEC void mpf_div_2exp (mpf_ptr, mpf_srcptr, mp_bitcnt_t);
1289 #define mpf_div_ui __gmpf_div_ui
1290 __GMP_DECLSPEC void mpf_div_ui (mpf_ptr, mpf_srcptr, unsigned long int);
1292 #define mpf_dump __gmpf_dump
1293 __GMP_DECLSPEC void mpf_dump (mpf_srcptr);
1295 #define mpf_eq __gmpf_eq
1296 __GMP_DECLSPEC int mpf_eq (mpf_srcptr, mpf_srcptr, mp_bitcnt_t) __GMP_ATTRIBUTE_PURE;
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;
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;
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;
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;
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;
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;
1316 #define mpf_floor __gmpf_floor
1317 __GMP_DECLSPEC void mpf_floor (mpf_ptr, mpf_srcptr);
1319 #define mpf_get_d __gmpf_get_d
1320 __GMP_DECLSPEC double mpf_get_d (mpf_srcptr) __GMP_ATTRIBUTE_PURE;
1322 #define mpf_get_d_2exp __gmpf_get_d_2exp
1323 __GMP_DECLSPEC double mpf_get_d_2exp (signed long int *, mpf_srcptr);
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;
1328 #define mpf_get_prec __gmpf_get_prec
1329 __GMP_DECLSPEC mp_bitcnt_t mpf_get_prec (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1331 #define mpf_get_si __gmpf_get_si
1332 __GMP_DECLSPEC long mpf_get_si (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1334 #define mpf_get_str __gmpf_get_str
1335 __GMP_DECLSPEC char *mpf_get_str (char *, mp_exp_t *, int, size_t, mpf_srcptr);
1337 #define mpf_get_ui __gmpf_get_ui
1338 __GMP_DECLSPEC unsigned long mpf_get_ui (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1340 #define mpf_init __gmpf_init
1341 __GMP_DECLSPEC void mpf_init (mpf_ptr);
1343 #define mpf_init2 __gmpf_init2
1344 __GMP_DECLSPEC void mpf_init2 (mpf_ptr, mp_bitcnt_t);
1346 #define mpf_inits __gmpf_inits
1347 __GMP_DECLSPEC void mpf_inits (mpf_ptr, ...);
1349 #define mpf_init_set __gmpf_init_set
1350 __GMP_DECLSPEC void mpf_init_set (mpf_ptr, mpf_srcptr);
1352 #define mpf_init_set_d __gmpf_init_set_d
1353 __GMP_DECLSPEC void mpf_init_set_d (mpf_ptr, double);
1355 #define mpf_init_set_si __gmpf_init_set_si
1356 __GMP_DECLSPEC void mpf_init_set_si (mpf_ptr, signed long int);
1358 #define mpf_init_set_str __gmpf_init_set_str
1359 __GMP_DECLSPEC int mpf_init_set_str (mpf_ptr, const char *, int);
1361 #define mpf_init_set_ui __gmpf_init_set_ui
1362 __GMP_DECLSPEC void mpf_init_set_ui (mpf_ptr, unsigned long int);
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);
1369 #define mpf_integer_p __gmpf_integer_p
1370 __GMP_DECLSPEC int mpf_integer_p (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1372 #define mpf_mul __gmpf_mul
1373 __GMP_DECLSPEC void mpf_mul (mpf_ptr, mpf_srcptr, mpf_srcptr);
1375 #define mpf_mul_2exp __gmpf_mul_2exp
1376 __GMP_DECLSPEC void mpf_mul_2exp (mpf_ptr, mpf_srcptr, mp_bitcnt_t);
1378 #define mpf_mul_ui __gmpf_mul_ui
1379 __GMP_DECLSPEC void mpf_mul_ui (mpf_ptr, mpf_srcptr, unsigned long int);
1381 #define mpf_neg __gmpf_neg
1382 __GMP_DECLSPEC void mpf_neg (mpf_ptr, mpf_srcptr);
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);
1389 #define mpf_pow_ui __gmpf_pow_ui
1390 __GMP_DECLSPEC void mpf_pow_ui (mpf_ptr, mpf_srcptr, unsigned long int);
1392 #define mpf_random2 __gmpf_random2
1393 __GMP_DECLSPEC void mpf_random2 (mpf_ptr, mp_size_t, mp_exp_t);
1395 #define mpf_reldiff __gmpf_reldiff
1396 __GMP_DECLSPEC void mpf_reldiff (mpf_ptr, mpf_srcptr, mpf_srcptr);
1398 #define mpf_set __gmpf_set
1399 __GMP_DECLSPEC void mpf_set (mpf_ptr, mpf_srcptr);
1401 #define mpf_set_d __gmpf_set_d
1402 __GMP_DECLSPEC void mpf_set_d (mpf_ptr, double);
1404 #define mpf_set_default_prec __gmpf_set_default_prec
1405 __GMP_DECLSPEC void mpf_set_default_prec (mp_bitcnt_t) __GMP_NOTHROW;
1407 #define mpf_set_prec __gmpf_set_prec
1408 __GMP_DECLSPEC void mpf_set_prec (mpf_ptr, mp_bitcnt_t);
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;
1413 #define mpf_set_q __gmpf_set_q
1414 __GMP_DECLSPEC void mpf_set_q (mpf_ptr, mpq_srcptr);
1416 #define mpf_set_si __gmpf_set_si
1417 __GMP_DECLSPEC void mpf_set_si (mpf_ptr, signed long int);
1419 #define mpf_set_str __gmpf_set_str
1420 __GMP_DECLSPEC int mpf_set_str (mpf_ptr, const char *, int);
1422 #define mpf_set_ui __gmpf_set_ui
1423 __GMP_DECLSPEC void mpf_set_ui (mpf_ptr, unsigned long int);
1425 #define mpf_set_z __gmpf_set_z
1426 __GMP_DECLSPEC void mpf_set_z (mpf_ptr, mpz_srcptr);
1428 #define mpf_size __gmpf_size
1429 __GMP_DECLSPEC size_t mpf_size (mpf_srcptr) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1431 #define mpf_sqrt __gmpf_sqrt
1432 __GMP_DECLSPEC void mpf_sqrt (mpf_ptr, mpf_srcptr);
1434 #define mpf_sqrt_ui __gmpf_sqrt_ui
1435 __GMP_DECLSPEC void mpf_sqrt_ui (mpf_ptr, unsigned long int);
1437 #define mpf_sub __gmpf_sub
1438 __GMP_DECLSPEC void mpf_sub (mpf_ptr, mpf_srcptr, mpf_srcptr);
1440 #define mpf_sub_ui __gmpf_sub_ui
1441 __GMP_DECLSPEC void mpf_sub_ui (mpf_ptr, mpf_srcptr, unsigned long int);
1443 #define mpf_swap __gmpf_swap
1444 __GMP_DECLSPEC void mpf_swap (mpf_ptr, mpf_ptr) __GMP_NOTHROW;
1446 #define mpf_trunc __gmpf_trunc
1447 __GMP_DECLSPEC void mpf_trunc (mpf_ptr, mpf_srcptr);
1449 #define mpf_ui_div __gmpf_ui_div
1450 __GMP_DECLSPEC void mpf_ui_div (mpf_ptr, unsigned long int, mpf_srcptr);
1452 #define mpf_ui_sub __gmpf_ui_sub
1453 __GMP_DECLSPEC void mpf_ui_sub (mpf_ptr, unsigned long int, mpf_srcptr);
1455 #define mpf_urandomb __gmpf_urandomb
1456 __GMP_DECLSPEC void mpf_urandomb (mpf_t, gmp_randstate_t, mp_bitcnt_t);
1459 /************ Low level positive-integer (i.e. N) routines. ************/
1461 /* This is ugly, but we need to make user calls reach the prefixed function. */
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);
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;
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);
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);
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;
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;
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);
1492 #define mpn_divexact_by3(dst,src,size) \
1493 mpn_divexact_by3c (dst, src, size, __GMP_CAST (mp_limb_t, 0))
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);
1498 #define mpn_divmod_1(qp,np,nsize,dlimb) \
1499 mpn_divrem_1 (qp, __GMP_CAST (mp_size_t, 0), np, nsize, dlimb)
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);
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);
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);
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);
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);
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);
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;
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;
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);
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);
1531 #define mpn_get_str __MPN(get_str)
1532 __GMP_DECLSPEC size_t mpn_get_str (unsigned char *, int, mp_ptr, mp_size_t);
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;
1537 #define mpn_lshift __MPN(lshift)
1538 __GMP_DECLSPEC mp_limb_t mpn_lshift (mp_ptr, mp_srcptr, mp_size_t, unsigned int);
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;
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);
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);
1549 #define mpn_mul_n __MPN(mul_n)
1550 __GMP_DECLSPEC void mpn_mul_n (mp_ptr, mp_srcptr, mp_srcptr, mp_size_t);
1552 #define mpn_sqr __MPN(sqr)
1553 __GMP_DECLSPEC void mpn_sqr (mp_ptr, mp_srcptr, mp_size_t);
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);
1560 #define mpn_com __MPN(com)
1561 __GMP_DECLSPEC void mpn_com (mp_ptr, mp_srcptr, mp_size_t);
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;
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;
1569 #define mpn_popcount __MPN(popcount)
1570 __GMP_DECLSPEC mp_bitcnt_t mpn_popcount (mp_srcptr, mp_size_t) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
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);
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;
1579 #define mpn_random __MPN(random)
1580 __GMP_DECLSPEC void mpn_random (mp_ptr, mp_size_t);
1582 #define mpn_random2 __MPN(random2)
1583 __GMP_DECLSPEC void mpn_random2 (mp_ptr, mp_size_t);
1585 #define mpn_rshift __MPN(rshift)
1586 __GMP_DECLSPEC mp_limb_t mpn_rshift (mp_ptr, mp_srcptr, mp_size_t, unsigned int);
1588 #define mpn_scan0 __MPN(scan0)
1589 __GMP_DECLSPEC mp_bitcnt_t mpn_scan0 (mp_srcptr, mp_bitcnt_t) __GMP_ATTRIBUTE_PURE;
1591 #define mpn_scan1 __MPN(scan1)
1592 __GMP_DECLSPEC mp_bitcnt_t mpn_scan1 (mp_srcptr, mp_bitcnt_t) __GMP_ATTRIBUTE_PURE;
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);
1597 #define mpn_sizeinbase __MPN(sizeinbase)
1598 __GMP_DECLSPEC size_t mpn_sizeinbase (mp_srcptr, mp_size_t, int);
1600 #define mpn_sqrtrem __MPN(sqrtrem)
1601 __GMP_DECLSPEC mp_size_t mpn_sqrtrem (mp_ptr, mp_ptr, mp_srcptr, mp_size_t);
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);
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;
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);
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);
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);
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);
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);
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);
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;
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;
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);
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;
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;
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;
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);
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;
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;
1697 /**************** mpz inlines ****************/
1699 /* The following are provided as inlines where possible, but always exist as
1700 library functions too, for binary compatibility.
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.
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. */
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)
1717 if (__gmp_w != __gmp_u)
1718 mpz_set (__gmp_w, __gmp_u);
1719 __gmp_w->_mp_size = __GMP_ABS (__gmp_w->_mp_size);
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));
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)));
1736 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_fits_uint_p)
1737 #if ! defined (__GMP_FORCE_mpz_fits_uint_p)
1741 mpz_fits_uint_p (mpz_srcptr __gmp_z) __GMP_NOTHROW
1743 __GMPZ_FITS_UTYPE_P (__gmp_z, UINT_MAX);
1747 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_fits_ulong_p)
1748 #if ! defined (__GMP_FORCE_mpz_fits_ulong_p)
1752 mpz_fits_ulong_p (mpz_srcptr __gmp_z) __GMP_NOTHROW
1754 __GMPZ_FITS_UTYPE_P (__gmp_z, ULONG_MAX);
1758 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_fits_ushort_p)
1759 #if ! defined (__GMP_FORCE_mpz_fits_ushort_p)
1763 mpz_fits_ushort_p (mpz_srcptr __gmp_z) __GMP_NOTHROW
1765 __GMPZ_FITS_UTYPE_P (__gmp_z, USHRT_MAX);
1769 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_get_ui)
1770 #if ! defined (__GMP_FORCE_mpz_get_ui)
1774 mpz_get_ui (mpz_srcptr __gmp_z) __GMP_NOTHROW
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);
1787 /* limb==long and nails, need two limbs when available */
1788 __gmp_n = __GMP_ABS (__gmp_n);
1790 return (__gmp_n != 0 ? __gmp_l : 0);
1792 return __gmp_l + (__gmp_p[1] << GMP_NUMB_BITS);
1797 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_getlimbn)
1798 #if ! defined (__GMP_FORCE_mpz_getlimbn)
1802 mpz_getlimbn (mpz_srcptr __gmp_z, mp_size_t __gmp_n) __GMP_NOTHROW
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;
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)
1815 if (__gmp_w != __gmp_u)
1816 mpz_set (__gmp_w, __gmp_u);
1817 __gmp_w->_mp_size = - __gmp_w->_mp_size;
1821 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_perfect_square_p)
1822 #if ! defined (__GMP_FORCE_mpz_perfect_square_p)
1826 mpz_perfect_square_p (mpz_srcptr __gmp_a)
1828 mp_size_t __gmp_asize;
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;
1839 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_popcount)
1840 #if ! defined (__GMP_FORCE_mpz_popcount)
1844 mpz_popcount (mpz_srcptr __gmp_u) __GMP_NOTHROW
1846 mp_size_t __gmp_usize;
1847 mp_bitcnt_t __gmp_result;
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;
1857 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_set_q)
1858 #if ! defined (__GMP_FORCE_mpz_set_q)
1862 mpz_set_q (mpz_ptr __gmp_w, mpq_srcptr __gmp_u)
1864 mpz_tdiv_q (__gmp_w, mpq_numref (__gmp_u), mpq_denref (__gmp_u));
1868 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_size)
1869 #if ! defined (__GMP_FORCE_mpz_size)
1873 mpz_size (mpz_srcptr __gmp_z) __GMP_NOTHROW
1875 return __GMP_ABS (__gmp_z->_mp_size);
1880 /**************** mpq inlines ****************/
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)
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);
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)
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;
1903 /**************** mpn inlines ****************/
1905 /* The comments with __GMPN_ADD_1 below apply here too.
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
1912 ysize==0 isn't a documented feature, but is used internally in a few
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. */
1920 #define __GMPN_AORS(cout, wp, xp, xsize, yp, ysize, FUNCTION, TEST) \
1922 mp_size_t __gmp_i; \
1923 mp_limb_t __gmp_x; \
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)); */ \
1930 __gmp_i = (ysize); \
1933 if (FUNCTION (wp, xp, yp, __gmp_i)) \
1937 if (__gmp_i >= (xsize)) \
1942 __gmp_x = (xp)[__gmp_i]; \
1948 __GMPN_COPY_REST (wp, xp, xsize, __gmp_i); \
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))
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.
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.
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.
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. */
1986 #if GMP_NAIL_BITS == 0
1987 #define __GMPN_AORS_1(cout, dst, src, n, v, OP, CB) \
1989 mp_size_t __gmp_i; \
1990 mp_limb_t __gmp_x, __gmp_r; \
1992 /* ASSERT ((n) >= 1); */ \
1993 /* ASSERT (MPN_SAME_OR_SEPARATE_P (dst, src, n)); */ \
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))) \
2001 for (__gmp_i = 1; __gmp_i < (n);) \
2003 __gmp_x = (src)[__gmp_i]; \
2004 __gmp_r = __gmp_x OP 1; \
2005 (dst)[__gmp_i] = __gmp_r; \
2007 if (!CB (__gmp_r, __gmp_x, 1)) \
2009 if ((src) != (dst)) \
2010 __GMPN_COPY_REST (dst, src, n, __gmp_i); \
2018 if ((src) != (dst)) \
2019 __GMPN_COPY_REST (dst, src, n, 1); \
2025 #if GMP_NAIL_BITS >= 1
2026 #define __GMPN_AORS_1(cout, dst, src, n, v, OP, CB) \
2028 mp_size_t __gmp_i; \
2029 mp_limb_t __gmp_x, __gmp_r; \
2031 /* ASSERT ((n) >= 1); */ \
2032 /* ASSERT (MPN_SAME_OR_SEPARATE_P (dst, src, n)); */ \
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) \
2040 for (__gmp_i = 1; __gmp_i < (n);) \
2042 __gmp_x = (src)[__gmp_i]; \
2043 __gmp_r = __gmp_x OP 1; \
2044 (dst)[__gmp_i] = __gmp_r & GMP_NUMB_MASK; \
2046 if (__gmp_r >> GMP_NUMB_BITS == 0) \
2048 if ((src) != (dst)) \
2049 __GMPN_COPY_REST (dst, src, n, __gmp_i); \
2057 if ((src) != (dst)) \
2058 __GMPN_COPY_REST (dst, src, n, 1); \
2064 #define __GMPN_ADDCB(r,x,y) ((r) < (y))
2065 #define __GMPN_SUBCB(r,x,y) ((x) < (y))
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)
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) \
2078 mp_size_t __gmp_i; \
2079 mp_limb_t __gmp_x, __gmp_y; \
2081 /* ASSERT ((size) >= 0); */ \
2085 while (--__gmp_i >= 0) \
2087 __gmp_x = (xp)[__gmp_i]; \
2088 __gmp_y = (yp)[__gmp_i]; \
2089 if (__gmp_x != __gmp_y) \
2091 /* Cannot use __gmp_x - __gmp_y, may overflow an "int" */ \
2092 (result) = (__gmp_x > __gmp_y ? 1 : -1); \
2099 #if defined (__GMPN_COPY) && ! defined (__GMPN_COPY_REST)
2100 #define __GMPN_COPY_REST(dst, src, size, start) \
2102 /* ASSERT ((start) >= 0); */ \
2103 /* ASSERT ((start) <= (size)); */ \
2104 __GMPN_COPY ((dst)+(start), (src)+(start), (size)-(start)); \
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,
2111 #if ! defined (__GMPN_COPY_REST)
2112 #define __GMPN_COPY_REST(dst, src, size, start) \
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]; \
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). */
2129 #if ! defined (__GMPN_COPY)
2130 #define __GMPN_COPY(dst, src, size) __GMPN_COPY_REST (dst, src, size, 0)
2134 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_add)
2135 #if ! defined (__GMP_FORCE_mpn_add)
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)
2142 __GMPN_ADD (__gmp_c, __gmp_wp, __gmp_xp, __gmp_xsize, __gmp_yp, __gmp_ysize);
2147 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_add_1)
2148 #if ! defined (__GMP_FORCE_mpn_add_1)
2152 mpn_add_1 (mp_ptr __gmp_dst, mp_srcptr __gmp_src, mp_size_t __gmp_size, mp_limb_t __gmp_n) __GMP_NOTHROW
2155 __GMPN_ADD_1 (__gmp_c, __gmp_dst, __gmp_src, __gmp_size, __gmp_n);
2160 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_cmp)
2161 #if ! defined (__GMP_FORCE_mpn_cmp)
2165 mpn_cmp (mp_srcptr __gmp_xp, mp_srcptr __gmp_yp, mp_size_t __gmp_size) __GMP_NOTHROW
2168 __GMPN_CMP (__gmp_result, __gmp_xp, __gmp_yp, __gmp_size);
2169 return __gmp_result;
2173 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_zero_p)
2174 #if ! defined (__GMP_FORCE_mpn_zero_p)
2178 mpn_zero_p (mp_srcptr __gmp_p, mp_size_t __gmp_n) __GMP_NOTHROW
2180 /* if (__GMP_LIKELY (__gmp_n > 0)) */
2182 if (__gmp_p[--__gmp_n] != 0)
2184 } while (__gmp_n != 0);
2189 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_sub)
2190 #if ! defined (__GMP_FORCE_mpn_sub)
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)
2197 __GMPN_SUB (__gmp_c, __gmp_wp, __gmp_xp, __gmp_xsize, __gmp_yp, __gmp_ysize);
2202 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_sub_1)
2203 #if ! defined (__GMP_FORCE_mpn_sub_1)
2207 mpn_sub_1 (mp_ptr __gmp_dst, mp_srcptr __gmp_src, mp_size_t __gmp_size, mp_limb_t __gmp_n) __GMP_NOTHROW
2210 __GMPN_SUB_1 (__gmp_c, __gmp_dst, __gmp_src, __gmp_size, __gmp_n);
2215 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_neg)
2216 #if ! defined (__GMP_FORCE_mpn_neg)
2220 mpn_neg (mp_ptr __gmp_rp, mp_srcptr __gmp_up, mp_size_t __gmp_n)
2222 while (*__gmp_up == 0) /* Low zero limbs are unchanged by negation. */
2225 if (!--__gmp_n) /* All zero */
2227 ++__gmp_up; ++__gmp_rp;
2230 *__gmp_rp = (- *__gmp_up) & GMP_NUMB_MASK;
2232 if (--__gmp_n) /* Higher limbs get complemented. */
2233 mpn_com (++__gmp_rp, ++__gmp_up, __gmp_n);
2239 #if defined (__cplusplus)
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)
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))
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)
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
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))
2282 /**************** C++ routines ****************/
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);
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)
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))
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
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
2325 /* Define CC and CFLAGS which were used to build this version of GMP */
2326 #define __GMP_CC "gcc"
2327 #define __GMP_CFLAGS "-O2 -pedantic -fomit-frame-pointer -m64 -mtune=znver2 -march=znver2"
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 1
2333 #define __GNU_MP_RELEASE (__GNU_MP_VERSION * 10000 + __GNU_MP_VERSION_MINOR * 100 + __GNU_MP_VERSION_PATCHLEVEL)
2336 #endif /* __GMP_H__ */