isl_printer_print_pw_qpolynomial_fold: use consistent names in C format
[platform/upstream/isl.git] / isl_convex_hull.c
index 663f142..1552492 100644 (file)
@@ -10,7 +10,7 @@
 #include "isl_lp.h"
 #include "isl_map.h"
 #include "isl_map_private.h"
-#include "isl_mat.h"
+#include <isl_mat_private.h>
 #include "isl_set.h"
 #include "isl_seq.h"
 #include "isl_equalities.h"
@@ -427,6 +427,7 @@ isl_int *isl_set_wrap_facet(__isl_keep isl_set *set,
        if (res == isl_lp_ok) {
                isl_int_neg(num, num);
                isl_seq_combine(facet, num, facet, den, ridge, dim);
+               isl_seq_normalize(ctx, facet, dim);
        }
        isl_int_clear(num);
        isl_int_clear(den);
@@ -789,7 +790,7 @@ error:
 static struct isl_set *set_project_out(struct isl_ctx *ctx,
        struct isl_set *set, unsigned n)
 {
-       return isl_set_remove_dims(set, isl_set_n_dim(set) - n, n);
+       return isl_set_remove_dims(set, isl_dim_set, isl_set_n_dim(set) - n, n);
 }
 
 static struct isl_basic_set *convex_hull_0d(struct isl_set *set)
@@ -868,7 +869,7 @@ static struct isl_basic_set *convex_hull_pair_elim(struct isl_basic_set *bset1,
                isl_int_set_si(hull->eq[k][2*(1+dim)+j], 1);
        }
        hull = isl_basic_set_set_rational(hull);
-       hull = isl_basic_set_remove_dims(hull, dim, 2*(1+dim));
+       hull = isl_basic_set_remove_dims(hull, isl_dim_set, dim, 2*(1+dim));
        hull = isl_basic_set_remove_redundancies(hull);
        isl_basic_set_free(bset1);
        isl_basic_set_free(bset2);
@@ -898,6 +899,25 @@ int isl_basic_set_is_bounded(__isl_keep isl_basic_set *bset)
        return bounded;
 }
 
+/* Is the image bounded for each value of the parameters and
+ * the domain variables?
+ */
+int isl_basic_map_image_is_bounded(__isl_keep isl_basic_map *bmap)
+{
+       unsigned nparam = isl_basic_map_dim(bmap, isl_dim_param);
+       unsigned n_in = isl_basic_map_dim(bmap, isl_dim_in);
+       int bounded;
+
+       bmap = isl_basic_map_copy(bmap);
+       bmap = isl_basic_map_cow(bmap);
+       bmap = isl_basic_map_move_dims(bmap, isl_dim_param, nparam,
+                                       isl_dim_in, 0, n_in);
+       bounded = isl_basic_set_is_bounded((isl_basic_set *)bmap);
+       isl_basic_map_free(bmap);
+
+       return bounded;
+}
+
 /* Is the set bounded for each value of the parameters?
  */
 int isl_set_is_bounded(__isl_keep isl_set *set)
@@ -1023,7 +1043,7 @@ static struct isl_basic_set *modulo_lineality(struct isl_set *set,
        Q = isl_mat_lin_to_aff(Q);
 
        set = isl_set_preimage(set, U);
-       set = isl_set_remove_dims(set, total - lin_dim, lin_dim);
+       set = isl_set_remove_dims(set, isl_dim_set, total - lin_dim, lin_dim);
        hull = uset_convex_hull(set);
        hull = isl_basic_set_preimage(hull, Q);
 
@@ -1936,6 +1956,19 @@ struct isl_basic_set *isl_set_convex_hull(struct isl_set *set)
                isl_map_convex_hull((struct isl_map *)set);
 }
 
+__isl_give isl_basic_map *isl_map_polyhedral_hull(__isl_take isl_map *map)
+{
+       isl_basic_map *hull;
+
+       hull = isl_map_convex_hull(map);
+       return isl_basic_map_remove_divs(hull);
+}
+
+__isl_give isl_basic_set *isl_set_polyhedral_hull(__isl_take isl_set *set)
+{
+       return (isl_basic_set *)isl_map_polyhedral_hull((isl_map *)set);
+}
+
 struct sh_data_entry {
        struct isl_hash_table   *table;
        struct isl_tab          *tab;