X-Git-Url: https://git.martlubbers.net/?a=blobdiff_plain;f=modelchecker%2Fmain.c;h=bcf8c7b6106bb1ebda1dc159d2722d3fe2452e57;hb=b3851097adbafa74d01e05b25ba550a814a5767d;hp=5e892113d2f8d5815b016859baac0ef67e96bb4e;hpb=7b58c6ef6de2382ed29f3e618d94b07a4b4f5791;p=mc1516pa.git diff --git a/modelchecker/main.c b/modelchecker/main.c index 5e89211..bcf8c7b 100644 --- a/modelchecker/main.c +++ b/modelchecker/main.c @@ -1,8 +1,196 @@ #include +#include +#include +#include +#include #include -int main(void){ - printf("Hello world!\n"); +#include "sokoban.h" +#include "coord.h" +#include "object.h" + +#define ERRPRINT(fmt, as...) fprintf(stderr, fmt, ## as); +#define DPRINT(fmt, as...) if(DEBUG) ERRPRINT(fmt, ## as); + +typedef enum {OBJECT, COORD, HYBRID} strategy; + +//Global variables +bool DEBUG = false; +strategy strat = HYBRID; + +void usage(char *prg) +{ + ERRPRINT("Usage:\n" + "\t%s [opts] [FILE [FILE [...]]]\n" + "\n" + "Options:\n" + "\tAll strategies are mutually exclusive\n" + "\t-c coordinate based strategy\n" + "\t-o object based strategy\n" + "\t-y hybrid strategy\n" +// "\t-l LURD lURD verification strategy\n" +// "\t-r show all positions that are a valid solution\n" + "\n" + "\t-d enable verbose debug output\n" + "\t-h show this help\n" + "\n" + "Positional arguments:\n" + "\tFILE zero or more sokoban screens\n" + "\t when no file is specified stdin will be used\n", prg); +} + +void solve(FILE *inputstream) +{ + clock_t time_start_read, time_end_read; + clock_t time_start_encode, time_end_encode; + + time_start_read = clock(); + sokoban_screen *screen = parse_screen(inputstream, false); + if (screen == NULL) printf("Something went wrong...\n"); + //sokoban_print(screen); + time_end_read = clock(); + + time_start_encode = clock(); + + lace_init(0, 1000000); + lace_startup(0, NULL, NULL); + LACE_ME; + sylvan_init_package(1LL<<21, 1LL<<27, 1LL<<20, 1LL<<26); + sylvan_init_bdd(6); + + state *init = encode_screen(screen); + state *goal = encode_goal(screen); + rels *rls = encode_rel(screen); + /* + printf("Doing left\n"); + test_trans(init, rls->rell); + printf("Doing up\n"); + test_trans(init, rls->relu); + printf("Doing right\n"); + test_trans(init, rls->relr); + printf("Doing down\n"); + test_trans(init, rls->reld); + */ + + BDD old = sylvan_false; + BDD new = init->bdd; + int iteration = 0; + if (goal != NULL){ + while(new != old){ + old = new; + ERRPRINT("Iteration %d\n", iteration++); + trans_t *t = rls->rell; + + //for testing + /* + BDDSET testvars = sylvan_set_fromarray(((BDDVAR[]){0,2,4,6,8,10,12,14,16}), 9); + BDD teststate = sylvan_cube(testvars, (uint8_t[]){1,0,1,0,0,1,0,1,0}); + BDD teststate1 = sylvan_or(teststate, sylvan_cube(testvars, (uint8_t[]){0,0,1,1,0,1,0,1,0})); + BDD teststate2 = sylvan_or(teststate1, sylvan_cube(testvars, (uint8_t[]){0,1,0,0,0,1,1,0,1})); + if (new == teststate2) printf("Wrong!\n"); + */ + while (t != NULL){ + new = sylvan_or(new, sylvan_relnext(new, t->bdd, t->varset.varset)); + t = t->next_rel; + //ERRPRINT("Satcount: %f\n", sylvan_satcount(new, init->vars.varset)); + } + + t = rls->relu; + while (t != NULL){ + new = sylvan_or(new, sylvan_relnext(new, t->bdd, t->varset.varset)); + t = t->next_rel; + //ERRPRINT("Satcount: %f\n", sylvan_satcount(new, init->vars.varset)); + } + + + t = rls->relr; + while (t != NULL){ + new = sylvan_or(new, sylvan_relnext(new, t->bdd, t->varset.varset)); + t = t->next_rel; + //ERRPRINT("Satcount: %f\n", sylvan_satcount(new, init->vars.varset)); + } + + + t = rls->reld; + while (t != NULL){ + new = sylvan_or(new, sylvan_relnext(new, t->bdd, t->varset.varset)); + t = t->next_rel; + //ERRPRINT("Satcount: %f\n", sylvan_satcount(new, init->vars.varset)); + } + double check = sylvan_satcount(sylvan_and(goal->bdd, new), init->vars.varset); + if (check > 0){ + printf("Solution is found after %d iteration(s)!\n", iteration); + break; + } + } + printf("Solution is not found after %d iteration(s)\n", iteration); + } + else printf("No solution found!\n"); + //ERRPRINT("Satcount: %f\n", sylvan_satcount(old, init->vars.varset)); +// sylvan_enum(old, init->vars.varset, (enum_cb)callback, NULL); +// sylvan_printdot_nc(old); + time_end_encode = clock(); + + //SOLVE??? + + sokoban_free(screen); + ERRPRINT("Reading: %fs\n", + ((double) (time_end_read-time_start_read))/CLOCKS_PER_SEC); + ERRPRINT("Encoding: %fs\n", + ((double) (time_end_encode-time_start_encode))/CLOCKS_PER_SEC); +} + +int main(int argc, char **argv) +{ + int optchar; + + while((optchar = getopt(argc, argv, "cdhoy")) != -1){ + switch(optchar){ + case 'c': + strat = COORD; + DPRINT("Strategy changed to Coordinate based\n"); + break; + case 'd': + DEBUG = true; + DPRINT("Debug enabled\n"); + break; + case 'h': + usage(argv[0]); + return 0; + case 'o': + strat = OBJECT; + DPRINT("Strategy changed to Object based\n"); + break; + case 'y': + strat = HYBRID; + DPRINT("Strategy changed to Hybrid\n"); + break; + case '?': + if(isprint(optopt)){ + ERRPRINT("Unknown option `-%c'.\n", optopt); + } else { + ERRPRINT("Unknown option char `-\\x%x'.\n", optopt); + } + return 2; + default: + break; + } + } + + if(optind == argc){ + ERRPRINT("You have not specified a file, reading from stdin\n"); + solve(stdin); + } + + for(int filepathindex = optind; filepathindex < argc; filepathindex++){ + char *currentfilepath = argv[filepathindex]; + ERRPRINT("Processing: %s\n", currentfilepath); + FILE *currentfile = fopen(currentfilepath, "r"); + DPRINT("Opening file\n"); + solve(currentfile); + DPRINT("Closing file\n"); + fclose(currentfile); + } return 0; }