some more trans var corrections
[mc1516pa.git] / modelchecker / coord.c
index b586629..3cd3a9a 100644 (file)
 #include "coord.h"
 #include "sokoban.h"
 
-/*
- * Global TODOs for now:
- * - update the helper maps so it is a single bimap
- * - make a data structure for relations
- * - encode relations as a set of four relations for each move
- */
 
 typedef struct {
        int x;
@@ -50,13 +44,6 @@ typedef struct {
        bddvar_xy_map *t;
 } bimap;
 
-typedef struct {
-       BDD *rell;
-       BDD *relu;
-       BDD *relr;
-       BDD *reld;
-} rels;
-
 xy_bddvar_map *getxy(int x, int y, xy_bddvar_map *map)
 {
        xy_bddvar_map k, *r = NULL;
@@ -76,51 +63,6 @@ bddvar_xy_map *getbdd(int key, bddvar_xy_map *map)
        return r;
 }
 
-/*
-xy_bddvar_map *create_xy_bddvar_map(sokoban_screen *screen)
-{
-       int varcount = 0;
-       sokoban_screen *r;
-       xy_bddvar_map *xybdd = NULL;
-       for(r=screen; r != NULL; r = (sokoban_screen *)(r->hh.next)){
-               xy_bddvar_map *f = NULL;
-               //bddvar_xy_map *t = NULL;
-
-               f = (xy_bddvar_map *)malloc(sizeof(xy_bddvar_map));
-               memset(f, 0, sizeof(xy_bddvar_map));
-               f->key.x = r->coord.x;
-               f->key.y = r->coord.y;
-               f->value.var[0] = varcount;
-               f->value.var[1] = varcount + 1;
-               f->value.var[2] = varcount + 2;
-               HASH_ADD(hh, xybdd, key, sizeof(xy), f);
-               varcount = varcount + 3;
-       }
-       return xybdd;
-}
-
-bddvar_xy_map *create_bddvar_xy_map(sokoban_screen *screen)
-{
-       int varcount = 0;
-       sokoban_screen *r;
-       bddvar_xy_map *bddxy = NULL;
-       for(r=screen; r != NULL; r = (sokoban_screen *)(r->hh.next)){
-               for (int i = 0; i <3; i++){
-                       bddvar_xy_map *f = NULL;
-                       //bddvar_xy_map *t = NULL;
-
-                       f = (bddvar_xy_map *)malloc(sizeof(bddvar_xy_map));
-                       memset(f, 0, sizeof(bddvar_xy_map));
-                       f->key = varcount + i;
-                       f->value.x = r->coord.x;
-                       f->value.y = r->coord.y;
-                       HASH_ADD(hh, bddxy, key, sizeof(int), f);
-               }
-               varcount = varcount + 3;
-       }
-       return bddxy;
-}
-*/
 
 bimap *create_bimap_helper(sokoban_screen *screen)
 {
@@ -276,227 +218,1104 @@ state *encode_screen(sokoban_screen *screen)
        return fullState;
 }
 
-int
-test_relprod()
-{
-    LACE_ME;
-
-    BDDVAR vars[] = {0,2,4};
-    BDDVAR all_vars[] = {4,5};
-
-    BDDSET vars_set = sylvan_set_fromarray(vars, 3);
-    BDDSET all_vars_set = sylvan_set_fromarray(all_vars, 2);
-
-    BDD s, t, next, prev;
-    BDD zeroes, ones;
-
-    // transition relation: 000 --> 111 and !000 --> 000
-    t = sylvan_false;
-    //t = sylvan_union_cube(t, all_vars_set, ((uint8_t[]){0,1,0,1,0,1}));
-    //t = sylvan_union_cube(t, all_vars_set, ((uint8_t[]){1,0,2,0,2,0}));
-    //t = sylvan_union_cube(t, all_vars_set, ((uint8_t[]){2,0,1,0,2,0}));
-    t = sylvan_union_cube(t, all_vars_set, ((uint8_t[]){1,0}));
-
-    s = sylvan_cube(vars_set, (uint8_t[]){0,0,1});
-    zeroes = sylvan_cube(vars_set, (uint8_t[]){0,0,0});
-    ones = sylvan_cube(vars_set, (uint8_t[]){0,0,1});
-
-    next = sylvan_relnext(s, t, all_vars_set);
-    prev = sylvan_relprev(t, next, all_vars_set);
-    if (next == zeroes) printf("Pass 1\n");
-    if (prev == ones || prev == zeroes) printf("Pass 2\n");
-
-    return 0;
-}
-
-BDD encode_rel(sokoban_screen *screen)
+trans_t *create_single_rel(sokoban_screen *screen, direction dir)
 {
        LACE_ME;
-       //      int tile_count = 0;
-       BDD t = sylvan_false;
-       //BDDVAR relvars[];
-       //BDDSET relvarset;
-       //uint8_t rel_enc[];
-
+       trans_t *trs, *trs_current;
        bimap *bm = create_bimap_helper(screen);
        int x = 0;
        int y = 0;
        bddvar_xy_map *bddxy = NULL;
+       int xdelta = 0;
+       int ydelta = 0;
+       int xgamma = 0;
+       int ygamma = 0;
+       trs = NULL;
+       switch(dir){
+       case LEFT:
+               xdelta = -1;
+               ydelta = 0;
+               xgamma = -2;
+               ygamma = 0;
+               break;
+       case UP:
+               xdelta = 0;
+               ydelta = -1;
+               xgamma = 0;
+               ygamma = -2;
+               break;
+       case RIGHT:
+               xdelta = 1;
+               ydelta = 0;
+               xgamma = 2;
+               ygamma = 0;
+               break;
+       case DOWN:
+               xdelta = 0;
+               ydelta = 1;
+               xgamma = 0;
+               ygamma = 2;
+               break;
+       }
 
-       //left relation
        for (unsigned int i = 0; i < HASH_COUNT(screen); i++){
-               printf("i:%d\n", i);
                bddxy = getbdd(i*3, bm->t);
                x = bddxy->value.x;
                y = bddxy->value.y;
-               if (check_space(x, y, LEFT, 1, bm) == 0){
-                       sokoban_screen *tmp_scr = get_coord(x, y, screen);
-                       if (tmp_scr->tile == AGENT){
-                               //1 1 0 0 1 1 (1 0 1 -> 1 0 1)
-                               BDDVAR relvars[6]  = {i * 3, i * 3 + 1, i * 3 + 2, i * 3 + 3, i * 3 + 4, i * 3 + 5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 6);
-                               uint8_t rel_enc[6] = {1, 1, 0, 0, 1, 1};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
-                       if (tmp_scr->tile == TARGAGENT){
-                               //1 1 1 1 0 0 (1 1 0 -> 1 1 0)
-                               BDDVAR relvars[6]  = {i * 3, i * 3 + 1, i * 3 + 2, i * 3 + 3, i * 3 + 4, i * 3 + 5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 6);
-                               uint8_t rel_enc[6] = {1, 1, 1, 1, 0, 0};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
+               if (check_space(x, y, dir, 1, bm) == 0){
+
+            //Agent -> Agent
+                       //1 1 0 0 1 1 (1 0 1 -> 1 0 1)
+                       BDDVAR relvars[6] = {i * 3, i * 3 + 1, i * 3 + 2, i * 3 + 3, i * 3 + 4, i * 3 + 5};
+                       BDDSET relvarset = sylvan_set_fromarray(relvars, 6);
+                       uint8_t rel_enc[6] = {1, 1, 0, 0, 1, 1};
+
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 6;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Targagent -> Targagent
+                       //1 1 1 1 0 0 (1 1 0 -> 1 1 0)
+                       relvars[0] = i * 3;
+                       relvars[1] = i * 3 + 1;
+                       relvars[2] = i * 3 + 2;
+                       relvars[3] = i * 3 + 3;
+                       relvars[4] = i * 3 + 4;
+                       relvars[5] = i * 3 + 5;
+                       relvarset = sylvan_set_fromarray(relvars, 6);
+                       rel_enc[0] = 1;
+                       rel_enc[1] = 1;
+                       rel_enc[2] = 1;
+                       rel_enc[3] = 1;
+                       rel_enc[4] = 0;
+                       rel_enc[5] = 0;
+
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 6;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
                }
-               else if (check_space(x, y, LEFT, 1, bm) == 1 && check_space(x, y, LEFT, 2, bm) == 0){
-                       sokoban_screen *tmp_scr = get_coord(x, y, screen);
-                       sokoban_screen *tmp_scr_d = get_coord(x - 1, y, screen);
-                       xy_bddvar_map *bddvar = getxy(x - 1, y, bm->f);
+               else if (check_space(x, y, dir, 1, bm) == 1 && check_space(x, y, dir, 2, bm) == 0){
+
+                       xy_bddvar_map *bddvar = getxy(x + xdelta, y + ydelta, bm->f);
                        int deltai = bddvar->value.var[0];
-                       if (tmp_scr->tile == AGENT && tmp_scr_d->tile == FREE){
-                               //(0 1 0 0 1 1 1 0 0 0 1 1)
-                               BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
-                               uint8_t rel_enc[12] = {0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
-                       if (tmp_scr->tile == AGENT && tmp_scr_d->tile == TARGET){
-                               //(0 1 1 0 1 1 1 0 0 1 1 1)
-                               BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
-                               uint8_t rel_enc[12] = {0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
-                       if (tmp_scr->tile == TARGAGENT && tmp_scr_d->tile == FREE){
-                               //(0 1 0 1 1 0 1 0 1 0 0 1)
-                               BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
-                               uint8_t rel_enc[12] = {0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0, 1};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
-                       if (tmp_scr->tile == TARGAGENT && tmp_scr_d->tile == TARGET){
-                               //(0 1 1 1 1 0 1 0 1 1 0 1)
-                               BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
-                               uint8_t rel_enc[12] = {0, 1, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
-                       if (tmp_scr->tile == AGENT && tmp_scr_d->tile == BOX){
-                               //(0 0 1 1 0 0 1 1 0 0 1 1)
-                               BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
-                               uint8_t rel_enc[12] = {0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
-                       if (tmp_scr->tile == AGENT && tmp_scr_d->tile == TARGBOX){
-                               //(1 1 0 0 0 0 1 1 0 0 1 1)
-                               BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
-                               uint8_t rel_enc[12] = {1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
-                       if (tmp_scr->tile == TARGAGENT && tmp_scr_d->tile == BOX){
-                               //(0 0 1 1 0 0 1 1 1 1 0 0)
-                               BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
-                               uint8_t rel_enc[12] = {0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 0, 0};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
-                       if (tmp_scr->tile == TARGAGENT && tmp_scr_d->tile == TARGBOX){
-                               //(1 1 0 0 0 0 1 1 1 1 0 0)
-                               BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
-                               uint8_t rel_enc[12] = {1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
+
+            //Free Agent -> Agent Free
+                       //(0 1 0 0 1 1 1 0 0 0 1 1)
+                       BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
+                       BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
+                       uint8_t rel_enc[12] = {0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1};
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 12;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Target Agent -> Targagent Free
+                       //(0 1 1 1 1 0 1 0 0 0 1 1)
+                       relvars[0] = deltai*3;
+                       relvars[1] = deltai*3+1;
+                       relvars[2] = deltai*3+2;
+                       relvars[3] = deltai*3+3;
+                       relvars[4] = deltai*3+4;
+                       relvars[5] = deltai*3+5;
+                       relvars[6] = i*3;
+                       relvars[7] = i*3+1;
+                       relvars[8] = i*3+2;
+                       relvars[9] = i*3+3;
+                       relvars[10] = i*3+4;
+                       relvars[11] = i*3+5;
+                       relvarset = sylvan_set_fromarray(relvars, 12);
+                       rel_enc[0] = 0;
+                       rel_enc[1] = 1;
+                       rel_enc[2] = 1;
+                       rel_enc[3] = 1;
+                       rel_enc[4] = 1;
+                       rel_enc[5] = 0;
+                       rel_enc[6] = 1;
+                       rel_enc[7] = 0;
+                       rel_enc[8] = 0;
+                       rel_enc[9] = 0;
+                       rel_enc[10] = 1;
+                       rel_enc[11] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 12;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Must be: Free Targagent -> Agent Target
+                       //(0 1 0 0 1 1 1 0 1 1 0 1)
+                       relvars[0] = deltai*3;
+                       relvars[1] = deltai*3+1;
+                       relvars[2] = deltai*3+2;
+                       relvars[3] = deltai*3+3;
+                       relvars[4] = deltai*3+4;
+                       relvars[5] = deltai*3+5;
+                       relvars[6] = i*3;
+                       relvars[7] = i*3+1;
+                       relvars[8] = i*3+2;
+                       relvars[9] = i*3+3;
+                       relvars[10] = i*3+4;
+                       relvars[11] = i*3+5;
+                       relvarset = sylvan_set_fromarray(relvars, 12);
+                       rel_enc[0] = 0;
+                       rel_enc[1] = 1;
+                       rel_enc[2] = 0;
+                       rel_enc[3] = 0;
+                       rel_enc[4] = 1;
+                       rel_enc[5] = 1;
+                       rel_enc[6] = 1;
+                       rel_enc[7] = 0;
+                       rel_enc[8] = 1;
+                       rel_enc[9] = 1;
+                       rel_enc[10] = 0;
+                       rel_enc[11] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 12;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Target Targagent -> Targagent Target
+                       //(0 1 1 1 1 0 1 0 1 1 0 1)
+                       relvars[0] = deltai*3;
+                       relvars[1] = deltai*3+1;
+                       relvars[2] = deltai*3+2;
+                       relvars[3] = deltai*3+3;
+                       relvars[4] = deltai*3+4;
+                       relvars[5] = deltai*3+5;
+                       relvars[6] = i*3;
+                       relvars[7] = i*3+1;
+                       relvars[8] = i*3+2;
+                       relvars[9] = i*3+3;
+                       relvars[10] = i*3+4;
+                       relvars[11] = i*3+5;
+                       relvarset = sylvan_set_fromarray(relvars, 12);
+                       rel_enc[0] = 0;
+                       rel_enc[1] = 1;
+                       rel_enc[2] = 1;
+                       rel_enc[3] = 1;
+                       rel_enc[4] = 1;
+                       rel_enc[5] = 0;
+                       rel_enc[6] = 1;
+                       rel_enc[7] = 0;
+                       rel_enc[8] = 1;
+                       rel_enc[9] = 1;
+                       rel_enc[10] = 0;
+                       rel_enc[11] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 12;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Box Agent -> Box Agent
+                       //(0 0 1 1 0 0 1 1 0 0 1 1)
+                       relvars[0] = deltai*3;
+                       relvars[1] = deltai*3+1;
+                       relvars[2] = deltai*3+2;
+                       relvars[3] = deltai*3+3;
+                       relvars[4] = deltai*3+4;
+                       relvars[5] = deltai*3+5;
+                       relvars[6] = i*3;
+                       relvars[7] = i*3+1;
+                       relvars[8] = i*3+2;
+                       relvars[9] = i*3+3;
+                       relvars[10] = i*3+4;
+                       relvars[11] = i*3+5;
+                       relvarset = sylvan_set_fromarray(relvars, 12);
+                       rel_enc[0] = 0;
+                       rel_enc[1] = 0;
+                       rel_enc[2] = 1;
+                       rel_enc[3] = 1;
+                       rel_enc[4] = 0;
+                       rel_enc[5] = 0;
+                       rel_enc[6] = 1;
+                       rel_enc[7] = 1;
+                       rel_enc[8] = 0;
+                       rel_enc[9] = 0;
+                       rel_enc[10] = 1;
+                       rel_enc[11] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 12;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Targbox Agent -> Targbox Agent
+                       //(1 1 0 0 0 0 1 1 0 0 1 1)
+                       relvars[0] = deltai*3;
+                       relvars[1] = deltai*3+1;
+                       relvars[2] = deltai*3+2;
+                       relvars[3] = deltai*3+3;
+                       relvars[4] = deltai*3+4;
+                       relvars[5] = deltai*3+5;
+                       relvars[6] = i*3;
+                       relvars[7] = i*3+1;
+                       relvars[8] = i*3+2;
+                       relvars[9] = i*3+3;
+                       relvars[10] = i*3+4;
+                       relvars[11] = i*3+5;
+                       relvarset = sylvan_set_fromarray(relvars, 12);
+                       rel_enc[0] = 1;
+                       rel_enc[1] = 1;
+                       rel_enc[2] = 0;
+                       rel_enc[3] = 0;
+                       rel_enc[4] = 0;
+                       rel_enc[5] = 0;
+                       rel_enc[6] = 1;
+                       rel_enc[7] = 1;
+                       rel_enc[8] = 0;
+                       rel_enc[9] = 0;
+                       rel_enc[10] = 1;
+                       rel_enc[11] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 12;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Box Targagent -> Box Targagent
+                       //(0 0 1 1 0 0 1 1 1 1 0 0)
+                       relvars[0] = deltai*3;
+                       relvars[1] = deltai*3+1;
+                       relvars[2] = deltai*3+2;
+                       relvars[3] = deltai*3+3;
+                       relvars[4] = deltai*3+4;
+                       relvars[5] = deltai*3+5;
+                       relvars[6] = i*3;
+                       relvars[7] = i*3+1;
+                       relvars[8] = i*3+2;
+                       relvars[9] = i*3+3;
+                       relvars[10] = i*3+4;
+                       relvars[11] = i*3+5;
+                       relvarset = sylvan_set_fromarray(relvars, 12);
+                       rel_enc[0] = 0;
+                       rel_enc[1] = 0;
+                       rel_enc[2] = 1;
+                       rel_enc[3] = 1;
+                       rel_enc[4] = 0;
+                       rel_enc[5] = 0;
+                       rel_enc[6] = 1;
+                       rel_enc[7] = 1;
+                       rel_enc[8] = 1;
+                       rel_enc[9] = 1;
+                       rel_enc[10] = 0;
+                       rel_enc[11] = 0;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 12;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Targbox Targagent -> Targbox Targagent
+                       //(1 1 0 0 0 0 1 1 1 1 0 0)
+                       relvars[0] = deltai*3;
+                       relvars[1] = deltai*3+1;
+                       relvars[2] = deltai*3+2;
+                       relvars[3] = deltai*3+3;
+                       relvars[4] = deltai*3+4;
+                       relvars[5] = deltai*3+5;
+                       relvars[6] = i*3;
+                       relvars[7] = i*3+1;
+                       relvars[8] = i*3+2;
+                       relvars[9] = i*3+3;
+                       relvars[10] = i*3+4;
+                       relvars[11] = i*3+5;
+                       relvarset = sylvan_set_fromarray(relvars, 12);
+                       rel_enc[0] = 1;
+                       rel_enc[1] = 1;
+                       rel_enc[2] = 0;
+                       rel_enc[3] = 0;
+                       rel_enc[4] = 0;
+                       rel_enc[5] = 0;
+                       rel_enc[6] = 1;
+                       rel_enc[7] = 1;
+                       rel_enc[8] = 1;
+                       rel_enc[9] = 1;
+                       rel_enc[10] = 0;
+                       rel_enc[11] = 0;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 12;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
                }
-               else if (check_space(x, y, LEFT, 1, bm) == 1 && check_space(x, y, LEFT, 2, bm) == 1){
-                       sokoban_screen *tmp_scr = get_coord(x, y, screen);
-                       sokoban_screen *tmp_scr_d = get_coord(x - 1, y, screen);
-                       sokoban_screen *tmp_scr_g = get_coord(x - 2, y, screen);
-                       xy_bddvar_map *bddvar = getxy(x - 1, y, bm->f);
+               else if (check_space(x, y, dir, 1, bm) == 1 && check_space(x, y, dir, 2, bm) == 1){
+
+                       xy_bddvar_map *bddvar = getxy(x + xdelta, y + ydelta, bm->f);
                        int deltai = bddvar->value.var[0];
-                       bddvar = getxy(x - 2, y, bm->f);
+                       bddvar = getxy(x + xgamma, y + ygamma, bm->f);
                        int gammai = bddvar->value.var[0];
-                       if (tmp_scr->tile == AGENT && tmp_scr_d->tile == FREE){
-                               //(0 1 0 0 1 1 1 0 0 0 1 1)
-                               BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
-                               uint8_t rel_enc[12] = {0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
-                       if (tmp_scr->tile == AGENT && tmp_scr_d->tile == TARGET){
-                               //(0 1 1 0 1 1 1 0 0 1 1 1)
-                               BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
-                               uint8_t rel_enc[12] = {0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
-                       if (tmp_scr->tile == TARGAGENT && tmp_scr_d->tile == FREE){
-                               //(0 1 0 1 1 0 1 0 1 0 0 1)
-                               BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
-                               uint8_t rel_enc[12] = {0, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0, 1};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
-                       if (tmp_scr->tile == TARGAGENT && tmp_scr_d->tile == TARGET){
-                               //(0 1 1 1 1 0 1 0 1 1 0 1)
-                               BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
-                               BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
-                               uint8_t rel_enc[12] = {0, 1, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1};
-                               t = sylvan_union_cube(t, relvarset, rel_enc);
-                       }
-                       if (tmp_scr->tile == AGENT && tmp_scr_d-> tile == BOX && tmp_scr_g->tile == BOX){
-                               //(0 0 1 1 0 0 0 0 1 1 0 0 1 1 0 0 1 1)
-                       }
-                       if (tmp_scr->tile == AGENT && tmp_scr_d-> tile == BOX && tmp_scr_g->tile == TARGBOX){
-                               //(1 1 0 0 0 0 0 0 1 1 0 0 1 1 0 0 1 1)
-                       }
-                       if (tmp_scr->tile == TARGAGENT && tmp_scr_d-> tile == BOX && tmp_scr_g->tile == BOX){
-                               //(0 0 1 1 0 0 0 0 1 1 0 0 1 1 1 1 0 0)
-                       }
-                       if (tmp_scr->tile == TARGAGENT && tmp_scr_d-> tile == BOX && tmp_scr_g->tile == TARGBOX){
-                               //(1 1 0 0 0 0 0 0 1 1 0 0 1 1 1 1 0 0)
-                       }
-                       if (tmp_scr->tile == AGENT && tmp_scr_d->tile == BOX && tmp_scr_g->tile == FREE){
-                               //free box agent -> box agent free
-                               //(0 0 0 1 1 0 0 1 1 0 0 1 1 0 0 0 1 1)
-                       }
-                       if (tmp_scr->tile == AGENT && tmp_scr_d->tile == TARGBOX && tmp_scr_g->tile == FREE){
-                               //(free targbox agent -> box targagent free)
-                               //(0 0 0 1 1 0 1 1 0 1 0 0 1 0 0 0 1 1)
-                       }
-                       if (tmp_scr->tile == AGENT && tmp_scr_d->tile == BOX && tmp_scr_g->tile == TARGET){
-                               //(target box agent -> targbox agent free)
-                               //(0 1 1 0 1 0 0 1 1 0 0 1 1 0 0 0 1 1)
-                       }
-                       if (tmp_scr->tile == AGENT && tmp_scr_d->tile == TARGBOX && tmp_scr_g->tile == TARGET){
-                               //(target targbox agent -> targbox targagent free)
-                               //(0 1 1 0 1 0 1 1 0 1 0 0 1 0 0 0 1 1)
-                       }
-                       if (tmp_scr->tile == TARGAGENT && tmp_scr_d->tile == BOX && tmp_scr_g->tile == FREE){
-                               //(free box targagent -> box agent target)
-                               //(0 0 0 1 1 0 0 1 1 0 0 1 1 0 1 1 0 1)
-                       }
-                       if (tmp_scr->tile == TARGAGENT && tmp_scr_d->tile == TARGBOX && tmp_scr_g->tile == FREE){
-                               //(free targbox targagent -> box targagent target)
-                               //(0 0 0 1 1 0 1 1 0 1 0 0 1 0 1 1 0 1)
-                       }
-                       if (tmp_scr->tile == TARGAGENT && tmp_scr_d->tile == BOX && tmp_scr_g->tile == TARGET){
-                               //(target box targagent -> targbox agent target)
-                               //(0 1 1 0 1 0 0 1 1 0 0 1 1 0 1 1 0 1)
-                       }
-                       if (tmp_scr->tile == TARGAGENT && tmp_scr_d->tile == TARGBOX && tmp_scr_g->tile == TARGET){
-                               //(target targbox targagent -> targbox targagent target)
-                               //(0 1 1 0 1 0 1 1 0 1 0 0 1 0 1 1 0 1)
-                       }
+
+            //Free Agent -> 1 0 1 0 0 1 Agent Free
+                       //(0 1 0 0 1 1 1 0 0 0 1 1)
+                       BDDVAR relvars[12]  = {deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
+                       BDDSET relvarset = sylvan_set_fromarray(relvars, 12);
+                       uint8_t rel_enc[12] = {0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1};
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 12;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Target Agent -> Targagent Free
+                       //(0 1 1 1 1 0 1 0 0 0 1 1)
+                       relvars[0] = deltai*3;
+                       relvars[1] = deltai*3+1;
+                       relvars[2] = deltai*3+2;
+                       relvars[3] = deltai*3+3;
+                       relvars[4] = deltai*3+4;
+                       relvars[5] = deltai*3+5;
+                       relvars[6] = i*3;
+                       relvars[7] = i*3+1;
+                       relvars[8] = i*3+2;
+                       relvars[9] = i*3+3;
+                       relvars[10] = i*3+4;
+                       relvars[11] = i*3+5;
+                       relvarset = sylvan_set_fromarray(relvars, 12);
+                       rel_enc[0] = 0;
+                       rel_enc[1] = 1;
+                       rel_enc[2] = 1;
+                       rel_enc[3] = 1;
+                       rel_enc[4] = 1;
+                       rel_enc[5] = 0;
+                       rel_enc[6] = 1;
+                       rel_enc[7] = 0;
+                       rel_enc[8] = 0;
+                       rel_enc[9] = 0;
+                       rel_enc[10] = 1;
+                       rel_enc[11] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 12;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Free Targagent -> Agent Target
+                       //(0 1 0 0 1 1 1 0 1 1 0 1)
+                       relvars[0] = deltai*3;
+                       relvars[1] = deltai*3+1;
+                       relvars[2] = deltai*3+2;
+                       relvars[3] = deltai*3+3;
+                       relvars[4] = deltai*3+4;
+                       relvars[5] = deltai*3+5;
+                       relvars[6] = i*3;
+                       relvars[7] = i*3+1;
+                       relvars[8] = i*3+2;
+                       relvars[9] = i*3+3;
+                       relvars[10] = i*3+4;
+                       relvars[11] = i*3+5;
+                       relvarset = sylvan_set_fromarray(relvars, 12);
+                       rel_enc[0] = 0;
+                       rel_enc[1] = 1;
+                       rel_enc[2] = 0;
+                       rel_enc[3] = 0;
+                       rel_enc[4] = 1;
+                       rel_enc[5] = 1;
+                       rel_enc[6] = 1;
+                       rel_enc[7] = 0;
+                       rel_enc[8] = 1;
+                       rel_enc[9] = 1;
+                       rel_enc[10] = 0;
+                       rel_enc[11] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 12;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Target Targagent -> Targagent Target
+                       //(0 1 1 1 1 0 1 0 1 1 0 1)
+                       relvars[0] = deltai*3;
+                       relvars[1] = deltai*3+1;
+                       relvars[2] = deltai*3+2;
+                       relvars[3] = deltai*3+3;
+                       relvars[4] = deltai*3+4;
+                       relvars[5] = deltai*3+5;
+                       relvars[6] = i*3;
+                       relvars[7] = i*3+1;
+                       relvars[8] = i*3+2;
+                       relvars[9] = i*3+3;
+                       relvars[10] = i*3+4;
+                       relvars[11] = i*3+5;
+                       relvarset = sylvan_set_fromarray(relvars, 12);
+                       rel_enc[0] = 0;
+                       rel_enc[1] = 1;
+                       rel_enc[2] = 1;
+                       rel_enc[3] = 1;
+                       rel_enc[4] = 1;
+                       rel_enc[5] = 0;
+                       rel_enc[6] = 1;
+                       rel_enc[7] = 0;
+                       rel_enc[8] = 1;
+                       rel_enc[9] = 1;
+                       rel_enc[10] = 0;
+                       rel_enc[11] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset, rel_enc);
+                       trs_current->varset.varset = relvarset;
+                       trs_current->varset.size = 12;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Box Box Agent -> Box Box Agent
+                       //(0 0 1 1 0 0 0 0 1 1 0 0 1 1 0 0 1 1)
+                       BDDVAR relvars1[18]  = {gammai*3, gammai*3+1, gammai*3+2, gammai*3+3, gammai*3+4, gammai*3+5, deltai*3, deltai*3+1, deltai*3+2, deltai*3+3, deltai*3+4, deltai*3+5, i*3, i*3+1, i*3+2, i*3+3, i*3+4, i*3+5};
+                       relvars1[0] = gammai*3;
+                       relvars1[1] = gammai*3+1;
+                       relvars1[2] = gammai*3+2;
+                       relvars1[3] = gammai*3+3;
+                       relvars1[4] = gammai*3+4;
+                       relvars1[5] = gammai*3+5;
+                       relvars1[6] = deltai*3;
+                       relvars1[7] = deltai*3+1;
+                       relvars1[8] = deltai*3+2;
+                       relvars1[9] = deltai*3+3;
+                       relvars1[10] = deltai*3+4;
+                       relvars1[11] = deltai*3+5;
+                       relvars1[12] = i*3;
+                       relvars1[13] = i*3+1;
+                       relvars1[14] = i*3+2;
+                       relvars1[15] = i*3+3;
+                       relvars1[16] = i*3+4;
+                       relvars1[17] = i*3+5;
+                       BDD relvarset1 = sylvan_set_fromarray(relvars1, 18);
+                       uint8_t rel_enc1[18];
+                       //(0 0 1 1 0 0 0 0 1 1 0 0 1 1 0 0 1 1)
+                       rel_enc1[0] = 0;
+                       rel_enc1[1] = 0;
+                       rel_enc1[2] = 1;
+                       rel_enc1[3] = 1;
+                       rel_enc1[4] = 0;
+                       rel_enc1[5] = 0;
+                       rel_enc1[6] = 0;
+                       rel_enc1[7] = 0;
+                       rel_enc1[8] = 1;
+                       rel_enc1[9] = 1;
+                       rel_enc1[10] = 0;
+                       rel_enc1[11] = 0;
+                       rel_enc1[12] = 1;
+                       rel_enc1[13] = 1;
+                       rel_enc1[14] = 0;
+                       rel_enc1[15] = 0;
+                       rel_enc1[16] = 1;
+                       rel_enc1[17] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset1, rel_enc1);
+                       trs_current->varset.varset = relvarset1;
+                       trs_current->varset.size = 18;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Targbox Box Agent -> Targbox Box Agent
+                       //(1 1 0 0 0 0 0 0 1 1 0 0 1 1 0 0 1 1)
+                       relvars1[0] = gammai*3;
+                       relvars1[1] = gammai*3+1;
+                       relvars1[2] = gammai*3+2;
+                       relvars1[3] = gammai*3+3;
+                       relvars1[4] = gammai*3+4;
+                       relvars1[5] = gammai*3+5;
+                       relvars1[6] = deltai*3;
+                       relvars1[7] = deltai*3+1;
+                       relvars1[8] = deltai*3+2;
+                       relvars1[9] = deltai*3+3;
+                       relvars1[10] = deltai*3+4;
+                       relvars1[11] = deltai*3+5;
+                       relvars1[12] = i*3;
+                       relvars1[13] = i*3+1;
+                       relvars1[14] = i*3+2;
+                       relvars1[15] = i*3+3;
+                       relvars1[16] = i*3+4;
+                       relvars1[17] = i*3+5;
+                       relvarset1 = sylvan_set_fromarray(relvars1, 18);
+                       //(1 1 0 0 0 0 0 0 1 1 0 0 1 1 0 0 1 1)
+                       rel_enc1[0] = 1;
+                       rel_enc1[1] = 1;
+                       rel_enc1[2] = 0;
+                       rel_enc1[3] = 0;
+                       rel_enc1[4] = 0;
+                       rel_enc1[5] = 0;
+                       rel_enc1[6] = 0;
+                       rel_enc1[7] = 0;
+                       rel_enc1[8] = 1;
+                       rel_enc1[9] = 1;
+                       rel_enc1[10] = 0;
+                       rel_enc1[11] = 0;
+                       rel_enc1[12] = 1;
+                       rel_enc1[13] = 1;
+                       rel_enc1[14] = 0;
+                       rel_enc1[15] = 0;
+                       rel_enc1[16] = 1;
+                       rel_enc1[17] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset1, rel_enc1);
+                       trs_current->varset.varset = relvarset1;
+                       trs_current->varset.size = 18;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Box Box Targagent -> Box Box Targagent
+                       //(0 0 1 1 0 0 0 0 1 1 0 0 1 1 1 1 0 0)
+                       relvars1[0] = gammai*3;
+                       relvars1[1] = gammai*3+1;
+                       relvars1[2] = gammai*3+2;
+                       relvars1[3] = gammai*3+3;
+                       relvars1[4] = gammai*3+4;
+                       relvars1[5] = gammai*3+5;
+                       relvars1[6] = deltai*3;
+                       relvars1[7] = deltai*3+1;
+                       relvars1[8] = deltai*3+2;
+                       relvars1[9] = deltai*3+3;
+                       relvars1[10] = deltai*3+4;
+                       relvars1[11] = deltai*3+5;
+                       relvars1[12] = i*3;
+                       relvars1[13] = i*3+1;
+                       relvars1[14] = i*3+2;
+                       relvars1[15] = i*3+3;
+                       relvars1[16] = i*3+4;
+                       relvars1[17] = i*3+5;
+                       relvarset1 = sylvan_set_fromarray(relvars1, 18);
+                       //(0 0 1 1 0 0 0 0 1 1 0 0 1 1 1 1 0 0)
+                       rel_enc1[0] = 0;
+                       rel_enc1[1] = 0;
+                       rel_enc1[2] = 1;
+                       rel_enc1[3] = 1;
+                       rel_enc1[4] = 0;
+                       rel_enc1[5] = 0;
+                       rel_enc1[6] = 0;
+                       rel_enc1[7] = 0;
+                       rel_enc1[8] = 1;
+                       rel_enc1[9] = 1;
+                       rel_enc1[10] = 0;
+                       rel_enc1[11] = 0;
+                       rel_enc1[12] = 1;
+                       rel_enc1[13] = 1;
+                       rel_enc1[14] = 1;
+                       rel_enc1[15] = 1;
+                       rel_enc1[16] = 0;
+                       rel_enc1[17] = 0;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset1, rel_enc1);
+                       trs_current->varset.varset = relvarset1;
+                       trs_current->varset.size = 18;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+            //Targbox Box Targagent -> Targbox Box Targagent
+                       //(1 1 0 0 0 0 0 0 1 1 0 0 1 1 1 1 0 0)
+                       relvars1[0] = gammai*3;
+                       relvars1[1] = gammai*3+1;
+                       relvars1[2] = gammai*3+2;
+                       relvars1[3] = gammai*3+3;
+                       relvars1[4] = gammai*3+4;
+                       relvars1[5] = gammai*3+5;
+                       relvars1[6] = deltai*3;
+                       relvars1[7] = deltai*3+1;
+                       relvars1[8] = deltai*3+2;
+                       relvars1[9] = deltai*3+3;
+                       relvars1[10] = deltai*3+4;
+                       relvars1[11] = deltai*3+5;
+                       relvars1[12] = i*3;
+                       relvars1[13] = i*3+1;
+                       relvars1[14] = i*3+2;
+                       relvars1[15] = i*3+3;
+                       relvars1[16] = i*3+4;
+                       relvars1[17] = i*3+5;
+                       relvarset1 = sylvan_set_fromarray(relvars1, 18);
+                       //(1 1 0 0 0 0 0 0 1 1 0 0 1 1 1 1 0 0)
+                       rel_enc1[0] = 1;
+                       rel_enc1[1] = 1;
+                       rel_enc1[2] = 0;
+                       rel_enc1[3] = 0;
+                       rel_enc1[4] = 0;
+                       rel_enc1[5] = 0;
+                       rel_enc1[6] = 0;
+                       rel_enc1[7] = 0;
+                       rel_enc1[8] = 1;
+                       rel_enc1[9] = 1;
+                       rel_enc1[10] = 0;
+                       rel_enc1[11] = 0;
+                       rel_enc1[12] = 1;
+                       rel_enc1[13] = 1;
+                       rel_enc1[14] = 1;
+                       rel_enc1[15] = 1;
+                       rel_enc1[16] = 0;
+                       rel_enc1[17] = 0;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset1, rel_enc1);
+                       trs_current->varset.varset = relvarset1;
+                       trs_current->varset.size = 18;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+                       //Free Box Agent -> Box Agent Free
+                       //(0 0 0 1 1 0 0 1 1 0 0 1 1 0 0 0 1 1)
+                       relvars1[0] = gammai*3;
+                       relvars1[1] = gammai*3+1;
+                       relvars1[2] = gammai*3+2;
+                       relvars1[3] = gammai*3+3;
+                       relvars1[4] = gammai*3+4;
+                       relvars1[5] = gammai*3+5;
+                       relvars1[6] = deltai*3;
+                       relvars1[7] = deltai*3+1;
+                       relvars1[8] = deltai*3+2;
+                       relvars1[9] = deltai*3+3;
+                       relvars1[10] = deltai*3+4;
+                       relvars1[11] = deltai*3+5;
+                       relvars1[12] = i*3;
+                       relvars1[13] = i*3+1;
+                       relvars1[14] = i*3+2;
+                       relvars1[15] = i*3+3;
+                       relvars1[16] = i*3+4;
+                       relvars1[17] = i*3+5;
+                       relvarset1 = sylvan_set_fromarray(relvars1, 18);
+                       //(0 0 0 1 1 0 0 1 1 0 0 1 1 0 0 0 1 1)
+                       rel_enc1[0] = 0;
+                       rel_enc1[1] = 0;
+                       rel_enc1[2] = 0;
+                       rel_enc1[3] = 1;
+                       rel_enc1[4] = 1;
+                       rel_enc1[5] = 0;
+                       rel_enc1[6] = 0;
+                       rel_enc1[7] = 1;
+                       rel_enc1[8] = 1;
+                       rel_enc1[9] = 0;
+                       rel_enc1[10] = 0;
+                       rel_enc1[11] = 1;
+                       rel_enc1[12] = 1;
+                       rel_enc1[13] = 0;
+                       rel_enc1[14] = 0;
+                       rel_enc1[15] = 0;
+                       rel_enc1[16] = 1;
+                       rel_enc1[17] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset1, rel_enc1);
+                       trs_current->varset.varset = relvarset1;
+                       trs_current->varset.size = 18;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+                       //Free Targbox Agent -> Box Targagent Free
+                       //(0 0 0 1 1 0 1 1 0 1 0 0 1 0 0 0 1 1)
+                       relvars1[0] = gammai*3;
+                       relvars1[1] = gammai*3+1;
+                       relvars1[2] = gammai*3+2;
+                       relvars1[3] = gammai*3+3;
+                       relvars1[4] = gammai*3+4;
+                       relvars1[5] = gammai*3+5;
+                       relvars1[6] = deltai*3;
+                       relvars1[7] = deltai*3+1;
+                       relvars1[8] = deltai*3+2;
+                       relvars1[9] = deltai*3+3;
+                       relvars1[10] = deltai*3+4;
+                       relvars1[11] = deltai*3+5;
+                       relvars1[12] = i*3;
+                       relvars1[13] = i*3+1;
+                       relvars1[14] = i*3+2;
+                       relvars1[15] = i*3+3;
+                       relvars1[16] = i*3+4;
+                       relvars1[17] = i*3+5;
+                       relvarset1 = sylvan_set_fromarray(relvars1, 18);
+                       //(0 0 0 1 1 0 1 1 0 1 0 0 1 0 0 0 1 1)
+                       rel_enc1[0] = 0;
+                       rel_enc1[1] = 0;
+                       rel_enc1[2] = 0;
+                       rel_enc1[3] = 1;
+                       rel_enc1[4] = 1;
+                       rel_enc1[5] = 0;
+                       rel_enc1[6] = 1;
+                       rel_enc1[7] = 1;
+                       rel_enc1[8] = 0;
+                       rel_enc1[9] = 1;
+                       rel_enc1[10] = 0;
+                       rel_enc1[11] = 0;
+                       rel_enc1[12] = 1;
+                       rel_enc1[13] = 0;
+                       rel_enc1[14] = 0;
+                       rel_enc1[15] = 0;
+                       rel_enc1[16] = 1;
+                       rel_enc1[17] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset1, rel_enc1);
+                       trs_current->varset.varset = relvarset1;
+                       trs_current->varset.size = 18;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+                       //Target Box Agent -> Targbox Agent Free
+                       //(0 1 1 0 1 0 0 1 1 0 0 1 1 0 0 0 1 1)
+                       relvars1[0] = gammai*3;
+                       relvars1[1] = gammai*3+1;
+                       relvars1[2] = gammai*3+2;
+                       relvars1[3] = gammai*3+3;
+                       relvars1[4] = gammai*3+4;
+                       relvars1[5] = gammai*3+5;
+                       relvars1[6] = deltai*3;
+                       relvars1[7] = deltai*3+1;
+                       relvars1[8] = deltai*3+2;
+                       relvars1[9] = deltai*3+3;
+                       relvars1[10] = deltai*3+4;
+                       relvars1[11] = deltai*3+5;
+                       relvars1[12] = i*3;
+                       relvars1[13] = i*3+1;
+                       relvars1[14] = i*3+2;
+                       relvars1[15] = i*3+3;
+                       relvars1[16] = i*3+4;
+                       relvars1[17] = i*3+5;
+                       relvarset1 = sylvan_set_fromarray(relvars1, 18);
+                       //(0 1 1 0 1 0 0 1 1 0 0 1 1 0 0 0 1 1)
+                       rel_enc1[0] = 0;
+                       rel_enc1[1] = 1;
+                       rel_enc1[2] = 1;
+                       rel_enc1[3] = 0;
+                       rel_enc1[4] = 1;
+                       rel_enc1[5] = 0;
+                       rel_enc1[6] = 0;
+                       rel_enc1[7] = 1;
+                       rel_enc1[8] = 1;
+                       rel_enc1[9] = 0;
+                       rel_enc1[10] = 0;
+                       rel_enc1[11] = 1;
+                       rel_enc1[12] = 1;
+                       rel_enc1[13] = 0;
+                       rel_enc1[14] = 0;
+                       rel_enc1[15] = 0;
+                       rel_enc1[16] = 1;
+                       rel_enc1[17] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset1, rel_enc1);
+                       trs_current->varset.varset = relvarset1;
+                       trs_current->varset.size = 18;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+                       //Target Targbox Agent -> Targbox Targagent Free
+                       //(0 1 1 0 1 0 1 1 0 1 0 0 1 0 0 0 1 1)
+                       relvars1[0] = gammai*3;
+                       relvars1[1] = gammai*3+1;
+                       relvars1[2] = gammai*3+2;
+                       relvars1[3] = gammai*3+3;
+                       relvars1[4] = gammai*3+4;
+                       relvars1[5] = gammai*3+5;
+                       relvars1[6] = deltai*3;
+                       relvars1[7] = deltai*3+1;
+                       relvars1[8] = deltai*3+2;
+                       relvars1[9] = deltai*3+3;
+                       relvars1[10] = deltai*3+4;
+                       relvars1[11] = deltai*3+5;
+                       relvars1[12] = i*3;
+                       relvars1[13] = i*3+1;
+                       relvars1[14] = i*3+2;
+                       relvars1[15] = i*3+3;
+                       relvars1[16] = i*3+4;
+                       relvars1[17] = i*3+5;
+                       relvarset1 = sylvan_set_fromarray(relvars1, 18);
+                       //(0 1 1 0 1 0 1 1 0 1 0 0 1 0 0 0 1 1)
+                       rel_enc1[0] = 0;
+                       rel_enc1[1] = 1;
+                       rel_enc1[2] = 1;
+                       rel_enc1[3] = 0;
+                       rel_enc1[4] = 1;
+                       rel_enc1[5] = 0;
+                       rel_enc1[6] = 1;
+                       rel_enc1[7] = 1;
+                       rel_enc1[8] = 0;
+                       rel_enc1[9] = 1;
+                       rel_enc1[10] = 0;
+                       rel_enc1[11] = 0;
+                       rel_enc1[12] = 1;
+                       rel_enc1[13] = 0;
+                       rel_enc1[14] = 0;
+                       rel_enc1[15] = 0;
+                       rel_enc1[16] = 1;
+                       rel_enc1[17] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset1, rel_enc1);
+                       trs_current->varset.varset = relvarset1;
+                       trs_current->varset.size = 18;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+                       //Free Box Targagent -> Box Agent Target
+                       //(0 0 0 1 1 0 0 1 1 0 0 1 1 0 1 1 0 1)
+                       relvars1[0] = gammai*3;
+                       relvars1[1] = gammai*3+1;
+                       relvars1[2] = gammai*3+2;
+                       relvars1[3] = gammai*3+3;
+                       relvars1[4] = gammai*3+4;
+                       relvars1[5] = gammai*3+5;
+                       relvars1[6] = deltai*3;
+                       relvars1[7] = deltai*3+1;
+                       relvars1[8] = deltai*3+2;
+                       relvars1[9] = deltai*3+3;
+                       relvars1[10] = deltai*3+4;
+                       relvars1[11] = deltai*3+5;
+                       relvars1[12] = i*3;
+                       relvars1[13] = i*3+1;
+                       relvars1[14] = i*3+2;
+                       relvars1[15] = i*3+3;
+                       relvars1[16] = i*3+4;
+                       relvars1[17] = i*3+5;
+                       relvarset1 = sylvan_set_fromarray(relvars1, 18);
+                       //(0 0 0 1 1 0 0 1 1 0 0 1 1 0 1 1 0 1)
+                       rel_enc1[0] = 0;
+                       rel_enc1[1] = 0;
+                       rel_enc1[2] = 0;
+                       rel_enc1[3] = 1;
+                       rel_enc1[4] = 1;
+                       rel_enc1[5] = 0;
+                       rel_enc1[6] = 0;
+                       rel_enc1[7] = 1;
+                       rel_enc1[8] = 1;
+                       rel_enc1[9] = 0;
+                       rel_enc1[10] = 0;
+                       rel_enc1[11] = 1;
+                       rel_enc1[12] = 1;
+                       rel_enc1[13] = 0;
+                       rel_enc1[14] = 1;
+                       rel_enc1[15] = 1;
+                       rel_enc1[16] = 0;
+                       rel_enc1[17] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset1, rel_enc1);
+                       trs_current->varset.varset = relvarset1;
+                       trs_current->varset.size = 18;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+                       //Free Targbox Targagent -> Box Targagent Target
+                       //(0 0 0 1 1 0 1 1 0 1 0 0 1 0 1 1 0 1)
+                       relvars1[0] = gammai*3;
+                       relvars1[1] = gammai*3+1;
+                       relvars1[2] = gammai*3+2;
+                       relvars1[3] = gammai*3+3;
+                       relvars1[4] = gammai*3+4;
+                       relvars1[5] = gammai*3+5;
+                       relvars1[6] = deltai*3;
+                       relvars1[7] = deltai*3+1;
+                       relvars1[8] = deltai*3+2;
+                       relvars1[9] = deltai*3+3;
+                       relvars1[10] = deltai*3+4;
+                       relvars1[11] = deltai*3+5;
+                       relvars1[12] = i*3;
+                       relvars1[13] = i*3+1;
+                       relvars1[14] = i*3+2;
+                       relvars1[15] = i*3+3;
+                       relvars1[16] = i*3+4;
+                       relvars1[17] = i*3+5;
+                       relvarset1 = sylvan_set_fromarray(relvars1, 18);
+                       //(0 0 0 1 1 0 1 1 0 1 0 0 1 0 1 1 0 1)
+                       rel_enc1[0] = 0;
+                       rel_enc1[1] = 0;
+                       rel_enc1[2] = 0;
+                       rel_enc1[3] = 1;
+                       rel_enc1[4] = 1;
+                       rel_enc1[5] = 0;
+                       rel_enc1[6] = 1;
+                       rel_enc1[7] = 1;
+                       rel_enc1[8] = 0;
+                       rel_enc1[9] = 1;
+                       rel_enc1[10] = 0;
+                       rel_enc1[11] = 0;
+                       rel_enc1[12] = 1;
+                       rel_enc1[13] = 0;
+                       rel_enc1[14] = 1;
+                       rel_enc1[15] = 1;
+                       rel_enc1[16] = 0;
+                       rel_enc1[17] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset1, rel_enc1);
+                       trs_current->varset.varset = relvarset1;
+                       trs_current->varset.size = 18;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+                       //Target Box Targagent -> Targbox Agent Target
+                       //(0 1 1 0 1 0 0 1 1 0 0 1 1 0 1 1 0 1)
+                       relvars1[0] = gammai*3;
+                       relvars1[1] = gammai*3+1;
+                       relvars1[2] = gammai*3+2;
+                       relvars1[3] = gammai*3+3;
+                       relvars1[4] = gammai*3+4;
+                       relvars1[5] = gammai*3+5;
+                       relvars1[6] = deltai*3;
+                       relvars1[7] = deltai*3+1;
+                       relvars1[8] = deltai*3+2;
+                       relvars1[9] = deltai*3+3;
+                       relvars1[10] = deltai*3+4;
+                       relvars1[11] = deltai*3+5;
+                       relvars1[12] = i*3;
+                       relvars1[13] = i*3+1;
+                       relvars1[14] = i*3+2;
+                       relvars1[15] = i*3+3;
+                       relvars1[16] = i*3+4;
+                       relvars1[17] = i*3+5;
+                       relvarset1 = sylvan_set_fromarray(relvars1, 18);
+                       //(0 1 1 0 1 0 0 1 1 0 0 1 1 0 1 1 0 1)
+                       rel_enc1[0] = 0;
+                       rel_enc1[1] = 1;
+                       rel_enc1[2] = 1;
+                       rel_enc1[3] = 0;
+                       rel_enc1[4] = 1;
+                       rel_enc1[5] = 0;
+                       rel_enc1[6] = 0;
+                       rel_enc1[7] = 1;
+                       rel_enc1[8] = 1;
+                       rel_enc1[9] = 0;
+                       rel_enc1[10] = 0;
+                       rel_enc1[11] = 1;
+                       rel_enc1[12] = 1;
+                       rel_enc1[13] = 0;
+                       rel_enc1[14] = 1;
+                       rel_enc1[15] = 1;
+                       rel_enc1[16] = 0;
+                       rel_enc1[17] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset1, rel_enc1);
+                       trs_current->varset.varset = relvarset1;
+                       trs_current->varset.size = 18;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
+
+                       //Target Targbox Targagent -> Targbox Targagent Target
+                       //(0 1 1 0 1 0 1 1 0 1 0 0 1 0 1 1 0 1)
+                       relvars1[0] = gammai*3;
+                       relvars1[1] = gammai*3+1;
+                       relvars1[2] = gammai*3+2;
+                       relvars1[3] = gammai*3+3;
+                       relvars1[4] = gammai*3+4;
+                       relvars1[5] = gammai*3+5;
+                       relvars1[6] = deltai*3;
+                       relvars1[7] = deltai*3+1;
+                       relvars1[8] = deltai*3+2;
+                       relvars1[9] = deltai*3+3;
+                       relvars1[10] = deltai*3+4;
+                       relvars1[11] = deltai*3+5;
+                       relvars1[12] = i*3;
+                       relvars1[13] = i*3+1;
+                       relvars1[14] = i*3+2;
+                       relvars1[15] = i*3+3;
+                       relvars1[16] = i*3+4;
+                       relvars1[17] = i*3+5;
+                       relvarset1 = sylvan_set_fromarray(relvars1, 18);
+                       //(0 1 1 0 1 0 1 1 0 1 0 0 1 0 1 1 0 1)
+                       rel_enc1[0] = 0;
+                       rel_enc1[1] = 1;
+                       rel_enc1[2] = 1;
+                       rel_enc1[3] = 0;
+                       rel_enc1[4] = 1;
+                       rel_enc1[5] = 0;
+                       rel_enc1[6] = 1;
+                       rel_enc1[7] = 1;
+                       rel_enc1[8] = 0;
+                       rel_enc1[9] = 1;
+                       rel_enc1[10] = 0;
+                       rel_enc1[11] = 0;
+                       rel_enc1[12] = 1;
+                       rel_enc1[13] = 0;
+                       rel_enc1[14] = 1;
+                       rel_enc1[15] = 1;
+                       rel_enc1[16] = 0;
+                       rel_enc1[17] = 1;
+                       trs_current = (trans_t *)malloc(sizeof(trans_t));
+                       trs_current->bdd = sylvan_cube(relvarset1, rel_enc1);
+                       trs_current->varset.varset = relvarset1;
+                       trs_current->varset.size = 18;
+                       trs_current->next_rel = trs;
+                       trs = trs_current;
                }
-               printf("x:%d y:%d\n", x, y);
        }
+       trs_current = trs;
+       return trs;
+}
+
+rels *encode_rel(sokoban_screen *screen)
+{
+       LACE_ME;
+
+       trans_t *tl = sylvan_false;
+
+       //left relation
+       tl = create_single_rel(screen, LEFT);
 
        //up relation
+       trans_t *tu = create_single_rel(screen, UP);
 
        //right relation
+       trans_t *tr = create_single_rel(screen, RIGHT);
 
        //down relation
+       trans_t *td = create_single_rel(screen, DOWN);
+
+       rels *rls = NULL;
+       rls = (rels *)malloc(sizeof(rels));
+       rls->rell = tl;
+       rls->relu = tu;
+       rls->relr = tr;
+       rls->reld = td;
+
+       return rls;
+}
+
+int
+test_relprod()
+{
+    LACE_ME;
+
+    BDDVAR vars[] = {0,2,4};
+       BDDVAR all_vars[] = {0,1,2,3,4,5};
+    BDDVAR all_vars2[] = {0,1};
+       //BDDVAR short_vars[] = {0,1};
+               /*BDDVAR short_vars2[] = {2,3};
+       BDDVAR short_vars3[] = {0,1,2,3};*/
 
-       return sylvan_true;
+    BDDSET vars_set = sylvan_set_fromarray(vars, 3);
+       BDDSET all_vars_set = sylvan_set_fromarray(all_vars, 6);
+    BDDSET all_vars_set2 = sylvan_set_fromarray(all_vars2, 2);
+       //BDDSET short_vars_set = sylvan_set_fromarray(short_vars, 2);
+               /*BDDSET short_vars_set2 = sylvan_set_fromarray(short_vars2, 2);
+       BDDSET short_vars_set3 = sylvan_set_fromarray(short_vars3, 4);*/
+
+    BDD s, t, next, prev;
+    BDD zeroes, ones;
+
+    // transition relation: 000 --> 111 and !000 --> 000
+    t = sylvan_false;
+       t = sylvan_union_cube(t, all_vars_set2, ((uint8_t[]){0,1}));
+    //t = sylvan_union_cube(t, all_vars_set, ((uint8_t[]){1,0,2,0,2,0}));
+    t = sylvan_union_cube(t, all_vars_set2, ((uint8_t[]){1,0}));
+       //    t = sylvan_union_cube(t, all_vars_set, ((uint8_t[]){0,1}));
+
+    s = sylvan_cube(vars_set, (uint8_t[]){0,0,0});
+    zeroes = sylvan_cube(vars_set, (uint8_t[]){1,0,0});
+    ones = sylvan_cube(vars_set, (uint8_t[]){0,0,0});
+
+    next = sylvan_relnext(s, t, all_vars_set);
+    prev = sylvan_relprev(t, next, all_vars_set);
+    if (next == zeroes) printf("Pass 1\n");
+    if (prev == ones) printf("Pass 2\n");
+       //trans *ts;
+       //ts = NULL;
+
+    return 0;
 }