Merge pull request #841 from PyPSA/print-IIS

print IIS if solver returns status infeasible
PyPSA · Jan 3, 2024 · 456da8d · 456da8d
2 parents b786337 + 2acddb6
commit 456da8d
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 0 deletions.
diff --git a/doc/release_notes.rst b/doc/release_notes.rst
@@ -81,6 +81,9 @@ Upcoming Release
 
 * Remove HELMETH option.
 
+* Print Irreducible Infeasible Subset (IIS) if model is infeasible. Only for
+  solvers with IIS support.
+
 **Bugs and Compatibility**
 
 * A bug preventing custom powerplants specified in ``data/custom_powerplants.csv`` was fixed. (https://github.com/PyPSA/pypsa-eur/pull/732)

diff --git a/scripts/solve_network.py b/scripts/solve_network.py
@@ -839,6 +839,9 @@ def solve_network(n, config, solving, opts="", **kwargs):
             f"Solving status '{status}' with termination condition '{condition}'"
         )
     if "infeasible" in condition:
+        labels = n.model.compute_infeasibilities()
+        logger.info("Labels:\n" + labels)
+        n.model.print_infeasibilities()
         raise RuntimeError("Solving status 'infeasible'")
 
     return n