|
|
|
@ -9,6 +9,7 @@ import sys
|
|
|
|
|
from io import StringIO
|
|
|
|
|
from os.path import exists
|
|
|
|
|
from typing import Callable, List, Any
|
|
|
|
|
import traceback
|
|
|
|
|
|
|
|
|
|
from ..h5 import H5File
|
|
|
|
|
from ..io import _RedirectOutput, gzip, _to_h5_filename
|
|
|
|
@ -29,52 +30,57 @@ class BasicCollector:
|
|
|
|
|
verbose: bool = False,
|
|
|
|
|
) -> None:
|
|
|
|
|
def _collect(data_filename: str) -> None:
|
|
|
|
|
h5_filename = _to_h5_filename(data_filename)
|
|
|
|
|
mps_filename = h5_filename.replace(".h5", ".mps")
|
|
|
|
|
|
|
|
|
|
if exists(h5_filename):
|
|
|
|
|
# Try to read optimal solution
|
|
|
|
|
mip_var_values = None
|
|
|
|
|
try:
|
|
|
|
|
with H5File(h5_filename, "r") as h5:
|
|
|
|
|
mip_var_values = h5.get_array("mip_var_values")
|
|
|
|
|
except:
|
|
|
|
|
pass
|
|
|
|
|
|
|
|
|
|
if mip_var_values is None:
|
|
|
|
|
print(f"Removing empty/corrupted h5 file: {h5_filename}")
|
|
|
|
|
os.remove(h5_filename)
|
|
|
|
|
else:
|
|
|
|
|
return
|
|
|
|
|
|
|
|
|
|
with H5File(h5_filename, "w") as h5:
|
|
|
|
|
streams: List[Any] = [StringIO()]
|
|
|
|
|
if verbose:
|
|
|
|
|
streams += [sys.stdout]
|
|
|
|
|
with _RedirectOutput(streams):
|
|
|
|
|
# Load and extract static features
|
|
|
|
|
model = build_model(data_filename)
|
|
|
|
|
model.extract_after_load(h5)
|
|
|
|
|
|
|
|
|
|
if not self.skip_lp:
|
|
|
|
|
# Solve LP relaxation
|
|
|
|
|
relaxed = model.relax()
|
|
|
|
|
relaxed.optimize()
|
|
|
|
|
relaxed.extract_after_lp(h5)
|
|
|
|
|
|
|
|
|
|
# Solve MIP
|
|
|
|
|
model.optimize()
|
|
|
|
|
model.extract_after_mip(h5)
|
|
|
|
|
|
|
|
|
|
if self.write_mps:
|
|
|
|
|
# Add lazy constraints to model
|
|
|
|
|
model._lazy_enforce_collected()
|
|
|
|
|
|
|
|
|
|
# Save MPS file
|
|
|
|
|
model.write(mps_filename)
|
|
|
|
|
gzip(mps_filename)
|
|
|
|
|
|
|
|
|
|
h5.put_scalar("mip_log", streams[0].getvalue())
|
|
|
|
|
try:
|
|
|
|
|
h5_filename = _to_h5_filename(data_filename)
|
|
|
|
|
mps_filename = h5_filename.replace(".h5", ".mps")
|
|
|
|
|
|
|
|
|
|
if exists(h5_filename):
|
|
|
|
|
# Try to read optimal solution
|
|
|
|
|
mip_var_values = None
|
|
|
|
|
try:
|
|
|
|
|
with H5File(h5_filename, "r") as h5:
|
|
|
|
|
mip_var_values = h5.get_array("mip_var_values")
|
|
|
|
|
except:
|
|
|
|
|
pass
|
|
|
|
|
|
|
|
|
|
if mip_var_values is None:
|
|
|
|
|
print(f"Removing empty/corrupted h5 file: {h5_filename}")
|
|
|
|
|
os.remove(h5_filename)
|
|
|
|
|
else:
|
|
|
|
|
return
|
|
|
|
|
|
|
|
|
|
with H5File(h5_filename, "w") as h5:
|
|
|
|
|
streams: List[Any] = [StringIO()]
|
|
|
|
|
if verbose:
|
|
|
|
|
streams += [sys.stdout]
|
|
|
|
|
with _RedirectOutput(streams):
|
|
|
|
|
# Load and extract static features
|
|
|
|
|
model = build_model(data_filename)
|
|
|
|
|
model.extract_after_load(h5)
|
|
|
|
|
|
|
|
|
|
if not self.skip_lp:
|
|
|
|
|
# Solve LP relaxation
|
|
|
|
|
relaxed = model.relax()
|
|
|
|
|
relaxed.optimize()
|
|
|
|
|
relaxed.extract_after_lp(h5)
|
|
|
|
|
|
|
|
|
|
# Solve MIP
|
|
|
|
|
model.optimize()
|
|
|
|
|
model.extract_after_mip(h5)
|
|
|
|
|
|
|
|
|
|
if self.write_mps:
|
|
|
|
|
# Add lazy constraints to model
|
|
|
|
|
model._lazy_enforce_collected()
|
|
|
|
|
|
|
|
|
|
# Save MPS file
|
|
|
|
|
model.write(mps_filename)
|
|
|
|
|
gzip(mps_filename)
|
|
|
|
|
|
|
|
|
|
h5.put_scalar("mip_log", streams[0].getvalue())
|
|
|
|
|
except:
|
|
|
|
|
print(f"Error processing: data_filename")
|
|
|
|
|
traceback.print_exc()
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
if n_jobs > 1:
|
|
|
|
|
p_umap(
|
|
|
|
|