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