1#ifndef CLOOG_INT_H
2#define CLOOG_INT_H
3
4#undef CLOOG_INT_INT
5#undef CLOOG_INT_LONG
6#undef CLOOG_INT_LONG_LONG
7#undef CLOOG_INT_GMP
8
9#define @CLOOG_INT_CONFIGURE@
10
11#include <assert.h>
12#include <stdio.h>
13#if defined(CLOOG_INT_GMP)
14#include <string.h>
15#include <gmp.h>
16#ifndef mp_get_memory_functions
17void mp_get_memory_functions(
18		void *(**alloc_func_ptr) (size_t),
19		void *(**realloc_func_ptr) (void *, size_t, size_t),
20		void (**free_func_ptr) (void *, size_t));
21#endif
22#endif
23
24#if defined(__cplusplus)
25extern "C"
26  {
27#endif
28
29#if defined(CLOOG_INT_INT)
30typedef int		cloog_int_t;
31#define CLOOG_INT_FORMAT	"%d"
32#elif defined(CLOOG_INT_LONG)
33typedef long		cloog_int_t;
34#define CLOOG_INT_FORMAT	"%ld"
35#elif defined(CLOOG_INT_LONG_LONG)
36typedef long long	cloog_int_t;
37#define CLOOG_INT_FORMAT	"%lld"
38#elif defined(CLOOG_INT_GMP)
39typedef mpz_t		cloog_int_t;
40#else
41#error "No integer type defined"
42#endif
43
44#if defined(CLOOG_INT_GMP)
45
46#define cloog_int_init(i)	mpz_init(i)
47#define cloog_int_clear(i)	mpz_clear(i);
48
49#define cloog_int_set(r,i)	mpz_set(r,i)
50#define cloog_int_set_si(r,i)	mpz_set_si(r,i)
51#define cloog_int_abs(r,i)	mpz_abs(r,i)
52#define cloog_int_neg(r,i)	mpz_neg(r,i)
53#define cloog_int_swap(i,j)	mpz_swap(i,j)
54
55#define cloog_int_add(r,i,j)	mpz_add(r,i,j)
56#define cloog_int_add_ui(r,i,j)	mpz_add_ui(r,i,j)
57#define cloog_int_sub(r,i,j)	mpz_sub(r,i,j)
58#define cloog_int_sub_ui(r,i,j)	mpz_sub_ui(r,i,j)
59#define cloog_int_mul(r,i,j)	mpz_mul(r,i,j)
60#define cloog_int_addmul(r,i,j)	mpz_addmul(r,i,j)
61#define cloog_int_divexact(r,i,j)	mpz_divexact(r,i,j)
62#define cloog_int_tdiv_q(r,i,j)	mpz_tdiv_q(r,i,j)
63#define cloog_int_fdiv_q(r,i,j)	mpz_fdiv_q(r,i,j)
64#define cloog_int_fdiv_r(r,i,j)	mpz_fdiv_r(r,i,j)
65#define cloog_int_cdiv_q(r,i,j)	mpz_cdiv_q(r,i,j)
66#define cloog_int_gcd(r,i,j)	mpz_gcd(r,i,j)
67
68#define cloog_int_sgn(i)	mpz_sgn(i)
69#define cloog_int_cmp(i,j)	mpz_cmp(i,j)
70#define cloog_int_abs_cmp(i,j)	mpz_cmpabs(i,j)
71#define cloog_int_cmp_si(i,si)	mpz_cmp_si(i,si)
72#define cloog_int_eq(i,j)	(mpz_cmp(i,j) == 0)
73#define cloog_int_ne(i,j)	(mpz_cmp(i,j) != 0)
74#define cloog_int_gt(i,j)	(mpz_cmp(i,j) > 0)
75
76#define cloog_int_is_divisible_by(i,j)	mpz_divisible_p(i,j)
77
78#define cloog_int_read(r,s)	mpz_set_str(r,s,10)
79typedef void (*cloog_int_print_gmp_free_t)(void *, size_t);
80#define cloog_int_print(out,i)						\
81	do {								\
82		char *s;						\
83		cloog_int_print_gmp_free_t gmp_free;			\
84		s = mpz_get_str(0, 10, i);				\
85		fprintf(out, "%s", s);					\
86		mp_get_memory_functions(NULL, NULL, &gmp_free);		\
87		(*gmp_free)(s, strlen(s)+1);				\
88	} while (0)
89
90#else
91
92#define cloog_int_init(i)	((i) = 0)
93#define cloog_int_clear(i)	do { } while (0)
94
95#define cloog_int_set(r,i)	((r) = (i))
96#define cloog_int_set_si(r,i)	((r) = (i))
97#define cloog_int_abs(r,i)	((r) = (i) > 0 ? (i) : -(i))
98#define cloog_int_neg(r,i)	((r) = -(i))
99#define cloog_int_swap(i,j)	do {					\
100					cloog_int_t _t = i;		\
101					i = j;				\
102					j = _t;				\
103				} while (0)
104
105#define cloog_int_add(r,i,j)	((r) = (i) + (j))
106#define cloog_int_add_ui(r,i,j)	((r) = (i) + (j))
107#define cloog_int_sub(r,i,j)	((r) = (i) - (j))
108#define cloog_int_sub_ui(r,i,j)	((r) = (i) - (j))
109#define cloog_int_mul(r,i,j)	((r) = (i) * (j))
110#define cloog_int_addmul(r,i,j)	((r) += (i) * (j))
111#define cloog_int_divexact(r,i,j)	((r) = (i) / (j))
112#define cloog_int_tdiv_q(r,i,j)	((r) = (i) / (j))
113#define cloog_int_fdiv_q(r,i,j)	do {					\
114					assert((j) > 0);		\
115					(r) = (i) >= 0 ? (i) / (j) :	\
116					    -((-(i)+(j)-1)/(j));	\
117				} while (0)
118#define cloog_int_fdiv_r(r,i,j)	do {					\
119					assert((j) > 0);		\
120					(r) = (i) >= 0 ? (i) % (j) :	\
121					    (j)-1 - ((-(i)+(j)-1)%(j));	\
122				} while (0)
123#define cloog_int_cdiv_q(r,i,j)	do {					\
124					assert((j) > 0);		\
125					(r) = (i) >= 0 ?		\
126					    ((i)+(j)-1) / (j) :		\
127					    -(-(i)/(j));		\
128				} while (0)
129cloog_int_t cloog_gcd(cloog_int_t a, cloog_int_t b);
130#define cloog_int_gcd(r,i,j)	(r) = cloog_gcd(i,j)
131
132#define cloog_int_sgn(i)	((i) > 0 ? 1 : (i) < 0 ? -1 : 0)
133#define cloog_int_cmp(i,j)	(i - j)
134#define cloog_int_abs_cmp(i,j)	(((i) > 0 ? (i) : -(i)) - ((j) > 0 ? (j) : -(j)))
135#define cloog_int_cmp_si(i,si)	(i - si)
136#define cloog_int_eq(i,j)	((i) == (j))
137#define cloog_int_ne(i,j)	((i) != (j))
138#define cloog_int_gt(i,j)	((i) > (j))
139
140#define cloog_int_is_divisible_by(i,j)	((i) % (j) == 0)
141
142#define cloog_int_read(i,s)	sscanf(s, CLOOG_INT_FORMAT, &i)
143#define cloog_int_print(out,i)	fprintf(out, CLOOG_INT_FORMAT, i)
144
145#endif
146
147#define cloog_int_is_pos(i)	(cloog_int_sgn(i) > 0)
148#define cloog_int_is_neg(i)	(cloog_int_sgn(i) < 0)
149#define cloog_int_is_zero(i)	(cloog_int_sgn(i) == 0)
150#define cloog_int_is_one(i)	(cloog_int_cmp_si(i,1) == 0)
151#define cloog_int_is_neg_one(i)	(cloog_int_cmp_si(i,-1) == 0)
152#define cloog_int_gt_si(i,si)	(cloog_int_cmp_si(i,si) > 0)
153#define cloog_int_ne_si(i,j)	(cloog_int_cmp_si(i,j) != 0)
154#define cloog_int_lt(i,j)	(cloog_int_cmp(i,j) < 0)
155#define cloog_int_le(i,j)	(cloog_int_cmp(i,j) <= 0)
156#define cloog_int_abs_ne(i,j)	(cloog_int_abs_cmp(i,j) != 0)
157#define cloog_int_abs_ge(i,j)	(cloog_int_abs_cmp(i,j) >= 0)
158#define cloog_int_abs_lt(i,j)	(cloog_int_abs_cmp(i,j) < 0)
159
160struct cloog_vec {
161	unsigned	size;
162	cloog_int_t	*p;
163};
164
165struct cloog_vec *cloog_vec_alloc(unsigned size);
166void cloog_vec_free(struct cloog_vec *vec);
167
168int cloog_seq_first_non_zero(cloog_int_t *p, unsigned len);
169void cloog_seq_cpy(cloog_int_t *dst, cloog_int_t *src, unsigned len);
170void cloog_seq_neg(cloog_int_t *dst, cloog_int_t *src, unsigned len);
171void cloog_seq_combine(cloog_int_t *dst, cloog_int_t m1, cloog_int_t *src1,
172			cloog_int_t m2, cloog_int_t *src2, unsigned len);
173void cloog_seq_gcd(cloog_int_t *p, unsigned len, cloog_int_t *gcd);
174int cloog_seq_is_neg(cloog_int_t *p1, cloog_int_t *p2, unsigned len);
175void cloog_seq_normalize(cloog_int_t *p, unsigned len);
176
177#if defined(__cplusplus)
178  }
179#endif
180
181#endif
182