isl_test.c

   1 /*
   2  * Copyright 2008-2009 Katholieke Universiteit Leuven
   3  *
   4  * Use of this software is governed by the GNU LGPLv2.1 license
   5  *
   6  * Written by Sven Verdoolaege, K.U.Leuven, Departement
   7  * Computerwetenschappen, Celestijnenlaan 200A, B-3001 Leuven, Belgium
   8  */
   9
  10 #include <assert.h>
  11 #include <stdio.h>
  12 #include <limits.h>
  13 #include <isl_ctx.h>
  14 #include <isl_set.h>
  15 #include <isl_constraint.h>
  16
  17 static char *srcdir;
  18
  19 void test_read(struct isl_ctx *ctx)
  20 {
  21         char filename[PATH_MAX];
  22         FILE *input;
  23         int n;
  24         struct isl_basic_set *bset1, *bset2;
  25         const char *str = "{[y]: Exists ( alpha : 2alpha = y)}";
  26
  27         n = snprintf(filename, sizeof(filename),
  28                         "%s/test_inputs/set.omega", srcdir);
  29         assert(n < sizeof(filename));
  30         input = fopen(filename, "r");
  31         assert(input);
  32
  33         bset1 = isl_basic_set_read_from_file(ctx, input, 0);
  34         bset2 = isl_basic_set_read_from_str(ctx, str, 0);
  35
  36         assert(isl_basic_set_is_equal(bset1, bset2) == 1);
  37
  38         isl_basic_set_free(bset1);
  39         isl_basic_set_free(bset2);
  40
  41         fclose(input);
  42 }
  43
  44 /* Construct the basic set { [i] : 5 <= i <= N } */
  45 void test_construction(struct isl_ctx *ctx)
  46 {
  47         isl_int v;
  48         struct isl_dim *dim;
  49         struct isl_basic_set *bset;
  50         struct isl_constraint *c;
  51
  52         isl_int_init(v);
  53
  54         dim = isl_dim_set_alloc(ctx, 1, 1);
  55         bset = isl_basic_set_universe(dim);
  56
  57         c = isl_inequality_alloc(isl_dim_copy(bset->dim));
  58         isl_int_set_si(v, -1);
  59         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
  60         isl_int_set_si(v, 1);
  61         isl_constraint_set_coefficient(c, isl_dim_param, 0, v);
  62         bset = isl_basic_set_add_constraint(bset, c);
  63
  64         c = isl_inequality_alloc(isl_dim_copy(bset->dim));
  65         isl_int_set_si(v, 1);
  66         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
  67         isl_int_set_si(v, -5);
  68         isl_constraint_set_constant(c, v);
  69         bset = isl_basic_set_add_constraint(bset, c);
  70
  71         isl_basic_set_free(bset);
  72
  73         isl_int_clear(v);
  74 }
  75
  76 void test_div(struct isl_ctx *ctx)
  77 {
  78         isl_int v;
  79         int pos;
  80         struct isl_dim *dim;
  81         struct isl_div *div;
  82         struct isl_basic_set *bset;
  83         struct isl_constraint *c;
  84
  85         isl_int_init(v);
  86
  87         /* test 1 */
  88         dim = isl_dim_set_alloc(ctx, 0, 1);
  89         bset = isl_basic_set_universe(dim);
  90
  91         c = isl_equality_alloc(isl_dim_copy(bset->dim));
  92         isl_int_set_si(v, -1);
  93         isl_constraint_set_constant(c, v);
  94         isl_int_set_si(v, 1);
  95         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
  96         div = isl_div_alloc(isl_dim_copy(bset->dim));
  97         c = isl_constraint_add_div(c, div, &pos);
  98         isl_int_set_si(v, 3);
  99         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 100         bset = isl_basic_set_add_constraint(bset, c);
 101
 102         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 103         isl_int_set_si(v, 1);
 104         isl_constraint_set_constant(c, v);
 105         isl_int_set_si(v, -1);
 106         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 107         div = isl_div_alloc(isl_dim_copy(bset->dim));
 108         c = isl_constraint_add_div(c, div, &pos);
 109         isl_int_set_si(v, 3);
 110         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 111         bset = isl_basic_set_add_constraint(bset, c);
 112
 113         assert(bset->n_div == 1);
 114         isl_basic_set_free(bset);
 115
 116         /* test 2 */
 117         dim = isl_dim_set_alloc(ctx, 0, 1);
 118         bset = isl_basic_set_universe(dim);
 119
 120         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 121         isl_int_set_si(v, 1);
 122         isl_constraint_set_constant(c, v);
 123         isl_int_set_si(v, -1);
 124         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 125         div = isl_div_alloc(isl_dim_copy(bset->dim));
 126         c = isl_constraint_add_div(c, div, &pos);
 127         isl_int_set_si(v, 3);
 128         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 129         bset = isl_basic_set_add_constraint(bset, c);
 130
 131         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 132         isl_int_set_si(v, -1);
 133         isl_constraint_set_constant(c, v);
 134         isl_int_set_si(v, 1);
 135         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 136         div = isl_div_alloc(isl_dim_copy(bset->dim));
 137         c = isl_constraint_add_div(c, div, &pos);
 138         isl_int_set_si(v, 3);
 139         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 140         bset = isl_basic_set_add_constraint(bset, c);
 141
 142         assert(bset->n_div == 1);
 143         isl_basic_set_free(bset);
 144
 145         /* test 3 */
 146         dim = isl_dim_set_alloc(ctx, 0, 1);
 147         bset = isl_basic_set_universe(dim);
 148
 149         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 150         isl_int_set_si(v, 1);
 151         isl_constraint_set_constant(c, v);
 152         isl_int_set_si(v, -1);
 153         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 154         div = isl_div_alloc(isl_dim_copy(bset->dim));
 155         c = isl_constraint_add_div(c, div, &pos);
 156         isl_int_set_si(v, 3);
 157         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 158         bset = isl_basic_set_add_constraint(bset, c);
 159
 160         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 161         isl_int_set_si(v, -3);
 162         isl_constraint_set_constant(c, v);
 163         isl_int_set_si(v, 1);
 164         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 165         div = isl_div_alloc(isl_dim_copy(bset->dim));
 166         c = isl_constraint_add_div(c, div, &pos);
 167         isl_int_set_si(v, 4);
 168         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 169         bset = isl_basic_set_add_constraint(bset, c);
 170
 171         assert(bset->n_div == 1);
 172         isl_basic_set_free(bset);
 173
 174         /* test 4 */
 175         dim = isl_dim_set_alloc(ctx, 0, 1);
 176         bset = isl_basic_set_universe(dim);
 177
 178         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 179         isl_int_set_si(v, 2);
 180         isl_constraint_set_constant(c, v);
 181         isl_int_set_si(v, -1);
 182         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 183         div = isl_div_alloc(isl_dim_copy(bset->dim));
 184         c = isl_constraint_add_div(c, div, &pos);
 185         isl_int_set_si(v, 3);
 186         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 187         bset = isl_basic_set_add_constraint(bset, c);
 188
 189         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 190         isl_int_set_si(v, -1);
 191         isl_constraint_set_constant(c, v);
 192         isl_int_set_si(v, 1);
 193         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 194         div = isl_div_alloc(isl_dim_copy(bset->dim));
 195         c = isl_constraint_add_div(c, div, &pos);
 196         isl_int_set_si(v, 6);
 197         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 198         bset = isl_basic_set_add_constraint(bset, c);
 199
 200         assert(isl_basic_set_is_empty(bset));
 201         isl_basic_set_free(bset);
 202
 203         /* test 5 */
 204         dim = isl_dim_set_alloc(ctx, 0, 2);
 205         bset = isl_basic_set_universe(dim);
 206
 207         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 208         isl_int_set_si(v, -1);
 209         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 210         div = isl_div_alloc(isl_dim_copy(bset->dim));
 211         c = isl_constraint_add_div(c, div, &pos);
 212         isl_int_set_si(v, 3);
 213         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 214         bset = isl_basic_set_add_constraint(bset, c);
 215
 216         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 217         isl_int_set_si(v, 1);
 218         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 219         isl_int_set_si(v, -3);
 220         isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
 221         bset = isl_basic_set_add_constraint(bset, c);
 222
 223         assert(bset->n_div == 0);
 224         isl_basic_set_free(bset);
 225
 226         /* test 6 */
 227         dim = isl_dim_set_alloc(ctx, 0, 2);
 228         bset = isl_basic_set_universe(dim);
 229
 230         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 231         isl_int_set_si(v, -1);
 232         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 233         div = isl_div_alloc(isl_dim_copy(bset->dim));
 234         c = isl_constraint_add_div(c, div, &pos);
 235         isl_int_set_si(v, 6);
 236         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 237         bset = isl_basic_set_add_constraint(bset, c);
 238
 239         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 240         isl_int_set_si(v, 1);
 241         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 242         isl_int_set_si(v, -3);
 243         isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
 244         bset = isl_basic_set_add_constraint(bset, c);
 245
 246         assert(bset->n_div == 1);
 247         isl_basic_set_free(bset);
 248
 249         /* test 7 */
 250         /* This test is a bit tricky.  We set up an equality
 251          *              a + 3b + 3c = 6 e0
 252          * Normalization of divs creates _two_ divs
 253          *              a = 3 e0
 254          *              c - b - e0 = 2 e1
 255          * Afterwards e0 is removed again because it has coefficient -1
 256          * and we end up with the original equality and div again.
 257          * Perhaps we can avoid the introduction of this temporary div.
 258          */
 259         dim = isl_dim_set_alloc(ctx, 0, 3);
 260         bset = isl_basic_set_universe(dim);
 261
 262         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 263         isl_int_set_si(v, -1);
 264         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 265         isl_int_set_si(v, -3);
 266         isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
 267         isl_int_set_si(v, -3);
 268         isl_constraint_set_coefficient(c, isl_dim_set, 2, v);
 269         div = isl_div_alloc(isl_dim_copy(bset->dim));
 270         c = isl_constraint_add_div(c, div, &pos);
 271         isl_int_set_si(v, 6);
 272         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 273         bset = isl_basic_set_add_constraint(bset, c);
 274
 275         assert(bset->n_div == 1);
 276         isl_basic_set_free(bset);
 277
 278         /* test 8 */
 279         dim = isl_dim_set_alloc(ctx, 0, 4);
 280         bset = isl_basic_set_universe(dim);
 281
 282         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 283         isl_int_set_si(v, -1);
 284         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 285         isl_int_set_si(v, -3);
 286         isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
 287         isl_int_set_si(v, -3);
 288         isl_constraint_set_coefficient(c, isl_dim_set, 3, v);
 289         div = isl_div_alloc(isl_dim_copy(bset->dim));
 290         c = isl_constraint_add_div(c, div, &pos);
 291         isl_int_set_si(v, 6);
 292         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 293         bset = isl_basic_set_add_constraint(bset, c);
 294
 295         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 296         isl_int_set_si(v, -1);
 297         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 298         isl_int_set_si(v, 1);
 299         isl_constraint_set_coefficient(c, isl_dim_set, 2, v);
 300         isl_int_set_si(v, 1);
 301         isl_constraint_set_constant(c, v);
 302         bset = isl_basic_set_add_constraint(bset, c);
 303
 304         assert(bset->n_div == 1);
 305         isl_basic_set_free(bset);
 306
 307         /* test 9 */
 308         dim = isl_dim_set_alloc(ctx, 0, 2);
 309         bset = isl_basic_set_universe(dim);
 310
 311         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 312         isl_int_set_si(v, 1);
 313         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 314         isl_int_set_si(v, -1);
 315         isl_constraint_set_coefficient(c, isl_dim_set, 1, v);
 316         div = isl_div_alloc(isl_dim_copy(bset->dim));
 317         c = isl_constraint_add_div(c, div, &pos);
 318         isl_int_set_si(v, -2);
 319         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 320         bset = isl_basic_set_add_constraint(bset, c);
 321
 322         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 323         isl_int_set_si(v, -1);
 324         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 325         div = isl_div_alloc(isl_dim_copy(bset->dim));
 326         c = isl_constraint_add_div(c, div, &pos);
 327         isl_int_set_si(v, 3);
 328         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 329         isl_int_set_si(v, 2);
 330         isl_constraint_set_constant(c, v);
 331         bset = isl_basic_set_add_constraint(bset, c);
 332
 333         bset = isl_basic_set_fix_si(bset, isl_dim_set, 0, 2);
 334
 335         assert(!isl_basic_set_is_empty(bset));
 336
 337         isl_basic_set_free(bset);
 338
 339         /* test 10 */
 340         dim = isl_dim_set_alloc(ctx, 0, 2);
 341         bset = isl_basic_set_universe(dim);
 342
 343         c = isl_equality_alloc(isl_dim_copy(bset->dim));
 344         isl_int_set_si(v, 1);
 345         isl_constraint_set_coefficient(c, isl_dim_set, 0, v);
 346         div = isl_div_alloc(isl_dim_copy(bset->dim));
 347         c = isl_constraint_add_div(c, div, &pos);
 348         isl_int_set_si(v, -2);
 349         isl_constraint_set_coefficient(c, isl_dim_div, pos, v);
 350         bset = isl_basic_set_add_constraint(bset, c);
 351
 352         bset = isl_basic_set_fix_si(bset, isl_dim_set, 0, 2);
 353
 354         isl_basic_set_free(bset);
 355
 356         isl_int_clear(v);
 357 }
 358
 359 void test_application_case(struct isl_ctx *ctx, const char *name)
 360 {
 361         char filename[PATH_MAX];
 362         FILE *input;
 363         int n;
 364         struct isl_basic_set *bset1, *bset2;
 365         struct isl_basic_map *bmap;
 366
 367         n = snprintf(filename, sizeof(filename),
 368                         "%s/test_inputs/%s.omega", srcdir, name);
 369         assert(n < sizeof(filename));
 370         input = fopen(filename, "r");
 371         assert(input);
 372
 373         bset1 = isl_basic_set_read_from_file(ctx, input, 0);
 374         bmap = isl_basic_map_read_from_file(ctx, input, 0);
 375
 376         bset1 = isl_basic_set_apply(bset1, bmap);
 377
 378         bset2 = isl_basic_set_read_from_file(ctx, input, 0);
 379
 380         assert(isl_basic_set_is_equal(bset1, bset2) == 1);
 381
 382         isl_basic_set_free(bset1);
 383         isl_basic_set_free(bset2);
 384
 385         fclose(input);
 386 }
 387
 388 void test_application(struct isl_ctx *ctx)
 389 {
 390         test_application_case(ctx, "application");
 391         test_application_case(ctx, "application2");
 392 }
 393
 394 void test_affine_hull_case(struct isl_ctx *ctx, const char *name)
 395 {
 396         char filename[PATH_MAX];
 397         FILE *input;
 398         int n;
 399         struct isl_basic_set *bset1, *bset2;
 400
 401         n = snprintf(filename, sizeof(filename),
 402                         "%s/test_inputs/%s.polylib", srcdir, name);
 403         assert(n < sizeof(filename));
 404         input = fopen(filename, "r");
 405         assert(input);
 406
 407         bset1 = isl_basic_set_read_from_file(ctx, input, 0);
 408         bset2 = isl_basic_set_read_from_file(ctx, input, 0);
 409
 410         bset1 = isl_basic_set_affine_hull(bset1);
 411
 412         assert(isl_basic_set_is_equal(bset1, bset2) == 1);
 413
 414         isl_basic_set_free(bset1);
 415         isl_basic_set_free(bset2);
 416
 417         fclose(input);
 418 }
 419
 420 void test_affine_hull(struct isl_ctx *ctx)
 421 {
 422         test_affine_hull_case(ctx, "affine2");
 423         test_affine_hull_case(ctx, "affine");
 424         test_affine_hull_case(ctx, "affine3");
 425 }
 426
 427 void test_convex_hull_case(struct isl_ctx *ctx, const char *name)
 428 {
 429         char filename[PATH_MAX];
 430         FILE *input;
 431         int n;
 432         struct isl_basic_set *bset1, *bset2;
 433         struct isl_set *set;
 434
 435         n = snprintf(filename, sizeof(filename),
 436                         "%s/test_inputs/%s.polylib", srcdir, name);
 437         assert(n < sizeof(filename));
 438         input = fopen(filename, "r");
 439         assert(input);
 440
 441         bset1 = isl_basic_set_read_from_file(ctx, input, 0);
 442         bset2 = isl_basic_set_read_from_file(ctx, input, 0);
 443
 444         set = isl_basic_set_union(bset1, bset2);
 445         bset1 = isl_set_convex_hull(set);
 446
 447         bset2 = isl_basic_set_read_from_file(ctx, input, 0);
 448
 449         assert(isl_basic_set_is_equal(bset1, bset2) == 1);
 450
 451         isl_basic_set_free(bset1);
 452         isl_basic_set_free(bset2);
 453
 454         fclose(input);
 455 }
 456
 457 void test_convex_hull(struct isl_ctx *ctx)
 458 {
 459         test_convex_hull_case(ctx, "convex0");
 460         test_convex_hull_case(ctx, "convex1");
 461         test_convex_hull_case(ctx, "convex2");
 462         test_convex_hull_case(ctx, "convex3");
 463         test_convex_hull_case(ctx, "convex4");
 464         test_convex_hull_case(ctx, "convex5");
 465         test_convex_hull_case(ctx, "convex6");
 466         test_convex_hull_case(ctx, "convex7");
 467         test_convex_hull_case(ctx, "convex8");
 468         test_convex_hull_case(ctx, "convex9");
 469         test_convex_hull_case(ctx, "convex10");
 470         test_convex_hull_case(ctx, "convex11");
 471         test_convex_hull_case(ctx, "convex12");
 472         test_convex_hull_case(ctx, "convex13");
 473         test_convex_hull_case(ctx, "convex14");
 474 }
 475
 476 void test_gist_case(struct isl_ctx *ctx, const char *name)
 477 {
 478         char filename[PATH_MAX];
 479         FILE *input;
 480         int n;
 481         struct isl_basic_set *bset1, *bset2;
 482
 483         n = snprintf(filename, sizeof(filename),
 484                         "%s/test_inputs/%s.polylib", srcdir, name);
 485         assert(n < sizeof(filename));
 486         input = fopen(filename, "r");
 487         assert(input);
 488
 489         bset1 = isl_basic_set_read_from_file(ctx, input, 0);
 490         bset2 = isl_basic_set_read_from_file(ctx, input, 0);
 491
 492         bset1 = isl_basic_set_gist(bset1, bset2);
 493
 494         bset2 = isl_basic_set_read_from_file(ctx, input, 0);
 495
 496         assert(isl_basic_set_is_equal(bset1, bset2) == 1);
 497
 498         isl_basic_set_free(bset1);
 499         isl_basic_set_free(bset2);
 500
 501         fclose(input);
 502 }
 503
 504 void test_gist(struct isl_ctx *ctx)
 505 {
 506         test_gist_case(ctx, "gist1");
 507 }
 508
 509 static struct isl_set *s_union(struct isl_ctx *ctx,
 510         const char *s1, const char *s2)
 511 {
 512         struct isl_basic_set *bset1;
 513         struct isl_basic_set *bset2;
 514         struct isl_set *set1, *set2;
 515
 516         bset1 = isl_basic_set_read_from_str(ctx, s1, 0);
 517         bset2 = isl_basic_set_read_from_str(ctx, s2, 0);
 518         set1 = isl_set_from_basic_set(bset1);
 519         set2 = isl_set_from_basic_set(bset2);
 520         return isl_set_union(set1, set2);
 521 }
 522
 523 void test_coalesce(struct isl_ctx *ctx)
 524 {
 525         struct isl_set *set;
 526
 527         set = s_union(ctx, "{[x,y]: x >= 0 & x <= 10 & y >= 0 & y <= 10}",
 528                            "{[x,y]: y >= x & x >= 2 & 5 >= y}");
 529         set = isl_set_coalesce(set);
 530         assert(set && set->n == 1);
 531         isl_set_free(set);
 532
 533         set = s_union(ctx, "{[x,y]: y >= 0 & 2x + y <= 30 & y <= 10 & x >= 0}",
 534                        "{[x,y]: x + y >= 10 & y <= x & x + y <= 20 & y >= 0}");
 535         set = isl_set_coalesce(set);
 536         assert(set && set->n == 1);
 537         isl_set_free(set);
 538
 539         set = s_union(ctx, "{[x,y]: y >= 0 & 2x + y <= 30 & y <= 10 & x >= 0}",
 540                        "{[x,y]: x + y >= 10 & y <= x & x + y <= 19 & y >= 0}");
 541         set = isl_set_coalesce(set);
 542         assert(set && set->n == 2);
 543         isl_set_free(set);
 544
 545         set = s_union(ctx, "{[x,y]: y >= 0 & x <= 5 & y <= x}",
 546                        "{[x,y]: y >= 0 & x >= 6 & x <= 10 & y <= x}");
 547         set = isl_set_coalesce(set);
 548         assert(set && set->n == 1);
 549         isl_set_free(set);
 550
 551         set = s_union(ctx, "{[x,y]: y >= 0 & x <= 5 & y <= x}",
 552                        "{[x,y]: y >= 0 & x >= 7 & x <= 10 & y <= x}");
 553         set = isl_set_coalesce(set);
 554         assert(set && set->n == 2);
 555         isl_set_free(set);
 556
 557         set = s_union(ctx, "{[x,y]: y >= 0 & x <= 5 & y <= x}",
 558                        "{[x,y]: y >= 0 & x >= 6 & x <= 10 & y + 1 <= x}");
 559         set = isl_set_coalesce(set);
 560         assert(set && set->n == 2);
 561         isl_set_free(set);
 562
 563         set = s_union(ctx, "{[x,y]: y >= 0 & x <= 5 & y <= x}",
 564                        "{[x,y]: y >= 0 & x = 6 & y <= 6}");
 565         set = isl_set_coalesce(set);
 566         assert(set && set->n == 1);
 567         isl_set_free(set);
 568
 569         set = s_union(ctx, "{[x,y]: y >= 0 & x <= 5 & y <= x}",
 570                        "{[x,y]: y >= 0 & x = 7 & y <= 6}");
 571         set = isl_set_coalesce(set);
 572         assert(set && set->n == 2);
 573         isl_set_free(set);
 574
 575         set = s_union(ctx, "{[x,y]: y >= 0 & x <= 5 & y <= x}",
 576                        "{[x,y]: y >= 0 & x = 6 & y <= 5}");
 577         set = isl_set_coalesce(set);
 578         assert(set && set->n == 2);
 579         isl_set_free(set);
 580
 581         set = s_union(ctx, "{[x,y]: y >= 0 & x <= 5 & y <= x}",
 582                        "{[x,y]: y >= 0 & x = 6 & y <= 7}");
 583         set = isl_set_coalesce(set);
 584         assert(set && set->n == 2);
 585         isl_set_free(set);
 586 }
 587
 588 int main()
 589 {
 590         struct isl_ctx *ctx;
 591
 592         srcdir = getenv("srcdir");
 593         assert(srcdir);
 594
 595         ctx = isl_ctx_alloc();
 596         test_read(ctx);
 597         test_construction(ctx);
 598         test_div(ctx);
 599         test_application(ctx);
 600         test_affine_hull(ctx);
 601         test_convex_hull(ctx);
 602         test_gist(ctx);
 603         test_coalesce(ctx);
 604         isl_ctx_free(ctx);
 605         return 0;
 606 }