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