gnu/usr.bin/gmp/gmp.h

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