Optional packages coinbrew, mibs, pao, pyomo

[mkoeppe]

• https://github.com/or-fusion/pao (as of 2023-10-06, no development since 2021-07)

@xuluze

📝 Checklist

• The title is concise, informative, and self-explanatory.
• The description explains in detail what this PR is about.
• I have linked a relevant issue or discussion.
• I have created tests covering the changes.
• I have updated the documentation accordingly.

⌛ Dependencies

[xuluze]

pao depends on module pyutilib.

sage: from pao.pyomo import *
...
ModuleNotFoundError: No module named 'pyutilib'

I can run the simple example with some deprecation warning:

sage: import pyomo.environ as pe
xxx/mambaforge/envs/sage-mibs/lib/python3.11/site-packages/pyutilib/misc/import_file.py:11: DeprecationWarning: the imp module is deprecated in favour of importlib and slated for removal in Python 3.12; see the module's documentation for alternative uses
  import imp
sage: from pao.pyomo import *
WARNING: DEPRECATED: add_docstring_list is deprecated.  Please use the
@document_kwargs_from_configdict() decorator.  (deprecated in 6.6.0) (called
from xxx/mambaforge/envs/sage-mibs/lib/python3.11/site-
packages/pao/common/solver.py:189)
sage: M = pe.ConcreteModel()
sage: M.x = pe.Var(bounds=(0,None))
....: M.y = pe.Var(bounds=(0,None))
sage: M.o = pe.Objective(expr=M.x - 4*M.y)
WARNING: Dynamically registering the following numeric type:
        sage.rings.integer.Integer
    Dynamic registration is supported for convenience, but there are known
    limitations to this approach.  We recommend explicitly registering numeric
    types using RegisterNumericType() or RegisterIntegerType().
sage: M.L = SubModel(fixed=M.x)
sage: M.L.o = pe.Objective(expr=M.y)
sage: M.L.c1 = pe.Constraint(expr=   -M.x -   M.y <= -3)
sage: M.L.c2 = pe.Constraint(expr= -2*M.x +   M.y <=  0)
sage: M.L.c3 = pe.Constraint(expr=  2*M.x +   M.y <= 12)
sage: M.L.c4 = pe.Constraint(expr=  3*M.x - 2*M.y <=  4)
sage: with Solver('pao.pyomo.FA') as solver: results = solver.solve(M)
sage: print(M.x.value)
4.0
sage: print(M.y.value)
4.0

If I replace it with Solver('pao.pyomo.MIBS'), it just output 0.

[mkoeppe]

I think this may be the same as:

• PAO not working with MibS or-fusion/pao#98

[github-actions]

Documentation preview for this PR (built with commit 4f76fb8; changes) is ready! 🎉