2 #include "isl_equalities.h"
6 /* The input of this program is the same as that of the "polytope_scan"
7 * program from the barvinok distribution.
9 * Constraints of set is PolyLib format.
11 * The input set is assumed to be bounded.
15 struct isl_scan_callback callback
;
16 struct isl_mat
*samples
;
19 static int scan_samples_add_sample(struct isl_scan_callback
*cb
,
20 __isl_take isl_vec
*sample
)
22 struct scan_samples
*ss
= (struct scan_samples
*)cb
;
24 ss
->samples
= isl_mat_extend(ss
->samples
, ss
->samples
->n_row
+ 1,
29 isl_seq_cpy(ss
->samples
->row
[ss
->samples
->n_row
- 1],
30 sample
->el
, sample
->size
);
39 static struct isl_mat
*isl_basic_set_scan_samples(struct isl_basic_set
*bset
)
42 struct scan_samples ss
;
44 dim
= isl_basic_set_total_dim(bset
);
45 ss
.callback
.add
= scan_samples_add_sample
;
46 ss
.samples
= isl_mat_alloc(bset
->ctx
, 0, 1 + dim
);
50 if (isl_basic_set_scan(bset
, &ss
.callback
) < 0) {
51 isl_mat_free(ss
.samples
);
57 isl_basic_set_free(bset
);
61 static struct isl_mat
*isl_basic_set_samples(struct isl_basic_set
*bset
)
64 struct isl_mat
*samples
;
67 return isl_basic_set_scan_samples(bset
);
69 bset
= isl_basic_set_remove_equalities(bset
, &T
, NULL
);
70 samples
= isl_basic_set_scan_samples(bset
);
71 return isl_mat_product(samples
, isl_mat_transpose(T
));
74 int main(int argc
, char **argv
)
76 struct isl_ctx
*ctx
= isl_ctx_alloc();
77 struct isl_basic_set
*bset
;
78 struct isl_mat
*samples
;
80 bset
= isl_basic_set_read_from_file(ctx
, stdin
, 0, ISL_FORMAT_POLYLIB
);
81 samples
= isl_basic_set_samples(bset
);
82 isl_mat_dump(samples
, stdout
, 0);
83 isl_mat_free(samples
);