glucose-syrup/incremental/Main.cc

   1
   2 #include "solver_interface.h"
   3 #include <errno.h>
   4
   5 #include <signal.h>
   6 #include <zlib.h>
   7 #include <sys/resource.h>
   8
   9 #include "utils/System.h"
  10 #include "utils/ParseUtils.h"
  11 #include "utils/Options.h"
  12 #include "core/Dimacs.h"
  13 #include "simp/SimpSolver.h"
  14
  15 using namespace Glucose;
  16
  17
  18 static const char* _certified = "CORE -- CERTIFIED UNSAT";
  19
  20 void printStats(Solver& solver)
  21 {
  22     double cpu_time = cpuTime();
  23     double mem_used = 0;//memUsedPeak();
  24     printf("c restarts              : %" PRIu64" (%" PRIu64" conflicts in avg)\n", solver.starts,(solver.starts>0 ?solver.conflicts/solver.starts : 0));
  25     printf("c blocked restarts      : %" PRIu64" (multiple: %" PRIu64") \n", solver.nbstopsrestarts,solver.nbstopsrestartssame);
  26     printf("c last block at restart : %" PRIu64"\n",solver.lastblockatrestart);
  27     printf("c nb ReduceDB           : %" PRIu64"\n", solver.nbReduceDB);
  28     printf("c nb removed Clauses    : %" PRIu64"\n",solver.nbRemovedClauses);
  29     printf("c nb learnts DL2        : %" PRIu64"\n", solver.nbDL2);
  30     printf("c nb learnts size 2     : %" PRIu64"\n", solver.nbBin);
  31     printf("c nb learnts size 1     : %" PRIu64"\n", solver.nbUn);
  32
  33     printf("c conflicts             : %-12" PRIu64"   (%.0f /sec)\n", solver.conflicts   , solver.conflicts   /cpu_time);
  34     printf("c decisions             : %-12" PRIu64"   (%4.2f %% random) (%.0f /sec)\n", solver.decisions, (float)solver.rnd_decisions*100 / (float)solver.decisions, solver.decisions   /cpu_time);
  35     printf("c propagations          : %-12" PRIu64"   (%.0f /sec)\n", solver.propagations, solver.propagations/cpu_time);
  36     printf("c conflict literals     : %-12" PRIu64"   (%4.2f %% deleted)\n", solver.tot_literals, (solver.max_literals - solver.tot_literals)*100 / (double)solver.max_literals);
  37     printf("c nb reduced Clauses    : %" PRIu64"\n",solver.nbReducedClauses);
  38
  39     if (mem_used != 0) printf("Memory used           : %.2f MB\n", mem_used);
  40     printf("c CPU time              : %g s\n", cpu_time);
  41 }
  42
  43
  44
  45 static Solver* solver;
  46 // Terminate by notifying the solver and back out gracefully. This is mainly to have a test-case
  47 // for this feature of the Solver as it may take longer than an immediate call to '_exit()'.
  48 static void SIGINT_interrupt(int signum) { solver->interrupt(); }
  49
  50 // Note that '_exit()' rather than 'exit()' has to be used. The reason is that 'exit()' calls
  51 // destructors and may cause deadlocks if a malloc/free function happens to be running (these
  52 // functions are guarded by locks for multithreaded use).
  53 static void SIGINT_exit(int signum) {
  54     printf("\n"); printf("*** INTERRUPTED ***\n");
  55     if (solver->verbosity > 0){
  56         printStats(*solver);
  57         printf("\n"); printf("*** INTERRUPTED ***\n"); }
  58     _exit(1); }
  59
  60
  61 int *buffer;
  62 int length;
  63 int offset;
  64
  65 int *outbuffer;
  66 int outoffset;
  67
  68 int getInt() {
  69   if (offset>=length) {
  70     offset = 0;
  71                 ssize_t ptr=read(0, buffer, sizeof(int)*IS_BUFFERSIZE);
  72                 if (ptr == -1 || ptr == 0)
  73                         exit(-1);
  74
  75     ssize_t bytestoread=(4-(ptr & 3)) & 3;
  76     while(bytestoread != 0) {
  77       ssize_t p=read(0, &((char *)buffer)[ptr], bytestoread);
  78       if (p == -1 || p == 0)
  79         exit(-1);
  80       bytestoread -= p;
  81       ptr += p;
  82     }
  83     length = ptr / 4;
  84     offset = 0;
  85   }
  86
  87   return buffer[offset++];
  88 }
  89
  90 void flushInts() {
  91   ssize_t bytestowrite=sizeof(int)*outoffset;
  92   ssize_t byteswritten=0;
  93   do {
  94     ssize_t n=write(IS_OUT_FD, &((char *)outbuffer)[byteswritten], bytestowrite);
  95     if (n == -1) {
  96       fprintf(stderr, "Write failure\n");
  97       exit(-1);
  98     }
  99     bytestowrite -= n;
 100     byteswritten += n;
 101   } while(bytestowrite != 0);
 102   outoffset = 0;
 103 }
 104
 105 void putInt(int value) {
 106   if (outoffset>=IS_BUFFERSIZE) {
 107     flushInts();
 108   }
 109   outbuffer[outoffset++]=value;
 110 }
 111
 112 void readClauses(Solver *solver) {
 113   vec<Lit> clause;
 114   int numvars = solver->nVars();
 115   bool haveClause = false;
 116   while(true) {
 117     int lit=getInt();
 118     if (lit!=0) {
 119       int var = abs(lit) - 1;
 120       while (var >= numvars) {
 121         numvars++;
 122         solver->newVar();
 123       }
 124       clause.push( (lit>0) ? mkLit(var) : ~mkLit(var));
 125       haveClause = true;
 126     } else {
 127       if (haveClause) {
 128         solver->addClause_(clause);
 129         haveClause = false;
 130         clause.clear();
 131       } else {
 132         //done with clauses
 133         return;
 134       }
 135     }
 136   }
 137 }
 138
 139 void processCommands(SimpSolver *solver) {
 140   while(true) {
 141     int command=getInt();
 142     switch(command) {
 143     case IS_FREEZE: {
 144       int var=getInt()-1;
 145       solver->setFrozen(var, true);
 146       break;
 147     }
 148     case IS_RUNSOLVER: {
 149       vec<Lit> dummy;
 150       lbool ret = solver->solveLimited(dummy);
 151       if (ret == l_True) {
 152         putInt(IS_SAT);
 153         putInt(solver->nVars());
 154         for(int i=0;i<solver->nVars();i++) {
 155           putInt(solver->model[i]==l_True);
 156         }
 157       } else if (ret == l_False) {
 158         putInt(IS_UNSAT);
 159       } else {
 160         putInt(IS_INDETER);
 161       }
 162       flushInts();
 163       return;
 164     }
 165     default:
 166       fprintf(stderr, "Unreconized command\n");
 167       exit(-1);
 168     }
 169   }
 170 }
 171
 172 void processSAT(SimpSolver *solver) {
 173   buffer=(int *) malloc(sizeof(int)*IS_BUFFERSIZE);
 174   offset=0;
 175   length=0;
 176   outbuffer=(int *) malloc(sizeof(int)*IS_BUFFERSIZE);
 177   outoffset=0;
 178
 179   while(true) {
 180     double initial_time = cpuTime();
 181     readClauses(solver);
 182     double parse_time = cpuTime();
 183     processCommands(solver);
 184     double finish_time = cpuTime();
 185     printf("Parse time: %12.2f s Solve time:%12.2f s\n", parse_time-initial_time, finish_time-parse_time);
 186   }
 187 }
 188
 189
 190
 191 //=================================================================================================
 192 // Main:
 193
 194 int main(int argc, char** argv) {
 195   try {
 196     printf("c\nc This is glucose 4.0 --  based on MiniSAT (Many thanks to MiniSAT team)\nc\n");
 197
 198
 199     setUsageHelp("c USAGE: %s [options] <input-file> <result-output-file>\n\n  where input may be either in plain or gzipped DIMACS.\n");
 200
 201
 202 #if defined(__linux__)
 203     fpu_control_t oldcw, newcw;
 204     _FPU_GETCW(oldcw); newcw = (oldcw & ~_FPU_EXTENDED) | _FPU_DOUBLE; _FPU_SETCW(newcw);
 205     //printf("c WARNING: for repeatability, setting FPU to use double precision\n");
 206 #endif
 207     // Extra options:
 208     //
 209     IntOption    verb   ("MAIN", "verb",   "Verbosity level (0=silent, 1=some, 2=more).", 0, IntRange(0, 2));
 210     BoolOption   mod   ("MAIN", "model",   "show model.", false);
 211     IntOption    vv  ("MAIN", "vv",   "Verbosity every vv conflicts", 10000, IntRange(1,INT32_MAX));
 212     BoolOption   pre    ("MAIN", "pre",    "Completely turn on/off any preprocessing.", true);
 213     StringOption dimacs ("MAIN", "dimacs", "If given, stop after preprocessing and write the result to this file.");
 214     IntOption    cpu_lim("MAIN", "cpu-lim","Limit on CPU time allowed in seconds.\n", INT32_MAX, IntRange(0, INT32_MAX));
 215     IntOption    mem_lim("MAIN", "mem-lim","Limit on memory usage in megabytes.\n", INT32_MAX, IntRange(0, INT32_MAX));
 216
 217     BoolOption    opt_certified      (_certified, "certified",    "Certified UNSAT using DRUP format", false);
 218     StringOption  opt_certified_file      (_certified, "certified-output",    "Certified UNSAT output file", "NULL");
 219
 220     parseOptions(argc, argv, true);
 221
 222     SimpSolver  S;
 223     S.parsing = 1;
 224     S.verbosity = verb;
 225     S.verbEveryConflicts = vv;
 226     S.showModel = mod;
 227     solver = &S;
 228     // Use signal handlers that forcibly quit until the solver will be
 229     // able to respond to interrupts:
 230     signal(SIGINT, SIGINT_exit);
 231     signal(SIGXCPU,SIGINT_exit);
 232
 233
 234     // Set limit on CPU-time:
 235     if (cpu_lim != INT32_MAX){
 236       rlimit rl;
 237       getrlimit(RLIMIT_CPU, &rl);
 238       if (rl.rlim_max == RLIM_INFINITY || (rlim_t)cpu_lim < rl.rlim_max){
 239         rl.rlim_cur = cpu_lim;
 240         if (setrlimit(RLIMIT_CPU, &rl) == -1)
 241           printf("c WARNING! Could not set resource limit: CPU-time.\n");
 242       } }
 243
 244     // Set limit on virtual memory:
 245     if (mem_lim != INT32_MAX){
 246       rlim_t new_mem_lim = (rlim_t)mem_lim * 1024*1024;
 247       rlimit rl;
 248       getrlimit(RLIMIT_AS, &rl);
 249       if (rl.rlim_max == RLIM_INFINITY || new_mem_lim < rl.rlim_max){
 250         rl.rlim_cur = new_mem_lim;
 251         if (setrlimit(RLIMIT_AS, &rl) == -1)
 252           printf("c WARNING! Could not set resource limit: Virtual memory.\n");
 253       } }
 254
 255     //do solver stuff here
 256     processSAT(&S);
 257
 258     printf("c |  Number of variables:  %12d                                                                   |\n", S.nVars());
 259     printf("c |  Number of clauses:    %12d                                                                   |\n", S.nClauses());
 260
 261
 262     // Change to signal-handlers that will only notify the solver and allow it to terminate
 263     // voluntarily:
 264     signal(SIGINT, SIGINT_interrupt);
 265     signal(SIGXCPU,SIGINT_interrupt);
 266
 267     S.parsing = 0;
 268
 269 #ifdef NDEBUG
 270     exit(0);
 271 #else
 272     return (0);
 273 #endif
 274   } catch (OutOfMemoryException&){
 275     printf("c =========================================================================================================\n");
 276     printf("INDETERMINATE\n");
 277     exit(0);
 278   }
 279 }