Add sb03md wrapper for slycot fallback

jamestjsp · claude · jamestjsp · commit ce850198a6f1 · 2026-02-02T06:35:31.000-06:00
Co-Authored-By: Claude Opus 4.5 &lt;noreply@anthropic.com&gt;
diff --git a/control/slicot_compat.py b/control/slicot_compat.py
@@ -1122,7 +1122,7 @@ def mb03rd(n, A, X, pmax=1.0, tol=0.0, ldwork=None):
 # If using slycot (not slicot), overwrite with direct imports from slycot
 if _use_slycot and not _use_slicot:
     from slycot import (  # noqa: F811
-        sb03md, sb03od, sb04md, sb04qd, sg03ad,
+        sb03od, sb04md, sb04qd, sg03ad,
         sb02md, sb02mt, sg02ad, sb01bd,
         sb10ad, sb10hd, ab08nd,
         ab09ad, ab09md, ab09nd,
@@ -1133,3 +1133,10 @@ def mb03rd(n, A, X, pmax=1.0, tol=0.0, ldwork=None):
         SlycotResultWarning as SlicotResultWarning,
         SlycotArithmeticError as SlicotArithmeticError,
     )
+
+    from slycot import sb03md57
+
+    def sb03md(n, C, A, U, dico, job='X', fact='N', trana='N', ldwork=None):  # noqa: F811
+        """Wrapper for slycot's sb03md57."""
+        ret = sb03md57(A, U, C, dico, job, fact, trana, ldwork)
+        return ret[2:]
