Merge branch 'master' of https://github.com/dopefishh/mc1516pa
[mc1516pa.git] / modelchecker / main.c
1 #include <stdio.h>
2 #include <stdbool.h>
3 #include <ctype.h>
4 #include <time.h>
5 #include <unistd.h>
6
7 #include <sylvan.h>
8
9 #include "sokoban.h"
10 #include "coord.h"
11 #include "object.h"
12
13 #define ERRPRINT(fmt, as...) fprintf(stderr, fmt, ## as);
14 #define DPRINT(fmt, as...) if(DEBUG) ERRPRINT(fmt, ## as);
15
16 typedef enum {OBJECT, COORD, HYBRID} strategy;
17
18 //Global variables
19 bool DEBUG = false;
20 strategy strat = HYBRID;
21
22 void usage(char *prg)
23 {
24 ERRPRINT("Usage:\n"
25 "\t%s [opts] [FILE [FILE [...]]]\n"
26 "\n"
27 "Options:\n"
28 "\tAll strategies are mutually exclusive\n"
29 "\t-c coordinate based strategy\n"
30 "\t-o object based strategy\n"
31 "\t-y hybrid strategy\n"
32 // "\t-l LURD lURD verification strategy\n"
33 // "\t-r show all positions that are a valid solution\n"
34 "\n"
35 "\t-d enable verbose debug output\n"
36 "\t-h show this help\n"
37 "\n"
38 "Positional arguments:\n"
39 "\tFILE zero or more sokoban screens\n"
40 "\t when no file is specified stdin will be used\n", prg);
41 }
42
43 void solve(FILE *inputstream)
44 {
45 clock_t time_start_read, time_end_read;
46 clock_t time_start_encode, time_end_encode;
47
48 time_start_read = clock();
49 sokoban_screen *screen = parse_screen(inputstream, false);
50 if (screen == NULL) printf("Something went wrong...\n");
51 //sokoban_print(screen);
52 time_end_read = clock();
53
54 time_start_encode = clock();
55
56 lace_init(0, 1000000);
57 lace_startup(0, NULL, NULL);
58 LACE_ME;
59 sylvan_init_package(1LL<<21, 1LL<<27, 1LL<<20, 1LL<<26);
60 sylvan_init_bdd(6);
61
62 state *init = encode_screen(screen);
63 rels *rls = encode_rel(screen);
64
65 BDD old = sylvan_false;
66 BDD new = init->bdd;
67 int iteration = 0;
68 while(new != old){
69 old = new;
70 ERRPRINT("Iteration %d\n", iteration++);
71 trans_t *t = rls->rell;
72
73 while (t != NULL){
74 new = sylvan_or(new, sylvan_relnext(new, t->bdd, t->varset.varset));
75 t = t->next_rel;
76 }
77 t = rls->relu;
78 while (t != NULL){
79 new = sylvan_or(new, sylvan_relnext(new, t->bdd, t->varset.varset));
80 t = t->next_rel;
81 }
82 t = rls->relr;
83 while (t != NULL){
84 new = sylvan_or(new, sylvan_relnext(new, t->bdd, t->varset.varset));
85 t = t->next_rel;
86 }
87 t = rls->reld;
88 while (t != NULL){
89 new = sylvan_or(new, sylvan_relnext(new, t->bdd, t->varset.varset));
90 t = t->next_rel;
91 }
92 }
93 ERRPRINT("Satcount: %f\n", sylvan_satcount(old, init->vars.varset));
94 // sylvan_enum(old, init->vars.varset, (enum_cb)callback, NULL);
95 // sylvan_printdot_nc(old);
96 time_end_encode = clock();
97
98 //SOLVE???
99
100 sokoban_free(screen);
101 ERRPRINT("Reading: %fs\n",
102 ((double) (time_end_read-time_start_read))/CLOCKS_PER_SEC);
103 ERRPRINT("Encoding: %fs\n",
104 ((double) (time_end_encode-time_start_encode))/CLOCKS_PER_SEC);
105 }
106
107 int main(int argc, char **argv)
108 {
109 int optchar;
110
111 while((optchar = getopt(argc, argv, "cdhoy")) != -1){
112 switch(optchar){
113 case 'c':
114 strat = COORD;
115 DPRINT("Strategy changed to Coordinate based\n");
116 break;
117 case 'd':
118 DEBUG = true;
119 DPRINT("Debug enabled\n");
120 break;
121 case 'h':
122 usage(argv[0]);
123 return 0;
124 case 'o':
125 strat = OBJECT;
126 DPRINT("Strategy changed to Object based\n");
127 break;
128 case 'y':
129 strat = HYBRID;
130 DPRINT("Strategy changed to Hybrid\n");
131 break;
132 case '?':
133 if(isprint(optopt)){
134 ERRPRINT("Unknown option `-%c'.\n", optopt);
135 } else {
136 ERRPRINT("Unknown option char `-\\x%x'.\n", optopt);
137 }
138 return 2;
139 default:
140 break;
141 }
142 }
143
144 if(optind == argc){
145 ERRPRINT("You have not specified a file, reading from stdin\n");
146 solve(stdin);
147 }
148
149 for(int filepathindex = optind; filepathindex < argc; filepathindex++){
150 char *currentfilepath = argv[filepathindex];
151 ERRPRINT("Processing: %s\n", currentfilepath);
152 FILE *currentfile = fopen(currentfilepath, "r");
153 DPRINT("Opening file\n");
154 solve(currentfile);
155 DPRINT("Closing file\n");
156 fclose(currentfile);
157 }
158 return 0;
159 }