1/*
2 * Copyright 2008-2009 Katholieke Universiteit Leuven
3 *
4 * Use of this software is governed by the MIT license
5 *
6 * Written by Sven Verdoolaege, K.U.Leuven, Departement
7 * Computerwetenschappen, Celestijnenlaan 200A, B-3001 Leuven, Belgium
8 */
9
10#ifndef ISL_INT_H
11#define ISL_INT_H
12
13#include <isl/hash.h>
14#include <string.h>
15#include <gmp.h>
16#include <isl_config.h>
17
18#ifndef mp_get_memory_functions
19void mp_get_memory_functions(
20		void *(**alloc_func_ptr) (size_t),
21		void *(**realloc_func_ptr) (void *, size_t, size_t),
22		void (**free_func_ptr) (void *, size_t));
23#endif
24
25/* isl_int is the basic integer type.  It currently always corresponds
26 * to a gmp mpz_t, but in the future, different types such as long long
27 * or cln::cl_I will be supported.
28 */
29typedef mpz_t	isl_int;
30
31#define isl_int_init(i)		mpz_init(i)
32#define isl_int_clear(i)	mpz_clear(i)
33
34#define isl_int_set(r,i)	mpz_set(r,i)
35#define isl_int_set_gmp(r,i)	mpz_set(r,i)
36#define isl_int_set_si(r,i)	mpz_set_si(r,i)
37#define isl_int_set_ui(r,i)	mpz_set_ui(r,i)
38#define isl_int_get_gmp(i,g)	mpz_set(g,i)
39#define isl_int_fits_slong(r)	mpz_fits_slong_p(r)
40#define isl_int_get_si(r)	mpz_get_si(r)
41#define isl_int_fits_ulong(r)	mpz_fits_ulong_p(r)
42#define isl_int_get_ui(r)	mpz_get_ui(r)
43#define isl_int_get_d(r)	mpz_get_d(r)
44#define isl_int_get_str(r)	mpz_get_str(0, 10, r)
45typedef void (*isl_int_print_gmp_free_t)(void *, size_t);
46#define isl_int_free_str(s)					\
47	do {								\
48		isl_int_print_gmp_free_t gmp_free;			\
49		mp_get_memory_functions(NULL, NULL, &gmp_free);		\
50		(*gmp_free)(s, strlen(s) + 1);				\
51	} while (0)
52#define isl_int_abs(r,i)	mpz_abs(r,i)
53#define isl_int_neg(r,i)	mpz_neg(r,i)
54#define isl_int_swap(i,j)	mpz_swap(i,j)
55#define isl_int_swap_or_set(i,j)	mpz_swap(i,j)
56#define isl_int_add_ui(r,i,j)	mpz_add_ui(r,i,j)
57#define isl_int_sub_ui(r,i,j)	mpz_sub_ui(r,i,j)
58
59#define isl_int_add(r,i,j)	mpz_add(r,i,j)
60#define isl_int_sub(r,i,j)	mpz_sub(r,i,j)
61#define isl_int_mul(r,i,j)	mpz_mul(r,i,j)
62#define isl_int_mul_2exp(r,i,j)	mpz_mul_2exp(r,i,j)
63#define isl_int_mul_si(r,i,j)	mpz_mul_si(r,i,j)
64#define isl_int_mul_ui(r,i,j)	mpz_mul_ui(r,i,j)
65#define isl_int_pow_ui(r,i,j)	mpz_pow_ui(r,i,j)
66#define isl_int_addmul(r,i,j)	mpz_addmul(r,i,j)
67#define isl_int_addmul_ui(r,i,j)	mpz_addmul_ui(r,i,j)
68#define isl_int_submul(r,i,j)	mpz_submul(r,i,j)
69#define isl_int_submul_ui(r,i,j)	mpz_submul_ui(r,i,j)
70
71#define isl_int_gcd(r,i,j)	mpz_gcd(r,i,j)
72#ifdef GMP_NORMALIZE_GCDEXT
73void isl_gmp_gcdext(mpz_t G, mpz_t S, mpz_t T, mpz_t A, mpz_t B);
74#define isl_int_gcdext(g,x,y,i,j)	isl_gmp_gcdext(g,x,y,i,j)
75#else
76#define isl_int_gcdext(g,x,y,i,j)	mpz_gcdext(g,x,y,i,j)
77#endif
78#define isl_int_lcm(r,i,j)	mpz_lcm(r,i,j)
79#define isl_int_divexact(r,i,j)	mpz_divexact(r,i,j)
80#define isl_int_divexact_ui(r,i,j)	mpz_divexact_ui(r,i,j)
81#define isl_int_tdiv_q(r,i,j)	mpz_tdiv_q(r,i,j)
82#define isl_int_cdiv_q(r,i,j)	mpz_cdiv_q(r,i,j)
83#define isl_int_fdiv_q(r,i,j)	mpz_fdiv_q(r,i,j)
84#define isl_int_fdiv_r(r,i,j)	mpz_fdiv_r(r,i,j)
85#define isl_int_fdiv_q_ui(r,i,j)	mpz_fdiv_q_ui(r,i,j)
86
87#define isl_int_read(r,s)	mpz_set_str(r,s,10)
88#define isl_int_print(out,i,width)					\
89	do {								\
90		char *s;						\
91		s = mpz_get_str(0, 10, i);				\
92		fprintf(out, "%*s", width, s);				\
93		isl_int_free_str(s);                                        \
94	} while (0)
95
96#define isl_int_sgn(i)		mpz_sgn(i)
97#define isl_int_cmp(i,j)	mpz_cmp(i,j)
98#define isl_int_cmp_si(i,si)	mpz_cmp_si(i,si)
99#define isl_int_eq(i,j)		(mpz_cmp(i,j) == 0)
100#define isl_int_ne(i,j)		(mpz_cmp(i,j) != 0)
101#define isl_int_lt(i,j)		(mpz_cmp(i,j) < 0)
102#define isl_int_le(i,j)		(mpz_cmp(i,j) <= 0)
103#define isl_int_gt(i,j)		(mpz_cmp(i,j) > 0)
104#define isl_int_ge(i,j)		(mpz_cmp(i,j) >= 0)
105#define isl_int_abs_eq(i,j)	(mpz_cmpabs(i,j) == 0)
106#define isl_int_abs_ne(i,j)	(mpz_cmpabs(i,j) != 0)
107#define isl_int_abs_lt(i,j)	(mpz_cmpabs(i,j) < 0)
108#define isl_int_abs_gt(i,j)	(mpz_cmpabs(i,j) > 0)
109#define isl_int_abs_ge(i,j)	(mpz_cmpabs(i,j) >= 0)
110
111
112#define isl_int_is_zero(i)	(isl_int_sgn(i) == 0)
113#define isl_int_is_one(i)	(isl_int_cmp_si(i,1) == 0)
114#define isl_int_is_negone(i)	(isl_int_cmp_si(i,-1) == 0)
115#define isl_int_is_pos(i)	(isl_int_sgn(i) > 0)
116#define isl_int_is_neg(i)	(isl_int_sgn(i) < 0)
117#define isl_int_is_nonpos(i)	(isl_int_sgn(i) <= 0)
118#define isl_int_is_nonneg(i)	(isl_int_sgn(i) >= 0)
119#define isl_int_is_divisible_by(i,j)	mpz_divisible_p(i,j)
120
121uint32_t isl_gmp_hash(mpz_t v, uint32_t hash);
122#define isl_int_hash(v,h)	isl_gmp_hash(v,h)
123
124#endif
125