1 /* Definitions for GNU multiple precision functions. -*- mode: c -*-
3 Copyright 1991, 1993, 1994, 1995, 1996, 1997, 1999, 2000, 2001, 2002, 2003,
4 2004, 2005, 2006, 2007, 2008, 2009 Free Software Foundation, Inc.
6 This file is part of the GNU MP Library.
8 The GNU MP Library is free software; you can redistribute it and/or modify
9 it under the terms of the GNU Lesser General Public License as published by
10 the Free Software Foundation; either version 3 of the License, or (at your
11 option) any later version.
13 The GNU MP Library is distributed in the hope that it will be useful, but
14 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
15 or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
16 License for more details.
18 You should have received a copy of the GNU Lesser General Public License
19 along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */
23 #if defined (__cplusplus)
24 #include <iosfwd> /* for std::istream, std::ostream, std::string */
29 /* Instantiated by configure. */
30 #if ! defined (__GMP_WITHIN_CONFIGURE)
32 #define __GMP_BITS_PER_MP_LIMB 32
33 #elif defined(__x86_64__)
34 #define __GMP_BITS_PER_MP_LIMB 64
38 #define __GMP_HAVE_HOST_CPU_FAMILY_power 0
39 #define __GMP_HAVE_HOST_CPU_FAMILY_powerpc 0
40 #define GMP_LIMB_BITS __GMP_BITS_PER_MP_LIMB
41 #define GMP_NAIL_BITS 0
43 #define GMP_NUMB_BITS (GMP_LIMB_BITS - GMP_NAIL_BITS)
44 #define GMP_NUMB_MASK ((~ __GMP_CAST (mp_limb_t, 0)) >> GMP_NAIL_BITS)
45 #define GMP_NUMB_MAX GMP_NUMB_MASK
46 #define GMP_NAIL_MASK (~ GMP_NUMB_MASK)
49 /* The following (everything under ifndef __GNU_MP__) must be identical in
50 gmp.h and mp.h to allow both to be included in an application or during
55 #define __need_size_t /* tell gcc stddef.h we only want size_t */
56 #if defined (__cplusplus)
57 #include <cstddef> /* for size_t */
59 #include <stddef.h> /* for size_t */
63 /* Instantiated by configure. */
64 #if ! defined (__GMP_WITHIN_CONFIGURE)
65 /* #undef _LONG_LONG_LIMB */
66 #define __GMP_LIBGMP_DLL 0
70 /* __STDC__ - some ANSI compilers define this only to 0, hence the use of
71 "defined" and not "__STDC__-0". In particular Sun workshop C 5.0
72 sets __STDC__ to 0, but requires "##" for token pasting.
74 _AIX - gnu ansidecl.h asserts that all known AIX compilers are ANSI but
75 don't always define __STDC__.
77 __DECC - current versions of DEC C (5.9 for instance) for alpha are ANSI,
78 but don't define __STDC__ in their default mode. Don't know if old
79 versions might have been K&R, but let's not worry about that unless
80 someone is still using one.
82 _mips - gnu ansidecl.h says the RISC/OS MIPS compiler is ANSI in SVR4
83 mode, but doesn't define __STDC__.
85 _MSC_VER - Microsoft C is ANSI, but __STDC__ is undefined unless the /Za
86 option is given (in which case it's 1).
88 _WIN32 - tested for by gnu ansidecl.h, no doubt on the assumption that
89 all w32 compilers are ansi.
91 Note: This same set of tests is used by gen-psqr.c and
92 demos/expr/expr-impl.h, so if anything needs adding, then be sure to
95 #if defined (__STDC__) \
96 || defined (__cplusplus) \
99 || (defined (__mips) && defined (_SYSTYPE_SVR4)) \
100 || defined (_MSC_VER) \
102 #define __GMP_HAVE_CONST 1
103 #define __GMP_HAVE_PROTOTYPES 1
104 #define __GMP_HAVE_TOKEN_PASTE 1
106 #define __GMP_HAVE_CONST 0
107 #define __GMP_HAVE_PROTOTYPES 0
108 #define __GMP_HAVE_TOKEN_PASTE 0
113 #define __gmp_const const
114 #define __gmp_signed signed
121 /* __GMP_DECLSPEC supports Windows DLL versions of libgmp, and is empty in
122 all other circumstances.
124 When compiling objects for libgmp, __GMP_DECLSPEC is an export directive,
125 or when compiling for an application it's an import directive. The two
126 cases are differentiated by __GMP_WITHIN_GMP defined by the GMP Makefiles
127 (and not defined from an application).
129 __GMP_DECLSPEC_XX is similarly used for libgmpxx. __GMP_WITHIN_GMPXX
130 indicates when building libgmpxx, and in that case libgmpxx functions are
131 exports, but libgmp functions which might get called are imports.
133 libmp.la uses __GMP_DECLSPEC, just as if it were libgmp.la. libgmp and
134 libmp don't call each other, so there's no conflict or confusion.
136 Libtool DLL_EXPORT define is not used.
138 There's no attempt to support GMP built both static and DLL. Doing so
139 would mean applications would have to tell us which of the two is going
140 to be used when linking, and that seems very tedious and error prone if
141 using GMP by hand, and equally tedious from a package since autoconf and
142 automake don't give much help.
144 __GMP_DECLSPEC is required on all documented global functions and
145 variables, the various internals in gmp-impl.h etc can be left unadorned.
146 But internals used by the test programs or speed measuring programs
147 should have __GMP_DECLSPEC, and certainly constants or variables must
148 have it or the wrong address will be resolved.
150 In gcc __declspec can go at either the start or end of a prototype.
152 In Microsoft C __declspec must go at the start, or after the type like
153 void __declspec(...) *foo()". There's no __dllexport or anything to
154 guard against someone foolish #defining dllexport. _export used to be
155 available, but no longer.
157 In Borland C _export still exists, but needs to go after the type, like
158 "void _export foo();". Would have to change the __GMP_DECLSPEC syntax to
159 make use of that. Probably more trouble than it's worth. */
161 #if defined (__GNUC__)
162 #define __GMP_DECLSPEC_EXPORT __declspec(__dllexport__)
163 #define __GMP_DECLSPEC_IMPORT __declspec(__dllimport__)
165 #if defined (_MSC_VER) || defined (__BORLANDC__)
166 #define __GMP_DECLSPEC_EXPORT __declspec(dllexport)
167 #define __GMP_DECLSPEC_IMPORT __declspec(dllimport)
170 #define __GMP_DECLSPEC_EXPORT __export
171 #define __GMP_DECLSPEC_IMPORT __import
174 #define __GMP_DECLSPEC_EXPORT _Export
175 #define __GMP_DECLSPEC_IMPORT _Import
180 /* compiling to go into a DLL libgmp */
181 #define __GMP_DECLSPEC __GMP_DECLSPEC_EXPORT
183 /* compiling to go into an application which will link to a DLL libgmp */
184 #define __GMP_DECLSPEC __GMP_DECLSPEC_IMPORT
187 /* all other cases */
188 #define __GMP_DECLSPEC
192 #ifdef __GMP_SHORT_LIMB
193 typedef unsigned int mp_limb_t;
194 typedef int mp_limb_signed_t;
196 #ifdef _LONG_LONG_LIMB
197 typedef unsigned long long int mp_limb_t;
198 typedef long long int mp_limb_signed_t;
200 typedef unsigned long int mp_limb_t;
201 typedef long int mp_limb_signed_t;
205 /* For reference, note that the name __mpz_struct gets into C++ mangled
206 function names, which means although the "__" suggests an internal, we
207 must leave this name for binary compatibility. */
210 int _mp_alloc; /* Number of *limbs* allocated and pointed
211 to by the _mp_d field. */
212 int _mp_size; /* abs(_mp_size) is the number of limbs the
213 last field points to. If _mp_size is
214 negative this is a negative number. */
215 mp_limb_t *_mp_d; /* Pointer to the limbs. */
218 #endif /* __GNU_MP__ */
221 typedef __mpz_struct MP_INT; /* gmp 1 source compatibility */
222 typedef __mpz_struct mpz_t[1];
224 typedef mp_limb_t * mp_ptr;
225 typedef __gmp_const mp_limb_t * mp_srcptr;
226 #if defined (_CRAY) && ! defined (_CRAYMPP)
227 /* plain `int' is much faster (48 bits) */
228 #define __GMP_MP_SIZE_T_INT 1
229 typedef int mp_size_t;
230 typedef int mp_exp_t;
232 #define __GMP_MP_SIZE_T_INT 0
233 typedef long int mp_size_t;
234 typedef long int mp_exp_t;
239 __mpz_struct _mp_num;
240 __mpz_struct _mp_den;
243 typedef __mpq_struct MP_RAT; /* gmp 1 source compatibility */
244 typedef __mpq_struct mpq_t[1];
248 int _mp_prec; /* Max precision, in number of `mp_limb_t's.
249 Set by mpf_init and modified by
250 mpf_set_prec. The area pointed to by the
251 _mp_d field contains `prec' + 1 limbs. */
252 int _mp_size; /* abs(_mp_size) is the number of limbs the
253 last field points to. If _mp_size is
254 negative this is a negative number. */
255 mp_exp_t _mp_exp; /* Exponent, in the base of `mp_limb_t'. */
256 mp_limb_t *_mp_d; /* Pointer to the limbs. */
259 /* typedef __mpf_struct MP_FLOAT; */
260 typedef __mpf_struct mpf_t[1];
262 /* Available random number generation algorithms. */
265 GMP_RAND_ALG_DEFAULT = 0,
266 GMP_RAND_ALG_LC = GMP_RAND_ALG_DEFAULT /* Linear congruential. */
269 /* Random state struct. */
272 mpz_t _mp_seed; /* _mp_d member points to state of the generator. */
273 gmp_randalg_t _mp_alg; /* Currently unused. */
275 void *_mp_lc; /* Pointer to function pointers structure. */
277 } __gmp_randstate_struct;
278 typedef __gmp_randstate_struct gmp_randstate_t[1];
280 /* Types for function declarations in gmp files. */
281 /* ??? Should not pollute user name space with these ??? */
282 typedef __gmp_const __mpz_struct *mpz_srcptr;
283 typedef __mpz_struct *mpz_ptr;
284 typedef __gmp_const __mpf_struct *mpf_srcptr;
285 typedef __mpf_struct *mpf_ptr;
286 typedef __gmp_const __mpq_struct *mpq_srcptr;
287 typedef __mpq_struct *mpq_ptr;
290 /* This is not wanted in mp.h, so put it outside the __GNU_MP__ common
293 #if __GMP_WITHIN_GMPXX
294 /* compiling to go into a DLL libgmpxx */
295 #define __GMP_DECLSPEC_XX __GMP_DECLSPEC_EXPORT
297 /* compiling to go into a application which will link to a DLL libgmpxx */
298 #define __GMP_DECLSPEC_XX __GMP_DECLSPEC_IMPORT
301 /* all other cases */
302 #define __GMP_DECLSPEC_XX
306 #if __GMP_HAVE_PROTOTYPES
307 #define __GMP_PROTO(x) x
309 #define __GMP_PROTO(x) ()
313 #if __GMP_HAVE_TOKEN_PASTE
314 #define __MPN(x) __gmpn_##x
316 #define __MPN(x) __gmpn_/**/x
320 /* For reference, "defined(EOF)" cannot be used here. In g++ 2.95.4,
321 <iostream> defines EOF but not FILE. */
323 || defined (H_STDIO) \
324 || defined (_H_STDIO) /* AIX */ \
325 || defined (_STDIO_H) /* glibc, Sun, SCO */ \
326 || defined (_STDIO_H_) /* BSD, OSF */ \
327 || defined (__STDIO_H) /* Borland */ \
328 || defined (__STDIO_H__) /* IRIX */ \
329 || defined (_STDIO_INCLUDED) /* HPUX */ \
330 || defined (__dj_include_stdio_h_) /* DJGPP */ \
331 || defined (_FILE_DEFINED) /* Microsoft */ \
332 || defined (__STDIO__) /* Apple MPW MrC */ \
333 || defined (_MSL_STDIO_H) /* Metrowerks */ \
334 || defined (_STDIO_H_INCLUDED) /* QNX4 */ \
335 || defined (_ISO_STDIO_ISO_H) /* Sun C++ */
336 #define _GMP_H_HAVE_FILE 1
339 /* In ISO C, if a prototype involving "struct obstack *" is given without
340 that structure defined, then the struct is scoped down to just the
341 prototype, causing a conflict if it's subsequently defined for real. So
342 only give prototypes if we've got obstack.h. */
343 #if defined (_OBSTACK_H) /* glibc <obstack.h> */
344 #define _GMP_H_HAVE_OBSTACK 1
347 /* The prototypes for gmp_vprintf etc are provided only if va_list is
348 available, via an application having included <stdarg.h> or <varargs.h>.
349 Usually va_list is a typedef so can't be tested directly, but C99
350 specifies that va_start is a macro (and it was normally a macro on past
351 systems too), so look for that.
353 <stdio.h> will define some sort of va_list for vprintf and vfprintf, but
354 let's not bother trying to use that since it's not standard and since
355 application uses for gmp_vprintf etc will almost certainly require the
356 whole <stdarg.h> or <varargs.h> anyway. */
359 #define _GMP_H_HAVE_VA_LIST 1
362 /* Test for gcc >= maj.min, as per __GNUC_PREREQ in glibc */
363 #if defined (__GNUC__) && defined (__GNUC_MINOR__)
364 #define __GMP_GNUC_PREREQ(maj, min) \
365 ((__GNUC__ << 16) + __GNUC_MINOR__ >= ((maj) << 16) + (min))
367 #define __GMP_GNUC_PREREQ(maj, min) 0
370 /* "pure" is in gcc 2.96 and up, see "(gcc)Function Attributes". Basically
371 it means a function does nothing but examine its arguments and memory
372 (global or via arguments) to generate a return value, but changes nothing
373 and has no side-effects. __GMP_NO_ATTRIBUTE_CONST_PURE lets
374 tune/common.c etc turn this off when trying to write timing loops. */
375 #if __GMP_GNUC_PREREQ (2,96) && ! defined (__GMP_NO_ATTRIBUTE_CONST_PURE)
376 #define __GMP_ATTRIBUTE_PURE __attribute__ ((__pure__))
378 #define __GMP_ATTRIBUTE_PURE
382 /* __GMP_CAST allows us to use static_cast in C++, so our macros are clean
383 to "g++ -Wold-style-cast".
385 Casts in "extern inline" code within an extern "C" block don't induce
386 these warnings, so __GMP_CAST only needs to be used on documented
390 #define __GMP_CAST(type, expr) (static_cast<type> (expr))
392 #define __GMP_CAST(type, expr) ((type) (expr))
396 /* An empty "throw ()" means the function doesn't throw any C++ exceptions,
397 this can save some stack frame info in applications.
399 Currently it's given only on functions which never divide-by-zero etc,
400 don't allocate memory, and are expected to never need to allocate memory.
401 This leaves open the possibility of a C++ throw from a future GMP
404 mpz_set_ui etc are omitted to leave open the lazy allocation scheme
405 described in doc/tasks.html. mpz_get_d etc are omitted to leave open
406 exceptions for float overflows.
408 Note that __GMP_NOTHROW must be given on any inlines the same as on their
409 prototypes (for g++ at least, where they're used together). Note also
410 that g++ 3.0 demands that __GMP_NOTHROW is before other attributes like
411 __GMP_ATTRIBUTE_PURE. */
413 #if defined (__cplusplus)
414 #define __GMP_NOTHROW throw ()
416 #define __GMP_NOTHROW
420 /* PORTME: What other compilers have a useful "extern inline"? "static
421 inline" would be an acceptable substitute if the compiler (or linker)
422 discards unused statics. */
424 /* gcc has __inline__ in all modes, including strict ansi. Give a prototype
425 for an inline too, so as to correctly specify "dllimport" on windows, in
426 case the function is called rather than inlined.
427 GCC 4.3 and above with -std=c99 or -std=gnu99 implements ISO C99
428 inline semantics, unless -fgnu89-inline is used. */
430 #if (defined __GNUC_STDC_INLINE__) || (__GNUC__ == 4 && __GNUC_MINOR__ == 2)
431 #define __GMP_EXTERN_INLINE extern __inline__ __attribute__ ((__gnu_inline__))
433 #define __GMP_EXTERN_INLINE extern __inline__
435 #define __GMP_INLINE_PROTOTYPES 1
438 /* DEC C (eg. version 5.9) supports "static __inline foo()", even in -std1
439 strict ANSI mode. Inlining is done even when not optimizing (ie. -O0
440 mode, which is the default), but an unnecessary local copy of foo is
441 emitted unless -O is used. "extern __inline" is accepted, but the
442 "extern" appears to be ignored, ie. it becomes a plain global function
443 but which is inlined within its file. Don't know if all old versions of
444 DEC C supported __inline, but as a start let's do the right thing for
447 #define __GMP_EXTERN_INLINE static __inline
450 /* SCO OpenUNIX 8 cc supports "static inline foo()" but not in -Xc strict
451 ANSI mode (__STDC__ is 1 in that mode). Inlining only actually takes
452 place under -O. Without -O "foo" seems to be emitted whether it's used
453 or not, which is wasteful. "extern inline foo()" isn't useful, the
454 "extern" is apparently ignored, so foo is inlined if possible but also
455 emitted as a global, which causes multiple definition errors when
456 building a shared libgmp. */
457 #ifdef __SCO_VERSION__
458 #if __SCO_VERSION__ > 400000000 && __STDC__ != 1 \
459 && ! defined (__GMP_EXTERN_INLINE)
460 #define __GMP_EXTERN_INLINE static inline
464 /* Microsoft's C compiler accepts __inline */
466 #define __GMP_EXTERN_INLINE __inline
469 /* Recent enough Sun C compilers accept "extern inline" */
470 #if defined (__SUNPRO_C) && __SUNPRO_C >= 0x560 \
471 && ! defined (__GMP_EXTERN_INLINE)
472 #define __GMP_EXTERN_INLINE extern inline
475 /* Somewhat older Sun C compilers accept "static inline" */
476 #if defined (__SUNPRO_C) && __SUNPRO_C >= 0x540 \
477 && ! defined (__GMP_EXTERN_INLINE)
478 #define __GMP_EXTERN_INLINE static inline
482 /* C++ always has "inline" and since it's a normal feature the linker should
483 discard duplicate non-inlined copies, or if it doesn't then that's a
484 problem for everyone, not just GMP. */
485 #if defined (__cplusplus) && ! defined (__GMP_EXTERN_INLINE)
486 #define __GMP_EXTERN_INLINE inline
489 /* Don't do any inlining within a configure run, since if the compiler ends
490 up emitting copies of the code into the object file it can end up
491 demanding the various support routines (like mpn_popcount) for linking,
492 making the "alloca" test and perhaps others fail. And on hppa ia64 a
493 pre-release gcc 3.2 was seen not respecting the "extern" in "extern
494 __inline__", triggering this problem too. */
495 #if defined (__GMP_WITHIN_CONFIGURE) && ! __GMP_WITHIN_CONFIGURE_INLINE
496 #undef __GMP_EXTERN_INLINE
499 /* By default, don't give a prototype when there's going to be an inline
500 version. Note in particular that Cray C++ objects to the combination of
501 prototype and inline. */
502 #ifdef __GMP_EXTERN_INLINE
503 #ifndef __GMP_INLINE_PROTOTYPES
504 #define __GMP_INLINE_PROTOTYPES 0
507 #define __GMP_INLINE_PROTOTYPES 1
511 #define __GMP_ABS(x) ((x) >= 0 ? (x) : -(x))
512 #define __GMP_MAX(h,i) ((h) > (i) ? (h) : (i))
514 /* __GMP_USHRT_MAX is not "~ (unsigned short) 0" because short is promoted
516 #define __GMP_UINT_MAX (~ (unsigned) 0)
517 #define __GMP_ULONG_MAX (~ (unsigned long) 0)
518 #define __GMP_USHRT_MAX ((unsigned short) ~0)
521 /* __builtin_expect is in gcc 3.0, and not in 2.95. */
522 #if __GMP_GNUC_PREREQ (3,0)
523 #define __GMP_LIKELY(cond) __builtin_expect ((cond) != 0, 1)
524 #define __GMP_UNLIKELY(cond) __builtin_expect ((cond) != 0, 0)
526 #define __GMP_LIKELY(cond) (cond)
527 #define __GMP_UNLIKELY(cond) (cond)
531 #define __GMP_CRAY_Pragma(str) _Pragma (str)
533 #define __GMP_CRAY_Pragma(str)
537 /* Allow direct user access to numerator and denominator of a mpq_t object. */
538 #define mpq_numref(Q) (&((Q)->_mp_num))
539 #define mpq_denref(Q) (&((Q)->_mp_den))
542 #if defined (__cplusplus)
547 #define mp_set_memory_functions __gmp_set_memory_functions
548 __GMP_DECLSPEC void mp_set_memory_functions __GMP_PROTO ((void *(*) (size_t),
549 void *(*) (void *, size_t, size_t),
550 void (*) (void *, size_t))) __GMP_NOTHROW;
552 #define mp_get_memory_functions __gmp_get_memory_functions
553 __GMP_DECLSPEC void mp_get_memory_functions __GMP_PROTO ((void *(**) (size_t),
554 void *(**) (void *, size_t, size_t),
555 void (**) (void *, size_t))) __GMP_NOTHROW;
557 #define mp_bits_per_limb __gmp_bits_per_limb
558 __GMP_DECLSPEC extern __gmp_const int mp_bits_per_limb;
560 #define gmp_errno __gmp_errno
561 __GMP_DECLSPEC extern int gmp_errno;
563 #define gmp_version __gmp_version
564 __GMP_DECLSPEC extern __gmp_const char * __gmp_const gmp_version;
567 /**************** Random number routines. ****************/
570 #define gmp_randinit __gmp_randinit
571 __GMP_DECLSPEC void gmp_randinit __GMP_PROTO ((gmp_randstate_t, gmp_randalg_t, ...));
573 #define gmp_randinit_default __gmp_randinit_default
574 __GMP_DECLSPEC void gmp_randinit_default __GMP_PROTO ((gmp_randstate_t));
576 #define gmp_randinit_lc_2exp __gmp_randinit_lc_2exp
577 __GMP_DECLSPEC void gmp_randinit_lc_2exp __GMP_PROTO ((gmp_randstate_t,
578 mpz_srcptr, unsigned long int,
581 #define gmp_randinit_lc_2exp_size __gmp_randinit_lc_2exp_size
582 __GMP_DECLSPEC int gmp_randinit_lc_2exp_size __GMP_PROTO ((gmp_randstate_t, unsigned long));
584 #define gmp_randinit_mt __gmp_randinit_mt
585 __GMP_DECLSPEC void gmp_randinit_mt __GMP_PROTO ((gmp_randstate_t));
587 #define gmp_randinit_set __gmp_randinit_set
588 __GMP_DECLSPEC void gmp_randinit_set __GMP_PROTO ((gmp_randstate_t, __gmp_const __gmp_randstate_struct *));
590 #define gmp_randseed __gmp_randseed
591 __GMP_DECLSPEC void gmp_randseed __GMP_PROTO ((gmp_randstate_t, mpz_srcptr));
593 #define gmp_randseed_ui __gmp_randseed_ui
594 __GMP_DECLSPEC void gmp_randseed_ui __GMP_PROTO ((gmp_randstate_t, unsigned long int));
596 #define gmp_randclear __gmp_randclear
597 __GMP_DECLSPEC void gmp_randclear __GMP_PROTO ((gmp_randstate_t));
599 #define gmp_urandomb_ui __gmp_urandomb_ui
600 __GMP_DECLSPEC unsigned long gmp_urandomb_ui __GMP_PROTO ((gmp_randstate_t, unsigned long));
602 #define gmp_urandomm_ui __gmp_urandomm_ui
603 __GMP_DECLSPEC unsigned long gmp_urandomm_ui __GMP_PROTO ((gmp_randstate_t, unsigned long));
606 /**************** Formatted output routines. ****************/
608 #define gmp_asprintf __gmp_asprintf
609 __GMP_DECLSPEC int gmp_asprintf __GMP_PROTO ((char **, __gmp_const char *, ...));
611 #define gmp_fprintf __gmp_fprintf
612 #ifdef _GMP_H_HAVE_FILE
613 __GMP_DECLSPEC int gmp_fprintf __GMP_PROTO ((FILE *, __gmp_const char *, ...));
616 #define gmp_obstack_printf __gmp_obstack_printf
617 #if defined (_GMP_H_HAVE_OBSTACK)
618 __GMP_DECLSPEC int gmp_obstack_printf __GMP_PROTO ((struct obstack *, __gmp_const char *, ...));
621 #define gmp_obstack_vprintf __gmp_obstack_vprintf
622 #if defined (_GMP_H_HAVE_OBSTACK) && defined (_GMP_H_HAVE_VA_LIST)
623 __GMP_DECLSPEC int gmp_obstack_vprintf __GMP_PROTO ((struct obstack *, __gmp_const char *, va_list));
626 #define gmp_printf __gmp_printf
627 __GMP_DECLSPEC int gmp_printf __GMP_PROTO ((__gmp_const char *, ...));
629 #define gmp_snprintf __gmp_snprintf
630 __GMP_DECLSPEC int gmp_snprintf __GMP_PROTO ((char *, size_t, __gmp_const char *, ...));
632 #define gmp_sprintf __gmp_sprintf
633 __GMP_DECLSPEC int gmp_sprintf __GMP_PROTO ((char *, __gmp_const char *, ...));
635 #define gmp_vasprintf __gmp_vasprintf
636 #if defined (_GMP_H_HAVE_VA_LIST)
637 __GMP_DECLSPEC int gmp_vasprintf __GMP_PROTO ((char **, __gmp_const char *, va_list));
640 #define gmp_vfprintf __gmp_vfprintf
641 #if defined (_GMP_H_HAVE_FILE) && defined (_GMP_H_HAVE_VA_LIST)
642 __GMP_DECLSPEC int gmp_vfprintf __GMP_PROTO ((FILE *, __gmp_const char *, va_list));
645 #define gmp_vprintf __gmp_vprintf
646 #if defined (_GMP_H_HAVE_VA_LIST)
647 __GMP_DECLSPEC int gmp_vprintf __GMP_PROTO ((__gmp_const char *, va_list));
650 #define gmp_vsnprintf __gmp_vsnprintf
651 #if defined (_GMP_H_HAVE_VA_LIST)
652 __GMP_DECLSPEC int gmp_vsnprintf __GMP_PROTO ((char *, size_t, __gmp_const char *, va_list));
655 #define gmp_vsprintf __gmp_vsprintf
656 #if defined (_GMP_H_HAVE_VA_LIST)
657 __GMP_DECLSPEC int gmp_vsprintf __GMP_PROTO ((char *, __gmp_const char *, va_list));
661 /**************** Formatted input routines. ****************/
663 #define gmp_fscanf __gmp_fscanf
664 #ifdef _GMP_H_HAVE_FILE
665 __GMP_DECLSPEC int gmp_fscanf __GMP_PROTO ((FILE *, __gmp_const char *, ...));
668 #define gmp_scanf __gmp_scanf
669 __GMP_DECLSPEC int gmp_scanf __GMP_PROTO ((__gmp_const char *, ...));
671 #define gmp_sscanf __gmp_sscanf
672 __GMP_DECLSPEC int gmp_sscanf __GMP_PROTO ((__gmp_const char *, __gmp_const char *, ...));
674 #define gmp_vfscanf __gmp_vfscanf
675 #if defined (_GMP_H_HAVE_FILE) && defined (_GMP_H_HAVE_VA_LIST)
676 __GMP_DECLSPEC int gmp_vfscanf __GMP_PROTO ((FILE *, __gmp_const char *, va_list));
679 #define gmp_vscanf __gmp_vscanf
680 #if defined (_GMP_H_HAVE_VA_LIST)
681 __GMP_DECLSPEC int gmp_vscanf __GMP_PROTO ((__gmp_const char *, va_list));
684 #define gmp_vsscanf __gmp_vsscanf
685 #if defined (_GMP_H_HAVE_VA_LIST)
686 __GMP_DECLSPEC int gmp_vsscanf __GMP_PROTO ((__gmp_const char *, __gmp_const char *, va_list));
690 /**************** Integer (i.e. Z) routines. ****************/
692 #define _mpz_realloc __gmpz_realloc
693 #define mpz_realloc __gmpz_realloc
694 __GMP_DECLSPEC void *_mpz_realloc __GMP_PROTO ((mpz_ptr, mp_size_t));
696 #define mpz_abs __gmpz_abs
697 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_abs)
698 __GMP_DECLSPEC void mpz_abs __GMP_PROTO ((mpz_ptr, mpz_srcptr));
701 #define mpz_add __gmpz_add
702 __GMP_DECLSPEC void mpz_add __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
704 #define mpz_add_ui __gmpz_add_ui
705 __GMP_DECLSPEC void mpz_add_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
707 #define mpz_addmul __gmpz_addmul
708 __GMP_DECLSPEC void mpz_addmul __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
710 #define mpz_addmul_ui __gmpz_addmul_ui
711 __GMP_DECLSPEC void mpz_addmul_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
713 #define mpz_and __gmpz_and
714 __GMP_DECLSPEC void mpz_and __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
716 #define mpz_array_init __gmpz_array_init
717 __GMP_DECLSPEC void mpz_array_init __GMP_PROTO ((mpz_ptr, mp_size_t, mp_size_t));
719 #define mpz_bin_ui __gmpz_bin_ui
720 __GMP_DECLSPEC void mpz_bin_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
722 #define mpz_bin_uiui __gmpz_bin_uiui
723 __GMP_DECLSPEC void mpz_bin_uiui __GMP_PROTO ((mpz_ptr, unsigned long int, unsigned long int));
725 #define mpz_cdiv_q __gmpz_cdiv_q
726 __GMP_DECLSPEC void mpz_cdiv_q __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
728 #define mpz_cdiv_q_2exp __gmpz_cdiv_q_2exp
729 __GMP_DECLSPEC void mpz_cdiv_q_2exp __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long));
731 #define mpz_cdiv_q_ui __gmpz_cdiv_q_ui
732 __GMP_DECLSPEC unsigned long int mpz_cdiv_q_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
734 #define mpz_cdiv_qr __gmpz_cdiv_qr
735 __GMP_DECLSPEC void mpz_cdiv_qr __GMP_PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
737 #define mpz_cdiv_qr_ui __gmpz_cdiv_qr_ui
738 __GMP_DECLSPEC unsigned long int mpz_cdiv_qr_ui __GMP_PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
740 #define mpz_cdiv_r __gmpz_cdiv_r
741 __GMP_DECLSPEC void mpz_cdiv_r __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
743 #define mpz_cdiv_r_2exp __gmpz_cdiv_r_2exp
744 __GMP_DECLSPEC void mpz_cdiv_r_2exp __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long));
746 #define mpz_cdiv_r_ui __gmpz_cdiv_r_ui
747 __GMP_DECLSPEC unsigned long int mpz_cdiv_r_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
749 #define mpz_cdiv_ui __gmpz_cdiv_ui
750 __GMP_DECLSPEC unsigned long int mpz_cdiv_ui __GMP_PROTO ((mpz_srcptr, unsigned long int)) __GMP_ATTRIBUTE_PURE;
752 #define mpz_clear __gmpz_clear
753 __GMP_DECLSPEC void mpz_clear __GMP_PROTO ((mpz_ptr));
755 #define mpz_clrbit __gmpz_clrbit
756 __GMP_DECLSPEC void mpz_clrbit __GMP_PROTO ((mpz_ptr, unsigned long int));
758 #define mpz_cmp __gmpz_cmp
759 __GMP_DECLSPEC int mpz_cmp __GMP_PROTO ((mpz_srcptr, mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
761 #define mpz_cmp_d __gmpz_cmp_d
762 __GMP_DECLSPEC int mpz_cmp_d __GMP_PROTO ((mpz_srcptr, double)) __GMP_ATTRIBUTE_PURE;
764 #define _mpz_cmp_si __gmpz_cmp_si
765 __GMP_DECLSPEC int _mpz_cmp_si __GMP_PROTO ((mpz_srcptr, signed long int)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
767 #define _mpz_cmp_ui __gmpz_cmp_ui
768 __GMP_DECLSPEC int _mpz_cmp_ui __GMP_PROTO ((mpz_srcptr, unsigned long int)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
770 #define mpz_cmpabs __gmpz_cmpabs
771 __GMP_DECLSPEC int mpz_cmpabs __GMP_PROTO ((mpz_srcptr, mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
773 #define mpz_cmpabs_d __gmpz_cmpabs_d
774 __GMP_DECLSPEC int mpz_cmpabs_d __GMP_PROTO ((mpz_srcptr, double)) __GMP_ATTRIBUTE_PURE;
776 #define mpz_cmpabs_ui __gmpz_cmpabs_ui
777 __GMP_DECLSPEC int mpz_cmpabs_ui __GMP_PROTO ((mpz_srcptr, unsigned long int)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
779 #define mpz_com __gmpz_com
780 __GMP_DECLSPEC void mpz_com __GMP_PROTO ((mpz_ptr, mpz_srcptr));
782 #define mpz_combit __gmpz_combit
783 __GMP_DECLSPEC void mpz_combit __GMP_PROTO ((mpz_ptr, unsigned long int));
785 #define mpz_congruent_p __gmpz_congruent_p
786 __GMP_DECLSPEC int mpz_congruent_p __GMP_PROTO ((mpz_srcptr, mpz_srcptr, mpz_srcptr)) __GMP_ATTRIBUTE_PURE;
788 #define mpz_congruent_2exp_p __gmpz_congruent_2exp_p
789 __GMP_DECLSPEC int mpz_congruent_2exp_p __GMP_PROTO ((mpz_srcptr, mpz_srcptr, unsigned long)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
791 #define mpz_congruent_ui_p __gmpz_congruent_ui_p
792 __GMP_DECLSPEC int mpz_congruent_ui_p __GMP_PROTO ((mpz_srcptr, unsigned long, unsigned long)) __GMP_ATTRIBUTE_PURE;
794 #define mpz_divexact __gmpz_divexact
795 __GMP_DECLSPEC void mpz_divexact __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
797 #define mpz_divexact_ui __gmpz_divexact_ui
798 __GMP_DECLSPEC void mpz_divexact_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long));
800 #define mpz_divisible_p __gmpz_divisible_p
801 __GMP_DECLSPEC int mpz_divisible_p __GMP_PROTO ((mpz_srcptr, mpz_srcptr)) __GMP_ATTRIBUTE_PURE;
803 #define mpz_divisible_ui_p __gmpz_divisible_ui_p
804 __GMP_DECLSPEC int mpz_divisible_ui_p __GMP_PROTO ((mpz_srcptr, unsigned long)) __GMP_ATTRIBUTE_PURE;
806 #define mpz_divisible_2exp_p __gmpz_divisible_2exp_p
807 __GMP_DECLSPEC int mpz_divisible_2exp_p __GMP_PROTO ((mpz_srcptr, unsigned long)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
809 #define mpz_dump __gmpz_dump
810 __GMP_DECLSPEC void mpz_dump __GMP_PROTO ((mpz_srcptr));
812 #define mpz_export __gmpz_export
813 __GMP_DECLSPEC void *mpz_export __GMP_PROTO ((void *, size_t *, int, size_t, int, size_t, mpz_srcptr));
815 #define mpz_fac_ui __gmpz_fac_ui
816 __GMP_DECLSPEC void mpz_fac_ui __GMP_PROTO ((mpz_ptr, unsigned long int));
818 #define mpz_fdiv_q __gmpz_fdiv_q
819 __GMP_DECLSPEC void mpz_fdiv_q __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
821 #define mpz_fdiv_q_2exp __gmpz_fdiv_q_2exp
822 __GMP_DECLSPEC void mpz_fdiv_q_2exp __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
824 #define mpz_fdiv_q_ui __gmpz_fdiv_q_ui
825 __GMP_DECLSPEC unsigned long int mpz_fdiv_q_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
827 #define mpz_fdiv_qr __gmpz_fdiv_qr
828 __GMP_DECLSPEC void mpz_fdiv_qr __GMP_PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
830 #define mpz_fdiv_qr_ui __gmpz_fdiv_qr_ui
831 __GMP_DECLSPEC unsigned long int mpz_fdiv_qr_ui __GMP_PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
833 #define mpz_fdiv_r __gmpz_fdiv_r
834 __GMP_DECLSPEC void mpz_fdiv_r __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
836 #define mpz_fdiv_r_2exp __gmpz_fdiv_r_2exp
837 __GMP_DECLSPEC void mpz_fdiv_r_2exp __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
839 #define mpz_fdiv_r_ui __gmpz_fdiv_r_ui
840 __GMP_DECLSPEC unsigned long int mpz_fdiv_r_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
842 #define mpz_fdiv_ui __gmpz_fdiv_ui
843 __GMP_DECLSPEC unsigned long int mpz_fdiv_ui __GMP_PROTO ((mpz_srcptr, unsigned long int)) __GMP_ATTRIBUTE_PURE;
845 #define mpz_fib_ui __gmpz_fib_ui
846 __GMP_DECLSPEC void mpz_fib_ui __GMP_PROTO ((mpz_ptr, unsigned long int));
848 #define mpz_fib2_ui __gmpz_fib2_ui
849 __GMP_DECLSPEC void mpz_fib2_ui __GMP_PROTO ((mpz_ptr, mpz_ptr, unsigned long int));
851 #define mpz_fits_sint_p __gmpz_fits_sint_p
852 __GMP_DECLSPEC int mpz_fits_sint_p __GMP_PROTO ((mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
854 #define mpz_fits_slong_p __gmpz_fits_slong_p
855 __GMP_DECLSPEC int mpz_fits_slong_p __GMP_PROTO ((mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
857 #define mpz_fits_sshort_p __gmpz_fits_sshort_p
858 __GMP_DECLSPEC int mpz_fits_sshort_p __GMP_PROTO ((mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
860 #define mpz_fits_uint_p __gmpz_fits_uint_p
861 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_fits_uint_p)
862 __GMP_DECLSPEC int mpz_fits_uint_p __GMP_PROTO ((mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
865 #define mpz_fits_ulong_p __gmpz_fits_ulong_p
866 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_fits_ulong_p)
867 __GMP_DECLSPEC int mpz_fits_ulong_p __GMP_PROTO ((mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
870 #define mpz_fits_ushort_p __gmpz_fits_ushort_p
871 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_fits_ushort_p)
872 __GMP_DECLSPEC int mpz_fits_ushort_p __GMP_PROTO ((mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
875 #define mpz_gcd __gmpz_gcd
876 __GMP_DECLSPEC void mpz_gcd __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
878 #define mpz_gcd_ui __gmpz_gcd_ui
879 __GMP_DECLSPEC unsigned long int mpz_gcd_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
881 #define mpz_gcdext __gmpz_gcdext
882 __GMP_DECLSPEC void mpz_gcdext __GMP_PROTO ((mpz_ptr, mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
884 #define mpz_get_d __gmpz_get_d
885 __GMP_DECLSPEC double mpz_get_d __GMP_PROTO ((mpz_srcptr)) __GMP_ATTRIBUTE_PURE;
887 #define mpz_get_d_2exp __gmpz_get_d_2exp
888 __GMP_DECLSPEC double mpz_get_d_2exp __GMP_PROTO ((signed long int *, mpz_srcptr));
890 #define mpz_get_si __gmpz_get_si
891 __GMP_DECLSPEC /* signed */ long int mpz_get_si __GMP_PROTO ((mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
893 #define mpz_get_str __gmpz_get_str
894 __GMP_DECLSPEC char *mpz_get_str __GMP_PROTO ((char *, int, mpz_srcptr));
896 #define mpz_get_ui __gmpz_get_ui
897 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_get_ui)
898 __GMP_DECLSPEC unsigned long int mpz_get_ui __GMP_PROTO ((mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
901 #define mpz_getlimbn __gmpz_getlimbn
902 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_getlimbn)
903 __GMP_DECLSPEC mp_limb_t mpz_getlimbn __GMP_PROTO ((mpz_srcptr, mp_size_t)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
906 #define mpz_hamdist __gmpz_hamdist
907 __GMP_DECLSPEC unsigned long int mpz_hamdist __GMP_PROTO ((mpz_srcptr, mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
909 #define mpz_import __gmpz_import
910 __GMP_DECLSPEC void mpz_import __GMP_PROTO ((mpz_ptr, size_t, int, size_t, int, size_t, __gmp_const void *));
912 #define mpz_init __gmpz_init
913 __GMP_DECLSPEC void mpz_init __GMP_PROTO ((mpz_ptr));
915 #define mpz_init2 __gmpz_init2
916 __GMP_DECLSPEC void mpz_init2 __GMP_PROTO ((mpz_ptr, unsigned long));
918 #define mpz_init_set __gmpz_init_set
919 __GMP_DECLSPEC void mpz_init_set __GMP_PROTO ((mpz_ptr, mpz_srcptr));
921 #define mpz_init_set_d __gmpz_init_set_d
922 __GMP_DECLSPEC void mpz_init_set_d __GMP_PROTO ((mpz_ptr, double));
924 #define mpz_init_set_si __gmpz_init_set_si
925 __GMP_DECLSPEC void mpz_init_set_si __GMP_PROTO ((mpz_ptr, signed long int));
927 #define mpz_init_set_str __gmpz_init_set_str
928 __GMP_DECLSPEC int mpz_init_set_str __GMP_PROTO ((mpz_ptr, __gmp_const char *, int));
930 #define mpz_init_set_ui __gmpz_init_set_ui
931 __GMP_DECLSPEC void mpz_init_set_ui __GMP_PROTO ((mpz_ptr, unsigned long int));
933 #define mpz_inp_raw __gmpz_inp_raw
934 #ifdef _GMP_H_HAVE_FILE
935 __GMP_DECLSPEC size_t mpz_inp_raw __GMP_PROTO ((mpz_ptr, FILE *));
938 #define mpz_inp_str __gmpz_inp_str
939 #ifdef _GMP_H_HAVE_FILE
940 __GMP_DECLSPEC size_t mpz_inp_str __GMP_PROTO ((mpz_ptr, FILE *, int));
943 #define mpz_invert __gmpz_invert
944 __GMP_DECLSPEC int mpz_invert __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
946 #define mpz_ior __gmpz_ior
947 __GMP_DECLSPEC void mpz_ior __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
949 #define mpz_jacobi __gmpz_jacobi
950 __GMP_DECLSPEC int mpz_jacobi __GMP_PROTO ((mpz_srcptr, mpz_srcptr)) __GMP_ATTRIBUTE_PURE;
952 #define mpz_kronecker mpz_jacobi /* alias */
954 #define mpz_kronecker_si __gmpz_kronecker_si
955 __GMP_DECLSPEC int mpz_kronecker_si __GMP_PROTO ((mpz_srcptr, long)) __GMP_ATTRIBUTE_PURE;
957 #define mpz_kronecker_ui __gmpz_kronecker_ui
958 __GMP_DECLSPEC int mpz_kronecker_ui __GMP_PROTO ((mpz_srcptr, unsigned long)) __GMP_ATTRIBUTE_PURE;
960 #define mpz_si_kronecker __gmpz_si_kronecker
961 __GMP_DECLSPEC int mpz_si_kronecker __GMP_PROTO ((long, mpz_srcptr)) __GMP_ATTRIBUTE_PURE;
963 #define mpz_ui_kronecker __gmpz_ui_kronecker
964 __GMP_DECLSPEC int mpz_ui_kronecker __GMP_PROTO ((unsigned long, mpz_srcptr)) __GMP_ATTRIBUTE_PURE;
966 #define mpz_lcm __gmpz_lcm
967 __GMP_DECLSPEC void mpz_lcm __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
969 #define mpz_lcm_ui __gmpz_lcm_ui
970 __GMP_DECLSPEC void mpz_lcm_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long));
972 #define mpz_legendre mpz_jacobi /* alias */
974 #define mpz_lucnum_ui __gmpz_lucnum_ui
975 __GMP_DECLSPEC void mpz_lucnum_ui __GMP_PROTO ((mpz_ptr, unsigned long int));
977 #define mpz_lucnum2_ui __gmpz_lucnum2_ui
978 __GMP_DECLSPEC void mpz_lucnum2_ui __GMP_PROTO ((mpz_ptr, mpz_ptr, unsigned long int));
980 #define mpz_millerrabin __gmpz_millerrabin
981 __GMP_DECLSPEC int mpz_millerrabin __GMP_PROTO ((mpz_srcptr, int)) __GMP_ATTRIBUTE_PURE;
983 #define mpz_mod __gmpz_mod
984 __GMP_DECLSPEC void mpz_mod __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
986 #define mpz_mod_ui mpz_fdiv_r_ui /* same as fdiv_r because divisor unsigned */
988 #define mpz_mul __gmpz_mul
989 __GMP_DECLSPEC void mpz_mul __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
991 #define mpz_mul_2exp __gmpz_mul_2exp
992 __GMP_DECLSPEC void mpz_mul_2exp __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
994 #define mpz_mul_si __gmpz_mul_si
995 __GMP_DECLSPEC void mpz_mul_si __GMP_PROTO ((mpz_ptr, mpz_srcptr, long int));
997 #define mpz_mul_ui __gmpz_mul_ui
998 __GMP_DECLSPEC void mpz_mul_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
1000 #define mpz_neg __gmpz_neg
1001 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_neg)
1002 __GMP_DECLSPEC void mpz_neg __GMP_PROTO ((mpz_ptr, mpz_srcptr));
1005 #define mpz_nextprime __gmpz_nextprime
1006 __GMP_DECLSPEC void mpz_nextprime __GMP_PROTO ((mpz_ptr, mpz_srcptr));
1008 #define mpz_out_raw __gmpz_out_raw
1009 #ifdef _GMP_H_HAVE_FILE
1010 __GMP_DECLSPEC size_t mpz_out_raw __GMP_PROTO ((FILE *, mpz_srcptr));
1013 #define mpz_out_str __gmpz_out_str
1014 #ifdef _GMP_H_HAVE_FILE
1015 __GMP_DECLSPEC size_t mpz_out_str __GMP_PROTO ((FILE *, int, mpz_srcptr));
1018 #define mpz_perfect_power_p __gmpz_perfect_power_p
1019 __GMP_DECLSPEC int mpz_perfect_power_p __GMP_PROTO ((mpz_srcptr)) __GMP_ATTRIBUTE_PURE;
1021 #define mpz_perfect_square_p __gmpz_perfect_square_p
1022 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_perfect_square_p)
1023 __GMP_DECLSPEC int mpz_perfect_square_p __GMP_PROTO ((mpz_srcptr)) __GMP_ATTRIBUTE_PURE;
1026 #define mpz_popcount __gmpz_popcount
1027 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_popcount)
1028 __GMP_DECLSPEC unsigned long int mpz_popcount __GMP_PROTO ((mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1031 #define mpz_pow_ui __gmpz_pow_ui
1032 __GMP_DECLSPEC void mpz_pow_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
1034 #define mpz_powm __gmpz_powm
1035 __GMP_DECLSPEC void mpz_powm __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr, mpz_srcptr));
1037 #define mpz_powm_ui __gmpz_powm_ui
1038 __GMP_DECLSPEC void mpz_powm_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int, mpz_srcptr));
1040 #define mpz_probab_prime_p __gmpz_probab_prime_p
1041 __GMP_DECLSPEC int mpz_probab_prime_p __GMP_PROTO ((mpz_srcptr, int)) __GMP_ATTRIBUTE_PURE;
1043 #define mpz_random __gmpz_random
1044 __GMP_DECLSPEC void mpz_random __GMP_PROTO ((mpz_ptr, mp_size_t));
1046 #define mpz_random2 __gmpz_random2
1047 __GMP_DECLSPEC void mpz_random2 __GMP_PROTO ((mpz_ptr, mp_size_t));
1049 #define mpz_realloc2 __gmpz_realloc2
1050 __GMP_DECLSPEC void mpz_realloc2 __GMP_PROTO ((mpz_ptr, unsigned long));
1052 #define mpz_remove __gmpz_remove
1053 __GMP_DECLSPEC unsigned long int mpz_remove __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
1055 #define mpz_root __gmpz_root
1056 __GMP_DECLSPEC int mpz_root __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
1058 #define mpz_rootrem __gmpz_rootrem
1059 __GMP_DECLSPEC void mpz_rootrem __GMP_PROTO ((mpz_ptr,mpz_ptr, mpz_srcptr, unsigned long int));
1061 #define mpz_rrandomb __gmpz_rrandomb
1062 __GMP_DECLSPEC void mpz_rrandomb __GMP_PROTO ((mpz_ptr, gmp_randstate_t, unsigned long int));
1064 #define mpz_scan0 __gmpz_scan0
1065 __GMP_DECLSPEC unsigned long int mpz_scan0 __GMP_PROTO ((mpz_srcptr, unsigned long int)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1067 #define mpz_scan1 __gmpz_scan1
1068 __GMP_DECLSPEC unsigned long int mpz_scan1 __GMP_PROTO ((mpz_srcptr, unsigned long int)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1070 #define mpz_set __gmpz_set
1071 __GMP_DECLSPEC void mpz_set __GMP_PROTO ((mpz_ptr, mpz_srcptr));
1073 #define mpz_set_d __gmpz_set_d
1074 __GMP_DECLSPEC void mpz_set_d __GMP_PROTO ((mpz_ptr, double));
1076 #define mpz_set_f __gmpz_set_f
1077 __GMP_DECLSPEC void mpz_set_f __GMP_PROTO ((mpz_ptr, mpf_srcptr));
1079 #define mpz_set_q __gmpz_set_q
1080 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_set_q)
1081 __GMP_DECLSPEC void mpz_set_q __GMP_PROTO ((mpz_ptr, mpq_srcptr));
1084 #define mpz_set_si __gmpz_set_si
1085 __GMP_DECLSPEC void mpz_set_si __GMP_PROTO ((mpz_ptr, signed long int));
1087 #define mpz_set_str __gmpz_set_str
1088 __GMP_DECLSPEC int mpz_set_str __GMP_PROTO ((mpz_ptr, __gmp_const char *, int));
1090 #define mpz_set_ui __gmpz_set_ui
1091 __GMP_DECLSPEC void mpz_set_ui __GMP_PROTO ((mpz_ptr, unsigned long int));
1093 #define mpz_setbit __gmpz_setbit
1094 __GMP_DECLSPEC void mpz_setbit __GMP_PROTO ((mpz_ptr, unsigned long int));
1096 #define mpz_size __gmpz_size
1097 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpz_size)
1098 __GMP_DECLSPEC size_t mpz_size __GMP_PROTO ((mpz_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1101 #define mpz_sizeinbase __gmpz_sizeinbase
1102 __GMP_DECLSPEC size_t mpz_sizeinbase __GMP_PROTO ((mpz_srcptr, int)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1104 #define mpz_sqrt __gmpz_sqrt
1105 __GMP_DECLSPEC void mpz_sqrt __GMP_PROTO ((mpz_ptr, mpz_srcptr));
1107 #define mpz_sqrtrem __gmpz_sqrtrem
1108 __GMP_DECLSPEC void mpz_sqrtrem __GMP_PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr));
1110 #define mpz_sub __gmpz_sub
1111 __GMP_DECLSPEC void mpz_sub __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
1113 #define mpz_sub_ui __gmpz_sub_ui
1114 __GMP_DECLSPEC void mpz_sub_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
1116 #define mpz_ui_sub __gmpz_ui_sub
1117 __GMP_DECLSPEC void mpz_ui_sub __GMP_PROTO ((mpz_ptr, unsigned long int, mpz_srcptr));
1119 #define mpz_submul __gmpz_submul
1120 __GMP_DECLSPEC void mpz_submul __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
1122 #define mpz_submul_ui __gmpz_submul_ui
1123 __GMP_DECLSPEC void mpz_submul_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
1125 #define mpz_swap __gmpz_swap
1126 __GMP_DECLSPEC void mpz_swap __GMP_PROTO ((mpz_ptr, mpz_ptr)) __GMP_NOTHROW;
1128 #define mpz_tdiv_ui __gmpz_tdiv_ui
1129 __GMP_DECLSPEC unsigned long int mpz_tdiv_ui __GMP_PROTO ((mpz_srcptr, unsigned long int)) __GMP_ATTRIBUTE_PURE;
1131 #define mpz_tdiv_q __gmpz_tdiv_q
1132 __GMP_DECLSPEC void mpz_tdiv_q __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
1134 #define mpz_tdiv_q_2exp __gmpz_tdiv_q_2exp
1135 __GMP_DECLSPEC void mpz_tdiv_q_2exp __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
1137 #define mpz_tdiv_q_ui __gmpz_tdiv_q_ui
1138 __GMP_DECLSPEC unsigned long int mpz_tdiv_q_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
1140 #define mpz_tdiv_qr __gmpz_tdiv_qr
1141 __GMP_DECLSPEC void mpz_tdiv_qr __GMP_PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
1143 #define mpz_tdiv_qr_ui __gmpz_tdiv_qr_ui
1144 __GMP_DECLSPEC unsigned long int mpz_tdiv_qr_ui __GMP_PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
1146 #define mpz_tdiv_r __gmpz_tdiv_r
1147 __GMP_DECLSPEC void mpz_tdiv_r __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
1149 #define mpz_tdiv_r_2exp __gmpz_tdiv_r_2exp
1150 __GMP_DECLSPEC void mpz_tdiv_r_2exp __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
1152 #define mpz_tdiv_r_ui __gmpz_tdiv_r_ui
1153 __GMP_DECLSPEC unsigned long int mpz_tdiv_r_ui __GMP_PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
1155 #define mpz_tstbit __gmpz_tstbit
1156 __GMP_DECLSPEC int mpz_tstbit __GMP_PROTO ((mpz_srcptr, unsigned long int)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1158 #define mpz_ui_pow_ui __gmpz_ui_pow_ui
1159 __GMP_DECLSPEC void mpz_ui_pow_ui __GMP_PROTO ((mpz_ptr, unsigned long int, unsigned long int));
1161 #define mpz_urandomb __gmpz_urandomb
1162 __GMP_DECLSPEC void mpz_urandomb __GMP_PROTO ((mpz_ptr, gmp_randstate_t, unsigned long int));
1164 #define mpz_urandomm __gmpz_urandomm
1165 __GMP_DECLSPEC void mpz_urandomm __GMP_PROTO ((mpz_ptr, gmp_randstate_t, mpz_srcptr));
1167 #define mpz_xor __gmpz_xor
1168 #define mpz_eor __gmpz_xor
1169 __GMP_DECLSPEC void mpz_xor __GMP_PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
1172 /**************** Rational (i.e. Q) routines. ****************/
1174 #define mpq_abs __gmpq_abs
1175 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpq_abs)
1176 __GMP_DECLSPEC void mpq_abs __GMP_PROTO ((mpq_ptr, mpq_srcptr));
1179 #define mpq_add __gmpq_add
1180 __GMP_DECLSPEC void mpq_add __GMP_PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
1182 #define mpq_canonicalize __gmpq_canonicalize
1183 __GMP_DECLSPEC void mpq_canonicalize __GMP_PROTO ((mpq_ptr));
1185 #define mpq_clear __gmpq_clear
1186 __GMP_DECLSPEC void mpq_clear __GMP_PROTO ((mpq_ptr));
1188 #define mpq_cmp __gmpq_cmp
1189 __GMP_DECLSPEC int mpq_cmp __GMP_PROTO ((mpq_srcptr, mpq_srcptr)) __GMP_ATTRIBUTE_PURE;
1191 #define _mpq_cmp_si __gmpq_cmp_si
1192 __GMP_DECLSPEC int _mpq_cmp_si __GMP_PROTO ((mpq_srcptr, long, unsigned long)) __GMP_ATTRIBUTE_PURE;
1194 #define _mpq_cmp_ui __gmpq_cmp_ui
1195 __GMP_DECLSPEC int _mpq_cmp_ui __GMP_PROTO ((mpq_srcptr, unsigned long int, unsigned long int)) __GMP_ATTRIBUTE_PURE;
1197 #define mpq_div __gmpq_div
1198 __GMP_DECLSPEC void mpq_div __GMP_PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
1200 #define mpq_div_2exp __gmpq_div_2exp
1201 __GMP_DECLSPEC void mpq_div_2exp __GMP_PROTO ((mpq_ptr, mpq_srcptr, unsigned long));
1203 #define mpq_equal __gmpq_equal
1204 __GMP_DECLSPEC int mpq_equal __GMP_PROTO ((mpq_srcptr, mpq_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1206 #define mpq_get_num __gmpq_get_num
1207 __GMP_DECLSPEC void mpq_get_num __GMP_PROTO ((mpz_ptr, mpq_srcptr));
1209 #define mpq_get_den __gmpq_get_den
1210 __GMP_DECLSPEC void mpq_get_den __GMP_PROTO ((mpz_ptr, mpq_srcptr));
1212 #define mpq_get_d __gmpq_get_d
1213 __GMP_DECLSPEC double mpq_get_d __GMP_PROTO ((mpq_srcptr)) __GMP_ATTRIBUTE_PURE;
1215 #define mpq_get_str __gmpq_get_str
1216 __GMP_DECLSPEC char *mpq_get_str __GMP_PROTO ((char *, int, mpq_srcptr));
1218 #define mpq_init __gmpq_init
1219 __GMP_DECLSPEC void mpq_init __GMP_PROTO ((mpq_ptr));
1221 #define mpq_inp_str __gmpq_inp_str
1222 #ifdef _GMP_H_HAVE_FILE
1223 __GMP_DECLSPEC size_t mpq_inp_str __GMP_PROTO ((mpq_ptr, FILE *, int));
1226 #define mpq_inv __gmpq_inv
1227 __GMP_DECLSPEC void mpq_inv __GMP_PROTO ((mpq_ptr, mpq_srcptr));
1229 #define mpq_mul __gmpq_mul
1230 __GMP_DECLSPEC void mpq_mul __GMP_PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
1232 #define mpq_mul_2exp __gmpq_mul_2exp
1233 __GMP_DECLSPEC void mpq_mul_2exp __GMP_PROTO ((mpq_ptr, mpq_srcptr, unsigned long));
1235 #define mpq_neg __gmpq_neg
1236 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpq_neg)
1237 __GMP_DECLSPEC void mpq_neg __GMP_PROTO ((mpq_ptr, mpq_srcptr));
1240 #define mpq_out_str __gmpq_out_str
1241 #ifdef _GMP_H_HAVE_FILE
1242 __GMP_DECLSPEC size_t mpq_out_str __GMP_PROTO ((FILE *, int, mpq_srcptr));
1245 #define mpq_set __gmpq_set
1246 __GMP_DECLSPEC void mpq_set __GMP_PROTO ((mpq_ptr, mpq_srcptr));
1248 #define mpq_set_d __gmpq_set_d
1249 __GMP_DECLSPEC void mpq_set_d __GMP_PROTO ((mpq_ptr, double));
1251 #define mpq_set_den __gmpq_set_den
1252 __GMP_DECLSPEC void mpq_set_den __GMP_PROTO ((mpq_ptr, mpz_srcptr));
1254 #define mpq_set_f __gmpq_set_f
1255 __GMP_DECLSPEC void mpq_set_f __GMP_PROTO ((mpq_ptr, mpf_srcptr));
1257 #define mpq_set_num __gmpq_set_num
1258 __GMP_DECLSPEC void mpq_set_num __GMP_PROTO ((mpq_ptr, mpz_srcptr));
1260 #define mpq_set_si __gmpq_set_si
1261 __GMP_DECLSPEC void mpq_set_si __GMP_PROTO ((mpq_ptr, signed long int, unsigned long int));
1263 #define mpq_set_str __gmpq_set_str
1264 __GMP_DECLSPEC int mpq_set_str __GMP_PROTO ((mpq_ptr, __gmp_const char *, int));
1266 #define mpq_set_ui __gmpq_set_ui
1267 __GMP_DECLSPEC void mpq_set_ui __GMP_PROTO ((mpq_ptr, unsigned long int, unsigned long int));
1269 #define mpq_set_z __gmpq_set_z
1270 __GMP_DECLSPEC void mpq_set_z __GMP_PROTO ((mpq_ptr, mpz_srcptr));
1272 #define mpq_sub __gmpq_sub
1273 __GMP_DECLSPEC void mpq_sub __GMP_PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
1275 #define mpq_swap __gmpq_swap
1276 __GMP_DECLSPEC void mpq_swap __GMP_PROTO ((mpq_ptr, mpq_ptr)) __GMP_NOTHROW;
1279 /**************** Float (i.e. F) routines. ****************/
1281 #define mpf_abs __gmpf_abs
1282 __GMP_DECLSPEC void mpf_abs __GMP_PROTO ((mpf_ptr, mpf_srcptr));
1284 #define mpf_add __gmpf_add
1285 __GMP_DECLSPEC void mpf_add __GMP_PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
1287 #define mpf_add_ui __gmpf_add_ui
1288 __GMP_DECLSPEC void mpf_add_ui __GMP_PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
1289 #define mpf_ceil __gmpf_ceil
1290 __GMP_DECLSPEC void mpf_ceil __GMP_PROTO ((mpf_ptr, mpf_srcptr));
1292 #define mpf_clear __gmpf_clear
1293 __GMP_DECLSPEC void mpf_clear __GMP_PROTO ((mpf_ptr));
1295 #define mpf_cmp __gmpf_cmp
1296 __GMP_DECLSPEC int mpf_cmp __GMP_PROTO ((mpf_srcptr, mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1298 #define mpf_cmp_d __gmpf_cmp_d
1299 __GMP_DECLSPEC int mpf_cmp_d __GMP_PROTO ((mpf_srcptr, double)) __GMP_ATTRIBUTE_PURE;
1301 #define mpf_cmp_si __gmpf_cmp_si
1302 __GMP_DECLSPEC int mpf_cmp_si __GMP_PROTO ((mpf_srcptr, signed long int)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1304 #define mpf_cmp_ui __gmpf_cmp_ui
1305 __GMP_DECLSPEC int mpf_cmp_ui __GMP_PROTO ((mpf_srcptr, unsigned long int)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1307 #define mpf_div __gmpf_div
1308 __GMP_DECLSPEC void mpf_div __GMP_PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
1310 #define mpf_div_2exp __gmpf_div_2exp
1311 __GMP_DECLSPEC void mpf_div_2exp __GMP_PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
1313 #define mpf_div_ui __gmpf_div_ui
1314 __GMP_DECLSPEC void mpf_div_ui __GMP_PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
1316 #define mpf_dump __gmpf_dump
1317 __GMP_DECLSPEC void mpf_dump __GMP_PROTO ((mpf_srcptr));
1319 #define mpf_eq __gmpf_eq
1320 __GMP_DECLSPEC int mpf_eq __GMP_PROTO ((mpf_srcptr, mpf_srcptr, unsigned long int)) __GMP_ATTRIBUTE_PURE;
1322 #define mpf_fits_sint_p __gmpf_fits_sint_p
1323 __GMP_DECLSPEC int mpf_fits_sint_p __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1325 #define mpf_fits_slong_p __gmpf_fits_slong_p
1326 __GMP_DECLSPEC int mpf_fits_slong_p __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1328 #define mpf_fits_sshort_p __gmpf_fits_sshort_p
1329 __GMP_DECLSPEC int mpf_fits_sshort_p __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1331 #define mpf_fits_uint_p __gmpf_fits_uint_p
1332 __GMP_DECLSPEC int mpf_fits_uint_p __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1334 #define mpf_fits_ulong_p __gmpf_fits_ulong_p
1335 __GMP_DECLSPEC int mpf_fits_ulong_p __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1337 #define mpf_fits_ushort_p __gmpf_fits_ushort_p
1338 __GMP_DECLSPEC int mpf_fits_ushort_p __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1340 #define mpf_floor __gmpf_floor
1341 __GMP_DECLSPEC void mpf_floor __GMP_PROTO ((mpf_ptr, mpf_srcptr));
1343 #define mpf_get_d __gmpf_get_d
1344 __GMP_DECLSPEC double mpf_get_d __GMP_PROTO ((mpf_srcptr)) __GMP_ATTRIBUTE_PURE;
1346 #define mpf_get_d_2exp __gmpf_get_d_2exp
1347 __GMP_DECLSPEC double mpf_get_d_2exp __GMP_PROTO ((signed long int *, mpf_srcptr));
1349 #define mpf_get_default_prec __gmpf_get_default_prec
1350 __GMP_DECLSPEC unsigned long int mpf_get_default_prec __GMP_PROTO ((void)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1352 #define mpf_get_prec __gmpf_get_prec
1353 __GMP_DECLSPEC unsigned long int mpf_get_prec __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1355 #define mpf_get_si __gmpf_get_si
1356 __GMP_DECLSPEC long mpf_get_si __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1358 #define mpf_get_str __gmpf_get_str
1359 __GMP_DECLSPEC char *mpf_get_str __GMP_PROTO ((char *, mp_exp_t *, int, size_t, mpf_srcptr));
1361 #define mpf_get_ui __gmpf_get_ui
1362 __GMP_DECLSPEC unsigned long mpf_get_ui __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1364 #define mpf_init __gmpf_init
1365 __GMP_DECLSPEC void mpf_init __GMP_PROTO ((mpf_ptr));
1367 #define mpf_init2 __gmpf_init2
1368 __GMP_DECLSPEC void mpf_init2 __GMP_PROTO ((mpf_ptr, unsigned long int));
1370 #define mpf_init_set __gmpf_init_set
1371 __GMP_DECLSPEC void mpf_init_set __GMP_PROTO ((mpf_ptr, mpf_srcptr));
1373 #define mpf_init_set_d __gmpf_init_set_d
1374 __GMP_DECLSPEC void mpf_init_set_d __GMP_PROTO ((mpf_ptr, double));
1376 #define mpf_init_set_si __gmpf_init_set_si
1377 __GMP_DECLSPEC void mpf_init_set_si __GMP_PROTO ((mpf_ptr, signed long int));
1379 #define mpf_init_set_str __gmpf_init_set_str
1380 __GMP_DECLSPEC int mpf_init_set_str __GMP_PROTO ((mpf_ptr, __gmp_const char *, int));
1382 #define mpf_init_set_ui __gmpf_init_set_ui
1383 __GMP_DECLSPEC void mpf_init_set_ui __GMP_PROTO ((mpf_ptr, unsigned long int));
1385 #define mpf_inp_str __gmpf_inp_str
1386 #ifdef _GMP_H_HAVE_FILE
1387 __GMP_DECLSPEC size_t mpf_inp_str __GMP_PROTO ((mpf_ptr, FILE *, int));
1390 #define mpf_integer_p __gmpf_integer_p
1391 __GMP_DECLSPEC int mpf_integer_p __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1393 #define mpf_mul __gmpf_mul
1394 __GMP_DECLSPEC void mpf_mul __GMP_PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
1396 #define mpf_mul_2exp __gmpf_mul_2exp
1397 __GMP_DECLSPEC void mpf_mul_2exp __GMP_PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
1399 #define mpf_mul_ui __gmpf_mul_ui
1400 __GMP_DECLSPEC void mpf_mul_ui __GMP_PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
1402 #define mpf_neg __gmpf_neg
1403 __GMP_DECLSPEC void mpf_neg __GMP_PROTO ((mpf_ptr, mpf_srcptr));
1405 #define mpf_out_str __gmpf_out_str
1406 #ifdef _GMP_H_HAVE_FILE
1407 __GMP_DECLSPEC size_t mpf_out_str __GMP_PROTO ((FILE *, int, size_t, mpf_srcptr));
1410 #define mpf_pow_ui __gmpf_pow_ui
1411 __GMP_DECLSPEC void mpf_pow_ui __GMP_PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
1413 #define mpf_random2 __gmpf_random2
1414 __GMP_DECLSPEC void mpf_random2 __GMP_PROTO ((mpf_ptr, mp_size_t, mp_exp_t));
1416 #define mpf_reldiff __gmpf_reldiff
1417 __GMP_DECLSPEC void mpf_reldiff __GMP_PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
1419 #define mpf_set __gmpf_set
1420 __GMP_DECLSPEC void mpf_set __GMP_PROTO ((mpf_ptr, mpf_srcptr));
1422 #define mpf_set_d __gmpf_set_d
1423 __GMP_DECLSPEC void mpf_set_d __GMP_PROTO ((mpf_ptr, double));
1425 #define mpf_set_default_prec __gmpf_set_default_prec
1426 __GMP_DECLSPEC void mpf_set_default_prec __GMP_PROTO ((unsigned long int)) __GMP_NOTHROW;
1428 #define mpf_set_prec __gmpf_set_prec
1429 __GMP_DECLSPEC void mpf_set_prec __GMP_PROTO ((mpf_ptr, unsigned long int));
1431 #define mpf_set_prec_raw __gmpf_set_prec_raw
1432 __GMP_DECLSPEC void mpf_set_prec_raw __GMP_PROTO ((mpf_ptr, unsigned long int)) __GMP_NOTHROW;
1434 #define mpf_set_q __gmpf_set_q
1435 __GMP_DECLSPEC void mpf_set_q __GMP_PROTO ((mpf_ptr, mpq_srcptr));
1437 #define mpf_set_si __gmpf_set_si
1438 __GMP_DECLSPEC void mpf_set_si __GMP_PROTO ((mpf_ptr, signed long int));
1440 #define mpf_set_str __gmpf_set_str
1441 __GMP_DECLSPEC int mpf_set_str __GMP_PROTO ((mpf_ptr, __gmp_const char *, int));
1443 #define mpf_set_ui __gmpf_set_ui
1444 __GMP_DECLSPEC void mpf_set_ui __GMP_PROTO ((mpf_ptr, unsigned long int));
1446 #define mpf_set_z __gmpf_set_z
1447 __GMP_DECLSPEC void mpf_set_z __GMP_PROTO ((mpf_ptr, mpz_srcptr));
1449 #define mpf_size __gmpf_size
1450 __GMP_DECLSPEC size_t mpf_size __GMP_PROTO ((mpf_srcptr)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1452 #define mpf_sqrt __gmpf_sqrt
1453 __GMP_DECLSPEC void mpf_sqrt __GMP_PROTO ((mpf_ptr, mpf_srcptr));
1455 #define mpf_sqrt_ui __gmpf_sqrt_ui
1456 __GMP_DECLSPEC void mpf_sqrt_ui __GMP_PROTO ((mpf_ptr, unsigned long int));
1458 #define mpf_sub __gmpf_sub
1459 __GMP_DECLSPEC void mpf_sub __GMP_PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
1461 #define mpf_sub_ui __gmpf_sub_ui
1462 __GMP_DECLSPEC void mpf_sub_ui __GMP_PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
1464 #define mpf_swap __gmpf_swap
1465 __GMP_DECLSPEC void mpf_swap __GMP_PROTO ((mpf_ptr, mpf_ptr)) __GMP_NOTHROW;
1467 #define mpf_trunc __gmpf_trunc
1468 __GMP_DECLSPEC void mpf_trunc __GMP_PROTO ((mpf_ptr, mpf_srcptr));
1470 #define mpf_ui_div __gmpf_ui_div
1471 __GMP_DECLSPEC void mpf_ui_div __GMP_PROTO ((mpf_ptr, unsigned long int, mpf_srcptr));
1473 #define mpf_ui_sub __gmpf_ui_sub
1474 __GMP_DECLSPEC void mpf_ui_sub __GMP_PROTO ((mpf_ptr, unsigned long int, mpf_srcptr));
1476 #define mpf_urandomb __gmpf_urandomb
1477 __GMP_DECLSPEC void mpf_urandomb __GMP_PROTO ((mpf_t, gmp_randstate_t, unsigned long int));
1480 /************ Low level positive-integer (i.e. N) routines. ************/
1482 /* This is ugly, but we need to make user calls reach the prefixed function. */
1484 #define mpn_add __MPN(add)
1485 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpn_add)
1486 __GMP_DECLSPEC mp_limb_t mpn_add __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr,mp_size_t));
1489 #define mpn_add_1 __MPN(add_1)
1490 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpn_add_1)
1491 __GMP_DECLSPEC mp_limb_t mpn_add_1 __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t)) __GMP_NOTHROW;
1494 #define mpn_add_n __MPN(add_n)
1495 __GMP_DECLSPEC mp_limb_t mpn_add_n __GMP_PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t));
1497 #define mpn_addmul_1 __MPN(addmul_1)
1498 __GMP_DECLSPEC mp_limb_t mpn_addmul_1 __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
1500 #define mpn_bdivmod __MPN(bdivmod)
1501 __GMP_DECLSPEC mp_limb_t mpn_bdivmod __GMP_PROTO ((mp_ptr, mp_ptr, mp_size_t, mp_srcptr, mp_size_t, unsigned long int));
1503 #define mpn_cmp __MPN(cmp)
1504 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpn_cmp)
1505 __GMP_DECLSPEC int mpn_cmp __GMP_PROTO ((mp_srcptr, mp_srcptr, mp_size_t)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1508 #define mpn_divexact_by3(dst,src,size) \
1509 mpn_divexact_by3c (dst, src, size, __GMP_CAST (mp_limb_t, 0))
1511 #define mpn_divexact_by3c __MPN(divexact_by3c)
1512 __GMP_DECLSPEC mp_limb_t mpn_divexact_by3c __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
1514 #define mpn_divmod_1(qp,np,nsize,dlimb) \
1515 mpn_divrem_1 (qp, __GMP_CAST (mp_size_t, 0), np, nsize, dlimb)
1517 #define mpn_divrem __MPN(divrem)
1518 __GMP_DECLSPEC mp_limb_t mpn_divrem __GMP_PROTO ((mp_ptr, mp_size_t, mp_ptr, mp_size_t, mp_srcptr, mp_size_t));
1520 #define mpn_divrem_1 __MPN(divrem_1)
1521 __GMP_DECLSPEC mp_limb_t mpn_divrem_1 __GMP_PROTO ((mp_ptr, mp_size_t, mp_srcptr, mp_size_t, mp_limb_t));
1523 #define mpn_divrem_2 __MPN(divrem_2)
1524 __GMP_DECLSPEC mp_limb_t mpn_divrem_2 __GMP_PROTO ((mp_ptr, mp_size_t, mp_ptr, mp_size_t, mp_srcptr));
1526 #define mpn_gcd __MPN(gcd)
1527 __GMP_DECLSPEC mp_size_t mpn_gcd __GMP_PROTO ((mp_ptr, mp_ptr, mp_size_t, mp_ptr, mp_size_t));
1529 #define mpn_gcd_1 __MPN(gcd_1)
1530 __GMP_DECLSPEC mp_limb_t mpn_gcd_1 __GMP_PROTO ((mp_srcptr, mp_size_t, mp_limb_t)) __GMP_ATTRIBUTE_PURE;
1532 #define mpn_gcdext_1 __MPN(gcdext_1)
1533 __GMP_DECLSPEC mp_limb_t mpn_gcdext_1 __GMP_PROTO ((mp_ptr, mp_ptr, mp_limb_t, mp_limb_t));
1535 #define mpn_gcdext __MPN(gcdext)
1536 __GMP_DECLSPEC mp_size_t mpn_gcdext __GMP_PROTO ((mp_ptr, mp_ptr, mp_size_t *, mp_ptr, mp_size_t, mp_ptr, mp_size_t));
1538 #define mpn_get_str __MPN(get_str)
1539 __GMP_DECLSPEC size_t mpn_get_str __GMP_PROTO ((unsigned char *, int, mp_ptr, mp_size_t));
1541 #define mpn_hamdist __MPN(hamdist)
1542 __GMP_DECLSPEC unsigned long int mpn_hamdist __GMP_PROTO ((mp_srcptr, mp_srcptr, mp_size_t)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1544 #define mpn_lshift __MPN(lshift)
1545 __GMP_DECLSPEC mp_limb_t mpn_lshift __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, unsigned int));
1547 #define mpn_mod_1 __MPN(mod_1)
1548 __GMP_DECLSPEC mp_limb_t mpn_mod_1 __GMP_PROTO ((mp_srcptr, mp_size_t, mp_limb_t)) __GMP_ATTRIBUTE_PURE;
1550 #define mpn_mul __MPN(mul)
1551 __GMP_DECLSPEC mp_limb_t mpn_mul __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t));
1553 #define mpn_mul_1 __MPN(mul_1)
1554 __GMP_DECLSPEC mp_limb_t mpn_mul_1 __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
1556 #define mpn_mul_n __MPN(mul_n)
1557 __GMP_DECLSPEC void mpn_mul_n __GMP_PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t));
1559 #define mpn_neg_n __MPN(neg_n)
1560 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpn_neg_n)
1561 __GMP_DECLSPEC mp_limb_t mpn_neg_n __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t));
1564 #define mpn_perfect_square_p __MPN(perfect_square_p)
1565 __GMP_DECLSPEC int mpn_perfect_square_p __GMP_PROTO ((mp_srcptr, mp_size_t)) __GMP_ATTRIBUTE_PURE;
1567 #define mpn_popcount __MPN(popcount)
1568 __GMP_DECLSPEC unsigned long int mpn_popcount __GMP_PROTO ((mp_srcptr, mp_size_t)) __GMP_NOTHROW __GMP_ATTRIBUTE_PURE;
1570 #define mpn_pow_1 __MPN(pow_1)
1571 __GMP_DECLSPEC mp_size_t mpn_pow_1 __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t, mp_ptr));
1573 /* undocumented now, but retained here for upward compatibility */
1574 #define mpn_preinv_mod_1 __MPN(preinv_mod_1)
1575 __GMP_DECLSPEC mp_limb_t mpn_preinv_mod_1 __GMP_PROTO ((mp_srcptr, mp_size_t, mp_limb_t, mp_limb_t)) __GMP_ATTRIBUTE_PURE;
1577 #define mpn_random __MPN(random)
1578 __GMP_DECLSPEC void mpn_random __GMP_PROTO ((mp_ptr, mp_size_t));
1580 #define mpn_random2 __MPN(random2)
1581 __GMP_DECLSPEC void mpn_random2 __GMP_PROTO ((mp_ptr, mp_size_t));
1583 #define mpn_rshift __MPN(rshift)
1584 __GMP_DECLSPEC mp_limb_t mpn_rshift __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, unsigned int));
1586 #define mpn_scan0 __MPN(scan0)
1587 __GMP_DECLSPEC unsigned long int mpn_scan0 __GMP_PROTO ((mp_srcptr, unsigned long int)) __GMP_ATTRIBUTE_PURE;
1589 #define mpn_scan1 __MPN(scan1)
1590 __GMP_DECLSPEC unsigned long int mpn_scan1 __GMP_PROTO ((mp_srcptr, unsigned long int)) __GMP_ATTRIBUTE_PURE;
1592 #define mpn_set_str __MPN(set_str)
1593 __GMP_DECLSPEC mp_size_t mpn_set_str __GMP_PROTO ((mp_ptr, __gmp_const unsigned char *, size_t, int));
1595 #define mpn_sqrtrem __MPN(sqrtrem)
1596 __GMP_DECLSPEC mp_size_t mpn_sqrtrem __GMP_PROTO ((mp_ptr, mp_ptr, mp_srcptr, mp_size_t));
1598 #define mpn_sub __MPN(sub)
1599 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpn_sub)
1600 __GMP_DECLSPEC mp_limb_t mpn_sub __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr,mp_size_t));
1603 #define mpn_sub_1 __MPN(sub_1)
1604 #if __GMP_INLINE_PROTOTYPES || defined (__GMP_FORCE_mpn_sub_1)
1605 __GMP_DECLSPEC mp_limb_t mpn_sub_1 __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t)) __GMP_NOTHROW;
1608 #define mpn_sub_n __MPN(sub_n)
1609 __GMP_DECLSPEC mp_limb_t mpn_sub_n __GMP_PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t));
1611 #define mpn_submul_1 __MPN(submul_1)
1612 __GMP_DECLSPEC mp_limb_t mpn_submul_1 __GMP_PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
1614 #define mpn_tdiv_qr __MPN(tdiv_qr)
1615 __GMP_DECLSPEC void mpn_tdiv_qr __GMP_PROTO ((mp_ptr, mp_ptr, mp_size_t, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t));
1618 /**************** mpz inlines ****************/
1620 /* The following are provided as inlines where possible, but always exist as
1621 library functions too, for binary compatibility.
1623 Within gmp itself this inlining generally isn't relied on, since it
1624 doesn't get done for all compilers, whereas if something is worth
1625 inlining then it's worth arranging always.
1627 There are two styles of inlining here. When the same bit of code is
1628 wanted for the inline as for the library version, then __GMP_FORCE_foo
1629 arranges for that code to be emitted and the __GMP_EXTERN_INLINE
1630 directive suppressed, eg. mpz_fits_uint_p. When a different bit of code
1631 is wanted for the inline than for the library version, then
1632 __GMP_FORCE_foo arranges the inline to be suppressed, eg. mpz_abs. */
1634 #if defined (__GMP_EXTERN_INLINE) && ! defined (__GMP_FORCE_mpz_abs)
1635 __GMP_EXTERN_INLINE void
1636 mpz_abs (mpz_ptr __gmp_w, mpz_srcptr __gmp_u)
1638 if (__gmp_w != __gmp_u)
1639 mpz_set (__gmp_w, __gmp_u);
1640 __gmp_w->_mp_size = __GMP_ABS (__gmp_w->_mp_size);
1644 #if GMP_NAIL_BITS == 0
1645 #define __GMPZ_FITS_UTYPE_P(z,maxval) \
1646 mp_size_t __gmp_n = z->_mp_size; \
1647 mp_ptr __gmp_p = z->_mp_d; \
1648 return (__gmp_n == 0 || (__gmp_n == 1 && __gmp_p[0] <= maxval));
1650 #define __GMPZ_FITS_UTYPE_P(z,maxval) \
1651 mp_size_t __gmp_n = z->_mp_size; \
1652 mp_ptr __gmp_p = z->_mp_d; \
1653 return (__gmp_n == 0 || (__gmp_n == 1 && __gmp_p[0] <= maxval) \
1654 || (__gmp_n == 2 && __gmp_p[1] <= ((mp_limb_t) maxval >> GMP_NUMB_BITS)));
1657 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_fits_uint_p)
1658 #if ! defined (__GMP_FORCE_mpz_fits_uint_p)
1662 mpz_fits_uint_p (mpz_srcptr __gmp_z) __GMP_NOTHROW
1664 __GMPZ_FITS_UTYPE_P (__gmp_z, __GMP_UINT_MAX);
1668 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_fits_ulong_p)
1669 #if ! defined (__GMP_FORCE_mpz_fits_ulong_p)
1673 mpz_fits_ulong_p (mpz_srcptr __gmp_z) __GMP_NOTHROW
1675 __GMPZ_FITS_UTYPE_P (__gmp_z, __GMP_ULONG_MAX);
1679 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_fits_ushort_p)
1680 #if ! defined (__GMP_FORCE_mpz_fits_ushort_p)
1684 mpz_fits_ushort_p (mpz_srcptr __gmp_z) __GMP_NOTHROW
1686 __GMPZ_FITS_UTYPE_P (__gmp_z, __GMP_USHRT_MAX);
1690 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_get_ui)
1691 #if ! defined (__GMP_FORCE_mpz_get_ui)
1695 mpz_get_ui (mpz_srcptr __gmp_z) __GMP_NOTHROW
1697 mp_ptr __gmp_p = __gmp_z->_mp_d;
1698 mp_size_t __gmp_n = __gmp_z->_mp_size;
1699 mp_limb_t __gmp_l = __gmp_p[0];
1700 /* This is a "#if" rather than a plain "if" so as to avoid gcc warnings
1701 about "<< GMP_NUMB_BITS" exceeding the type size, and to avoid Borland
1702 C++ 6.0 warnings about condition always true for something like
1703 "__GMP_ULONG_MAX < GMP_NUMB_MASK". */
1704 #if GMP_NAIL_BITS == 0 || defined (_LONG_LONG_LIMB)
1705 /* limb==long and no nails, or limb==longlong, one limb is enough */
1706 return (__gmp_n != 0 ? __gmp_l : 0);
1708 /* limb==long and nails, need two limbs when available */
1709 __gmp_n = __GMP_ABS (__gmp_n);
1711 return (__gmp_n != 0 ? __gmp_l : 0);
1713 return __gmp_l + (__gmp_p[1] << GMP_NUMB_BITS);
1718 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_getlimbn)
1719 #if ! defined (__GMP_FORCE_mpz_getlimbn)
1723 mpz_getlimbn (mpz_srcptr __gmp_z, mp_size_t __gmp_n) __GMP_NOTHROW
1725 mp_limb_t __gmp_result = 0;
1726 if (__GMP_LIKELY (__gmp_n >= 0 && __gmp_n < __GMP_ABS (__gmp_z->_mp_size)))
1727 __gmp_result = __gmp_z->_mp_d[__gmp_n];
1728 return __gmp_result;
1732 #if defined (__GMP_EXTERN_INLINE) && ! defined (__GMP_FORCE_mpz_neg)
1733 __GMP_EXTERN_INLINE void
1734 mpz_neg (mpz_ptr __gmp_w, mpz_srcptr __gmp_u)
1736 if (__gmp_w != __gmp_u)
1737 mpz_set (__gmp_w, __gmp_u);
1738 __gmp_w->_mp_size = - __gmp_w->_mp_size;
1742 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_perfect_square_p)
1743 #if ! defined (__GMP_FORCE_mpz_perfect_square_p)
1747 mpz_perfect_square_p (mpz_srcptr __gmp_a)
1749 mp_size_t __gmp_asize;
1752 __gmp_asize = __gmp_a->_mp_size;
1753 __gmp_result = (__gmp_asize >= 0); /* zero is a square, negatives are not */
1754 if (__GMP_LIKELY (__gmp_asize > 0))
1755 __gmp_result = mpn_perfect_square_p (__gmp_a->_mp_d, __gmp_asize);
1756 return __gmp_result;
1760 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_popcount)
1761 #if ! defined (__GMP_FORCE_mpz_popcount)
1765 mpz_popcount (mpz_srcptr __gmp_u) __GMP_NOTHROW
1767 mp_size_t __gmp_usize;
1768 unsigned long __gmp_result;
1770 __gmp_usize = __gmp_u->_mp_size;
1771 __gmp_result = (__gmp_usize < 0 ? __GMP_ULONG_MAX : 0);
1772 if (__GMP_LIKELY (__gmp_usize > 0))
1773 __gmp_result = mpn_popcount (__gmp_u->_mp_d, __gmp_usize);
1774 return __gmp_result;
1778 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_set_q)
1779 #if ! defined (__GMP_FORCE_mpz_set_q)
1783 mpz_set_q (mpz_ptr __gmp_w, mpq_srcptr __gmp_u)
1785 mpz_tdiv_q (__gmp_w, mpq_numref (__gmp_u), mpq_denref (__gmp_u));
1789 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpz_size)
1790 #if ! defined (__GMP_FORCE_mpz_size)
1794 mpz_size (mpz_srcptr __gmp_z) __GMP_NOTHROW
1796 return __GMP_ABS (__gmp_z->_mp_size);
1801 /**************** mpq inlines ****************/
1803 #if defined (__GMP_EXTERN_INLINE) && ! defined (__GMP_FORCE_mpq_abs)
1804 __GMP_EXTERN_INLINE void
1805 mpq_abs (mpq_ptr __gmp_w, mpq_srcptr __gmp_u)
1807 if (__gmp_w != __gmp_u)
1808 mpq_set (__gmp_w, __gmp_u);
1809 __gmp_w->_mp_num._mp_size = __GMP_ABS (__gmp_w->_mp_num._mp_size);
1813 #if defined (__GMP_EXTERN_INLINE) && ! defined (__GMP_FORCE_mpq_neg)
1814 __GMP_EXTERN_INLINE void
1815 mpq_neg (mpq_ptr __gmp_w, mpq_srcptr __gmp_u)
1817 if (__gmp_w != __gmp_u)
1818 mpq_set (__gmp_w, __gmp_u);
1819 __gmp_w->_mp_num._mp_size = - __gmp_w->_mp_num._mp_size;
1824 /**************** mpn inlines ****************/
1826 /* The comments with __GMPN_ADD_1 below apply here too.
1828 The test for FUNCTION returning 0 should predict well. If it's assumed
1829 {yp,ysize} will usually have a random number of bits then the high limb
1830 won't be full and a carry out will occur a good deal less than 50% of the
1833 ysize==0 isn't a documented feature, but is used internally in a few
1836 Producing cout last stops it using up a register during the main part of
1837 the calculation, though gcc (as of 3.0) on an "if (mpn_add (...))"
1838 doesn't seem able to move the true and false legs of the conditional up
1839 to the two places cout is generated. */
1841 #define __GMPN_AORS(cout, wp, xp, xsize, yp, ysize, FUNCTION, TEST) \
1843 mp_size_t __gmp_i; \
1844 mp_limb_t __gmp_x; \
1846 /* ASSERT ((ysize) >= 0); */ \
1847 /* ASSERT ((xsize) >= (ysize)); */ \
1848 /* ASSERT (MPN_SAME_OR_SEPARATE2_P (wp, xsize, xp, xsize)); */ \
1849 /* ASSERT (MPN_SAME_OR_SEPARATE2_P (wp, xsize, yp, ysize)); */ \
1851 __gmp_i = (ysize); \
1854 if (FUNCTION (wp, xp, yp, __gmp_i)) \
1858 if (__gmp_i >= (xsize)) \
1863 __gmp_x = (xp)[__gmp_i]; \
1869 __GMPN_COPY_REST (wp, xp, xsize, __gmp_i); \
1875 #define __GMPN_ADD(cout, wp, xp, xsize, yp, ysize) \
1876 __GMPN_AORS (cout, wp, xp, xsize, yp, ysize, mpn_add_n, \
1877 (((wp)[__gmp_i++] = (__gmp_x + 1) & GMP_NUMB_MASK) == 0))
1878 #define __GMPN_SUB(cout, wp, xp, xsize, yp, ysize) \
1879 __GMPN_AORS (cout, wp, xp, xsize, yp, ysize, mpn_sub_n, \
1880 (((wp)[__gmp_i++] = (__gmp_x - 1) & GMP_NUMB_MASK), __gmp_x == 0))
1883 /* The use of __gmp_i indexing is designed to ensure a compile time src==dst
1884 remains nice and clear to the compiler, so that __GMPN_COPY_REST can
1885 disappear, and the load/add/store gets a chance to become a
1886 read-modify-write on CISC CPUs.
1890 Using a pair of pointers instead of indexing would be possible, but gcc
1891 isn't able to recognise compile-time src==dst in that case, even when the
1892 pointers are incremented more or less together. Other compilers would
1893 very likely have similar difficulty.
1895 gcc could use "if (__builtin_constant_p(src==dst) && src==dst)" or
1896 similar to detect a compile-time src==dst. This works nicely on gcc
1897 2.95.x, it's not good on gcc 3.0 where __builtin_constant_p(p==p) seems
1898 to be always false, for a pointer p. But the current code form seems
1899 good enough for src==dst anyway.
1901 gcc on x86 as usual doesn't give particularly good flags handling for the
1902 carry/borrow detection. It's tempting to want some multi instruction asm
1903 blocks to help it, and this was tried, but in truth there's only a few
1904 instructions to save and any gain is all too easily lost by register
1905 juggling setting up for the asm. */
1907 #if GMP_NAIL_BITS == 0
1908 #define __GMPN_AORS_1(cout, dst, src, n, v, OP, CB) \
1910 mp_size_t __gmp_i; \
1911 mp_limb_t __gmp_x, __gmp_r; \
1913 /* ASSERT ((n) >= 1); */ \
1914 /* ASSERT (MPN_SAME_OR_SEPARATE_P (dst, src, n)); */ \
1916 __gmp_x = (src)[0]; \
1917 __gmp_r = __gmp_x OP (v); \
1918 (dst)[0] = __gmp_r; \
1919 if (CB (__gmp_r, __gmp_x, (v))) \
1922 for (__gmp_i = 1; __gmp_i < (n);) \
1924 __gmp_x = (src)[__gmp_i]; \
1925 __gmp_r = __gmp_x OP 1; \
1926 (dst)[__gmp_i] = __gmp_r; \
1928 if (!CB (__gmp_r, __gmp_x, 1)) \
1930 if ((src) != (dst)) \
1931 __GMPN_COPY_REST (dst, src, n, __gmp_i); \
1939 if ((src) != (dst)) \
1940 __GMPN_COPY_REST (dst, src, n, 1); \
1946 #if GMP_NAIL_BITS >= 1
1947 #define __GMPN_AORS_1(cout, dst, src, n, v, OP, CB) \
1949 mp_size_t __gmp_i; \
1950 mp_limb_t __gmp_x, __gmp_r; \
1952 /* ASSERT ((n) >= 1); */ \
1953 /* ASSERT (MPN_SAME_OR_SEPARATE_P (dst, src, n)); */ \
1955 __gmp_x = (src)[0]; \
1956 __gmp_r = __gmp_x OP (v); \
1957 (dst)[0] = __gmp_r & GMP_NUMB_MASK; \
1958 if (__gmp_r >> GMP_NUMB_BITS != 0) \
1961 for (__gmp_i = 1; __gmp_i < (n);) \
1963 __gmp_x = (src)[__gmp_i]; \
1964 __gmp_r = __gmp_x OP 1; \
1965 (dst)[__gmp_i] = __gmp_r & GMP_NUMB_MASK; \
1967 if (__gmp_r >> GMP_NUMB_BITS == 0) \
1969 if ((src) != (dst)) \
1970 __GMPN_COPY_REST (dst, src, n, __gmp_i); \
1978 if ((src) != (dst)) \
1979 __GMPN_COPY_REST (dst, src, n, 1); \
1985 #define __GMPN_ADDCB(r,x,y) ((r) < (y))
1986 #define __GMPN_SUBCB(r,x,y) ((x) < (y))
1988 #define __GMPN_ADD_1(cout, dst, src, n, v) \
1989 __GMPN_AORS_1(cout, dst, src, n, v, +, __GMPN_ADDCB)
1990 #define __GMPN_SUB_1(cout, dst, src, n, v) \
1991 __GMPN_AORS_1(cout, dst, src, n, v, -, __GMPN_SUBCB)
1994 /* Compare {xp,size} and {yp,size}, setting "result" to positive, zero or
1995 negative. size==0 is allowed. On random data usually only one limb will
1996 need to be examined to get a result, so it's worth having it inline. */
1997 #define __GMPN_CMP(result, xp, yp, size) \
1999 mp_size_t __gmp_i; \
2000 mp_limb_t __gmp_x, __gmp_y; \
2002 /* ASSERT ((size) >= 0); */ \
2006 while (--__gmp_i >= 0) \
2008 __gmp_x = (xp)[__gmp_i]; \
2009 __gmp_y = (yp)[__gmp_i]; \
2010 if (__gmp_x != __gmp_y) \
2012 /* Cannot use __gmp_x - __gmp_y, may overflow an "int" */ \
2013 (result) = (__gmp_x > __gmp_y ? 1 : -1); \
2020 #if defined (__GMPN_COPY) && ! defined (__GMPN_COPY_REST)
2021 #define __GMPN_COPY_REST(dst, src, size, start) \
2023 /* ASSERT ((start) >= 0); */ \
2024 /* ASSERT ((start) <= (size)); */ \
2025 __GMPN_COPY ((dst)+(start), (src)+(start), (size)-(start)); \
2029 /* Copy {src,size} to {dst,size}, starting at "start". This is designed to
2030 keep the indexing dst[j] and src[j] nice and simple for __GMPN_ADD_1,
2032 #if ! defined (__GMPN_COPY_REST)
2033 #define __GMPN_COPY_REST(dst, src, size, start) \
2035 mp_size_t __gmp_j; \
2036 /* ASSERT ((size) >= 0); */ \
2037 /* ASSERT ((start) >= 0); */ \
2038 /* ASSERT ((start) <= (size)); */ \
2039 /* ASSERT (MPN_SAME_OR_SEPARATE_P (dst, src, size)); */ \
2040 __GMP_CRAY_Pragma ("_CRI ivdep"); \
2041 for (__gmp_j = (start); __gmp_j < (size); __gmp_j++) \
2042 (dst)[__gmp_j] = (src)[__gmp_j]; \
2046 /* Enhancement: Use some of the smarter code from gmp-impl.h. Maybe use
2047 mpn_copyi if there's a native version, and if we don't mind demanding
2048 binary compatibility for it (on targets which use it). */
2050 #if ! defined (__GMPN_COPY)
2051 #define __GMPN_COPY(dst, src, size) __GMPN_COPY_REST (dst, src, size, 0)
2055 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_add)
2056 #if ! defined (__GMP_FORCE_mpn_add)
2060 mpn_add (mp_ptr __gmp_wp, mp_srcptr __gmp_xp, mp_size_t __gmp_xsize, mp_srcptr __gmp_yp, mp_size_t __gmp_ysize)
2063 __GMPN_ADD (__gmp_c, __gmp_wp, __gmp_xp, __gmp_xsize, __gmp_yp, __gmp_ysize);
2068 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_add_1)
2069 #if ! defined (__GMP_FORCE_mpn_add_1)
2073 mpn_add_1 (mp_ptr __gmp_dst, mp_srcptr __gmp_src, mp_size_t __gmp_size, mp_limb_t __gmp_n) __GMP_NOTHROW
2076 __GMPN_ADD_1 (__gmp_c, __gmp_dst, __gmp_src, __gmp_size, __gmp_n);
2081 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_cmp)
2082 #if ! defined (__GMP_FORCE_mpn_cmp)
2086 mpn_cmp (mp_srcptr __gmp_xp, mp_srcptr __gmp_yp, mp_size_t __gmp_size) __GMP_NOTHROW
2089 __GMPN_CMP (__gmp_result, __gmp_xp, __gmp_yp, __gmp_size);
2090 return __gmp_result;
2094 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_sub)
2095 #if ! defined (__GMP_FORCE_mpn_sub)
2099 mpn_sub (mp_ptr __gmp_wp, mp_srcptr __gmp_xp, mp_size_t __gmp_xsize, mp_srcptr __gmp_yp, mp_size_t __gmp_ysize)
2102 __GMPN_SUB (__gmp_c, __gmp_wp, __gmp_xp, __gmp_xsize, __gmp_yp, __gmp_ysize);
2107 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_sub_1)
2108 #if ! defined (__GMP_FORCE_mpn_sub_1)
2112 mpn_sub_1 (mp_ptr __gmp_dst, mp_srcptr __gmp_src, mp_size_t __gmp_size, mp_limb_t __gmp_n) __GMP_NOTHROW
2115 __GMPN_SUB_1 (__gmp_c, __gmp_dst, __gmp_src, __gmp_size, __gmp_n);
2120 #if defined (__GMP_EXTERN_INLINE) || defined (__GMP_FORCE_mpn_neg_n)
2121 #if ! defined (__GMP_FORCE_mpn_neg_n)
2125 mpn_neg_n (mp_ptr __gmp_rp, mp_srcptr __gmp_up, mp_size_t __gmp_n)
2127 mp_limb_t __gmp_ul, __gmp_cy;
2130 __gmp_ul = *__gmp_up++;
2131 *__gmp_rp++ = -__gmp_ul - __gmp_cy;
2132 __gmp_cy |= __gmp_ul != 0;
2133 } while (--__gmp_n != 0);
2138 #if defined (__cplusplus)
2143 /* Allow faster testing for negative, zero, and positive. */
2144 #define mpz_sgn(Z) ((Z)->_mp_size < 0 ? -1 : (Z)->_mp_size > 0)
2145 #define mpf_sgn(F) ((F)->_mp_size < 0 ? -1 : (F)->_mp_size > 0)
2146 #define mpq_sgn(Q) ((Q)->_mp_num._mp_size < 0 ? -1 : (Q)->_mp_num._mp_size > 0)
2148 /* When using GCC, optimize certain common comparisons. */
2149 #if defined (__GNUC__) && __GNUC__ >= 2
2150 #define mpz_cmp_ui(Z,UI) \
2151 (__builtin_constant_p (UI) && (UI) == 0 \
2152 ? mpz_sgn (Z) : _mpz_cmp_ui (Z,UI))
2153 #define mpz_cmp_si(Z,SI) \
2154 (__builtin_constant_p (SI) && (SI) == 0 ? mpz_sgn (Z) \
2155 : __builtin_constant_p (SI) && (SI) > 0 \
2156 ? _mpz_cmp_ui (Z, __GMP_CAST (unsigned long int, SI)) \
2157 : _mpz_cmp_si (Z,SI))
2158 #define mpq_cmp_ui(Q,NUI,DUI) \
2159 (__builtin_constant_p (NUI) && (NUI) == 0 \
2160 ? mpq_sgn (Q) : _mpq_cmp_ui (Q,NUI,DUI))
2161 #define mpq_cmp_si(q,n,d) \
2162 (__builtin_constant_p ((n) >= 0) && (n) >= 0 \
2163 ? mpq_cmp_ui (q, __GMP_CAST (unsigned long, n), d) \
2164 : _mpq_cmp_si (q, n, d))
2166 #define mpz_cmp_ui(Z,UI) _mpz_cmp_ui (Z,UI)
2167 #define mpz_cmp_si(Z,UI) _mpz_cmp_si (Z,UI)
2168 #define mpq_cmp_ui(Q,NUI,DUI) _mpq_cmp_ui (Q,NUI,DUI)
2169 #define mpq_cmp_si(q,n,d) _mpq_cmp_si(q,n,d)
2173 /* Using "&" rather than "&&" means these can come out branch-free. Every
2174 mpz_t has at least one limb allocated, so fetching the low limb is always
2176 #define mpz_odd_p(z) (((z)->_mp_size != 0) & __GMP_CAST (int, (z)->_mp_d[0]))
2177 #define mpz_even_p(z) (! mpz_odd_p (z))
2180 /**************** C++ routines ****************/
2183 __GMP_DECLSPEC_XX std::ostream& operator<< (std::ostream &, mpz_srcptr);
2184 __GMP_DECLSPEC_XX std::ostream& operator<< (std::ostream &, mpq_srcptr);
2185 __GMP_DECLSPEC_XX std::ostream& operator<< (std::ostream &, mpf_srcptr);
2186 __GMP_DECLSPEC_XX std::istream& operator>> (std::istream &, mpz_ptr);
2187 __GMP_DECLSPEC_XX std::istream& operator>> (std::istream &, mpq_ptr);
2188 __GMP_DECLSPEC_XX std::istream& operator>> (std::istream &, mpf_ptr);
2192 /* Source-level compatibility with GMP 2 and earlier. */
2193 #define mpn_divmod(qp,np,nsize,dp,dsize) \
2194 mpn_divrem (qp, __GMP_CAST (mp_size_t, 0), np, nsize, dp, dsize)
2196 /* Source-level compatibility with GMP 1. */
2197 #define mpz_mdiv mpz_fdiv_q
2198 #define mpz_mdivmod mpz_fdiv_qr
2199 #define mpz_mmod mpz_fdiv_r
2200 #define mpz_mdiv_ui mpz_fdiv_q_ui
2201 #define mpz_mdivmod_ui(q,r,n,d) \
2202 (((r) == 0) ? mpz_fdiv_q_ui (q,n,d) : mpz_fdiv_qr_ui (q,r,n,d))
2203 #define mpz_mmod_ui(r,n,d) \
2204 (((r) == 0) ? mpz_fdiv_ui (n,d) : mpz_fdiv_r_ui (r,n,d))
2206 /* Useful synonyms, but not quite compatible with GMP 1. */
2207 #define mpz_div mpz_fdiv_q
2208 #define mpz_divmod mpz_fdiv_qr
2209 #define mpz_div_ui mpz_fdiv_q_ui
2210 #define mpz_divmod_ui mpz_fdiv_qr_ui
2211 #define mpz_div_2exp mpz_fdiv_q_2exp
2212 #define mpz_mod_2exp mpz_fdiv_r_2exp
2217 GMP_ERROR_UNSUPPORTED_ARGUMENT = 1,
2218 GMP_ERROR_DIVISION_BY_ZERO = 2,
2219 GMP_ERROR_SQRT_OF_NEGATIVE = 4,
2220 GMP_ERROR_INVALID_ARGUMENT = 8
2223 /* Define CC and CFLAGS which were used to build this version of GMP */
2224 #define __GMP_CC "gcc"
2225 #define __GMP_CFLAGS "-m32 -O2 -pedantic -fomit-frame-pointer -mtune=k8 -march=k8"
2227 /* Major version number is the value of __GNU_MP__ too, above and in mp.h. */
2228 #define __GNU_MP_VERSION 4
2229 #define __GNU_MP_VERSION_MINOR 3
2230 #define __GNU_MP_VERSION_PATCHLEVEL 1
2233 #endif /* __GMP_H__ */