Update

Author: odow

papers/Dowson_et_al_2026/Project.toml

@@ -1,6 +1,7 @@
 [deps]
+DataFrames = "a93c6f00-e57d-5684-b7b6-d8193f3e46c0"
 Gurobi = "2e9cd046-0924-5485-92f1-d5272153d98b"
 HiGHS = "87dc4568-4c63-4d18-b0c0-bb2238e4078b"
 JuMP = "4076af6c-e467-56ae-b986-b466b2749572"
+MathOptInterface = "b8f27783-ece8-5eb3-8dc8-9495eed66fee"
 MultiObjectiveAlgorithms = "0327d340-17cd-11ea-3e99-2fd5d98cecda"
-Polyhedra = "67491407-f73d-577b-9b50-8179a7c68029"

papers/Dowson_et_al_2026/main.jl

@@ -4,9 +4,12 @@
 #  obtain one at http://mozilla.org/MPL/2.0/.
 
 using JuMP
+import DataFrames
 import Gurobi
 import HiGHS
+import MathOptInterface
 import MultiObjectiveAlgorithms as MOA
+import Statistics
 
 function solve_tamby_mokp(;
     p::Int,
@@ -15,7 +18,7 @@ function solve_tamby_mokp(;
     optimizer::Any,
     algorithm::MOA.AbstractAlgorithm,
     silent::Bool = false,
-    time_limit::Float64 = 600.0,
+    time_limit::Float64 = 1200.0,
 )
     filename = joinpath(@__DIR__, "MOKP", "MOKP_p-$(p)_n-$(n)_$(instance).dat")
     lines = readlines(filename)
@@ -65,7 +68,7 @@ TESTSETS = Dict(
 
 results_filename = "results.log"
 results = Any[]
-for ((p, n), (algorithms, solvers)) in TESTSETS, instance in 1:1 #10
+for ((p, n), (algorithms, solvers)) in TESTSETS, instance in 1:10
     for alg in algorithms, solver in solvers
         algorithm, optimizer = ALGORITHMS[alg], OPTIMIZERS[solver]
         println("Running: p-$(p)_n-$(n)_$(instance) :: $alg :: $solver")
@@ -76,3 +79,21 @@ for ((p, n), (algorithms, solvers)) in TESTSETS, instance in 1:1 #10
         end
     end
 end
+
+l = eval.(Meta.parse.(readlines(results_filename)))
+df = DataFrames.DataFrame(l)
+df.percent_time = round.(Int, 100 .* df.solve_time_inner ./ df.solve_time)
+ret = DataFrames.combine(
+	DataFrames.groupby(df, [:p, :n, :alg, :solver]),
+	:result_count => Statistics.mean => :result_count,
+	:solve_time => Statistics.mean => :solve_time,
+	:subproblem_count => Statistics.mean => :subproblem_count,
+	:percent_time => Statistics.mean => :percent_time,
+)
+ret.solve_time .= max.(1, round.(Int, ret.solve_time))
+sort!(ret, [:p, :n, :alg, :solver])
+open("final_results.out", "w") do io
+    print(io, ret)
+end
+display(ret)
+