1 /* gmp.h -- Definitions for GNU multiple precision functions.
2
3 Copyright (C) 1991, 1993, 1994, 1995, 1996 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 the GNU Lesser General Public License as published by
9 the Free Software Foundation; either version 2.1 of the License, or (at your
10 option) any later version.
11
12 The GNU MP Library is distributed in the hope that it will be useful, but
13 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
14 or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
15 License for more details.
16
17 You should have received a copy of the GNU Lesser General Public License
18 along with the GNU MP Library; see the file COPYING.LIB. If not, write to
19 the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
20 MA 02111-1307, USA. */
21
22 #ifndef __GMP_H__
23
24 #include <features.h>
25
26 #ifndef __GNU_MP__
27 #define __GNU_MP__ 2
28 #define __need_size_t
29 #include <stddef.h>
30 #undef __need_size_t
31
32 #if defined (__STDC__) || defined (__cplusplus)
33 #define __gmp_const const
34 #else
35 #define __gmp_const
36 #endif
37
38 #if defined (__GNUC__)
39 #define __gmp_inline __inline__
40 #else
41 #define __gmp_inline
42 #endif
43
44 #ifndef _EXTERN_INLINE
45 #ifdef __GNUC__
46 #define _EXTERN_INLINE __extern_inline
47 #else
48 #define _EXTERN_INLINE static
49 #endif
50 #endif
51
52 #ifdef _SHORT_LIMB
53 typedef unsigned int mp_limb_t;
54 typedef int mp_limb_signed_t;
55 #else
56 #ifdef _LONG_LONG_LIMB
57 typedef unsigned long long int mp_limb_t;
58 typedef long long int mp_limb_signed_t;
59 #else
60 typedef unsigned long int mp_limb_t;
61 typedef long int mp_limb_signed_t;
62 #endif
63 #endif
64
65 typedef mp_limb_t * mp_ptr;
66 typedef __gmp_const mp_limb_t * mp_srcptr;
67 typedef long int mp_size_t;
68 typedef long int mp_exp_t;
69
70 #ifndef __MP_SMALL__
71 typedef struct
72 {
73 int _mp_alloc; /* Number of *limbs* allocated and pointed
74 to by the D field. */
75 int _mp_size; /* abs(SIZE) is the number of limbs
76 the last field points to. If SIZE
77 is negative this is a negative
78 number. */
79 mp_limb_t *_mp_d; /* Pointer to the limbs. */
80 } __mpz_struct;
81 #else
82 typedef struct
83 {
84 short int _mp_alloc; /* Number of *limbs* allocated and pointed
85 to by the D field. */
86 short int _mp_size; /* abs(SIZE) is the number of limbs
87 the last field points to. If SIZE
88 is negative this is a negative
89 number. */
90 mp_limb_t *_mp_d; /* Pointer to the limbs. */
91 } __mpz_struct;
92 #endif
93 #endif /* __GNU_MP__ */
94
95 /* User-visible types. */
96 typedef __mpz_struct MP_INT;
97 typedef __mpz_struct mpz_t[1];
98
99 /* Structure for rational numbers. Zero is represented as 0/any, i.e.
100 the denominator is ignored. Negative numbers have the sign in
101 the numerator. */
102 typedef struct
103 {
104 __mpz_struct _mp_num;
105 __mpz_struct _mp_den;
106 #if 0
107 int _mp_num_alloc; /* Number of limbs allocated
108 for the numerator. */
109 int _mp_num_size; /* The absolute value of this field is the
110 length of the numerator; the sign is the
111 sign of the entire rational number. */
112 mp_ptr _mp_num; /* Pointer to the numerator limbs. */
113 int _mp_den_alloc; /* Number of limbs allocated
114 for the denominator. */
115 int _mp_den_size; /* Length of the denominator. (This field
116 should always be positive.) */
117 mp_ptr _mp_den; /* Pointer to the denominator limbs. */
118 #endif
119 } __mpq_struct;
120
121 typedef __mpq_struct MP_RAT;
122 typedef __mpq_struct mpq_t[1];
123
124 typedef struct
125 {
126 int _mp_prec; /* Max precision, in number of `mp_limb_t's.
127 Set by mpf_init and modified by
128 mpf_set_prec. The area pointed to
129 by the `d' field contains `prec' + 1
130 limbs. */
131 int _mp_size; /* abs(SIZE) is the number of limbs
132 the last field points to. If SIZE
133 is negative this is a negative
134 number. */
135 mp_exp_t _mp_exp; /* Exponent, in the base of `mp_limb_t'. */
136 mp_limb_t *_mp_d; /* Pointer to the limbs. */
137 } __mpf_struct;
138
139 /* typedef __mpf_struct MP_FLOAT; */
140 typedef __mpf_struct mpf_t[1];
141
142 /* Types for function declarations in gmp files. */
143 /* ??? Should not pollute user name space with these ??? */
144 typedef __gmp_const __mpz_struct *mpz_srcptr;
145 typedef __mpz_struct *mpz_ptr;
146 typedef __gmp_const __mpf_struct *mpf_srcptr;
147 typedef __mpf_struct *mpf_ptr;
148 typedef __gmp_const __mpq_struct *mpq_srcptr;
149 typedef __mpq_struct *mpq_ptr;
150
151 #ifndef _PROTO
152 #if defined (__STDC__) || defined (__cplusplus)
153 #define _PROTO(x) x
154 #else
155 #define _PROTO(x) ()
156 #endif
157 #endif
158
159 #ifndef __MPN
160 #if defined (__STDC__) || defined (__cplusplus)
161 #define __MPN(x) __mpn_##x
162 #else
163 #define __MPN(x) __mpn_/**/x
164 #endif
165 #endif
166
167 #if defined (FILE) || defined (_STDIO_H_) || defined (__STDIO_H__) || defined (H_STDIO)
168 #define _GMP_H_HAVE_FILE 1
169 #endif
170
171 void mp_set_memory_functions _PROTO ((void *(*) (size_t),
172 void *(*) (void *, size_t, size_t),
173 void (*) (void *, size_t)));
174 extern const int mp_bits_per_limb;
175
176 /**************** Integer (i.e. Z) routines. ****************/
177
178 #if defined (__cplusplus)
179 extern "C" {
180 #endif
181 void *_mpz_realloc _PROTO ((mpz_ptr, mp_size_t));
182
183 void mpz_abs _PROTO ((mpz_ptr, mpz_srcptr));
184 void mpz_add _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
185 void mpz_add_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
186 void mpz_and _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
187 void mpz_array_init _PROTO ((mpz_ptr, mp_size_t, mp_size_t));
188 void mpz_cdiv_q _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
189 unsigned long int mpz_cdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
190 void mpz_cdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
191 unsigned long int mpz_cdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
192 void mpz_cdiv_r _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
193 unsigned long int mpz_cdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
194 unsigned long int mpz_cdiv_ui _PROTO ((mpz_srcptr, unsigned long int));
195 void mpz_clear _PROTO ((mpz_ptr));
196 void mpz_clrbit _PROTO ((mpz_ptr, unsigned long int));
197 int mpz_cmp _PROTO ((mpz_srcptr, mpz_srcptr));
198 int mpz_cmp_si _PROTO ((mpz_srcptr, signed long int));
199 int mpz_cmp_ui _PROTO ((mpz_srcptr, unsigned long int));
200 void mpz_com _PROTO ((mpz_ptr, mpz_srcptr));
201 void mpz_divexact _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
202 void mpz_fac_ui _PROTO ((mpz_ptr, unsigned long int));
203 void mpz_fdiv_q _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
204 void mpz_fdiv_q_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
205 unsigned long int mpz_fdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
206 void mpz_fdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
207 unsigned long int mpz_fdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
208 void mpz_fdiv_r _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
209 void mpz_fdiv_r_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
210 unsigned long int mpz_fdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
211 unsigned long int mpz_fdiv_ui _PROTO ((mpz_srcptr, unsigned long int));
212 void mpz_gcd _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
213 unsigned long int mpz_gcd_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
214 void mpz_gcdext _PROTO ((mpz_ptr, mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
215 /* signed */ long int mpz_get_si _PROTO ((mpz_srcptr));
216 char *mpz_get_str _PROTO ((char *, int, mpz_srcptr));
217 unsigned long int mpz_get_ui _PROTO ((mpz_srcptr));
218 mp_limb_t mpz_getlimbn _PROTO ((mpz_srcptr, mp_size_t));
219 unsigned long int mpz_hamdist _PROTO ((mpz_srcptr, mpz_srcptr));
220 void mpz_init _PROTO ((mpz_ptr));
221 #ifdef _GMP_H_HAVE_FILE
222 size_t mpz_inp_binary _PROTO ((mpz_ptr, FILE *));
223 size_t mpz_inp_raw _PROTO ((mpz_ptr, FILE *));
224 size_t mpz_inp_str _PROTO ((mpz_ptr, FILE *, int));
225 #endif
226 void mpz_init_set _PROTO ((mpz_ptr, mpz_srcptr));
227 void mpz_init_set_d _PROTO ((mpz_ptr, double));
228 void mpz_init_set_si _PROTO ((mpz_ptr, signed long int));
229 int mpz_init_set_str _PROTO ((mpz_ptr, const char *, int));
230 void mpz_init_set_ui _PROTO ((mpz_ptr, unsigned long int));
231 int mpz_invert _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
232 void mpz_ior _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
233 int mpz_jacobi _PROTO ((mpz_srcptr, mpz_srcptr));
234 int mpz_legendre _PROTO ((mpz_srcptr, mpz_srcptr));
235 void mpz_mod _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
236 void mpz_mul _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
237 void mpz_mul_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
238 void mpz_mul_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
239 void mpz_neg _PROTO ((mpz_ptr, mpz_srcptr));
240 #ifdef _GMP_H_HAVE_FILE
241 size_t mpz_out_binary _PROTO ((FILE *, mpz_srcptr));
242 size_t mpz_out_raw _PROTO ((FILE *, mpz_srcptr));
243 size_t mpz_out_str _PROTO ((FILE *, int, mpz_srcptr));
244 #endif
245 int mpz_perfect_square_p _PROTO ((mpz_srcptr));
246 unsigned long int mpz_popcount _PROTO ((mpz_srcptr));
247 void mpz_pow_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
248 void mpz_powm _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr, mpz_srcptr));
249 void mpz_powm_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int, mpz_srcptr));
250 int mpz_probab_prime_p _PROTO ((mpz_srcptr, int));
251 void mpz_random _PROTO ((mpz_ptr, mp_size_t));
252 void mpz_random2 _PROTO ((mpz_ptr, mp_size_t));
253 unsigned long int mpz_scan0 _PROTO ((mpz_srcptr, unsigned long int));
254 unsigned long int mpz_scan1 _PROTO ((mpz_srcptr, unsigned long int));
255 void mpz_set _PROTO ((mpz_ptr, mpz_srcptr));
256 void mpz_set_d _PROTO ((mpz_ptr, double));
257 void mpz_set_si _PROTO ((mpz_ptr, signed long int));
258 int mpz_set_str _PROTO ((mpz_ptr, const char *, int));
259 void mpz_set_ui _PROTO ((mpz_ptr, unsigned long int));
260 void mpz_setbit _PROTO ((mpz_ptr, unsigned long int));
261 size_t mpz_size _PROTO ((mpz_srcptr));
262 size_t mpz_sizeinbase _PROTO ((mpz_srcptr, int));
263 void mpz_sqrt _PROTO ((mpz_ptr, mpz_srcptr));
264 void mpz_sqrtrem _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr));
265 void mpz_sub _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
266 void mpz_sub_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
267 void mpz_tdiv_q _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
268 void mpz_tdiv_q_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
269 void mpz_tdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
270 void mpz_tdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
271 void mpz_tdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
272 void mpz_tdiv_r _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
273 void mpz_tdiv_r_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
274 void mpz_tdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
275 void mpz_ui_pow_ui _PROTO ((mpz_ptr, unsigned long int, unsigned long int));
276
277 /**************** Rational (i.e. Q) routines. ****************/
278
279 void mpq_init _PROTO ((mpq_ptr));
280 void mpq_clear _PROTO ((mpq_ptr));
281 void mpq_set _PROTO ((mpq_ptr, mpq_srcptr));
282 void mpq_set_ui _PROTO ((mpq_ptr, unsigned long int, unsigned long int));
283 void mpq_set_si _PROTO ((mpq_ptr, signed long int, unsigned long int));
284 void mpq_add _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
285 void mpq_sub _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
286 void mpq_mul _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
287 void mpq_div _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
288 void mpq_neg _PROTO ((mpq_ptr, mpq_srcptr));
289 int mpq_cmp _PROTO ((mpq_srcptr, mpq_srcptr));
290 int mpq_cmp_ui _PROTO ((mpq_srcptr, unsigned long int, unsigned long int));
291 void mpq_inv _PROTO ((mpq_ptr, mpq_srcptr));
292 void mpq_set_num _PROTO ((mpq_ptr, mpz_srcptr));
293 void mpq_set_den _PROTO ((mpq_ptr, mpz_srcptr));
294 void mpq_get_num _PROTO ((mpz_ptr, mpq_srcptr));
295 void mpq_get_den _PROTO ((mpz_ptr, mpq_srcptr));
296 double mpq_get_d _PROTO ((mpq_srcptr));
297 void mpq_canonicalize _PROTO ((mpq_ptr));
298
299 /**************** Float (i.e. F) routines. ****************/
300
301 void mpf_abs _PROTO ((mpf_ptr, mpf_srcptr));
302 void mpf_add _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
303 void mpf_add_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
304 void mpf_clear _PROTO ((mpf_ptr));
305 int mpf_cmp _PROTO ((mpf_srcptr, mpf_srcptr));
306 int mpf_cmp_si _PROTO ((mpf_srcptr, signed long int));
307 int mpf_cmp_ui _PROTO ((mpf_srcptr, unsigned long int));
308 void mpf_div _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
309 void mpf_div_2exp _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
310 void mpf_div_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
311 void mpf_dump _PROTO ((mpf_srcptr));
312 int mpf_eq _PROTO ((mpf_srcptr, mpf_srcptr, unsigned long int));
313 unsigned long int mpf_get_prec _PROTO ((mpf_srcptr));
314 char *mpf_get_str _PROTO ((char *, mp_exp_t *, int, size_t, mpf_srcptr));
315 void mpf_init _PROTO ((mpf_ptr));
316 void mpf_init2 _PROTO ((mpf_ptr, unsigned long int));
317 #ifdef _GMP_H_HAVE_FILE
318 size_t mpf_inp_str _PROTO ((mpf_ptr, FILE *, int));
319 #endif
320 void mpf_init_set _PROTO ((mpf_ptr, mpf_srcptr));
321 void mpf_init_set_d _PROTO ((mpf_ptr, double));
322 void mpf_init_set_si _PROTO ((mpf_ptr, signed long int));
323 int mpf_init_set_str _PROTO ((mpf_ptr, char *, int));
324 void mpf_init_set_ui _PROTO ((mpf_ptr, unsigned long int));
325 void mpf_mul _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
326 void mpf_mul_2exp _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
327 void mpf_mul_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
328 void mpf_neg _PROTO ((mpf_ptr, mpf_srcptr));
329 #ifdef _GMP_H_HAVE_FILE
330 size_t mpf_out_str _PROTO ((FILE *, int, size_t, mpf_srcptr));
331 #endif
332 void mpf_random2 _PROTO ((mpf_ptr, mp_size_t, mp_exp_t));
333 void mpf_reldiff _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
334 void mpf_set _PROTO ((mpf_ptr, mpf_srcptr));
335 void mpf_set_d _PROTO ((mpf_ptr, double));
336 void mpf_set_default_prec _PROTO ((unsigned long int));
337 void mpf_set_prec _PROTO ((mpf_ptr, unsigned long int));
338 void mpf_set_prec_raw _PROTO ((mpf_ptr, unsigned long int));
339 void mpf_set_si _PROTO ((mpf_ptr, signed long int));
340 int mpf_set_str _PROTO ((mpf_ptr, const char *, int));
341 void mpf_set_ui _PROTO ((mpf_ptr, unsigned long int));
342 size_t mpf_size _PROTO ((mpf_srcptr));
343 void mpf_sqrt _PROTO ((mpf_ptr, mpf_srcptr));
344 void mpf_sqrt_ui _PROTO ((mpf_ptr, unsigned long int));
345 void mpf_sub _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
346 void mpf_sub_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
347 void mpf_ui_div _PROTO ((mpf_ptr, unsigned long int, mpf_srcptr));
348 void mpf_ui_sub _PROTO ((mpf_ptr, unsigned long int, mpf_srcptr));
349 #if defined (__cplusplus)
350 }
351 #endif
352 /************ Low level positive-integer (i.e. N) routines. ************/
353
354 /* This is ugly, but we need to make usr calls reach the prefixed function. */
355 #define mpn_add __MPN(add)
356 #define mpn_add_1 __MPN(add_1)
357 #define mpn_add_n __MPN(add_n)
358 #define mpn_addmul_1 __MPN(addmul_1)
359 #define mpn_bdivmod __MPN(bdivmod)
360 #define mpn_cmp __MPN(cmp)
361 #define mpn_divmod_1 __MPN(divmod_1)
362 #define mpn_divrem __MPN(divrem)
363 #define mpn_divrem_1 __MPN(divrem_1)
364 #define mpn_dump __MPN(dump)
365 #define mpn_gcd __MPN(gcd)
366 #define mpn_gcd_1 __MPN(gcd_1)
367 #define mpn_gcdext __MPN(gcdext)
368 #define mpn_get_str __MPN(get_str)
369 #define mpn_hamdist __MPN(hamdist)
370 #define mpn_lshift __MPN(lshift)
371 #define mpn_mod_1 __MPN(mod_1)
372 #define mpn_mul __MPN(mul)
373 #define mpn_mul_1 __MPN(mul_1)
374 #define mpn_mul_n __MPN(mul_n)
375 #define mpn_perfect_square_p __MPN(perfect_square_p)
376 #define mpn_popcount __MPN(popcount)
377 #define mpn_preinv_mod_1 __MPN(preinv_mod_1)
378 #define mpn_random2 __MPN(random2)
379 #define mpn_rshift __MPN(rshift)
380 #define mpn_scan0 __MPN(scan0)
381 #define mpn_scan1 __MPN(scan1)
382 #define mpn_set_str __MPN(set_str)
383 #define mpn_sqrtrem __MPN(sqrtrem)
384 #define mpn_sub __MPN(sub)
385 #define mpn_sub_1 __MPN(sub_1)
386 #define mpn_sub_n __MPN(sub_n)
387 #define mpn_submul_1 __MPN(submul_1)
388 #define mpn_udiv_w_sdiv __MPN(udiv_w_sdiv)
389
390 #if defined (__cplusplus)
391 extern "C" {
392 #endif
393 mp_limb_t mpn_add _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr,mp_size_t));
394 mp_limb_t mpn_add_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
395 mp_limb_t mpn_add_n _PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t));
396 mp_limb_t mpn_addmul_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
397 mp_limb_t mpn_bdivmod _PROTO ((mp_ptr, mp_ptr, mp_size_t, mp_srcptr, mp_size_t, unsigned long int));
398 int mpn_cmp _PROTO ((mp_srcptr, mp_srcptr, mp_size_t));
399 mp_limb_t mpn_divmod_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
400 mp_limb_t mpn_divrem _PROTO ((mp_ptr, mp_size_t, mp_ptr, mp_size_t, mp_srcptr, mp_size_t));
401 mp_limb_t mpn_divrem_1 _PROTO ((mp_ptr, mp_size_t, mp_srcptr, mp_size_t, mp_limb_t));
402 void mpn_dump _PROTO ((mp_srcptr, mp_size_t));
403 mp_size_t mpn_gcd _PROTO ((mp_ptr, mp_ptr, mp_size_t, mp_ptr, mp_size_t));
404 mp_limb_t mpn_gcd_1 _PROTO ((mp_srcptr, mp_size_t, mp_limb_t));
405 mp_size_t mpn_gcdext _PROTO ((mp_ptr, mp_ptr, mp_ptr, mp_size_t, mp_ptr, mp_size_t));
406 size_t mpn_get_str _PROTO ((unsigned char *, int, mp_ptr, mp_size_t));
407 unsigned long int mpn_hamdist _PROTO ((mp_srcptr, mp_srcptr, mp_size_t));
408 mp_limb_t mpn_lshift _PROTO ((mp_ptr, mp_srcptr, mp_size_t, unsigned int));
409 mp_limb_t mpn_mod_1 _PROTO ((mp_srcptr, mp_size_t, mp_limb_t));
410 mp_limb_t mpn_mul _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t));
411 mp_limb_t mpn_mul_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
412 void mpn_mul_n _PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t));
413 int mpn_perfect_square_p _PROTO ((mp_srcptr, mp_size_t));
414 unsigned long int mpn_popcount _PROTO ((mp_srcptr, mp_size_t));
415 mp_limb_t mpn_preinv_mod_1 _PROTO ((mp_srcptr, mp_size_t, mp_limb_t, mp_limb_t));
416 void mpn_random2 _PROTO ((mp_ptr, mp_size_t));
417 mp_limb_t mpn_rshift _PROTO ((mp_ptr, mp_srcptr, mp_size_t, unsigned int));
418 unsigned long int mpn_scan0 _PROTO ((mp_srcptr, unsigned long int));
419 unsigned long int mpn_scan1 _PROTO ((mp_srcptr, unsigned long int));
420 mp_size_t mpn_set_str _PROTO ((mp_ptr, const unsigned char *, size_t, int));
421 mp_size_t mpn_sqrtrem _PROTO ((mp_ptr, mp_ptr, mp_srcptr, mp_size_t));
422 mp_limb_t mpn_sub _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr,mp_size_t));
423 mp_limb_t mpn_sub_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
424 mp_limb_t mpn_sub_n _PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t));
425 mp_limb_t mpn_submul_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
426 #if defined (__cplusplus)
427 }
428 #endif
429
430 #if defined (__GNUC__) || defined (_FORCE_INLINES)
431 _EXTERN_INLINE mp_limb_t
432 #if defined (__STDC__) || defined (__cplusplus)
mpn_add_1(register mp_ptr res_ptr,register mp_srcptr s1_ptr,register mp_size_t s1_size,register mp_limb_t s2_limb)433 mpn_add_1 (register mp_ptr res_ptr,
434 register mp_srcptr s1_ptr,
435 register mp_size_t s1_size,
436 register mp_limb_t s2_limb)
437 #else
438 mpn_add_1 (res_ptr, s1_ptr, s1_size, s2_limb)
439 register mp_ptr res_ptr;
440 register mp_srcptr s1_ptr;
441 register mp_size_t s1_size;
442 register mp_limb_t s2_limb;
443 #endif
444 {
445 register mp_limb_t x;
446
447 x = *s1_ptr++;
448 s2_limb = x + s2_limb;
449 *res_ptr++ = s2_limb;
450 if (s2_limb < x)
451 {
452 while (--s1_size != 0)
453 {
454 x = *s1_ptr++ + 1;
455 *res_ptr++ = x;
456 if (x != 0)
457 goto fin;
458 }
459
460 return 1;
461 }
462
463 fin:
464 if (res_ptr != s1_ptr)
465 {
466 mp_size_t i;
467 for (i = 0; i < s1_size - 1; i++)
468 res_ptr[i] = s1_ptr[i];
469 }
470 return 0;
471 }
472
473 _EXTERN_INLINE mp_limb_t
474 #if defined (__STDC__) || defined (__cplusplus)
mpn_add(register mp_ptr res_ptr,register mp_srcptr s1_ptr,register mp_size_t s1_size,register mp_srcptr s2_ptr,register mp_size_t s2_size)475 mpn_add (register mp_ptr res_ptr,
476 register mp_srcptr s1_ptr,
477 register mp_size_t s1_size,
478 register mp_srcptr s2_ptr,
479 register mp_size_t s2_size)
480 #else
481 mpn_add (res_ptr, s1_ptr, s1_size, s2_ptr, s2_size)
482 register mp_ptr res_ptr;
483 register mp_srcptr s1_ptr;
484 register mp_size_t s1_size;
485 register mp_srcptr s2_ptr;
486 register mp_size_t s2_size;
487 #endif
488 {
489 mp_limb_t cy_limb = 0;
490
491 if (s2_size != 0)
492 cy_limb = mpn_add_n (res_ptr, s1_ptr, s2_ptr, s2_size);
493
494 if (s1_size - s2_size != 0)
495 cy_limb = mpn_add_1 (res_ptr + s2_size,
496 s1_ptr + s2_size,
497 s1_size - s2_size,
498 cy_limb);
499 return cy_limb;
500 }
501
502 _EXTERN_INLINE mp_limb_t
503 #if defined (__STDC__) || defined (__cplusplus)
mpn_sub_1(register mp_ptr res_ptr,register mp_srcptr s1_ptr,register mp_size_t s1_size,register mp_limb_t s2_limb)504 mpn_sub_1 (register mp_ptr res_ptr,
505 register mp_srcptr s1_ptr,
506 register mp_size_t s1_size,
507 register mp_limb_t s2_limb)
508 #else
509 mpn_sub_1 (res_ptr, s1_ptr, s1_size, s2_limb)
510 register mp_ptr res_ptr;
511 register mp_srcptr s1_ptr;
512 register mp_size_t s1_size;
513 register mp_limb_t s2_limb;
514 #endif
515 {
516 register mp_limb_t x;
517
518 x = *s1_ptr++;
519 s2_limb = x - s2_limb;
520 *res_ptr++ = s2_limb;
521 if (s2_limb > x)
522 {
523 while (--s1_size != 0)
524 {
525 x = *s1_ptr++;
526 *res_ptr++ = x - 1;
527 if (x != 0)
528 goto fin;
529 }
530
531 return 1;
532 }
533
534 fin:
535 if (res_ptr != s1_ptr)
536 {
537 mp_size_t i;
538 for (i = 0; i < s1_size - 1; i++)
539 res_ptr[i] = s1_ptr[i];
540 }
541 return 0;
542 }
543
544 _EXTERN_INLINE mp_limb_t
545 #if defined (__STDC__) || defined (__cplusplus)
mpn_sub(register mp_ptr res_ptr,register mp_srcptr s1_ptr,register mp_size_t s1_size,register mp_srcptr s2_ptr,register mp_size_t s2_size)546 mpn_sub (register mp_ptr res_ptr,
547 register mp_srcptr s1_ptr,
548 register mp_size_t s1_size,
549 register mp_srcptr s2_ptr,
550 register mp_size_t s2_size)
551 #else
552 mpn_sub (res_ptr, s1_ptr, s1_size, s2_ptr, s2_size)
553 register mp_ptr res_ptr;
554 register mp_srcptr s1_ptr;
555 register mp_size_t s1_size;
556 register mp_srcptr s2_ptr;
557 register mp_size_t s2_size;
558 #endif
559 {
560 mp_limb_t cy_limb = 0;
561
562 if (s2_size != 0)
563 cy_limb = mpn_sub_n (res_ptr, s1_ptr, s2_ptr, s2_size);
564
565 if (s1_size - s2_size != 0)
566 cy_limb = mpn_sub_1 (res_ptr + s2_size,
567 s1_ptr + s2_size,
568 s1_size - s2_size,
569 cy_limb);
570 return cy_limb;
571 }
572 #endif /* __GNUC__ */
573
574 /* Allow faster testing for negative, zero, and positive. */
575 #define mpz_sgn(Z) ((Z)->_mp_size < 0 ? -1 : (Z)->_mp_size > 0)
576 #define mpf_sgn(F) ((F)->_mp_size < 0 ? -1 : (F)->_mp_size > 0)
577 #define mpq_sgn(Q) ((Q)->_mp_num._mp_size < 0 ? -1 : (Q)->_mp_num._mp_size > 0)
578
579 /* Allow direct user access to numerator and denominator of a mpq_t object. */
580 #define mpq_numref(Q) (&((Q)->_mp_num))
581 #define mpq_denref(Q) (&((Q)->_mp_den))
582
583 /* When using GCC, optimize certain common comparisons. */
584 #if defined (__GNUC__)
585 #define mpz_cmp_ui(Z,UI) \
586 (__builtin_constant_p (UI) && (UI) == 0 \
587 ? mpz_sgn (Z) : mpz_cmp_ui (Z,UI))
588 #define mpz_cmp_si(Z,UI) \
589 (__builtin_constant_p (UI) && (UI) == 0 ? mpz_sgn (Z) \
590 : __builtin_constant_p (UI) && (UI) > 0 ? mpz_cmp_ui (Z,UI) \
591 : mpz_cmp_si (Z,UI))
592 #define mpq_cmp_ui(Q,NUI,DUI) \
593 (__builtin_constant_p (NUI) && (NUI) == 0 \
594 ? mpq_sgn (Q) : mpq_cmp_ui (Q,NUI,DUI))
595 #endif
596
597 #define mpn_divmod(qp,np,nsize,dp,dsize) mpn_divrem (qp,0,np,nsize,dp,dsize)
598 #if 0
599 #define mpn_divmod_1(qp,np,nsize,dlimb) mpn_divrem_1 (qp,0,np,nsize,dlimb)
600 #endif
601
602 /* Compatibility with GMP 1. */
603 #define mpz_mdiv mpz_fdiv_q
604 #define mpz_mdivmod mpz_fdiv_qr
605 #define mpz_mmod mpz_fdiv_r
606 #define mpz_mdiv_ui mpz_fdiv_q_ui
607 #define mpz_mdivmod_ui(q,r,n,d) \
608 ((r == 0) ? mpz_fdiv_q_ui (q,n,d) : mpz_fdiv_qr_ui (q,r,n,d))
609 #define mpz_mmod_ui(r,n,d) \
610 ((r == 0) ? mpz_fdiv_ui (n,d) : mpz_fdiv_r_ui (r,n,d))
611
612 /* Useful synonyms, but not quite compatible with GMP 1. */
613 #define mpz_div mpz_fdiv_q
614 #define mpz_divmod mpz_fdiv_qr
615 #define mpz_div_ui mpz_fdiv_q_ui
616 #define mpz_divmod_ui mpz_fdiv_qr_ui
617 #define mpz_mod_ui mpz_fdiv_r_ui
618 #define mpz_div_2exp mpz_fdiv_q_2exp
619 #define mpz_mod_2exp mpz_fdiv_r_2exp
620
621 #define __GNU_MP_VERSION 2
622 #define __GNU_MP_VERSION_MINOR 0
623 #define __GMP_H__
624 #endif /* __GMP_H__ */
625