commit 9071b2aa0fc0c1cd3e390c6259c6733fb2e4223e · dekker.one/on-restart-benchmarks

+105 -59
analyse_gecode.py
···

       13
       13
        
       import re

     

       14
       14
        
       import sys

     

       15
       15
        
       

     

       16
       16
       +
       from statistics import stdev

     

       16
       17
        
       

     

       17
       17
       -
       def compute_area(file):

     

       18
       18
       -
           area = -1

     

       19
       18
        
       

     

       19
       19
       +
       def compute_area(file):

     

       20
       20
        
           objectives = []

     

       21
       21
        
           times = [0]

     

       22
       22
       -
           timeout = -1

     

       23
       22
        
           objectives.append(0)

     

       24
       24
       -
           for line in contents:

     

       23
       23
       +
           for line in file:

     

       24
       24
       +
               # match = re.match(r'%\sinit_area\s=\s(\d+)', line)

     

       25
       25
        
               # match = re.match(r'%\sinit_area\s=\s(\d+)', line)

     

       26
       26
        
               # if match:

     

       27
       27
        
               #     objectives.append(int(match[1]))

     
···

       30
       30
        
               if match:

     

       31
       31
        
                   objectives.append(int(match.group(1)))

     

       32
       32
        
                   continue

     

       33
       33
       -
               match = re.match(r"%\stime elapsed:\s(\d\.\d+)\ss", line)

     

       33
       33
       +
               match = re.match(r"%\stime elapsed:\s(\d+\.\d+)\ss", line)

     

       34
       34
        
               if match:

     

       35
       35
        
                   times.append(float(match.group(1)))

     

       36
       36
        
                   continue

     
···

       43
       43
        
           assert len(objectives) + 1 == len(times)

     

       44
       44
        
           area = 0

     

       45
       45
        
           for i in range(len(objectives)):

     

       46
       46
       -
               area += ((times[i + 1] - times[i]) / 1000) * objectives[i]

     

       46
       46
       +
               area += (times[i + 1] - times[i]) * objectives[i]

     

       47
       47
        
           return int(area)

     

       48
       48
        
       

     

       49
       49
        
       

     

       50
       50
        
       folder = sys.argv[1]

     

       51
       51
       -
       statistics = {}

     

       52
       52
       -
       for root, dirs, files in os.walk(folder):

     

       53
       53
       -
           for name in files:

     

       54
       54
       -
               if name.endswith(".sol"):

     

       55
       55
       -
                   seed = 1

     

       56
       56
       -
                   match = re.search(r"\.(\d+)\.sol", name)

     

       57
       57
       -
                   if match:

     

       58
       58
       -
                       seed = int(match.group(1))

     

       59
       59
       -
                   with open(os.path.join(root, name)) as f:

     

       60
       60
       -
                       contents = f.readlines()

     

       61
       61
       -
                   # Area

     

       62
       62
       -
                   area = compute_area(contents)

     

       51
       51
       +
       statistics = dict()

     

       52
       52
       +
       instances = set()

     

       53
       53
       +
       for config in ["original", "restart", "replay"]:

     

       54
       54
       +
           for root, dirs, files in os.walk(folder + "/" + config):

     

       55
       55
       +
               for name in files:

     

       56
       56
       +
                   if name.endswith(".sol"):

     

       57
       57
       +
                       components = name[:-(4)].split(".")

     

       58
       58
       +
                       data = components[0]

     

       59
       59
       +
                       instances.add(data)

     

       60
       60
       +
                       seed = 1

     

       61
       61
       +
                       if len(components) > 1:

     

       62
       62
       +
                           assert len(components) == 2

     

       63
       63
       +
                           seed = components[1]

     

       63
       64
        
       

     

       64
       64
       -
                   objective = "UNSAT"

     

       65
       65
       -
                   for line in contents[::-1]:

     

       66
       66
       -
                       # Best objective

     

       67
       67
       -
                       match = re.match(r"objective\s=\s(\d+)", line)

     

       68
       68
       -
                       if match:

     

       69
       69
       -
                           objective = int(match.group(1))

     

       70
       70
       -
                           break

     

       65
       65
       +
                       if data not in statistics:

     

       66
       66
       +
                           statistics[data] = dict()

     

       67
       67
       +
                       if config not in statistics[data]:

     

       68
       68
       +
                           statistics[data][config] = []

     

       71
       69
        
       

     

       72
       72
       -
                   nodes = -1

     

       73
       73
       -
                   solvetime = -1

     

       74
       74
       -
                   restarts = -1

     

       75
       75
       -
                   for line in contents:

     

       76
       76
       -
                       # Evaluation time

     

       77
       77
       -
                       match = re.search(r"copies:\s+(\d+)", line)

     

       78
       78
       -
                       if match:

     

       79
       79
       -
                           nodes = int(match.group(1))

     

       80
       80
       -
                           continue

     

       81
       81
       -
                       # Solve time

     

       82
       82
       -
                       match = re.search(r"solveTime=(\d+(.\d+)?)", line)

     

       83
       83
       -
                       if match:

     

       84
       84
       -
                           solvetime = float(match.group(1))

     

       85
       85
       -
                           continue

     

       86
       86
       -
                       # Restarts

     

       87
       87
       -
                       match = re.search(r"restarts=(\d+)", line)

     

       88
       88
       -
                       if match:

     

       89
       89
       -
                           restarts = int(match.group(1))

     

       90
       90
       -
                           continue

     

       91
       91
       -
                   statistics[name[:-(4)].replace(".", ",")] = (

     

       92
       92
       -
                       area,

     

       93
       93
       -
                       objective,

     

       94
       94
       -
                       solvetime,

     

       95
       95
       -
                       restarts,

     

       96
       96
       -
                       nodes,

     

       97
       97
       -
                   )

     

       70
       70
       +
                       with open(os.path.join(root, name)) as f:

     

       71
       71
       +
                           contents = f.readlines()

     

       72
       72
       +
                       # Area

     

       73
       73
       +
                       area = compute_area(contents)

     

       98
       74
        
       

     

       99
       99
       -
       sorted_stats = sorted(statistics.items())

     

       100
       100
       -
       a = sorted_stats[0][0][: sorted_stats[0][0].find(",")]

     

       101
       101
       -
       for key, val in sorted_stats:

     

       102
       102
       -
           if key[: key.find(",")] != a:

     

       103
       103
       -
               print("\n\n")

     

       104
       104
       -
               a = key[: key.find(",")]

     

       105
       105
       -
           print("%s,%s" % (key, ",".join([v.__str__() for v in val])))

     

       75
       75
       +
                       objective = "UNSAT"

     

       76
       76
       +
                       for line in contents[::-1]:

     

       77
       77
       +
                           # Best objective

     

       78
       78
       +
                           match = re.match(r"objective\s=\s(\d+)", line)

     

       79
       79
       +
                           if match:

     

       80
       80
       +
                               objective = int(match.group(1))

     

       81
       81
       +
                               break

     

       82
       82
       +
       

     

       83
       83
       +
                       nodes = -1

     

       84
       84
       +
                       solvetime = -1

     

       85
       85
       +
                       restarts = -1

     

       86
       86
       +
                       for line in contents:

     

       87
       87
       +
                           # Evaluation time

     

       88
       88
       +
                           match = re.search(r"copies:\s+(\d+)", line)

     

       89
       89
       +
                           if match:

     

       90
       90
       +
                               nodes = int(match.group(1))

     

       91
       91
       +
                               continue

     

       92
       92
       +
                           # Solve time

     

       93
       93
       +
                           match = re.search(r"solveTime=(\d+(.\d+)?)", line)

     

       94
       94
       +
                           if match:

     

       95
       95
       +
                               solvetime = float(match.group(1))

     

       96
       96
       +
                               continue

     

       97
       97
       +
                           # Restarts

     

       98
       98
       +
                           match = re.search(r"restarts=(\d+)", line)

     

       99
       99
       +
                           if match:

     

       100
       100
       +
                               restarts = int(match.group(1))

     

       101
       101
       +
                               continue

     

       102
       102
       +
                       statistics[data][config].append(

     

       103
       103
       +
                           (

     

       104
       104
       +
                               area,

     

       105
       105
       +
                               objective,

     

       106
       106
       +
                               solvetime,

     

       107
       107
       +
                               restarts,

     

       108
       108
       +
                               nodes,

     

       109
       109
       +
                           )

     

       110
       110
       +
                       )

     

       106
       111
        
       

     

       107
       107
       -
       exit(1)

     

       112
       112
       +
       for data in instances:

     

       113
       113
       +
           for config in ["original", "restart", "replay"]:

     

       114
       114
       +
               stats = statistics[data][config]

     

       115
       115
       +
               cumulative = stats[0]

     

       116
       116
       +
               for i in range(1, len(stats)):

     

       117
       117
       +
                   cumulative = tuple(map(sum, zip(cumulative, stats[i])))

     

       118
       118
       +
               avg = tuple([x / len(stats) for x in cumulative])

     

       119
       119
       +
               dev = stdev([x[1] for x in stats]) if len(stats) > 1 else 0

     

       120
       120
       +
               # (avg area, avg objective, stdev objective)

     

       121
       121
       +
               statistics[data][config] = (avg[0], avg[1], dev)

     

       122
       122
       +
       

     

       123
       123
       +
       # Print header

     

       124
       124
       +
       print(

     

       125
       125
       +
           """

     

       126
       126
       +
       \\begin{tabular}{l|rr|rr|rr}

     

       127
       127
       +
       \\toprule

     

       128
       128
       +
       & \multicolumn{2}{c|}{Gecode} & \multicolumn{2}{c|}{Gecode Restart} & \multicolumn{2}{c}{Gecode Replay}\\\\

     

       129
       129
       +
       Instance & $\intobj$ & $\minobj$ & $\intobj$ & $\minobj$ & $\intobj$ & $\minobj$ \\\\

     

       130
       130
       +
       \midrule

     

       131
       131
       +
       """

     

       132
       132
       +
       )

     

       133
       133
       +
       

     

       134
       134
       +
       sorted_instances = sorted(instances)

     

       135
       135
       +
       for data in sorted_instances:

     

       136
       136
       +
           print(f"{data}", end="")

     

       137
       137
       +
           for config in ["original", "restart", "replay"]:

     

       138
       138
       +
               print(

     

       139
       139
       +
                   f" & {int(statistics[data][config][0] / 1000) }k & {int(statistics[data][config][1])}",

     

       140
       140
       +
                   end="",

     

       141
       141
       +
               )

     

       142
       142
       +
               if statistics[data][config][2] != 0:

     

       143
       143
       +
                   print("^{", end="")

     

       144
       144
       +
                   print(

     

       145
       145
       +
                       int(statistics[data][config][2] / statistics[data][config][1] * 100),

     

       146
       146
       +
                       # int(statistics[data][config][2]),

     

       147
       147
       +
                       end="",

     

       148
       148
       +
                   )

     

       149
       149
       +
                   print("}", end="")

     

       150
       150
       +
           print(" \\\\")

     

       151
       151
       +
       

     

       152
       152
       +
       # Print footer

     

       153
       153
       +
       print("\n\\bottomrule\n\end{tabular}")