add some isl_map_plain_is_fixed tests
[platform/upstream/isl.git] / isl_test.c
index 55a50d5..cf4fafb 100644 (file)
@@ -12,6 +12,7 @@
 #include <limits.h>
 #include <isl_ctx_private.h>
 #include <isl_map_private.h>
+#include <isl/aff.h>
 #include <isl/set.h>
 #include <isl/flow.h>
 #include <isl/constraint.h>
 
 static char *srcdir;
 
+static char *get_filename(isl_ctx *ctx, const char *name, const char *suffix) {
+       char *filename;
+       int length;
+       char *pattern = "%s/test_inputs/%s.%s";
+
+       length = strlen(pattern) - 6 + strlen(srcdir) + strlen(name)
+               + strlen(suffix) + 1;
+       filename = isl_alloc_array(ctx, char, length);
+
+       if (!filename)
+               return NULL;
+
+       sprintf(filename, pattern, srcdir, name, suffix);
+
+       return filename;
+}
+
 void test_parse_map(isl_ctx *ctx, const char *str)
 {
        isl_map *map;
@@ -108,15 +126,13 @@ void test_parse(struct isl_ctx *ctx)
 
 void test_read(struct isl_ctx *ctx)
 {
-       char filename[PATH_MAX];
+       char *filename;
        FILE *input;
-       int n;
        struct isl_basic_set *bset1, *bset2;
        const char *str = "{[y]: Exists ( alpha : 2alpha = y)}";
 
-       n = snprintf(filename, sizeof(filename),
-                       "%s/test_inputs/set.omega", srcdir);
-       assert(n < sizeof(filename));
+       filename = get_filename(ctx, "set", "omega");
+       assert(filename);
        input = fopen(filename, "r");
        assert(input);
 
@@ -127,6 +143,7 @@ void test_read(struct isl_ctx *ctx)
 
        isl_basic_set_free(bset1);
        isl_basic_set_free(bset2);
+       free(filename);
 
        fclose(input);
 }
@@ -218,16 +235,14 @@ void test_dim(struct isl_ctx *ctx)
 void test_div(struct isl_ctx *ctx)
 {
        isl_int v;
-       int pos;
        struct isl_dim *dim;
-       struct isl_div *div;
        struct isl_basic_set *bset;
        struct isl_constraint *c;
 
        isl_int_init(v);
 
        /* test 1 */
-       dim = isl_dim_set_alloc(ctx, 0, 1);
+       dim = isl_dim_set_alloc(ctx, 0, 3);
        bset = isl_basic_set_universe(dim);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -235,10 +250,8 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_constant(c, v);
        isl_int_set_si(v, 1);
        isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, 3);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -246,17 +259,17 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_constant(c, v);
        isl_int_set_si(v, -1);
        isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, 3);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 2, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
+       bset = isl_basic_set_project_out(bset, isl_dim_set, 1, 2);
+
        assert(bset && bset->n_div == 1);
        isl_basic_set_free(bset);
 
        /* test 2 */
-       dim = isl_dim_set_alloc(ctx, 0, 1);
+       dim = isl_dim_set_alloc(ctx, 0, 3);
        bset = isl_basic_set_universe(dim);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -264,10 +277,8 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_constant(c, v);
        isl_int_set_si(v, -1);
        isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, 3);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -275,17 +286,17 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_constant(c, v);
        isl_int_set_si(v, 1);
        isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, 3);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 2, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
+       bset = isl_basic_set_project_out(bset, isl_dim_set, 1, 2);
+
        assert(bset && bset->n_div == 1);
        isl_basic_set_free(bset);
 
        /* test 3 */
-       dim = isl_dim_set_alloc(ctx, 0, 1);
+       dim = isl_dim_set_alloc(ctx, 0, 3);
        bset = isl_basic_set_universe(dim);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -293,10 +304,8 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_constant(c, v);
        isl_int_set_si(v, -1);
        isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, 3);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -304,17 +313,17 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_constant(c, v);
        isl_int_set_si(v, 1);
        isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, 4);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 2, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
+       bset = isl_basic_set_project_out(bset, isl_dim_set, 1, 2);
+
        assert(bset && bset->n_div == 1);
        isl_basic_set_free(bset);
 
        /* test 4 */
-       dim = isl_dim_set_alloc(ctx, 0, 1);
+       dim = isl_dim_set_alloc(ctx, 0, 3);
        bset = isl_basic_set_universe(dim);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -322,10 +331,8 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_constant(c, v);
        isl_int_set_si(v, -1);
        isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, 3);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -333,26 +340,24 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_constant(c, v);
        isl_int_set_si(v, 1);
        isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, 6);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 2, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
+       bset = isl_basic_set_project_out(bset, isl_dim_set, 1, 2);
+
        assert(isl_basic_set_is_empty(bset));
        isl_basic_set_free(bset);
 
        /* test 5 */
-       dim = isl_dim_set_alloc(ctx, 0, 2);
+       dim = isl_dim_set_alloc(ctx, 0, 3);
        bset = isl_basic_set_universe(dim);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
        isl_int_set_si(v, -1);
        isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, 3);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 2, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -362,20 +367,20 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
+       bset = isl_basic_set_project_out(bset, isl_dim_set, 2, 1);
+
        assert(bset && bset->n_div == 0);
        isl_basic_set_free(bset);
 
        /* test 6 */
-       dim = isl_dim_set_alloc(ctx, 0, 2);
+       dim = isl_dim_set_alloc(ctx, 0, 3);
        bset = isl_basic_set_universe(dim);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
        isl_int_set_si(v, -1);
        isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, 6);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 2, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -385,6 +390,8 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
+       bset = isl_basic_set_project_out(bset, isl_dim_set, 2, 1);
+
        assert(bset && bset->n_div == 1);
        isl_basic_set_free(bset);
 
@@ -398,7 +405,7 @@ void test_div(struct isl_ctx *ctx)
         * and we end up with the original equality and div again.
         * Perhaps we can avoid the introduction of this temporary div.
         */
-       dim = isl_dim_set_alloc(ctx, 0, 3);
+       dim = isl_dim_set_alloc(ctx, 0, 4);
        bset = isl_basic_set_universe(dim);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -408,12 +415,12 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
        isl_int_set_si(v, -3);
        isl_constraint_set_coefficient(c, isl_dim_set, 2, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, 6);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 3, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
+       bset = isl_basic_set_project_out(bset, isl_dim_set, 3, 1);
+
        /* Test disabled for now */
        /*
        assert(bset && bset->n_div == 1);
@@ -421,7 +428,7 @@ void test_div(struct isl_ctx *ctx)
        isl_basic_set_free(bset);
 
        /* test 8 */
-       dim = isl_dim_set_alloc(ctx, 0, 4);
+       dim = isl_dim_set_alloc(ctx, 0, 5);
        bset = isl_basic_set_universe(dim);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -431,10 +438,8 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
        isl_int_set_si(v, -3);
        isl_constraint_set_coefficient(c, isl_dim_set, 3, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, 6);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 4, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -446,6 +451,8 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_constant(c, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
+       bset = isl_basic_set_project_out(bset, isl_dim_set, 4, 1);
+
        /* Test disabled for now */
        /*
        assert(bset && bset->n_div == 1);
@@ -453,7 +460,7 @@ void test_div(struct isl_ctx *ctx)
        isl_basic_set_free(bset);
 
        /* test 9 */
-       dim = isl_dim_set_alloc(ctx, 0, 2);
+       dim = isl_dim_set_alloc(ctx, 0, 4);
        bset = isl_basic_set_universe(dim);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
@@ -461,23 +468,21 @@ void test_div(struct isl_ctx *ctx)
        isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
        isl_int_set_si(v, -1);
        isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, -2);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 2, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
        isl_int_set_si(v, -1);
        isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, 3);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 3, v);
        isl_int_set_si(v, 2);
        isl_constraint_set_constant(c, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
+       bset = isl_basic_set_project_out(bset, isl_dim_set, 2, 2);
+
        bset = isl_basic_set_fix_si(bset, isl_dim_set, 0, 2);
 
        assert(!isl_basic_set_is_empty(bset));
@@ -485,18 +490,18 @@ void test_div(struct isl_ctx *ctx)
        isl_basic_set_free(bset);
 
        /* test 10 */
-       dim = isl_dim_set_alloc(ctx, 0, 2);
+       dim = isl_dim_set_alloc(ctx, 0, 3);
        bset = isl_basic_set_universe(dim);
 
        c = isl_equality_alloc(isl_basic_set_get_dim(bset));
        isl_int_set_si(v, 1);
        isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
-       div = isl_div_alloc(isl_basic_set_get_dim(bset));
-       c = isl_constraint_add_div(c, div, &pos);
        isl_int_set_si(v, -2);
-       isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
+       isl_constraint_set_coefficient(c, isl_dim_set, 2, v);
        bset = isl_basic_set_add_constraint(bset, c);
 
+       bset = isl_basic_set_project_out(bset, isl_dim_set, 2, 1);
+
        bset = isl_basic_set_fix_si(bset, isl_dim_set, 0, 2);
 
        isl_basic_set_free(bset);
@@ -506,15 +511,13 @@ void test_div(struct isl_ctx *ctx)
 
 void test_application_case(struct isl_ctx *ctx, const char *name)
 {
-       char filename[PATH_MAX];
+       char *filename;
        FILE *input;
-       int n;
        struct isl_basic_set *bset1, *bset2;
        struct isl_basic_map *bmap;
 
-       n = snprintf(filename, sizeof(filename),
-                       "%s/test_inputs/%s.omega", srcdir, name);
-       assert(n < sizeof(filename));
+       filename = get_filename(ctx, name, "omega");
+       assert(filename);
        input = fopen(filename, "r");
        assert(input);
 
@@ -529,6 +532,7 @@ void test_application_case(struct isl_ctx *ctx, const char *name)
 
        isl_basic_set_free(bset1);
        isl_basic_set_free(bset2);
+       free(filename);
 
        fclose(input);
 }
@@ -541,14 +545,12 @@ void test_application(struct isl_ctx *ctx)
 
 void test_affine_hull_case(struct isl_ctx *ctx, const char *name)
 {
-       char filename[PATH_MAX];
+       char *filename;
        FILE *input;
-       int n;
        struct isl_basic_set *bset1, *bset2;
 
-       n = snprintf(filename, sizeof(filename),
-                       "%s/test_inputs/%s.polylib", srcdir, name);
-       assert(n < sizeof(filename));
+       filename = get_filename(ctx, name, "polylib");
+       assert(filename);
        input = fopen(filename, "r");
        assert(input);
 
@@ -561,6 +563,7 @@ void test_affine_hull_case(struct isl_ctx *ctx, const char *name)
 
        isl_basic_set_free(bset1);
        isl_basic_set_free(bset2);
+       free(filename);
 
        fclose(input);
 }
@@ -574,15 +577,13 @@ void test_affine_hull(struct isl_ctx *ctx)
 
 void test_convex_hull_case(struct isl_ctx *ctx, const char *name)
 {
-       char filename[PATH_MAX];
+       char *filename;
        FILE *input;
-       int n;
        struct isl_basic_set *bset1, *bset2;
        struct isl_set *set;
 
-       n = snprintf(filename, sizeof(filename),
-                       "%s/test_inputs/%s.polylib", srcdir, name);
-       assert(n < sizeof(filename));
+       filename = get_filename(ctx, name, "polylib");
+       assert(filename);
        input = fopen(filename, "r");
        assert(input);
 
@@ -598,6 +599,7 @@ void test_convex_hull_case(struct isl_ctx *ctx, const char *name)
 
        isl_basic_set_free(bset1);
        isl_basic_set_free(bset2);
+       free(filename);
 
        fclose(input);
 }
@@ -648,14 +650,12 @@ void test_convex_hull(struct isl_ctx *ctx)
 
 void test_gist_case(struct isl_ctx *ctx, const char *name)
 {
-       char filename[PATH_MAX];
+       char *filename;
        FILE *input;
-       int n;
        struct isl_basic_set *bset1, *bset2;
 
-       n = snprintf(filename, sizeof(filename),
-                       "%s/test_inputs/%s.polylib", srcdir, name);
-       assert(n < sizeof(filename));
+       filename = get_filename(ctx, name, "polylib");
+       assert(filename);
        input = fopen(filename, "r");
        assert(input);
 
@@ -670,6 +670,7 @@ void test_gist_case(struct isl_ctx *ctx, const char *name)
 
        isl_basic_set_free(bset1);
        isl_basic_set_free(bset2);
+       free(filename);
 
        fclose(input);
 }
@@ -1216,9 +1217,26 @@ void test_lex(struct isl_ctx *ctx)
        isl_map_free(map);
 }
 
+static int consume_lexmin(__isl_take isl_basic_set *dom,
+       __isl_take isl_aff_list *list, void *user)
+{
+       isl_dim *dim;
+       isl_basic_map *bmap;
+       isl_map **map = user;
+
+       dim = isl_basic_set_get_dim(dom);
+       bmap = isl_basic_map_from_aff_list(dim, list);
+       bmap = isl_basic_map_intersect_domain(bmap, dom);
+
+       *map = isl_map_union(*map, isl_map_from_basic_map(bmap));
+
+       return 0;
+}
+
 void test_lexmin(struct isl_ctx *ctx)
 {
        const char *str;
+       isl_basic_map *bmap;
        isl_map *map, *map2;
        isl_set *set;
        isl_set *set2;
@@ -1278,6 +1296,16 @@ void test_lexmin(struct isl_ctx *ctx)
        assert(isl_map_is_equal(map, map2));
        isl_map_free(map);
        isl_map_free(map2);
+
+       str = "{ [i] -> [i', j] : j = i - 8i' and i' >= 0 and i' <= 7 and "
+                               " 8i' <= i and 8i' >= -7 + i }";
+       bmap = isl_basic_map_read_from_str(ctx, str, -1);
+       map2 = isl_map_empty(isl_basic_map_get_dim(bmap));
+       isl_basic_map_foreach_lexmin(bmap, &consume_lexmin, &map2);
+       map = isl_map_from_basic_map(bmap);
+       assert(isl_map_is_equal(map, map2));
+       isl_map_free(map);
+       isl_map_free(map2);
 }
 
 struct must_may {
@@ -1670,6 +1698,14 @@ void test_pwqp(struct isl_ctx *ctx)
        assert(isl_pw_qpolynomial_is_zero(pwqp1));
 
        isl_pw_qpolynomial_free(pwqp1);
+
+       str = "{ [i] -> ([i/2]) : i >= 0; [i] -> ([i/3]) : i < 0 }";
+       pwqp1 = isl_pw_qpolynomial_read_from_str(ctx, str);
+       pwqp2 = isl_pw_qpolynomial_read_from_str(ctx, str);
+       pwqp1 = isl_pw_qpolynomial_coalesce(pwqp1);
+       pwqp1 = isl_pw_qpolynomial_sub(pwqp1, pwqp2);
+       assert(isl_pw_qpolynomial_is_zero(pwqp1));
+       isl_pw_qpolynomial_free(pwqp1);
 }
 
 void test_split_periods(isl_ctx *ctx)
@@ -1884,8 +1920,9 @@ int test_one_schedule(isl_ctx *ctx, const char *d, const char *w,
                is_tilable = 1;
                is_parallel = 1;
                is_nonneg = 1;
+               isl_union_set_free(delta);
        } else {
-               delta_set = isl_union_set_copy_set(delta);
+               delta_set = isl_set_from_union_set(delta);
 
                slice = isl_set_universe(isl_set_get_dim(delta_set));
                for (i = 0; i < tilable; ++i)
@@ -1911,7 +1948,6 @@ int test_one_schedule(isl_ctx *ctx, const char *d, const char *w,
                isl_set_free(origin);
                isl_set_free(delta_set);
        }
-       isl_union_set_free(delta);
 
        if (is_nonneg < 0 || is_parallel < 0 || is_tilable < 0 ||
            is_injection < 0 || is_complete < 0)
@@ -2182,6 +2218,209 @@ int test_injective(isl_ctx *ctx)
        return 0;
 }
 
+int test_aff(isl_ctx *ctx)
+{
+       const char *str;
+       isl_set *set;
+       isl_dim *dim;
+       isl_local_space *ls;
+       isl_aff *aff;
+       int zero;
+
+       dim = isl_dim_set_alloc(ctx, 0, 1);
+       ls = isl_local_space_from_dim(dim);
+       aff = isl_aff_zero(ls);
+
+       aff = isl_aff_add_coefficient_si(aff, isl_dim_set, 0, 1);
+       aff = isl_aff_scale_down_ui(aff, 3);
+       aff = isl_aff_floor(aff);
+       aff = isl_aff_add_coefficient_si(aff, isl_dim_set, 0, 1);
+       aff = isl_aff_scale_down_ui(aff, 2);
+       aff = isl_aff_floor(aff);
+       aff = isl_aff_add_coefficient_si(aff, isl_dim_set, 0, 1);
+
+       str = "{ [10] }";
+       set = isl_set_read_from_str(ctx, str, 0);
+       aff = isl_aff_gist(aff, set);
+
+       aff = isl_aff_add_constant_si(aff, -16);
+       zero = isl_aff_plain_is_zero(aff);
+       isl_aff_free(aff);
+
+       if (zero < 0)
+               return -1;
+       if (!zero)
+               isl_die(ctx, isl_error_unknown, "unexpected result", return -1);
+
+       return 0;
+}
+
+int test_dim_max(isl_ctx *ctx)
+{
+       int equal;
+       const char *str;
+       isl_map *map, *map2;
+       isl_set *set;
+       isl_pw_aff *pwaff;
+
+       str = "[N] -> { [i] : 0 <= i <= min(N,10) }";
+       set = isl_set_read_from_str(ctx, str, -1);
+       pwaff = isl_set_dim_max(set, 0);
+       map = isl_map_from_pw_aff(pwaff);
+       str = "[N] -> { [] -> [10] : N >= 10; [] -> [N] : N <= 9 and N >= 0 }";
+       map2 = isl_map_read_from_str(ctx, str, -1);
+       equal = isl_map_is_equal(map, map2);
+       isl_map_free(map);
+       isl_map_free(map2);
+       if (equal < 0)
+               return -1;
+       if (!equal)
+               isl_die(ctx, isl_error_unknown, "unexpected result", return -1);
+
+       str = "[N] -> { [i] : 0 <= i <= max(2N,N+6) }";
+       set = isl_set_read_from_str(ctx, str, -1);
+       pwaff = isl_set_dim_max(set, 0);
+       map = isl_map_from_pw_aff(pwaff);
+       str = "[N] -> { [] -> [6 + N] : -6 <= N <= 5; [] -> [2N] : N >= 6 }";
+       map2 = isl_map_read_from_str(ctx, str, -1);
+       equal = isl_map_is_equal(map, map2);
+       isl_map_free(map);
+       isl_map_free(map2);
+       if (equal < 0)
+               return -1;
+       if (!equal)
+               isl_die(ctx, isl_error_unknown, "unexpected result", return -1);
+
+       str = "[N] -> { [i] : 0 <= i <= 2N or 0 <= i <= N+6 }";
+       set = isl_set_read_from_str(ctx, str, -1);
+       pwaff = isl_set_dim_max(set, 0);
+       map = isl_map_from_pw_aff(pwaff);
+       str = "[N] -> { [] -> [6 + N] : -6 <= N <= 5; [] -> [2N] : N >= 6 }";
+       map2 = isl_map_read_from_str(ctx, str, -1);
+       equal = isl_map_is_equal(map, map2);
+       isl_map_free(map);
+       isl_map_free(map2);
+       if (equal < 0)
+               return -1;
+       if (!equal)
+               isl_die(ctx, isl_error_unknown, "unexpected result", return -1);
+
+       str = "[N,M] -> { [i,j] -> [([i/16]), i%16, ([j/16]), j%16] : "
+                       "0 <= i < N and 0 <= j < M }";
+       map = isl_map_read_from_str(ctx, str, -1);
+       set = isl_map_range(map);
+
+       pwaff = isl_set_dim_max(isl_set_copy(set), 0);
+       map = isl_map_from_pw_aff(pwaff);
+       str = "[N,M] -> { [] -> [([(N-1)/16])] : M,N > 0 }";
+       map2 = isl_map_read_from_str(ctx, str, -1);
+       equal = isl_map_is_equal(map, map2);
+       isl_map_free(map);
+       isl_map_free(map2);
+
+       pwaff = isl_set_dim_max(isl_set_copy(set), 3);
+       map = isl_map_from_pw_aff(pwaff);
+       str = "[N,M] -> { [] -> [t] : t = min(M-1,15) and M,N > 0 }";
+       map2 = isl_map_read_from_str(ctx, str, -1);
+       if (equal >= 0 && equal)
+               equal = isl_map_is_equal(map, map2);
+       isl_map_free(map);
+       isl_map_free(map2);
+
+       isl_set_free(set);
+
+       if (equal < 0)
+               return -1;
+       if (!equal)
+               isl_die(ctx, isl_error_unknown, "unexpected result", return -1);
+
+       return 0;
+}
+
+int test_product(isl_ctx *ctx)
+{
+       const char *str;
+       isl_set *set;
+       int ok;
+
+       str = "{ A[i] }";
+       set = isl_set_read_from_str(ctx, str, -1);
+       set = isl_set_product(set, isl_set_copy(set));
+       ok = isl_set_is_wrapping(set);
+       isl_set_free(set);
+       if (ok < 0)
+               return -1;
+       if (!ok)
+               isl_die(ctx, isl_error_unknown, "unexpected result", return -1);
+
+       return 0;
+}
+
+int test_equal(isl_ctx *ctx)
+{
+       const char *str;
+       isl_set *set, *set2;
+       int equal;
+
+       str = "{ S_6[i] }";
+       set = isl_set_read_from_str(ctx, str, -1);
+       str = "{ S_7[i] }";
+       set2 = isl_set_read_from_str(ctx, str, -1);
+       equal = isl_set_is_equal(set, set2);
+       isl_set_free(set);
+       isl_set_free(set2);
+       if (equal < 0)
+               return -1;
+       if (equal)
+               isl_die(ctx, isl_error_unknown, "unexpected result", return -1);
+
+       return 0;
+}
+
+static int test_plain_fixed(isl_ctx *ctx, __isl_take isl_map *map,
+       enum isl_dim_type type, unsigned pos, int fixed)
+{
+       int test;
+
+       test = isl_map_plain_is_fixed(map, type, pos, NULL);
+       isl_map_free(map);
+       if (test < 0)
+               return -1;
+       if (test == fixed)
+               return 0;
+       if (fixed)
+               isl_die(ctx, isl_error_unknown,
+                       "map not detected as fixed", return -1);
+       else
+               isl_die(ctx, isl_error_unknown,
+                       "map detected as fixed", return -1);
+}
+
+int test_fixed(isl_ctx *ctx)
+{
+       const char *str;
+       isl_map *map;
+
+       str = "{ [i] -> [i] }";
+       map = isl_map_read_from_str(ctx, str, -1);
+       if (test_plain_fixed(ctx, map, isl_dim_out, 0, 0))
+               return -1;
+       str = "{ [i] -> [1] }";
+       map = isl_map_read_from_str(ctx, str, -1);
+       if (test_plain_fixed(ctx, map, isl_dim_out, 0, 1))
+               return -1;
+       str = "{ S_1[p1] -> [o0] : o0 = -2 and p1 >= 1 and p1 <= 7 }";
+       map = isl_map_read_from_str(ctx, str, -1);
+       if (test_plain_fixed(ctx, map, isl_dim_out, 0, 1))
+               return -1;
+       map = isl_map_read_from_str(ctx, str, -1);
+       map = isl_map_neg(map);
+       if (test_plain_fixed(ctx, map, isl_dim_out, 0, 1))
+               return -1;
+
+       return 0;
+}
+
 int main()
 {
        struct isl_ctx *ctx;
@@ -2190,6 +2429,16 @@ int main()
        assert(srcdir);
 
        ctx = isl_ctx_alloc();
+       if (test_fixed(ctx) < 0)
+               goto error;
+       if (test_equal(ctx) < 0)
+               goto error;
+       if (test_product(ctx) < 0)
+               goto error;
+       if (test_dim_max(ctx) < 0)
+               goto error;
+       if (test_aff(ctx) < 0)
+               goto error;
        if (test_injective(ctx) < 0)
                goto error;
        if (test_schedule(ctx) < 0)