1 /* gmp.h -- Definitions for GNU multiple precision functions.
3 Copyright (C) 1991, 1993, 1994, 1995, 1996 Free Software Foundation, Inc.
5 This file is part of the GNU MP Library.
7 The GNU MP Library is free software; you can redistribute it and/or modify
8 it under the terms of the GNU Library General Public License as published by
9 the Free Software Foundation; either version 2 of the License, or (at your
10 option) any later version.
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 Library General Public
15 License for more details.
17 You should have received a copy of the GNU Library 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. */
25 /* Get DLL_IMPORT_RTS */
26 #include "../../includes/config.h"
27 #include "../../includes/StgDLL.h"
36 #if defined (__STDC__) || defined (__cplusplus)
37 #define __gmp_const const
42 #if defined (__GNUC__)
43 #define __gmp_inline __inline__
48 #ifndef _EXTERN_INLINE
50 #define _EXTERN_INLINE extern __inline__
52 #define _EXTERN_INLINE static
57 typedef unsigned int mp_limb_t;
58 typedef int mp_limb_signed_t;
60 #ifdef _LONG_LONG_LIMB
61 typedef unsigned long long int mp_limb_t;
62 typedef long long int mp_limb_signed_t;
64 typedef unsigned long int mp_limb_t;
65 typedef long int mp_limb_signed_t;
69 typedef mp_limb_t * mp_ptr;
70 typedef __gmp_const mp_limb_t * mp_srcptr;
71 typedef long int mp_size_t;
72 typedef long int mp_exp_t;
77 int _mp_alloc; /* Number of *limbs* allocated and pointed
79 int _mp_size; /* abs(SIZE) is the number of limbs
80 the last field points to. If SIZE
81 is negative this is a negative
83 mp_limb_t *_mp_d; /* Pointer to the limbs. */
88 short int _mp_alloc; /* Number of *limbs* allocated and pointed
90 short int _mp_size; /* abs(SIZE) is the number of limbs
91 the last field points to. If SIZE
92 is negative this is a negative
94 mp_limb_t *_mp_d; /* Pointer to the limbs. */
97 #endif /* __GNU_MP__ */
99 /* User-visible types. */
100 typedef __mpz_struct MP_INT;
101 typedef __mpz_struct mpz_t[1];
103 /* Structure for rational numbers. Zero is represented as 0/any, i.e.
104 the denominator is ignored. Negative numbers have the sign in
108 __mpz_struct _mp_num;
109 __mpz_struct _mp_den;
111 int _mp_num_alloc; /* Number of limbs allocated
112 for the numerator. */
113 int _mp_num_size; /* The absolute value of this field is the
114 length of the numerator; the sign is the
115 sign of the entire rational number. */
116 mp_ptr _mp_num; /* Pointer to the numerator limbs. */
117 int _mp_den_alloc; /* Number of limbs allocated
118 for the denominator. */
119 int _mp_den_size; /* Length of the denominator. (This field
120 should always be positive.) */
121 mp_ptr _mp_den; /* Pointer to the denominator limbs. */
125 typedef __mpq_struct MP_RAT;
126 typedef __mpq_struct mpq_t[1];
130 int _mp_prec; /* Max precision, in number of `mp_limb_t's.
131 Set by mpf_init and modified by
132 mpf_set_prec. The area pointed to
133 by the `d' field contains `prec' + 1
135 int _mp_size; /* abs(SIZE) is the number of limbs
136 the last field points to. If SIZE
137 is negative this is a negative
139 mp_exp_t _mp_exp; /* Exponent, in the base of `mp_limb_t'. */
140 mp_limb_t *_mp_d; /* Pointer to the limbs. */
143 /* typedef __mpf_struct MP_FLOAT; */
144 typedef __mpf_struct mpf_t[1];
146 /* Types for function declarations in gmp files. */
147 /* ??? Should not pollute user name space with these ??? */
148 typedef __gmp_const __mpz_struct *mpz_srcptr;
149 typedef __mpz_struct *mpz_ptr;
150 typedef __gmp_const __mpf_struct *mpf_srcptr;
151 typedef __mpf_struct *mpf_ptr;
152 typedef __gmp_const __mpq_struct *mpq_srcptr;
153 typedef __mpq_struct *mpq_ptr;
156 #if defined (__STDC__) || defined (__cplusplus)
164 #if defined (__STDC__) || defined (__cplusplus)
165 #define __MPN(x) __mpn_##x
167 #define __MPN(x) __mpn_/**/x
171 #if defined (FILE) || defined (H_STDIO) || defined (_H_STDIO) \
172 || defined (_STDIO_H) || defined (_STDIO_H_) || defined (__STDIO_H__) \
173 || defined (_STDIO_INCLUDED)
174 #define _GMP_H_HAVE_FILE 1
177 DLL_IMPORT_RTS void mp_set_memory_functions _PROTO ((void *(*) (size_t),
178 void *(*) (void *, size_t, size_t),
179 void (*) (void *, size_t)));
180 DLL_IMPORT_RTS extern __gmp_const int mp_bits_per_limb;
182 /**************** Integer (i.e. Z) routines. ****************/
184 #if defined (__cplusplus)
187 DLL_IMPORT_RTS void *_mpz_realloc _PROTO ((mpz_ptr, mp_size_t));
189 DLL_IMPORT_RTS void mpz_abs _PROTO ((mpz_ptr, mpz_srcptr));
190 DLL_IMPORT_RTS void mpz_add _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
191 DLL_IMPORT_RTS void mpz_add_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
192 DLL_IMPORT_RTS void mpz_and _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
193 DLL_IMPORT_RTS void mpz_array_init _PROTO ((mpz_ptr, mp_size_t, mp_size_t));
194 DLL_IMPORT_RTS void mpz_cdiv_q _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
195 DLL_IMPORT_RTS unsigned long int mpz_cdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
196 DLL_IMPORT_RTS void mpz_cdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
197 DLL_IMPORT_RTS unsigned long int mpz_cdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
198 DLL_IMPORT_RTS void mpz_cdiv_r _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
199 DLL_IMPORT_RTS unsigned long int mpz_cdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
200 DLL_IMPORT_RTS unsigned long int mpz_cdiv_ui _PROTO ((mpz_srcptr, unsigned long int));
201 DLL_IMPORT_RTS void mpz_clear _PROTO ((mpz_ptr));
202 DLL_IMPORT_RTS void mpz_clrbit _PROTO ((mpz_ptr, unsigned long int));
203 DLL_IMPORT_RTS int mpz_cmp _PROTO ((mpz_srcptr, mpz_srcptr));
204 DLL_IMPORT_RTS int mpz_cmp_si _PROTO ((mpz_srcptr, signed long int));
205 DLL_IMPORT_RTS int mpz_cmp_ui _PROTO ((mpz_srcptr, unsigned long int));
206 DLL_IMPORT_RTS void mpz_com _PROTO ((mpz_ptr, mpz_srcptr));
207 DLL_IMPORT_RTS void mpz_divexact _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
208 DLL_IMPORT_RTS void mpz_fac_ui _PROTO ((mpz_ptr, unsigned long int));
209 DLL_IMPORT_RTS void mpz_fdiv_q _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
210 DLL_IMPORT_RTS void mpz_fdiv_q_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
211 DLL_IMPORT_RTS unsigned long int mpz_fdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
212 DLL_IMPORT_RTS void mpz_fdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
213 DLL_IMPORT_RTS unsigned long int mpz_fdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
214 DLL_IMPORT_RTS void mpz_fdiv_r _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
215 DLL_IMPORT_RTS void mpz_fdiv_r_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
216 DLL_IMPORT_RTS unsigned long int mpz_fdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
217 DLL_IMPORT_RTS unsigned long int mpz_fdiv_ui _PROTO ((mpz_srcptr, unsigned long int));
218 DLL_IMPORT_RTS void mpz_gcd _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
219 DLL_IMPORT_RTS unsigned long int mpz_gcd_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
220 DLL_IMPORT_RTS void mpz_gcdext _PROTO ((mpz_ptr, mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
221 DLL_IMPORT_RTS double mpz_get_d _PROTO ((mpz_srcptr));
222 /* signed */ long int mpz_get_si _PROTO ((mpz_srcptr));
223 DLL_IMPORT_RTS char *mpz_get_str _PROTO ((char *, int, mpz_srcptr));
224 DLL_IMPORT_RTS unsigned long int mpz_get_ui _PROTO ((mpz_srcptr));
225 DLL_IMPORT_RTS mp_limb_t mpz_getlimbn _PROTO ((mpz_srcptr, mp_size_t));
226 DLL_IMPORT_RTS unsigned long int mpz_hamdist _PROTO ((mpz_srcptr, mpz_srcptr));
227 DLL_IMPORT_RTS void mpz_init _PROTO ((mpz_ptr));
228 #ifdef _GMP_H_HAVE_FILE
229 DLL_IMPORT_RTS size_t mpz_inp_binary _PROTO ((mpz_ptr, FILE *));
230 DLL_IMPORT_RTS size_t mpz_inp_raw _PROTO ((mpz_ptr, FILE *));
231 DLL_IMPORT_RTS size_t mpz_inp_str _PROTO ((mpz_ptr, FILE *, int));
233 DLL_IMPORT_RTS void mpz_init_set _PROTO ((mpz_ptr, mpz_srcptr));
234 DLL_IMPORT_RTS void mpz_init_set_d _PROTO ((mpz_ptr, double));
235 DLL_IMPORT_RTS void mpz_init_set_si _PROTO ((mpz_ptr, signed long int));
236 DLL_IMPORT_RTS int mpz_init_set_str _PROTO ((mpz_ptr, const char *, int));
237 DLL_IMPORT_RTS void mpz_init_set_ui _PROTO ((mpz_ptr, unsigned long int));
238 DLL_IMPORT_RTS int mpz_invert _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
239 DLL_IMPORT_RTS void mpz_ior _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
240 DLL_IMPORT_RTS int mpz_jacobi _PROTO ((mpz_srcptr, mpz_srcptr));
241 DLL_IMPORT_RTS int mpz_legendre _PROTO ((mpz_srcptr, mpz_srcptr));
242 DLL_IMPORT_RTS void mpz_mod _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
243 DLL_IMPORT_RTS void mpz_mul _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
244 DLL_IMPORT_RTS void mpz_mul_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
245 DLL_IMPORT_RTS void mpz_mul_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
246 DLL_IMPORT_RTS void mpz_neg _PROTO ((mpz_ptr, mpz_srcptr));
247 #ifdef _GMP_H_HAVE_FILE
248 DLL_IMPORT_RTS size_t mpz_out_binary _PROTO ((FILE *, mpz_srcptr));
249 DLL_IMPORT_RTS size_t mpz_out_raw _PROTO ((FILE *, mpz_srcptr));
250 DLL_IMPORT_RTS size_t mpz_out_str _PROTO ((FILE *, int, mpz_srcptr));
252 DLL_IMPORT_RTS int mpz_perfect_square_p _PROTO ((mpz_srcptr));
253 DLL_IMPORT_RTS unsigned long int mpz_popcount _PROTO ((mpz_srcptr));
254 DLL_IMPORT_RTS void mpz_pow_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
255 DLL_IMPORT_RTS void mpz_powm _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr, mpz_srcptr));
256 DLL_IMPORT_RTS void mpz_powm_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int, mpz_srcptr));
257 DLL_IMPORT_RTS int mpz_probab_prime_p _PROTO ((mpz_srcptr, int));
258 DLL_IMPORT_RTS void mpz_random _PROTO ((mpz_ptr, mp_size_t));
259 DLL_IMPORT_RTS void mpz_random2 _PROTO ((mpz_ptr, mp_size_t));
260 DLL_IMPORT_RTS unsigned long int mpz_scan0 _PROTO ((mpz_srcptr, unsigned long int));
261 DLL_IMPORT_RTS unsigned long int mpz_scan1 _PROTO ((mpz_srcptr, unsigned long int));
262 DLL_IMPORT_RTS void mpz_set _PROTO ((mpz_ptr, mpz_srcptr));
263 DLL_IMPORT_RTS void mpz_set_d _PROTO ((mpz_ptr, double));
264 DLL_IMPORT_RTS void mpz_set_f _PROTO ((mpz_ptr, mpf_srcptr));
265 DLL_IMPORT_RTS void mpz_set_q _PROTO ((mpz_ptr, mpq_srcptr));
266 DLL_IMPORT_RTS void mpz_set_si _PROTO ((mpz_ptr, signed long int));
267 DLL_IMPORT_RTS int mpz_set_str _PROTO ((mpz_ptr, const char *, int));
268 DLL_IMPORT_RTS void mpz_set_ui _PROTO ((mpz_ptr, unsigned long int));
269 DLL_IMPORT_RTS void mpz_setbit _PROTO ((mpz_ptr, unsigned long int));
270 DLL_IMPORT_RTS size_t mpz_size _PROTO ((mpz_srcptr));
271 DLL_IMPORT_RTS size_t mpz_sizeinbase _PROTO ((mpz_srcptr, int));
272 DLL_IMPORT_RTS void mpz_sqrt _PROTO ((mpz_ptr, mpz_srcptr));
273 DLL_IMPORT_RTS void mpz_sqrtrem _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr));
274 DLL_IMPORT_RTS void mpz_sub _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
275 DLL_IMPORT_RTS void mpz_sub_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
276 DLL_IMPORT_RTS void mpz_tdiv_q _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
277 DLL_IMPORT_RTS void mpz_tdiv_q_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
278 DLL_IMPORT_RTS void mpz_tdiv_q_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
279 DLL_IMPORT_RTS void mpz_tdiv_qr _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, mpz_srcptr));
280 DLL_IMPORT_RTS void mpz_tdiv_qr_ui _PROTO ((mpz_ptr, mpz_ptr, mpz_srcptr, unsigned long int));
281 DLL_IMPORT_RTS void mpz_tdiv_r _PROTO ((mpz_ptr, mpz_srcptr, mpz_srcptr));
282 DLL_IMPORT_RTS void mpz_tdiv_r_2exp _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
283 DLL_IMPORT_RTS void mpz_tdiv_r_ui _PROTO ((mpz_ptr, mpz_srcptr, unsigned long int));
284 DLL_IMPORT_RTS void mpz_ui_pow_ui _PROTO ((mpz_ptr, unsigned long int, unsigned long int));
286 /**************** Rational (i.e. Q) routines. ****************/
288 DLL_IMPORT_RTS void mpq_init _PROTO ((mpq_ptr));
289 DLL_IMPORT_RTS void mpq_clear _PROTO ((mpq_ptr));
290 DLL_IMPORT_RTS void mpq_set _PROTO ((mpq_ptr, mpq_srcptr));
291 DLL_IMPORT_RTS void mpq_set_ui _PROTO ((mpq_ptr, unsigned long int, unsigned long int));
292 DLL_IMPORT_RTS void mpq_set_si _PROTO ((mpq_ptr, signed long int, unsigned long int));
293 DLL_IMPORT_RTS void mpq_set_z _PROTO ((mpq_ptr, mpz_srcptr));
294 DLL_IMPORT_RTS void mpq_add _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
295 DLL_IMPORT_RTS void mpq_sub _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
296 DLL_IMPORT_RTS void mpq_mul _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
297 DLL_IMPORT_RTS void mpq_div _PROTO ((mpq_ptr, mpq_srcptr, mpq_srcptr));
298 DLL_IMPORT_RTS void mpq_neg _PROTO ((mpq_ptr, mpq_srcptr));
299 DLL_IMPORT_RTS int mpq_cmp _PROTO ((mpq_srcptr, mpq_srcptr));
300 DLL_IMPORT_RTS int mpq_cmp_ui _PROTO ((mpq_srcptr, unsigned long int, unsigned long int));
301 DLL_IMPORT_RTS int mpq_equal _PROTO ((mpq_srcptr, mpq_srcptr));
302 DLL_IMPORT_RTS void mpq_inv _PROTO ((mpq_ptr, mpq_srcptr));
303 DLL_IMPORT_RTS void mpq_set_num _PROTO ((mpq_ptr, mpz_srcptr));
304 DLL_IMPORT_RTS void mpq_set_den _PROTO ((mpq_ptr, mpz_srcptr));
305 DLL_IMPORT_RTS void mpq_get_num _PROTO ((mpz_ptr, mpq_srcptr));
306 DLL_IMPORT_RTS void mpq_get_den _PROTO ((mpz_ptr, mpq_srcptr));
307 DLL_IMPORT_RTS double mpq_get_d _PROTO ((mpq_srcptr));
308 DLL_IMPORT_RTS void mpq_canonicalize _PROTO ((mpq_ptr));
310 /**************** Float (i.e. F) routines. ****************/
312 DLL_IMPORT_RTS void mpf_abs _PROTO ((mpf_ptr, mpf_srcptr));
313 DLL_IMPORT_RTS void mpf_add _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
314 DLL_IMPORT_RTS void mpf_add_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
315 DLL_IMPORT_RTS void mpf_clear _PROTO ((mpf_ptr));
316 DLL_IMPORT_RTS int mpf_cmp _PROTO ((mpf_srcptr, mpf_srcptr));
317 DLL_IMPORT_RTS int mpf_cmp_si _PROTO ((mpf_srcptr, signed long int));
318 DLL_IMPORT_RTS int mpf_cmp_ui _PROTO ((mpf_srcptr, unsigned long int));
319 DLL_IMPORT_RTS void mpf_div _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
320 DLL_IMPORT_RTS void mpf_div_2exp _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
321 DLL_IMPORT_RTS void mpf_div_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
322 DLL_IMPORT_RTS void mpf_dump _PROTO ((mpf_srcptr));
323 DLL_IMPORT_RTS int mpf_eq _PROTO ((mpf_srcptr, mpf_srcptr, unsigned long int));
324 DLL_IMPORT_RTS double mpf_get_d _PROTO ((mpf_srcptr));
325 DLL_IMPORT_RTS unsigned long int mpf_get_prec _PROTO ((mpf_srcptr));
326 DLL_IMPORT_RTS char *mpf_get_str _PROTO ((char *, mp_exp_t *, int, size_t, mpf_srcptr));
327 DLL_IMPORT_RTS void mpf_init _PROTO ((mpf_ptr));
328 DLL_IMPORT_RTS void mpf_init2 _PROTO ((mpf_ptr, unsigned long int));
329 #ifdef _GMP_H_HAVE_FILE
330 DLL_IMPORT_RTS size_t mpf_inp_str _PROTO ((mpf_ptr, FILE *, int));
332 DLL_IMPORT_RTS void mpf_init_set _PROTO ((mpf_ptr, mpf_srcptr));
333 DLL_IMPORT_RTS void mpf_init_set_d _PROTO ((mpf_ptr, double));
334 DLL_IMPORT_RTS void mpf_init_set_si _PROTO ((mpf_ptr, signed long int));
335 DLL_IMPORT_RTS int mpf_init_set_str _PROTO ((mpf_ptr, const char *, int));
336 DLL_IMPORT_RTS void mpf_init_set_ui _PROTO ((mpf_ptr, unsigned long int));
337 DLL_IMPORT_RTS void mpf_mul _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
338 DLL_IMPORT_RTS void mpf_mul_2exp _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
339 DLL_IMPORT_RTS void mpf_mul_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
340 DLL_IMPORT_RTS void mpf_neg _PROTO ((mpf_ptr, mpf_srcptr));
341 #ifdef _GMP_H_HAVE_FILE
342 DLL_IMPORT_RTS size_t mpf_out_str _PROTO ((FILE *, int, size_t, mpf_srcptr));
344 DLL_IMPORT_RTS void mpf_random2 _PROTO ((mpf_ptr, mp_size_t, mp_exp_t));
345 DLL_IMPORT_RTS void mpf_reldiff _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
346 DLL_IMPORT_RTS void mpf_set _PROTO ((mpf_ptr, mpf_srcptr));
347 DLL_IMPORT_RTS void mpf_set_d _PROTO ((mpf_ptr, double));
348 DLL_IMPORT_RTS void mpf_set_default_prec _PROTO ((unsigned long int));
349 DLL_IMPORT_RTS void mpf_set_prec _PROTO ((mpf_ptr, unsigned long int));
350 DLL_IMPORT_RTS void mpf_set_prec_raw _PROTO ((mpf_ptr, unsigned long int));
351 DLL_IMPORT_RTS void mpf_set_q _PROTO ((mpf_ptr, mpq_srcptr));
352 DLL_IMPORT_RTS void mpf_set_si _PROTO ((mpf_ptr, signed long int));
353 DLL_IMPORT_RTS int mpf_set_str _PROTO ((mpf_ptr, const char *, int));
354 DLL_IMPORT_RTS void mpf_set_ui _PROTO ((mpf_ptr, unsigned long int));
355 DLL_IMPORT_RTS void mpf_set_z _PROTO ((mpf_ptr, mpz_srcptr));
356 DLL_IMPORT_RTS size_t mpf_size _PROTO ((mpf_srcptr));
357 DLL_IMPORT_RTS void mpf_sqrt _PROTO ((mpf_ptr, mpf_srcptr));
358 DLL_IMPORT_RTS void mpf_sqrt_ui _PROTO ((mpf_ptr, unsigned long int));
359 DLL_IMPORT_RTS void mpf_sub _PROTO ((mpf_ptr, mpf_srcptr, mpf_srcptr));
360 DLL_IMPORT_RTS void mpf_sub_ui _PROTO ((mpf_ptr, mpf_srcptr, unsigned long int));
361 DLL_IMPORT_RTS void mpf_ui_div _PROTO ((mpf_ptr, unsigned long int, mpf_srcptr));
362 DLL_IMPORT_RTS void mpf_ui_sub _PROTO ((mpf_ptr, unsigned long int, mpf_srcptr));
363 #if defined (__cplusplus)
366 /************ Low level positive-integer (i.e. N) routines. ************/
368 /* This is ugly, but we need to make usr calls reach the prefixed function. */
369 #define mpn_add __MPN(add)
370 #define mpn_add_1 __MPN(add_1)
371 #define mpn_add_n __MPN(add_n)
372 #define mpn_addmul_1 __MPN(addmul_1)
373 #define mpn_bdivmod __MPN(bdivmod)
374 #define mpn_cmp __MPN(cmp)
375 #define mpn_divmod_1 __MPN(divmod_1)
376 #define mpn_divrem __MPN(divrem)
377 #define mpn_divrem_1 __MPN(divrem_1)
378 #define mpn_dump __MPN(dump)
379 #define mpn_gcd __MPN(gcd)
380 #define mpn_gcd_1 __MPN(gcd_1)
381 #define mpn_gcdext __MPN(gcdext)
382 #define mpn_get_str __MPN(get_str)
383 #define mpn_hamdist __MPN(hamdist)
384 #define mpn_lshift __MPN(lshift)
385 #define mpn_mod_1 __MPN(mod_1)
386 #define mpn_mul __MPN(mul)
387 #define mpn_mul_1 __MPN(mul_1)
388 #define mpn_mul_n __MPN(mul_n)
389 #define mpn_perfect_square_p __MPN(perfect_square_p)
390 #define mpn_popcount __MPN(popcount)
391 #define mpn_preinv_mod_1 __MPN(preinv_mod_1)
392 #define mpn_random2 __MPN(random2)
393 #define mpn_rshift __MPN(rshift)
394 #define mpn_scan0 __MPN(scan0)
395 #define mpn_scan1 __MPN(scan1)
396 #define mpn_set_str __MPN(set_str)
397 #define mpn_sqrtrem __MPN(sqrtrem)
398 #define mpn_sub __MPN(sub)
399 #define mpn_sub_1 __MPN(sub_1)
400 #define mpn_sub_n __MPN(sub_n)
401 #define mpn_submul_1 __MPN(submul_1)
402 #define mpn_udiv_w_sdiv __MPN(udiv_w_sdiv)
404 #if defined (__cplusplus)
407 DLL_IMPORT_RTS mp_limb_t mpn_add _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr,mp_size_t));
408 DLL_IMPORT_RTS mp_limb_t mpn_add_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
409 DLL_IMPORT_RTS mp_limb_t mpn_add_n _PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t));
410 DLL_IMPORT_RTS mp_limb_t mpn_addmul_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
411 DLL_IMPORT_RTS mp_limb_t mpn_bdivmod _PROTO ((mp_ptr, mp_ptr, mp_size_t, mp_srcptr, mp_size_t, unsigned long int));
412 DLL_IMPORT_RTS int mpn_cmp _PROTO ((mp_srcptr, mp_srcptr, mp_size_t));
413 DLL_IMPORT_RTS mp_limb_t mpn_divmod_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
414 DLL_IMPORT_RTS mp_limb_t mpn_divrem _PROTO ((mp_ptr, mp_size_t, mp_ptr, mp_size_t, mp_srcptr, mp_size_t));
415 DLL_IMPORT_RTS mp_limb_t mpn_divrem_1 _PROTO ((mp_ptr, mp_size_t, mp_srcptr, mp_size_t, mp_limb_t));
416 DLL_IMPORT_RTS void mpn_dump _PROTO ((mp_srcptr, mp_size_t));
417 DLL_IMPORT_RTS mp_size_t mpn_gcd _PROTO ((mp_ptr, mp_ptr, mp_size_t, mp_ptr, mp_size_t));
418 DLL_IMPORT_RTS mp_limb_t mpn_gcd_1 _PROTO ((mp_srcptr, mp_size_t, mp_limb_t));
419 DLL_IMPORT_RTS mp_size_t mpn_gcdext _PROTO ((mp_ptr, mp_ptr, mp_ptr, mp_size_t, mp_ptr, mp_size_t));
420 DLL_IMPORT_RTS size_t mpn_get_str _PROTO ((unsigned char *, int, mp_ptr, mp_size_t));
421 DLL_IMPORT_RTS unsigned long int mpn_hamdist _PROTO ((mp_srcptr, mp_srcptr, mp_size_t));
422 DLL_IMPORT_RTS mp_limb_t mpn_lshift _PROTO ((mp_ptr, mp_srcptr, mp_size_t, unsigned int));
423 DLL_IMPORT_RTS mp_limb_t mpn_mod_1 _PROTO ((mp_srcptr, mp_size_t, mp_limb_t));
424 DLL_IMPORT_RTS mp_limb_t mpn_mul _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr, mp_size_t));
425 DLL_IMPORT_RTS mp_limb_t mpn_mul_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
426 DLL_IMPORT_RTS void mpn_mul_n _PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t));
427 DLL_IMPORT_RTS int mpn_perfect_square_p _PROTO ((mp_srcptr, mp_size_t));
428 DLL_IMPORT_RTS unsigned long int mpn_popcount _PROTO ((mp_srcptr, mp_size_t));
429 DLL_IMPORT_RTS mp_limb_t mpn_preinv_mod_1 _PROTO ((mp_srcptr, mp_size_t, mp_limb_t, mp_limb_t));
430 DLL_IMPORT_RTS void mpn_random2 _PROTO ((mp_ptr, mp_size_t));
431 DLL_IMPORT_RTS mp_limb_t mpn_rshift _PROTO ((mp_ptr, mp_srcptr, mp_size_t, unsigned int));
432 DLL_IMPORT_RTS unsigned long int mpn_scan0 _PROTO ((mp_srcptr, unsigned long int));
433 DLL_IMPORT_RTS unsigned long int mpn_scan1 _PROTO ((mp_srcptr, unsigned long int));
434 DLL_IMPORT_RTS mp_size_t mpn_set_str _PROTO ((mp_ptr, const unsigned char *, size_t, int));
435 DLL_IMPORT_RTS mp_size_t mpn_sqrtrem _PROTO ((mp_ptr, mp_ptr, mp_srcptr, mp_size_t));
436 DLL_IMPORT_RTS mp_limb_t mpn_sub _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_srcptr,mp_size_t));
437 DLL_IMPORT_RTS mp_limb_t mpn_sub_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
438 DLL_IMPORT_RTS mp_limb_t mpn_sub_n _PROTO ((mp_ptr, mp_srcptr, mp_srcptr, mp_size_t));
439 DLL_IMPORT_RTS mp_limb_t mpn_submul_1 _PROTO ((mp_ptr, mp_srcptr, mp_size_t, mp_limb_t));
440 #if defined (__cplusplus)
444 #if defined (__GNUC__) || defined (_FORCE_INLINES)
445 _EXTERN_INLINE mp_limb_t
446 #if defined (__STDC__) || defined (__cplusplus)
447 mpn_add_1 (register mp_ptr res_ptr,
448 register mp_srcptr s1_ptr,
449 register mp_size_t s1_size,
450 register mp_limb_t s2_limb)
452 mpn_add_1 (res_ptr, s1_ptr, s1_size, s2_limb)
453 register mp_ptr res_ptr;
454 register mp_srcptr s1_ptr;
455 register mp_size_t s1_size;
456 register mp_limb_t s2_limb;
459 register mp_limb_t x;
462 s2_limb = x + s2_limb;
463 *res_ptr++ = s2_limb;
466 while (--s1_size != 0)
478 if (res_ptr != s1_ptr)
481 for (i = 0; i < s1_size - 1; i++)
482 res_ptr[i] = s1_ptr[i];
487 _EXTERN_INLINE mp_limb_t
488 #if defined (__STDC__) || defined (__cplusplus)
489 mpn_add (register mp_ptr res_ptr,
490 register mp_srcptr s1_ptr,
491 register mp_size_t s1_size,
492 register mp_srcptr s2_ptr,
493 register mp_size_t s2_size)
495 mpn_add (res_ptr, s1_ptr, s1_size, s2_ptr, s2_size)
496 register mp_ptr res_ptr;
497 register mp_srcptr s1_ptr;
498 register mp_size_t s1_size;
499 register mp_srcptr s2_ptr;
500 register mp_size_t s2_size;
503 mp_limb_t cy_limb = 0;
506 cy_limb = mpn_add_n (res_ptr, s1_ptr, s2_ptr, s2_size);
508 if (s1_size - s2_size != 0)
509 cy_limb = mpn_add_1 (res_ptr + s2_size,
516 _EXTERN_INLINE mp_limb_t
517 #if defined (__STDC__) || defined (__cplusplus)
518 mpn_sub_1 (register mp_ptr res_ptr,
519 register mp_srcptr s1_ptr,
520 register mp_size_t s1_size,
521 register mp_limb_t s2_limb)
523 mpn_sub_1 (res_ptr, s1_ptr, s1_size, s2_limb)
524 register mp_ptr res_ptr;
525 register mp_srcptr s1_ptr;
526 register mp_size_t s1_size;
527 register mp_limb_t s2_limb;
530 register mp_limb_t x;
533 s2_limb = x - s2_limb;
534 *res_ptr++ = s2_limb;
537 while (--s1_size != 0)
549 if (res_ptr != s1_ptr)
552 for (i = 0; i < s1_size - 1; i++)
553 res_ptr[i] = s1_ptr[i];
558 _EXTERN_INLINE mp_limb_t
559 #if defined (__STDC__) || defined (__cplusplus)
560 mpn_sub (register mp_ptr res_ptr,
561 register mp_srcptr s1_ptr,
562 register mp_size_t s1_size,
563 register mp_srcptr s2_ptr,
564 register mp_size_t s2_size)
566 mpn_sub (res_ptr, s1_ptr, s1_size, s2_ptr, s2_size)
567 register mp_ptr res_ptr;
568 register mp_srcptr s1_ptr;
569 register mp_size_t s1_size;
570 register mp_srcptr s2_ptr;
571 register mp_size_t s2_size;
574 mp_limb_t cy_limb = 0;
577 cy_limb = mpn_sub_n (res_ptr, s1_ptr, s2_ptr, s2_size);
579 if (s1_size - s2_size != 0)
580 cy_limb = mpn_sub_1 (res_ptr + s2_size,
586 #endif /* __GNUC__ */
588 /* Allow faster testing for negative, zero, and positive. */
589 #define mpz_sgn(Z) ((Z)->_mp_size < 0 ? -1 : (Z)->_mp_size > 0)
590 #define mpf_sgn(F) ((F)->_mp_size < 0 ? -1 : (F)->_mp_size > 0)
591 #define mpq_sgn(Q) ((Q)->_mp_num._mp_size < 0 ? -1 : (Q)->_mp_num._mp_size > 0)
593 /* Allow direct user access to numerator and denominator of a mpq_t object. */
594 #define mpq_numref(Q) (&((Q)->_mp_num))
595 #define mpq_denref(Q) (&((Q)->_mp_den))
597 /* When using GCC, optimize certain common comparisons. */
598 #if defined (__GNUC__)
599 #define mpz_cmp_ui(Z,UI) \
600 (__builtin_constant_p (UI) && (UI) == 0 \
601 ? mpz_sgn (Z) : mpz_cmp_ui (Z,UI))
602 #define mpz_cmp_si(Z,UI) \
603 (__builtin_constant_p (UI) && (UI) == 0 ? mpz_sgn (Z) \
604 : __builtin_constant_p (UI) && (UI) > 0 ? mpz_cmp_ui (Z,UI) \
606 #define mpq_cmp_ui(Q,NUI,DUI) \
607 (__builtin_constant_p (NUI) && (NUI) == 0 \
608 ? mpq_sgn (Q) : mpq_cmp_ui (Q,NUI,DUI))
611 #define mpn_divmod(qp,np,nsize,dp,dsize) mpn_divrem (qp,0,np,nsize,dp,dsize)
613 #define mpn_divmod_1(qp,np,nsize,dlimb) mpn_divrem_1 (qp,0,np,nsize,dlimb)
616 /* Compatibility with GMP 1. */
617 #define mpz_mdiv mpz_fdiv_q
618 #define mpz_mdivmod mpz_fdiv_qr
619 #define mpz_mmod mpz_fdiv_r
620 #define mpz_mdiv_ui mpz_fdiv_q_ui
621 #define mpz_mdivmod_ui(q,r,n,d) \
622 ((r == 0) ? mpz_fdiv_q_ui (q,n,d) : mpz_fdiv_qr_ui (q,r,n,d))
623 #define mpz_mmod_ui(r,n,d) \
624 ((r == 0) ? mpz_fdiv_ui (n,d) : mpz_fdiv_r_ui (r,n,d))
626 /* Useful synonyms, but not quite compatible with GMP 1. */
627 #define mpz_div mpz_fdiv_q
628 #define mpz_divmod mpz_fdiv_qr
629 #define mpz_div_ui mpz_fdiv_q_ui
630 #define mpz_divmod_ui mpz_fdiv_qr_ui
631 #define mpz_mod_ui mpz_fdiv_r_ui
632 #define mpz_div_2exp mpz_fdiv_q_2exp
633 #define mpz_mod_2exp mpz_fdiv_r_2exp
635 #define __GNU_MP_VERSION 2
636 #define __GNU_MP_VERSION_MINOR 0
638 #endif /* __GMP_H__ */