Refactoring ...
[satune.git] / src / Tuner / autotuner.cc
1 #include "autotuner.h"
2 #include "csolver.h"
3 #include "searchtuner.h"
4 #include <math.h>
5 #include <stdlib.h>
6 #include <float.h>
7
8 #define UNSETVALUE -1
9
10 AutoTuner::AutoTuner(uint _budget) :
11         budget(_budget), result(UNSETVALUE) {
12 }
13
14 void AutoTuner::addProblem(CSolver *solver) {
15         solvers.push(solver);
16 }
17
18 long long AutoTuner::evaluate(CSolver *problem, SearchTuner *tuner) {
19         CSolver *copy = problem->clone();
20         copy->setTuner(tuner);
21         model_print("**********************\n");
22         int sat = copy->solve();
23         if (result == UNSETVALUE)
24                 result = sat;
25         else if (result != sat) {
26                 model_print("&&&&&&&&&&&&&&&&&& Result has changed &&&&&&&&&&&&&\n");
27                 copy->printConstraints();
28         }
29         //model_print("SAT %d\n", result);
30         long long elapsedTime = copy->getElapsedTime();
31 //      long long encodeTime = copy->getEncodeTime();
32 //      long long solveTime = copy->getSolveTime();
33         long long metric = elapsedTime;
34 //      model_print("Elapsed Time: %llu\n", elapsedTime);
35 //      model_print("Encode Time: %llu\n", encodeTime);
36 //      model_print("Solve Time: %llu\n", solveTime);
37         delete copy;
38         return metric;
39 }
40
41 double AutoTuner::evaluateAll(SearchTuner *tuner) {
42         double product = 1;
43         for (uint i = 0; i < solvers.getSize(); i++) {
44                 CSolver *problem = solvers.get(i);
45                 double score = evaluate(problem, tuner);
46                 product *= score;
47         }
48         return pow(product, 1 / ((double)solvers.getSize()));
49 }
50
51 SearchTuner *AutoTuner::mutateTuner(SearchTuner *oldTuner, uint k) {
52         SearchTuner *newTuner = oldTuner->copyUsed();
53         uint numSettings = oldTuner->getSize();
54         uint settingsToMutate = (uint)(AUTOTUNERFACTOR * (((double)numSettings) * (budget - k)) / (budget));
55         if (settingsToMutate < 1)
56                 settingsToMutate = 1;
57         model_print("Mutating %u settings\n", settingsToMutate);
58         while (settingsToMutate-- != 0) {
59                 newTuner->randomMutate();
60         }
61         return newTuner;
62 }
63
64
65 void AutoTuner::tune() {
66         SearchTuner *bestTuner = NULL;
67         double bestScore = DBL_MAX;
68
69         SearchTuner *oldTuner = new SearchTuner();
70         double base_temperature = evaluateAll(oldTuner);
71         double oldScore = base_temperature;
72
73         for (uint i = 0; i < budget; i++) {
74                 SearchTuner *newTuner = mutateTuner(oldTuner, i);
75                 double newScore = evaluateAll(newTuner);
76                 newTuner->printUsed();
77                 model_print("Received score %f\n", newScore);
78                 double scoreDiff = newScore - oldScore; //smaller is better
79                 if (newScore < bestScore) {
80                         if (bestTuner != NULL)
81                                 delete bestTuner;
82                         bestScore = newScore;
83                         bestTuner = newTuner->copyUsed();
84                 }
85
86                 double acceptanceP;
87                 if (scoreDiff < 0) {
88                         acceptanceP = 1;
89                 } else {
90                         double currTemp = base_temperature * (((double)budget - i) / budget);
91                         acceptanceP = exp(-scoreDiff / currTemp);
92                 }
93                 double ran = ((double)random()) / RAND_MAX;
94                 if (ran <= acceptanceP) {
95                         delete oldTuner;
96                         oldScore = newScore;
97                         oldTuner = newTuner;
98                 } else {
99                         delete newTuner;
100                 }
101         }
102         model_print("Best tuner:\n");
103         bestTuner->print();
104         bestTuner->serialize();
105         model_print("Received score %f\n", bestScore);
106         if (bestTuner != NULL)
107                 delete bestTuner;
108         delete oldTuner;
109 }