add isl_set_lower_bound_dim

diff --git a/include/isl_set.h b/include/isl_set.h
index a4262d0..c1bad36 100644 (file)
--- a/include/isl_set.h
+++ b/include/isl_set.h
@@ -132,6 +132,8 @@ struct isl_set *isl_set_subtract(struct isl_set *set1, struct isl_set *set2);
 struct isl_set *isl_set_apply(struct isl_set *set, struct isl_map *map);
 struct isl_set *isl_set_fix_dim_si(struct isl_set *set,
                unsigned dim, int value);
+struct isl_set *isl_set_lower_bound_dim(struct isl_set *set,
+               unsigned dim, isl_int value);
 struct isl_basic_set *isl_basic_set_remove_dims(struct isl_basic_set *bset,
                unsigned first, unsigned n);
 struct isl_set *isl_set_remove_dims(struct isl_set *set,

diff --git a/isl_map.c b/isl_map.c
index 47a1218..f405b96 100644 (file)
--- a/isl_map.c
+++ b/isl_map.c
@@ -2557,6 +2557,47 @@ error:
        return NULL;
 }
 
+struct isl_basic_set *isl_basic_set_lower_bound_dim(struct isl_basic_set *bset,
+       unsigned dim, isl_int value)
+{
+       int j;
+
+       bset = isl_basic_set_cow(bset);
+       bset = isl_basic_set_extend(bset, bset->nparam, bset->dim, 0, 0, 1);
+       j = isl_basic_set_alloc_inequality(bset);
+       if (j < 0)
+               goto error;
+       isl_seq_clr(bset->ineq[j], 1 + bset->nparam + bset->dim + bset->extra);
+       isl_int_set_si(bset->ineq[j][1 + bset->nparam + dim], 1);
+       isl_int_neg(bset->ineq[j][0], value);
+       bset = isl_basic_set_simplify(bset);
+       return isl_basic_set_finalize(bset);
+error:
+       isl_basic_set_free(bset);
+       return NULL;
+}
+
+struct isl_set *isl_set_lower_bound_dim(struct isl_set *set, unsigned dim,
+                                       isl_int value)
+{
+       int i;
+
+       set = isl_set_cow(set);
+       if (!set)
+               return NULL;
+
+       isl_assert(set->ctx, dim < set->dim, goto error);
+       for (i = 0; i < set->n; ++i) {
+               set->p[i] = isl_basic_set_lower_bound_dim(set->p[i], dim, value);
+               if (!set->p[i])
+                       goto error;
+       }
+       return set;
+error:
+       isl_set_free(set);
+       return NULL;
+}
+
 struct isl_map *isl_map_reverse(struct isl_map *map)
 {
        int i;