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