int main(int argc, char **argv)
{
struct isl_ctx *ctx = isl_ctx_alloc();
struct isl_basic_set *bset;
struct isl_vec *sample;
isl_printer *p;
bset = isl_basic_set_read_from_file(ctx, stdin);
sample = isl_basic_set_sample_vec(isl_basic_set_copy(bset));
p = isl_printer_to_file(ctx, stdout);
p = isl_printer_print_vec(p, sample);
p = isl_printer_end_line(p);
isl_printer_free(p);
assert(sample);
if (sample->size > 0)
assert(isl_basic_set_contains(bset, sample));
isl_basic_set_free(bset);
isl_vec_free(sample);
isl_ctx_free(ctx);
return 0;
}
/* Look for all equalities satisfied by the integer points in bset,
* which is assumed to be bounded.
*
* The equalities are obtained by successively looking for
* a point that is affinely independent of the points found so far.
* In particular, for each equality satisfied by the points so far,
* we check if there is any point on a hyperplane parallel to the
* corresponding hyperplane shifted by at least one (in either direction).
*/
static struct isl_basic_set *uset_affine_hull_bounded(struct isl_basic_set *bset)
{
struct isl_vec *sample = NULL;
struct isl_basic_set *hull;
struct isl_tab *tab = NULL;
unsigned dim;
if (isl_basic_set_plain_is_empty(bset))
return bset;
dim = isl_basic_set_n_dim(bset);
if (bset->sample && bset->sample->size == 1 + dim) {
int contains = isl_basic_set_contains(bset, bset->sample);
if (contains < 0)
goto error;
if (contains) {
if (dim == 0)
return bset;
sample = isl_vec_copy(bset->sample);
} else {
isl_vec_free(bset->sample);
bset->sample = NULL;
}
}
tab = isl_tab_from_basic_set(bset);
if (!tab)
goto error;
if (tab->empty) {
isl_tab_free(tab);
isl_vec_free(sample);
return isl_basic_set_set_to_empty(bset);
}
if (isl_tab_track_bset(tab, isl_basic_set_copy(bset)) < 0)
goto error;
if (!sample) {
struct isl_tab_undo *snap;
snap = isl_tab_snap(tab);
sample = isl_tab_sample(tab);
if (isl_tab_rollback(tab, snap) < 0)
goto error;
isl_vec_free(tab->bmap->sample);
tab->bmap->sample = isl_vec_copy(sample);
}
if (!sample)
goto error;
if (sample->size == 0) {
isl_tab_free(tab);
isl_vec_free(sample);
return isl_basic_set_set_to_empty(bset);
}
hull = isl_basic_set_from_vec(sample);
isl_basic_set_free(bset);
hull = extend_affine_hull(tab, hull);
isl_tab_free(tab);
return hull;
error:
isl_vec_free(sample);
isl_tab_free(tab);
isl_basic_set_free(bset);
return NULL;
}
