Home | History | Annotate | Line # | Download | only in sh3el
      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                      32
     44 #define GMP_NAIL_BITS                      0
     45 #endif
     46 #define GMP_NUMB_BITS     (GMP_LIMB_BITS - GMP_NAIL_BITS)
     47 #define GMP_NUMB_MASK     ((~ __GMP_CAST (mp_limb_t, 0)) >> GMP_NAIL_BITS)
     48 #define GMP_NUMB_MAX      GMP_NUMB_MASK
     49 #define GMP_NAIL_MASK     (~ GMP_NUMB_MASK)
     50 
     51 
     52 #ifndef __GNU_MP__
     53 #define __GNU_MP__ 6
     54 
     55 #include <stddef.h>    /* for size_t */
     56 #include <limits.h>
     57 
     58 /* Instantiated by configure. */
     59 #if ! defined (__GMP_WITHIN_CONFIGURE)
     60 /* #undef _LONG_LONG_LIMB */
     61 #define __GMP_LIBGMP_DLL  0
     62 #endif
     63 
     64 
     65 /* __GMP_DECLSPEC supports Windows DLL versions of libgmp, and is empty in
     66    all other circumstances.
     67 
     68    When compiling objects for libgmp, __GMP_DECLSPEC is an export directive,
     69    or when compiling for an application it's an import directive.  The two
     70    cases are differentiated by __GMP_WITHIN_GMP defined by the GMP Makefiles
     71    (and not defined from an application).
     72 
     73    __GMP_DECLSPEC_XX is similarly used for libgmpxx.  __GMP_WITHIN_GMPXX
     74    indicates when building libgmpxx, and in that case libgmpxx functions are
     75    exports, but libgmp functions which might get called are imports.
     76 
     77    Libtool DLL_EXPORT define is not used.
     78 
     79    There's no attempt to support GMP built both static and DLL.  Doing so
     80    would mean applications would have to tell us which of the two is going
     81    to be used when linking, and that seems very tedious and error prone if
     82    using GMP by hand, and equally tedious from a package since autoconf and
     83    automake don't give much help.
     84 
     85    __GMP_DECLSPEC is required on all documented global functions and
     86    variables, the various internals in gmp-impl.h etc can be left unadorned.
     87    But internals used by the test programs or speed measuring programs
     88    should have __GMP_DECLSPEC, and certainly constants or variables must
     89    have it or the wrong address will be resolved.
     90 
     91    In gcc __declspec can go at either the start or end of a prototype.
     92 
     93    In Microsoft C __declspec must go at the start, or after the type like
     94    void __declspec(...) *foo()".  There's no __dllexport or anything to
     95    guard against someone foolish #defining dllexport.  _export used to be
     96    available, but no longer.
     97 
     98    In Borland C _export still exists, but needs to go after the type, like
     99    "void _export foo();".  Would have to change the __GMP_DECLSPEC syntax to
    100    make use of that.  Probably more trouble than it's worth.  */
    101 
    102 #if defined (__GNUC__)
    103 #define __GMP_DECLSPEC_EXPORT  __declspec(__dllexport__)
    104 #define __GMP_DECLSPEC_IMPORT  __declspec(__dllimport__)
    105 #endif
    106 #if defined (_MSC_VER) || defined (__BORLANDC__)
    107 #define __GMP_DECLSPEC_EXPORT  __declspec(dllexport)
    108 #define __GMP_DECLSPEC_IMPORT  __declspec(dllimport)
    109 #endif
    110 #ifdef __WATCOMC__
    111 #define __GMP_DECLSPEC_EXPORT  __export
    112 #define __GMP_DECLSPEC_IMPORT  __import
    113 #endif
    114 #ifdef __IBMC__
    115 #define __GMP_DECLSPEC_EXPORT  _Export
    116 #define __GMP_DECLSPEC_IMPORT  _Import
    117 #endif
    118 
    119 #if __GMP_LIBGMP_DLL
    120 #ifdef __GMP_WITHIN_GMP
    121 /* compiling to go into a DLL libgmp */
    122 #define __GMP_DECLSPEC  __GMP_DECLSPEC_EXPORT
    123 #else
    124 /* compiling to go into an application which will link to a DLL libgmp */
    125 #define __GMP_DECLSPEC  __GMP_DECLSPEC_IMPORT
    126 #endif
    127 #else
    128 /* all other cases */
    129 #define __GMP_DECLSPEC
    130 #endif
    131 
    132 
    133 #ifdef __GMP_SHORT_LIMB
    134 typedef unsigned int		mp_limb_t;
    135 typedef int			mp_limb_signed_t;
    136 #else
    137 #ifdef _LONG_LONG_LIMB
    138 typedef unsigned long long int	mp_limb_t;
    139 typedef long long int		mp_limb_signed_t;
    140 #else
    141 typedef unsigned long int	mp_limb_t;
    142 typedef long int		mp_limb_signed_t;
    143 #endif
    144 #endif
    145 typedef unsigned long int	mp_bitcnt_t;
    146 
    147 /* For reference, note that the name __mpz_struct gets into C++ mangled
    148    function names, which means although the "__" suggests an internal, we
    149    must leave this name for binary compatibility.  */
    150 typedef struct
    151 {
    152   int _mp_alloc;		/* Number of *limbs* allocated and pointed
    153 				   to by the _mp_d field.  */
    154   int _mp_size;			/* abs(_mp_size) is the number of limbs the
    155 				   last field points to.  If _mp_size is
    156 				   negative this is a negative number.  */
    157   mp_limb_t *_mp_d;		/* Pointer to the limbs.  */
    158 } __mpz_struct;
    159 
    160 #endif /* __GNU_MP__ */
    161 
    162 
    163 typedef __mpz_struct MP_INT;    /* gmp 1 source compatibility */
    164 typedef __mpz_struct mpz_t[1];
    165 
    166 typedef mp_limb_t *		mp_ptr;
    167 typedef const mp_limb_t *	mp_srcptr;
    168 #if defined (_CRAY) && ! defined (_CRAYMPP)
    169 /* plain `int' is much faster (48 bits) */
    170 #define __GMP_MP_SIZE_T_INT     1
    171 typedef int			mp_size_t;
    172 typedef int			mp_exp_t;
    173 #else
    174 #define __GMP_MP_SIZE_T_INT     0
    175 typedef long int		mp_size_t;
    176 typedef long int		mp_exp_t;
    177 #endif
    178 
    179 typedef struct
    180 {
    181   __mpz_struct _mp_num;
    182   __mpz_struct _mp_den;
    183 } __mpq_struct;
    184 
    185 typedef __mpq_struct MP_RAT;    /* gmp 1 source compatibility */
    186 typedef __mpq_struct mpq_t[1];
    187 
    188 typedef struct
    189 {
    190   int _mp_prec;			/* Max precision, in number of `mp_limb_t's.
    191 				   Set by mpf_init and modified by
    192 				   mpf_set_prec.  The area pointed to by the
    193 				   _mp_d field contains `prec' + 1 limbs.  */
    194   int _mp_size;			/* abs(_mp_size) is the number of limbs the
    195 				   last field points to.  If _mp_size is
    196 				   negative this is a negative number.  */
    197   mp_exp_t _mp_exp;		/* Exponent, in the base of `mp_limb_t'.  */
    198   mp_limb_t *_mp_d;		/* Pointer to the limbs.  */
    199 } __mpf_struct;
    200 
    201 /* typedef __mpf_struct MP_FLOAT; */
    202 typedef __mpf_struct mpf_t[1];
    203 
    204 /* Available random number generation algorithms.  */
    205 typedef enum
    206 {
    207   GMP_RAND_ALG_DEFAULT = 0,
    208   GMP_RAND_ALG_LC = GMP_RAND_ALG_DEFAULT /* Linear congruential.  */
    209 } gmp_randalg_t;
    210 
    211 /* Random state struct.  */
    212 typedef struct
    213 {
    214   mpz_t _mp_seed;	  /* _mp_d member points to state of the generator. */
    215   gmp_randalg_t _mp_alg;  /* Currently unused. */
    216   union {
    217     void *_mp_lc;         /* Pointer to function pointers structure.  */
    218   } _mp_algdata;
    219 } __gmp_randstate_struct;
    220 typedef __gmp_randstate_struct gmp_randstate_t[1];
    221 
    222 /* Types for function declarations in gmp files.  */
    223 /* ??? Should not pollute user name space with these ??? */
    224 typedef const __mpz_struct *mpz_srcptr;
    225 typedef __mpz_struct *mpz_ptr;
    226 typedef const __mpf_struct *mpf_srcptr;
    227 typedef __mpf_struct *mpf_ptr;
    228 typedef const __mpq_struct *mpq_srcptr;
    229 typedef __mpq_struct *mpq_ptr;
    230 
    231 
    232 #if __GMP_LIBGMP_DLL
    233 #ifdef __GMP_WITHIN_GMPXX
    234 /* compiling to go into a DLL libgmpxx */
    235 #define __GMP_DECLSPEC_XX  __GMP_DECLSPEC_EXPORT
    236 #else
    237 /* compiling to go into a application which will link to a DLL libgmpxx */
    238 #define __GMP_DECLSPEC_XX  __GMP_DECLSPEC_IMPORT
    239 #endif
    240 #else
    241 /* all other cases */
    242 #define __GMP_DECLSPEC_XX
    243 #endif
    244 
    245 
    246 #ifndef __MPN
    247 #define __MPN(x) __gmpn_##x
    248 #endif
    249 
    250 /* For reference, "defined(EOF)" cannot be used here.  In g++ 2.95.4,
    251    <iostream> defines EOF but not FILE.  */
    252 #if defined (FILE)                                              \
    253   || defined (H_STDIO)                                          \
    254   || defined (_H_STDIO)               /* AIX */                 \
    255   || defined (_STDIO_H)               /* glibc, Sun, SCO */     \
    256   || defined (_STDIO_H_)              /* BSD, OSF */            \
    257   || defined (__STDIO_H)              /* Borland */             \
    258   || defined (__STDIO_H__)            /* IRIX */                \
    259   || defined (_STDIO_INCLUDED)        /* HPUX */                \
    260   || defined (__dj_include_stdio_h_)  /* DJGPP */               \
    261   || defined (_FILE_DEFINED)          /* Microsoft */           \
    262   || defined (__STDIO__)              /* Apple MPW MrC */       \
    263   || defined (_MSL_STDIO_H)           /* Metrowerks */          \
    264   || defined (_STDIO_H_INCLUDED)      /* QNX4 */		\
    265   || defined (_ISO_STDIO_ISO_H)       /* Sun C++ */		\
    266   || defined (__STDIO_LOADED)         /* VMS */			\
    267   || defined (__DEFINED_FILE)         /* musl */
    268 #define _GMP_H_HAVE_FILE 1
    269 #endif
    270 
    271 /* In ISO C, if a prototype involving "struct obstack *" is given without
    272    that structure defined, then the struct is scoped down to just the
    273    prototype, causing a conflict if it's subsequently defined for real.  So
    274    only give prototypes if we've got obstack.h.  */
    275 #if defined (_OBSTACK_H)   /* glibc <obstack.h> */
    276 #define _GMP_H_HAVE_OBSTACK 1
    277 #endif
    278 
    279 /* The prototypes for gmp_vprintf etc are provided only if va_list is defined,
    280    via an application having included <stdarg.h>.  Usually va_list is a typedef
    281    so can't be tested directly, but C99 specifies that va_start is a macro.
    282 
    283    <stdio.h> will define some sort of va_list for vprintf and vfprintf, but
    284    let's not bother trying to use that since it's not standard and since
    285    application uses for gmp_vprintf etc will almost certainly require the
    286    whole <stdarg.h> anyway.  */
    287 
    288 #ifdef va_start
    289 #define _GMP_H_HAVE_VA_LIST 1
    290 #endif
    291 
    292 /* Test for gcc >= maj.min, as per __GNUC_PREREQ in glibc */
    293 #if defined (__GNUC__) && defined (__GNUC_MINOR__)
    294 #define __GMP_GNUC_PREREQ(maj, min) \
    295   ((__GNUC__ << 16) + __GNUC_MINOR__ >= ((maj) << 16) + (min))
    296 #else
    297 #define __GMP_GNUC_PREREQ(maj, min)  0
    298 #endif
    299 
    300 /* "pure" is in gcc 2.96 and up, see "(gcc)Function Attributes".  Basically
    301    it means a function does nothing but examine its arguments and memory
    302    (global or via arguments) to generate a return value, but changes nothing
    303    and has no side-effects.  __GMP_NO_ATTRIBUTE_CONST_PURE lets
    304    tune/common.c etc turn this off when trying to write timing loops.  */
    305 #if __GMP_GNUC_PREREQ (2,96) && ! defined (__GMP_NO_ATTRIBUTE_CONST_PURE)
    306 #define __GMP_ATTRIBUTE_PURE   __attribute__ ((__pure__))
    307 #else
    308 #define __GMP_ATTRIBUTE_PURE
    309 #endif
    310 
    311 
    312 /* __GMP_CAST allows us to use static_cast in C++, so our macros are clean
    313    to "g++ -Wold-style-cast".
    314 
    315    Casts in "extern inline" code within an extern "C" block don't induce
    316    these warnings, so __GMP_CAST only needs to be used on documented
    317    macros.  */
    318 
    319 #ifdef __cplusplus
    320 #define __GMP_CAST(type, expr)  (static_cast<type> (expr))
    321 #else
    322 #define __GMP_CAST(type, expr)  ((type) (expr))
    323 #endif
    324 
    325 
    326 /* An empty "throw ()" means the function doesn't throw any C++ exceptions,
    327    this can save some stack frame info in applications.
    328 
    329    Currently it's given only on functions which never divide-by-zero etc,
    330    don't allocate memory, and are expected to never need to allocate memory.
    331    This leaves open the possibility of a C++ throw from a future GMP
    332    exceptions scheme.
    333 
    334    mpz_set_ui etc are omitted to leave open the lazy allocation scheme
    335    described in doc/tasks.html.  mpz_get_d etc are omitted to leave open
    336    exceptions for float overflows.
    337 
    338    Note that __GMP_NOTHROW must be given on any inlines the same as on their
    339    prototypes (for g++ at least, where they're used together).  Note also
    340    that g++ 3.0 demands that __GMP_NOTHROW is before other attributes like
    341    __GMP_ATTRIBUTE_PURE.  */
    342 
    343 #if defined (__cplusplus)
    344 #if __cplusplus >= 201103L
    345 #define __GMP_NOTHROW  noexcept
    346 #else
    347 #define __GMP_NOTHROW  throw ()
    348 #endif
    349 #else
    350 #define __GMP_NOTHROW
    351 #endif
    352 
    353 
    354 /* PORTME: What other compilers have a useful "extern inline"?  "static
    355    inline" would be an acceptable substitute if the compiler (or linker)
    356    discards unused statics.  */
    357 
    358  /* gcc has __inline__ in all modes, including strict ansi.  Give a prototype
    359     for an inline too, so as to correctly specify "dllimport" on windows, in
    360     case the function is called rather than inlined.
    361     GCC 4.3 and above with -std=c99 or -std=gnu99 implements ISO C99
    362     inline semantics, unless -fgnu89-inline is used.  */
    363 #ifdef __GNUC__
    364 #if (defined __GNUC_STDC_INLINE__) || \
    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
   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
   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
   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
   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
   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
   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
   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
   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
   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
   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
   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
   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
   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
   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
   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
   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
   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
   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
   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
   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"
   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