1 #include <isl_point_private.h>
3 #include <isl_map_private.h>
4 #include <isl_sample.h>
8 __isl_give isl_point *isl_point_alloc(__isl_take isl_dim *dim,
9 __isl_take isl_vec *vec)
11 struct isl_point *pnt;
16 pnt = isl_alloc_type(dim->ctx, struct isl_point);
20 if (vec->size > 1 + isl_dim_total(dim)) {
21 vec = isl_vec_cow(vec);
24 vec->size = 1 + isl_dim_total(dim);
38 __isl_give isl_point *isl_point_zero(__isl_take isl_dim *dim)
44 vec = isl_vec_alloc(dim->ctx, 1 + isl_dim_total(dim));
47 isl_int_set_si(vec->el[0], 1);
48 isl_seq_clr(vec->el + 1, vec->size - 1);
49 return isl_point_alloc(dim, vec);
55 __isl_give isl_point *isl_point_dup(__isl_keep isl_point *pnt)
57 struct isl_point *pnt2;
61 pnt2 = isl_point_alloc(isl_dim_copy(pnt->dim), isl_vec_copy(pnt->vec));
65 __isl_give isl_point *isl_point_cow(__isl_take isl_point *pnt)
67 struct isl_point *pnt2;
74 pnt2 = isl_point_dup(pnt);
79 __isl_give isl_point *isl_point_copy(__isl_keep isl_point *pnt)
88 void isl_point_free(__isl_take isl_point *pnt)
96 isl_dim_free(pnt->dim);
97 isl_vec_free(pnt->vec);
101 __isl_give isl_point *isl_point_void(__isl_take isl_dim *dim)
106 return isl_point_alloc(dim, isl_vec_alloc(dim->ctx, 0));
109 int isl_point_is_void(__isl_keep isl_point *pnt)
114 return pnt->vec->size == 0;
117 void isl_point_get_coordinate(__isl_keep isl_point *pnt,
118 enum isl_dim_type type, int pos, isl_int *v)
120 if (!pnt || isl_point_is_void(pnt))
122 if (type == isl_dim_set)
123 pos += isl_dim_size(pnt->dim, isl_dim_param);
124 isl_int_set(*v, pnt->vec->el[1 + pos]);
127 __isl_give isl_point *isl_point_set_coordinate(__isl_take isl_point *pnt,
128 enum isl_dim_type type, int pos, isl_int v)
130 if (!pnt || isl_point_is_void(pnt))
133 pnt = isl_point_cow(pnt);
136 pnt->vec = isl_vec_cow(pnt->vec);
140 if (type == isl_dim_set)
141 pos += isl_dim_size(pnt->dim, isl_dim_param);
143 isl_int_set(pnt->vec->el[1 + pos], v);
151 __isl_give isl_point *isl_point_add_ui(__isl_take isl_point *pnt,
152 enum isl_dim_type type, int pos, unsigned val)
154 if (!pnt || isl_point_is_void(pnt))
157 pnt = isl_point_cow(pnt);
160 pnt->vec = isl_vec_cow(pnt->vec);
164 if (type == isl_dim_set)
165 pos += isl_dim_size(pnt->dim, isl_dim_param);
167 isl_int_add_ui(pnt->vec->el[1 + pos], pnt->vec->el[1 + pos], val);
175 __isl_give isl_point *isl_point_sub_ui(__isl_take isl_point *pnt,
176 enum isl_dim_type type, int pos, unsigned val)
178 if (!pnt || isl_point_is_void(pnt))
181 pnt = isl_point_cow(pnt);
184 pnt->vec = isl_vec_cow(pnt->vec);
188 if (type == isl_dim_set)
189 pos += isl_dim_size(pnt->dim, isl_dim_param);
191 isl_int_sub_ui(pnt->vec->el[1 + pos], pnt->vec->el[1 + pos], val);
199 struct isl_foreach_point {
200 struct isl_scan_callback callback;
201 int (*fn)(__isl_take isl_point *pnt, void *user);
206 static int foreach_point(struct isl_scan_callback *cb, __isl_take isl_vec *sample)
208 struct isl_foreach_point *fp = (struct isl_foreach_point *)cb;
211 pnt = isl_point_alloc(isl_dim_copy(fp->dim), sample);
213 return fp->fn(pnt, fp->user);
216 int isl_set_foreach_point(__isl_keep isl_set *set,
217 int (*fn)(__isl_take isl_point *pnt, void *user), void *user)
219 struct isl_foreach_point fp = { { &foreach_point }, fn, user };
225 fp.dim = isl_set_get_dim(set);
229 set = isl_set_copy(set);
230 set = isl_set_cow(set);
231 set = isl_set_make_disjoint(set);
232 set = isl_set_compute_divs(set);
236 for (i = 0; i < set->n; ++i)
237 if (isl_basic_set_scan(isl_basic_set_copy(set->p[i]),
242 isl_dim_free(fp.dim);
247 isl_dim_free(fp.dim);
251 __isl_give isl_set *isl_set_box_from_points(__isl_take isl_point *pnt1,
252 __isl_take isl_point *pnt2)
265 isl_assert(pnt1->dim->ctx,
266 isl_dim_equal(pnt1->dim, pnt2->dim), goto error);
268 if (isl_point_is_void(pnt1) && isl_point_is_void(pnt2)) {
269 isl_dim *dim = isl_dim_copy(pnt1->dim);
270 isl_point_free(pnt1);
271 isl_point_free(pnt2);
273 return isl_set_empty(dim);
275 if (isl_point_is_void(pnt1)) {
276 isl_basic_set *model;
277 model = isl_basic_set_empty(isl_dim_copy(pnt2->dim));
278 bset = isl_basic_set_from_vec(isl_vec_copy(pnt2->vec));
279 bset = isl_basic_set_from_underlying_set(bset, model);
280 isl_point_free(pnt1);
281 isl_point_free(pnt2);
283 return isl_set_from_basic_set(bset);
285 if (isl_point_is_void(pnt2)) {
286 isl_basic_set *model;
287 model = isl_basic_set_empty(isl_dim_copy(pnt1->dim));
288 bset = isl_basic_set_from_vec(isl_vec_copy(pnt1->vec));
289 bset = isl_basic_set_from_underlying_set(bset, model);
290 isl_point_free(pnt1);
291 isl_point_free(pnt2);
293 return isl_set_from_basic_set(bset);
296 total = isl_dim_total(pnt1->dim);
297 bset = isl_basic_set_alloc_dim(isl_dim_copy(pnt1->dim), 0, 0, 2 * total);
299 for (i = 0; i < total; ++i) {
300 isl_int_mul(t, pnt1->vec->el[1 + i], pnt2->vec->el[0]);
301 isl_int_submul(t, pnt2->vec->el[1 + i], pnt1->vec->el[0]);
303 k = isl_basic_set_alloc_inequality(bset);
306 isl_seq_clr(bset->ineq[k] + 1, total);
307 if (isl_int_is_pos(t)) {
308 isl_int_set_si(bset->ineq[k][1 + i], -1);
309 isl_int_set(bset->ineq[k][0], pnt1->vec->el[1 + i]);
311 isl_int_set_si(bset->ineq[k][1 + i], 1);
312 isl_int_neg(bset->ineq[k][0], pnt1->vec->el[1 + i]);
314 isl_int_fdiv_q(bset->ineq[k][0], bset->ineq[k][0], pnt1->vec->el[0]);
316 k = isl_basic_set_alloc_inequality(bset);
319 isl_seq_clr(bset->ineq[k] + 1, total);
320 if (isl_int_is_pos(t)) {
321 isl_int_set_si(bset->ineq[k][1 + i], 1);
322 isl_int_neg(bset->ineq[k][0], pnt2->vec->el[1 + i]);
324 isl_int_set_si(bset->ineq[k][1 + i], -1);
325 isl_int_set(bset->ineq[k][0], pnt2->vec->el[1 + i]);
327 isl_int_fdiv_q(bset->ineq[k][0], bset->ineq[k][0], pnt2->vec->el[0]);
330 bset = isl_basic_set_finalize(bset);
332 isl_point_free(pnt1);
333 isl_point_free(pnt2);
337 return isl_set_from_basic_set(bset);
339 isl_point_free(pnt1);
340 isl_point_free(pnt2);
345 void isl_point_print(__isl_keep isl_point *pnt, FILE *out)
353 if (isl_point_is_void(pnt)) {
354 fprintf(out, "void\n");
358 nparam = isl_dim_size(pnt->dim, isl_dim_param);
359 dim = isl_dim_size(pnt->dim, isl_dim_set);
362 for (i = 0; i < nparam; ++i) {
366 name = isl_dim_get_name(pnt->dim, isl_dim_param, i);
368 fprintf(out, "%s = ", name);
369 isl_int_print(out, pnt->vec->el[1 + i], 0);
370 if (!isl_int_is_one(pnt->vec->el[0])) {
372 isl_int_print(out, pnt->vec->el[0], 0);
375 fprintf(out, "] -> ");
378 for (i = 0; i < dim; ++i) {
381 isl_int_print(out, pnt->vec->el[1 + nparam + i], 0);
382 if (!isl_int_is_one(pnt->vec->el[0])) {
384 isl_int_print(out, pnt->vec->el[0], 0);