|
|
@ -8,6 +8,7 @@ import sys
|
|
|
|
from abc import ABC, abstractmethod
|
|
|
|
from abc import ABC, abstractmethod
|
|
|
|
from copy import deepcopy
|
|
|
|
from copy import deepcopy
|
|
|
|
from io import StringIO
|
|
|
|
from io import StringIO
|
|
|
|
|
|
|
|
from typing import Optional, List
|
|
|
|
|
|
|
|
|
|
|
|
import pyomo.core.kernel.objective
|
|
|
|
import pyomo.core.kernel.objective
|
|
|
|
import pyomo.environ as pe
|
|
|
|
import pyomo.environ as pe
|
|
|
@ -24,8 +25,8 @@ logger = logging.getLogger(__name__)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
# Global memory for multiprocessing
|
|
|
|
# Global memory for multiprocessing
|
|
|
|
SOLVER = [None]
|
|
|
|
SOLVER = [None] # type: List[Optional[LearningSolver]]
|
|
|
|
INSTANCES = [None]
|
|
|
|
INSTANCES = [None] # type: List[Optional[dict]]
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def _parallel_solve(instance_idx):
|
|
|
|
def _parallel_solve(instance_idx):
|
|
|
@ -298,14 +299,39 @@ class InternalSolver(ABC):
|
|
|
|
ws_value = float(matches[0])
|
|
|
|
ws_value = float(matches[0])
|
|
|
|
return ws_value
|
|
|
|
return ws_value
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def set_threads(self, threads):
|
|
|
|
|
|
|
|
key = self._get_threads_option_name()
|
|
|
|
|
|
|
|
self._pyomo_solver.options[key] = threads
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def set_time_limit(self, time_limit):
|
|
|
|
|
|
|
|
key = self._get_time_limit_option_name()
|
|
|
|
|
|
|
|
self._pyomo_solver.options[key] = time_limit
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def set_gap_tolerance(self, gap_tolerance):
|
|
|
|
|
|
|
|
key = self._get_gap_tolerance_option_name()
|
|
|
|
|
|
|
|
self._pyomo_solver.options[key] = gap_tolerance
|
|
|
|
|
|
|
|
|
|
|
|
@abstractmethod
|
|
|
|
@abstractmethod
|
|
|
|
def _get_warm_start_regexp(self):
|
|
|
|
def _get_warm_start_regexp(self):
|
|
|
|
pass
|
|
|
|
pass
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
@abstractmethod
|
|
|
|
|
|
|
|
def _get_threads_option_name(self):
|
|
|
|
|
|
|
|
pass
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
@abstractmethod
|
|
|
|
|
|
|
|
def _get_time_limit_option_name(self):
|
|
|
|
|
|
|
|
pass
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
@abstractmethod
|
|
|
|
|
|
|
|
def _get_gap_tolerance_option_name(self):
|
|
|
|
|
|
|
|
pass
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
class GurobiSolver(InternalSolver):
|
|
|
|
class GurobiSolver(InternalSolver):
|
|
|
|
def __init__(self,
|
|
|
|
def __init__(self,
|
|
|
|
use_lazy_callbacks=False):
|
|
|
|
use_lazy_callbacks=False,
|
|
|
|
|
|
|
|
options=None):
|
|
|
|
"""
|
|
|
|
"""
|
|
|
|
Creates a new GurobiSolver.
|
|
|
|
Creates a new GurobiSolver.
|
|
|
|
|
|
|
|
|
|
|
@ -314,21 +340,18 @@ class GurobiSolver(InternalSolver):
|
|
|
|
use_lazy_callbacks: bool
|
|
|
|
use_lazy_callbacks: bool
|
|
|
|
If true, lazy constraints will be enforced via lazy callbacks.
|
|
|
|
If true, lazy constraints will be enforced via lazy callbacks.
|
|
|
|
Otherwise, they will be enforced via a simple solve-check loop.
|
|
|
|
Otherwise, they will be enforced via a simple solve-check loop.
|
|
|
|
|
|
|
|
options: dict
|
|
|
|
|
|
|
|
Dictionary of options to pass to the Pyomo solver. For example,
|
|
|
|
|
|
|
|
{"Threads": 4} to set the number of threads.
|
|
|
|
"""
|
|
|
|
"""
|
|
|
|
super().__init__()
|
|
|
|
super().__init__()
|
|
|
|
self._use_lazy_callbacks = use_lazy_callbacks
|
|
|
|
self._use_lazy_callbacks = use_lazy_callbacks
|
|
|
|
self._pyomo_solver = pe.SolverFactory('gurobi_persistent')
|
|
|
|
self._pyomo_solver = pe.SolverFactory('gurobi_persistent')
|
|
|
|
self._pyomo_solver.options["Seed"] = randint(low=0, high=1000).rvs()
|
|
|
|
self._pyomo_solver.options["Seed"] = randint(low=0, high=1000).rvs()
|
|
|
|
|
|
|
|
if options is not None:
|
|
|
|
|
|
|
|
for (key, value) in options.items():
|
|
|
|
|
|
|
|
self._pyomo_solver.options[key] = value
|
|
|
|
|
|
|
|
|
|
|
|
def set_threads(self, threads):
|
|
|
|
|
|
|
|
self._pyomo_solver.options["Threads"] = threads
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def set_time_limit(self, time_limit):
|
|
|
|
|
|
|
|
self._pyomo_solver.options["TimeLimit"] = time_limit
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def set_gap_tolerance(self, gap_tolerance):
|
|
|
|
|
|
|
|
self._pyomo_solver.options["MIPGap"] = gap_tolerance
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def solve(self, tee=False):
|
|
|
|
def solve(self, tee=False):
|
|
|
|
if self._use_lazy_callbacks:
|
|
|
|
if self._use_lazy_callbacks:
|
|
|
|
return self._solve_with_callbacks(tee)
|
|
|
|
return self._solve_with_callbacks(tee)
|
|
|
@ -376,26 +399,35 @@ class GurobiSolver(InternalSolver):
|
|
|
|
def _get_warm_start_regexp(self):
|
|
|
|
def _get_warm_start_regexp(self):
|
|
|
|
return "MIP start with objective ([0-9.e+-]*)"
|
|
|
|
return "MIP start with objective ([0-9.e+-]*)"
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def _get_threads_option_name(self):
|
|
|
|
|
|
|
|
return "Threads"
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def _get_time_limit_option_name(self):
|
|
|
|
|
|
|
|
return "TimeLimit"
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def _get_gap_tolerance_option_name(self):
|
|
|
|
|
|
|
|
return "MIPGap"
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
class CPLEXSolver(InternalSolver):
|
|
|
|
class CPLEXSolver(InternalSolver):
|
|
|
|
def __init__(self,
|
|
|
|
def __init__(self, options=None):
|
|
|
|
presolve=1,
|
|
|
|
"""
|
|
|
|
mip_display=4,
|
|
|
|
Creates a new CPLEXSolver.
|
|
|
|
threads=None,
|
|
|
|
|
|
|
|
time_limit=None,
|
|
|
|
Parameters
|
|
|
|
gap_tolerance=None):
|
|
|
|
----------
|
|
|
|
|
|
|
|
options: dict
|
|
|
|
|
|
|
|
Dictionary of options to pass to the Pyomo solver. For example,
|
|
|
|
|
|
|
|
{"mip_display": 5} to increase the log verbosity.
|
|
|
|
|
|
|
|
"""
|
|
|
|
super().__init__()
|
|
|
|
super().__init__()
|
|
|
|
self._pyomo_solver = pe.SolverFactory('cplex_persistent')
|
|
|
|
self._pyomo_solver = pe.SolverFactory('cplex_persistent')
|
|
|
|
self._pyomo_solver.options["randomseed"] = randint(low=0, high=1000).rvs()
|
|
|
|
self._pyomo_solver.options["randomseed"] = randint(low=0, high=1000).rvs()
|
|
|
|
self._pyomo_solver.options["preprocessing_presolve"] = presolve
|
|
|
|
self._pyomo_solver.options["mip_display"] = 4
|
|
|
|
self._pyomo_solver.options["mip_display"] = mip_display
|
|
|
|
if options is not None:
|
|
|
|
if threads is not None:
|
|
|
|
for (key, value) in options.items():
|
|
|
|
self.set_threads(threads)
|
|
|
|
self._pyomo_solver.options[key] = value
|
|
|
|
if time_limit is not None:
|
|
|
|
|
|
|
|
self.set_time_limit(time_limit)
|
|
|
|
|
|
|
|
if gap_tolerance is not None:
|
|
|
|
|
|
|
|
self.set_gap_tolerance(gap_tolerance)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def set_threads(self, threads):
|
|
|
|
def set_threads(self, threads):
|
|
|
|
self._pyomo_solver.options["threads"] = threads
|
|
|
|
self._pyomo_solver.options["threads"] = threads
|
|
|
|
|
|
|
|
|
|
|
@ -420,6 +452,15 @@ class CPLEXSolver(InternalSolver):
|
|
|
|
def _get_warm_start_regexp(self):
|
|
|
|
def _get_warm_start_regexp(self):
|
|
|
|
return "MIP start .* with objective ([0-9.e+-]*)\\."
|
|
|
|
return "MIP start .* with objective ([0-9.e+-]*)\\."
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def _get_threads_option_name(self):
|
|
|
|
|
|
|
|
return "threads"
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def _get_time_limit_option_name(self):
|
|
|
|
|
|
|
|
return "timelimit"
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def _get_gap_tolerance_option_name(self):
|
|
|
|
|
|
|
|
return "mip_gap_tolerances_mipgap"
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
class LearningSolver:
|
|
|
|
class LearningSolver:
|
|
|
|
"""
|
|
|
|
"""
|
|
|
|