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