isl_basic_map_contains_point: use isl_local_extend_point_vec
[isl.git] / polytope_scan.c
blobfa53cfa65c7a2f00c8e639fbbde63a3fa5ad06d3
1 /*
2 * Copyright 2008-2009 Katholieke Universiteit Leuven
4 * Use of this software is governed by the MIT license
6 * Written by Sven Verdoolaege, K.U.Leuven, Departement
7 * Computerwetenschappen, Celestijnenlaan 200A, B-3001 Leuven, Belgium
8 */
10 #include <assert.h>
11 #include <isl_map_private.h>
12 #include "isl_equalities.h"
13 #include <isl_seq.h>
14 #include "isl_scan.h"
15 #include <isl_mat_private.h>
16 #include <isl_vec_private.h>
18 /* The input of this program is the same as that of the "polytope_scan"
19 * program from the barvinok distribution.
21 * Constraints of set is PolyLib format.
23 * The input set is assumed to be bounded.
26 struct scan_samples {
27 struct isl_scan_callback callback;
28 struct isl_mat *samples;
31 static isl_stat scan_samples_add_sample(struct isl_scan_callback *cb,
32 __isl_take isl_vec *sample)
34 struct scan_samples *ss = (struct scan_samples *)cb;
36 ss->samples = isl_mat_extend(ss->samples, ss->samples->n_row + 1,
37 ss->samples->n_col);
38 if (!ss->samples)
39 goto error;
41 isl_seq_cpy(ss->samples->row[ss->samples->n_row - 1],
42 sample->el, sample->size);
44 isl_vec_free(sample);
45 return isl_stat_ok;
46 error:
47 isl_vec_free(sample);
48 return isl_stat_error;
51 static __isl_give isl_mat *isl_basic_set_scan_samples(
52 __isl_take isl_basic_set *bset)
54 isl_ctx *ctx;
55 unsigned dim;
56 struct scan_samples ss;
58 ctx = isl_basic_set_get_ctx(bset);
59 dim = isl_basic_set_total_dim(bset);
60 ss.callback.add = scan_samples_add_sample;
61 ss.samples = isl_mat_alloc(ctx, 0, 1 + dim);
62 if (!ss.samples)
63 goto error;
65 if (isl_basic_set_scan(bset, &ss.callback) < 0) {
66 isl_mat_free(ss.samples);
67 return NULL;
70 return ss.samples;
71 error:
72 isl_basic_set_free(bset);
73 return NULL;
76 static struct isl_mat *isl_basic_set_samples(struct isl_basic_set *bset)
78 struct isl_mat *T;
79 struct isl_mat *samples;
81 if (!bset)
82 return NULL;
84 if (bset->n_eq == 0)
85 return isl_basic_set_scan_samples(bset);
87 bset = isl_basic_set_remove_equalities(bset, &T, NULL);
88 samples = isl_basic_set_scan_samples(bset);
89 return isl_mat_product(samples, isl_mat_transpose(T));
92 int main(int argc, char **argv)
94 struct isl_ctx *ctx = isl_ctx_alloc();
95 struct isl_basic_set *bset;
96 struct isl_mat *samples;
98 bset = isl_basic_set_read_from_file(ctx, stdin);
99 samples = isl_basic_set_samples(bset);
100 isl_mat_print_internal(samples, stdout, 0);
101 isl_mat_free(samples);
102 isl_ctx_free(ctx);
104 return 0;