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);