src/csolver.cc

   1 #include "csolver.h"
   2 #include "set.h"
   3 #include "mutableset.h"
   4 #include "element.h"
   5 #include "boolean.h"
   6 #include "predicate.h"
   7 #include "order.h"
   8 #include "table.h"
   9 #include "function.h"
  10 #include "satencoder.h"
  11 #include "sattranslator.h"
  12 #include "tunable.h"
  13 #include "polarityassignment.h"
  14 #include "transformer.h"
  15 #include "autotuner.h"
  16 #include "astops.h"
  17 #include "structs.h"
  18
  19 CSolver::CSolver() :
  20         boolTrue(new BooleanConst(true)),
  21         boolFalse(new BooleanConst(false)),
  22         unsat(false),
  23         tuner(NULL),
  24         elapsedTime(0)
  25 {
  26         satEncoder = new SATEncoder(this);
  27         transformer = new Transformer(this);
  28 }
  29
  30 /** This function tears down the solver and the entire AST */
  31
  32 CSolver::~CSolver() {
  33         uint size = allBooleans.getSize();
  34         for (uint i = 0; i < size; i++) {
  35                 delete allBooleans.get(i);
  36         }
  37
  38         size = allSets.getSize();
  39         for (uint i = 0; i < size; i++) {
  40                 delete allSets.get(i);
  41         }
  42
  43         size = allElements.getSize();
  44         for (uint i = 0; i < size; i++) {
  45                 delete allElements.get(i);
  46         }
  47
  48         size = allTables.getSize();
  49         for (uint i = 0; i < size; i++) {
  50                 delete allTables.get(i);
  51         }
  52
  53         size = allPredicates.getSize();
  54         for (uint i = 0; i < size; i++) {
  55                 delete allPredicates.get(i);
  56         }
  57
  58         size = allOrders.getSize();
  59         for (uint i = 0; i < size; i++) {
  60                 delete allOrders.get(i);
  61         }
  62
  63         size = allFunctions.getSize();
  64         for (uint i = 0; i < size; i++) {
  65                 delete allFunctions.get(i);
  66         }
  67
  68         delete boolTrue;
  69         delete boolFalse;
  70         delete satEncoder;
  71         delete transformer;
  72 }
  73
  74 CSolver *CSolver::clone() {
  75         CSolver *copy = new CSolver();
  76         CloneMap map;
  77         SetIteratorBoolean *it = getConstraints();
  78         while (it->hasNext()) {
  79                 Boolean *b = it->next();
  80                 copy->addConstraint(b->clone(copy, &map));
  81         }
  82         delete it;
  83         return copy;
  84 }
  85
  86 Set *CSolver::createSet(VarType type, uint64_t *elements, uint numelements) {
  87         Set *set = new Set(type, elements, numelements);
  88         allSets.push(set);
  89         return set;
  90 }
  91
  92 Set *CSolver::createRangeSet(VarType type, uint64_t lowrange, uint64_t highrange) {
  93         Set *set = new Set(type, lowrange, highrange);
  94         allSets.push(set);
  95         return set;
  96 }
  97
  98 Element *CSolver::createRangeVar(VarType type, uint64_t lowrange, uint64_t highrange) {
  99         Set *s = createRangeSet(type, lowrange, highrange);
 100         return getElementVar(s);
 101 }
 102
 103 MutableSet *CSolver::createMutableSet(VarType type) {
 104         MutableSet *set = new MutableSet(type);
 105         allSets.push(set);
 106         return set;
 107 }
 108
 109 void CSolver::addItem(MutableSet *set, uint64_t element) {
 110         set->addElementMSet(element);
 111 }
 112
 113 uint64_t CSolver::createUniqueItem(MutableSet *set) {
 114         uint64_t element = set->getNewUniqueItem();
 115         set->addElementMSet(element);
 116         return element;
 117 }
 118
 119 Element *CSolver::getElementVar(Set *set) {
 120         Element *element = new ElementSet(set);
 121         allElements.push(element);
 122         return element;
 123 }
 124
 125 Element *CSolver::getElementConst(VarType type, uint64_t value) {
 126         uint64_t array[] = {value};
 127         Set *set = new Set(type, array, 1);
 128         Element *element = new ElementConst(value, type, set);
 129         Element *e = elemMap.get(element);
 130         if (e == NULL) {
 131                 allSets.push(set);
 132                 allElements.push(element);
 133                 elemMap.put(element, element);
 134                 return element;
 135         } else {
 136                 delete set;
 137                 delete element;
 138                 return e;
 139         }
 140 }
 141
 142 Element *CSolver::applyFunction(Function *function, Element **array, uint numArrays, Boolean *overflowstatus) {
 143         Element *element = new ElementFunction(function,array,numArrays,overflowstatus);
 144         Element *e = elemMap.get(element);
 145         if (e == NULL) {
 146                 allElements.push(element);
 147                 elemMap.put(element, element);
 148                 return element;
 149         } else {
 150                 delete element;
 151                 return e;
 152         }
 153 }
 154
 155 Function *CSolver::createFunctionOperator(ArithOp op, Set **domain, uint numDomain, Set *range,OverFlowBehavior overflowbehavior) {
 156         Function *function = new FunctionOperator(op, domain, numDomain, range, overflowbehavior);
 157         allFunctions.push(function);
 158         return function;
 159 }
 160
 161 Predicate *CSolver::createPredicateOperator(CompOp op, Set **domain, uint numDomain) {
 162         Predicate *predicate = new PredicateOperator(op, domain,numDomain);
 163         allPredicates.push(predicate);
 164         return predicate;
 165 }
 166
 167 Predicate *CSolver::createPredicateTable(Table *table, UndefinedBehavior behavior) {
 168         Predicate *predicate = new PredicateTable(table, behavior);
 169         allPredicates.push(predicate);
 170         return predicate;
 171 }
 172
 173 Table *CSolver::createTable(Set **domains, uint numDomain, Set *range) {
 174         Table *table = new Table(domains,numDomain,range);
 175         allTables.push(table);
 176         return table;
 177 }
 178
 179 Table *CSolver::createTableForPredicate(Set **domains, uint numDomain) {
 180         return createTable(domains, numDomain, NULL);
 181 }
 182
 183 void CSolver::addTableEntry(Table *table, uint64_t *inputs, uint inputSize, uint64_t result) {
 184         table->addNewTableEntry(inputs, inputSize, result);
 185 }
 186
 187 Function *CSolver::completeTable(Table *table, UndefinedBehavior behavior) {
 188         Function *function = new FunctionTable(table, behavior);
 189         allFunctions.push(function);
 190         return function;
 191 }
 192
 193 Boolean *CSolver::getBooleanVar(VarType type) {
 194         Boolean *boolean = new BooleanVar(type);
 195         allBooleans.push(boolean);
 196         return boolean;
 197 }
 198
 199 Boolean *CSolver::getBooleanTrue() {
 200         return boolTrue;
 201 }
 202
 203 Boolean *CSolver::getBooleanFalse() {
 204         return boolFalse;
 205 }
 206
 207 Boolean *CSolver::applyPredicate(Predicate *predicate, Element **inputs, uint numInputs) {
 208         return applyPredicateTable(predicate, inputs, numInputs, NULL);
 209 }
 210
 211 Boolean *CSolver::applyPredicateTable(Predicate *predicate, Element **inputs, uint numInputs, Boolean *undefinedStatus) {
 212         BooleanPredicate *boolean = new BooleanPredicate(predicate, inputs, numInputs, undefinedStatus);
 213         Boolean * b = boolMap.get(boolean);
 214         if (b == NULL) {
 215                 boolMap.put(boolean, boolean);
 216                 allBooleans.push(boolean);
 217                 return boolean;
 218         } else {
 219                 delete boolean;
 220                 return b;
 221         }
 222 }
 223
 224 bool CSolver::isTrue(Boolean *b) {
 225         return b->isTrue();
 226 }
 227
 228 bool CSolver::isFalse(Boolean *b) {
 229         return b->isFalse();
 230 }
 231
 232 Boolean *CSolver::applyLogicalOperation(LogicOp op, Boolean * arg1, Boolean * arg2) {
 233         Boolean * array[] = {arg1, arg2};
 234         return applyLogicalOperation(op, array, 2);
 235 }
 236
 237 Boolean *CSolver::applyLogicalOperation(LogicOp op, Boolean *arg) {
 238         Boolean * array[] = {arg};
 239         return applyLogicalOperation(op, array, 1);
 240 }
 241
 242
 243 Boolean *CSolver::applyLogicalOperation(LogicOp op, Boolean **array, uint asize) {
 244         Boolean * newarray[asize];
 245         switch(op) {
 246         case SATC_NOT: {
 247                 if (array[0]->type == LOGICOP && ((BooleanLogic *)array[0])->op==SATC_NOT) {
 248                         return ((BooleanLogic *) array[0])->inputs.get(0);
 249                 } else if (array[0]->type == BOOLCONST) {
 250                         return array[0]->isTrue() ? boolFalse : boolTrue;
 251                 }
 252                 break;
 253         }
 254         case SATC_IFF: {
 255                 for(uint i=0;i<2;i++) {
 256                         if (array[i]->type == BOOLCONST) {
 257                                 if (array[i]->isTrue()) {
 258                                         return array[1-i];
 259                                 } else {
 260                                         newarray[0]=array[1-i];
 261                                         return applyLogicalOperation(SATC_NOT, newarray, 1);
 262                                 }
 263                         }
 264                 }
 265                 break;
 266         }
 267         case SATC_XOR: {
 268                 for(uint i=0;i<2;i++) {
 269                         if (array[i]->type == BOOLCONST) {
 270                                 if (array[i]->isTrue()) {
 271                                         newarray[0]=array[1-i];
 272                                         return applyLogicalOperation(SATC_NOT, newarray, 1);
 273                                 } else
 274                                         return array[1-i];
 275                         }
 276                 }
 277                 break;
 278         }
 279         case SATC_OR: {
 280                 uint newindex=0;
 281                 for(uint i=0;i<asize;i++) {
 282                         Boolean *b=array[i];
 283                         if (b->type == BOOLCONST) {
 284                                 if (b->isTrue())
 285                                         return b;
 286                                 else
 287                                         continue;
 288                         } else
 289                                 newarray[newindex++]=b;
 290                 }
 291                 if (newindex==1)
 292                         return newarray[0];
 293                 else if (newindex == 2) {
 294                         bool isNot0 = (newarray[0]->type==BOOLCONST) && ((BooleanLogic *)newarray[0])->op == SATC_NOT;
 295                         bool isNot1 = (newarray[1]->type==BOOLCONST) && ((BooleanLogic *)newarray[1])->op == SATC_NOT;
 296
 297                         if (isNot0 != isNot1) {
 298                                 if (isNot0) {
 299                                         newarray[0] = ((BooleanLogic *) newarray[0])->inputs.get(0);
 300                                 } else {
 301                                         Boolean *tmp =  ((BooleanLogic *) array[1])->inputs.get(0);
 302                                         array[1] = array[0];
 303                                         array[0] = tmp;
 304                                 }
 305                                 return applyLogicalOperation(SATC_IMPLIES, newarray, 2);
 306                         }
 307                 } else {
 308                         array = newarray;
 309                         asize = newindex;
 310                 }
 311                 break;
 312         }
 313         case SATC_AND: {
 314                 uint newindex=0;
 315                 for(uint i=0;i<asize;i++) {
 316                         Boolean *b=array[i];
 317                         if (b->type == BOOLCONST) {
 318                                 if (b->isTrue())
 319                                         continue;
 320                                 else
 321                                         return b;
 322                         } else
 323                                 newarray[newindex++]=b;
 324                 }
 325                 if(newindex==1) {
 326                         return newarray[0];
 327                 } else {
 328                         array = newarray;
 329                         asize = newindex;
 330                 }
 331                 break;
 332         }
 333         case SATC_IMPLIES: {
 334                 if (array[0]->type == BOOLCONST) {
 335                         if (array[0]->isTrue()) {
 336                                 return array[1];
 337                         } else {
 338                                 return boolTrue;
 339                         }
 340                 } else if (array[1]->type == BOOLCONST) {
 341                         if (array[1]->isTrue()) {
 342                                 return array[1];
 343                         } else {
 344                                 return applyLogicalOperation(SATC_NOT, array, 1);
 345                         }
 346                 }
 347                 break;
 348         }
 349         }
 350
 351         Boolean *boolean = new BooleanLogic(this, op, array, asize);
 352         Boolean *b = boolMap.get(boolean);
 353         if (b == NULL) {
 354                 boolMap.put(boolean, boolean);
 355                 allBooleans.push(boolean);
 356                 return boolean;
 357         } else {
 358                 delete boolean;
 359                 return b;
 360         }
 361 }
 362
 363 Boolean *CSolver::orderConstraint(Order *order, uint64_t first, uint64_t second) {
 364         Boolean *constraint = new BooleanOrder(order, first, second);
 365         allBooleans.push(constraint);
 366         return constraint;
 367 }
 368
 369 void CSolver::addConstraint(Boolean *constraint) {
 370         if (constraint == boolTrue)
 371                 return;
 372         else if (constraint == boolFalse)
 373                 setUnSAT();
 374         else
 375                 constraints.add(constraint);
 376 }
 377
 378 Order *CSolver::createOrder(OrderType type, Set *set) {
 379         Order *order = new Order(type, set);
 380         allOrders.push(order);
 381         return order;
 382 }
 383
 384 int CSolver::solve() {
 385         bool deleteTuner = false;
 386         if (tuner == NULL) {
 387                 tuner = new DefaultTuner();
 388                 deleteTuner = true;
 389         }
 390
 391         long long startTime = getTimeNano();
 392         computePolarities(this);
 393         transformer->orderAnalysis();
 394         naiveEncodingDecision(this);
 395         satEncoder->encodeAllSATEncoder(this);
 396         int result = unsat ? IS_UNSAT : satEncoder->solve();
 397         long long finishTime = getTimeNano();
 398         elapsedTime = finishTime - startTime;
 399         if (deleteTuner) {
 400                 delete tuner;
 401                 tuner = NULL;
 402         }
 403         return result;
 404 }
 405
 406 uint64_t CSolver::getElementValue(Element *element) {
 407         switch (element->type) {
 408         case ELEMSET:
 409         case ELEMCONST:
 410         case ELEMFUNCRETURN:
 411                 return getElementValueSATTranslator(this, element);
 412         default:
 413                 ASSERT(0);
 414         }
 415         exit(-1);
 416 }
 417
 418 bool CSolver::getBooleanValue(Boolean *boolean) {
 419         switch (boolean->type) {
 420         case BOOLEANVAR:
 421                 return getBooleanVariableValueSATTranslator(this, boolean);
 422         default:
 423                 ASSERT(0);
 424         }
 425         exit(-1);
 426 }
 427
 428 HappenedBefore CSolver::getOrderConstraintValue(Order *order, uint64_t first, uint64_t second) {
 429         return getOrderConstraintValueSATTranslator(this, order, first, second);
 430 }
 431
 432 long long CSolver::getEncodeTime() { return satEncoder->getEncodeTime(); }
 433
 434 long long CSolver::getSolveTime() { return satEncoder->getSolveTime(); }
 435
 436 void CSolver::autoTune(uint budget) {
 437         AutoTuner * autotuner=new AutoTuner(budget);
 438         autotuner->addProblem(this);
 439         autotuner->tune();
 440         delete autotuner;
 441 }