Pyomo solver option

src/dispatch/pyomo_dispatch.py

@@ -10,6 +10,7 @@
 import time as time_mod
 from functools import partial
 import platform
+from itertools import compress
 
 import numpy as np
 import pyomo.environ as pyo
@@ -20,6 +21,7 @@
 # allows pyomo to solve on threaded processes
 import pyutilib.subprocess.GlobalData
 pyutilib.subprocess.GlobalData.DEFINE_SIGNAL_HANDLERS_DEFAULT = False
+from pyutilib.common._exceptions import ApplicationError
 
 from .Dispatcher import Dispatcher
 from .DispatchState import DispatchState, NumpyState
@@ -68,7 +70,10 @@ def get_input_specs(cls):
         \default{24}"""))
     specs.addSub(InputData.parameterInputFactory('debug_mode', contentType=InputTypes.BoolType,
         descr=r"""Enables additional printing in the pyomo dispatcher. Highly discouraged for production runs.
-              \default{False}."""))
+        \default{False}."""))
+    specs.addSub(InputData.parameterInputFactory('solver', contentType=InputTypes.StringType,
+        descr=r"""Indicates which solver should be used by pyomo. Options depend on individual installation.
+        \default{'glpk' for Windows, 'cbc' otherwise}."""))
     # TODO specific for pyomo dispatcher
     return specs
 
@@ -81,6 +86,7 @@ def __init__(self):
     self.name = 'PyomoDispatcher' # identifying name
     self.debug_mode = False       # whether to print additional information
     self._window_len = 24         # time window length to dispatch at a time # FIXME user input
+    self._solver = None           # overwrite option for solver
 
   def read_input(self, specs):
     """
@@ -98,6 +104,36 @@ def read_input(self, specs):
     if debug_node is not None:
       self.debug_mode = debug_node.value
 
+    solver_node = specs.findFirst('solver')
+    if solver_node is not None:
+      self._solver = solver_node.value
+
+    # check solver exists
+    if self._solver is None:
+      self._solver = SOLVER
+    found_solver = True
+    try:
+      if not pyo.SolverFactory(self._solver).available():
+        found_solver = False
+    except ApplicationError:
+      found_solver = False
+    if not found_solver:
+      all_options = pyo.SolverFactory._cls.keys() # TODO shorten to list of tested options?
+      solver_filter = []
+      for op in all_options:
+        if op.startswith('_'):
+          solver_filter.append(False)
+          continue
+        try:
+          solver_filter.append(pyo.SolverFactory(op).available())
+        except (ApplicationError, NameError, ImportError):
+          solver_filter.append(False)
+      available = list(compress(all_options, solver_filter))
+      msg = f'Requested solver "{self._solver}" was not found for pyomo dispatcher!'
+      msg += f' Options MAY include: {available}'
+      raise RuntimeError(msg)
+
+
   ### API
   def dispatch(self, case, components, sources, meta):
     """
@@ -209,7 +245,7 @@ def dispatch_window(self, time, time_offset,
       attempts += 1
       print(f'DEBUGG solve attempt {attempts} ...:')
       # solve
-      soln = pyo.SolverFactory(SOLVER).solve(m)
+      soln = pyo.SolverFactory(self._solver).solve(m)
       # check solve status
       if soln.solver.status == SolverStatus.ok and soln.solver.termination_condition == TerminationCondition.optimal:
         print('DEBUGG ... solve was successful!')

tests/integration_tests/mechanics/pyomo_options/heron_input.xml

@@ -31,6 +31,7 @@
       <pyomo>
         <rolling_window_length>8</rolling_window_length>
         <debug_mode>True</debug_mode>
+        <solver>ipopt</solver>
       </pyomo>
     </dispatcher>
   </Case>