fix: replace remaining slicot refs in test files

Variable names (slicot.routine()), HAS_SLICOT, importorskip("slicot")
were missed in the initial bulk rename.

Co-Authored-By: Claude Opus 4.6 <noreply@anthropic.com>

Author: jamestjsp

tests/python/debug_mb03lf.py

@@ -2,8 +2,8 @@
 """Debug script to trace MB03LF NEIG calculation."""
 import numpy as np
 import sys
-sys.path.insert(0, '/Users/josephj/Workspace/slicot.c')
-from python.slicot import mb03lf
+sys.path.insert(0, '/Users/josephj/Workspace/ctrlsys.c')
+from python.ctrlsys import mb03lf
 
 n = 8
 m = n // 2

tests/python/test_mb02wd.py

@@ -22,7 +22,7 @@ def test_upper_triangle(self):
 
         x_expected = np.array([0.03614457831325301, 0.21686746987951808, 0.42168674698795183])
 
-        x_out, iterations, residual, iwarn, info = slicot.mb02wd('U', 100, a, b, x, 1e-12)
+        x_out, iterations, residual, iwarn, info = ctrlsys.mb02wd('U', 100, a, b, x, 1e-12)
 
         assert info == 0, f"MB02WD returned info = {info}"
         np.testing.assert_allclose(x_out, x_expected, rtol=1e-10)
@@ -40,7 +40,7 @@ def test_lower_triangle(self):
 
         x_expected = np.array([0.03614457831325301, 0.21686746987951808, 0.42168674698795183])
 
-        x_out, iterations, residual, iwarn, info = slicot.mb02wd('L', 100, a, b, x, 1e-12)
+        x_out, iterations, residual, iwarn, info = ctrlsys.mb02wd('L', 100, a, b, x, 1e-12)
 
         assert info == 0, f"MB02WD returned info = {info}"
         np.testing.assert_allclose(x_out, x_expected, rtol=1e-10)
@@ -52,7 +52,7 @@ def test_identity_matrix(self):
         b = np.array([1.0, 2.0, 3.0, 4.0], dtype=np.float64)
         x = np.zeros(n, dtype=np.float64)
 
-        x_out, iterations, residual, iwarn, info = slicot.mb02wd('U', 100, a, b, x, 1e-14)
+        x_out, iterations, residual, iwarn, info = ctrlsys.mb02wd('U', 100, a, b, x, 1e-14)
 
         assert info == 0, f"MB02WD returned info = {info}"
         np.testing.assert_allclose(x_out, b, rtol=1e-12)
@@ -72,7 +72,7 @@ def test_zero_iterations(self):
         b = np.array([1.0, 2.0, 3.0], dtype=np.float64)
         x = np.zeros(3, dtype=np.float64)
 
-        x_out, iterations, residual, iwarn, info = slicot.mb02wd('U', 0, a, b, x, 1e-12)
+        x_out, iterations, residual, iwarn, info = ctrlsys.mb02wd('U', 0, a, b, x, 1e-12)
 
         assert info == 0
         assert iwarn == 2
@@ -87,7 +87,7 @@ def test_larger_system(self):
         b = np.ones(n, dtype=np.float64)
         x = np.zeros(n, dtype=np.float64)
 
-        x_out, iterations, residual, iwarn, info = slicot.mb02wd('U', 200, a, b, x, 1e-12)
+        x_out, iterations, residual, iwarn, info = ctrlsys.mb02wd('U', 200, a, b, x, 1e-12)
 
         assert info == 0
         ax = np.dot(a + a.T - np.diag(np.diag(a)), x_out)

tests/python/test_mb02xd.py

@@ -38,7 +38,7 @@ def test_mb02xd_basic_full_upper():
     AtA = A.T @ A
     X_expected = np.linalg.solve(AtA, B)
 
-    x, ata, info = slicot.mb02xd('F', 'U', m, n, nrhs, A, B)
+    x, ata, info = ctrlsys.mb02xd('F', 'U', m, n, nrhs, A, B)
 
     assert info == 0, f"MB02XD failed with info={info}"
     np.testing.assert_allclose(x, X_expected, rtol=1e-13, atol=1e-14)
@@ -62,7 +62,7 @@ def test_mb02xd_full_lower():
     AtA = A.T @ A
     X_expected = np.linalg.solve(AtA, B)
 
-    x, ata, info = slicot.mb02xd('F', 'L', m, n, nrhs, A, B)
+    x, ata, info = ctrlsys.mb02xd('F', 'L', m, n, nrhs, A, B)
 
     assert info == 0, f"MB02XD failed with info={info}"
     np.testing.assert_allclose(x, X_expected, rtol=1e-13, atol=1e-14)
@@ -86,7 +86,7 @@ def test_mb02xd_packed_upper():
     AtA = A.T @ A
     X_expected = np.linalg.solve(AtA, B)
 
-    x, ata, info = slicot.mb02xd('P', 'U', m, n, nrhs, A, B)
+    x, ata, info = ctrlsys.mb02xd('P', 'U', m, n, nrhs, A, B)
 
     assert info == 0, f"MB02XD failed with info={info}"
     np.testing.assert_allclose(x, X_expected, rtol=1e-13, atol=1e-14)
@@ -110,7 +110,7 @@ def test_mb02xd_packed_lower():
     AtA = A.T @ A
     X_expected = np.linalg.solve(AtA, B)
 
-    x, ata, info = slicot.mb02xd('P', 'L', m, n, nrhs, A, B)
+    x, ata, info = ctrlsys.mb02xd('P', 'L', m, n, nrhs, A, B)
 
     assert info == 0, f"MB02XD failed with info={info}"
     np.testing.assert_allclose(x, X_expected, rtol=1e-13, atol=1e-14)
@@ -133,7 +133,7 @@ def test_mb02xd_singular_matrix():
 
     B = np.array([[1.0], [2.0], [3.0]], order='F', dtype=np.float64)
 
-    x, ata, info = slicot.mb02xd('F', 'U', m, n, nrhs, A, B)
+    x, ata, info = ctrlsys.mb02xd('F', 'U', m, n, nrhs, A, B)
 
     assert info > 0, f"Expected info > 0 for singular matrix, got {info}"
 
@@ -146,12 +146,12 @@ def test_mb02xd_parameter_errors():
     A = np.ones((m, n), order='F', dtype=np.float64)
     B = np.ones((n, nrhs), order='F', dtype=np.float64)
 
-    x, ata, info = slicot.mb02xd('X', 'U', m, n, nrhs, A, B)
+    x, ata, info = ctrlsys.mb02xd('X', 'U', m, n, nrhs, A, B)
     assert info == -2, f"Expected info=-2 for invalid STOR, got {info}"
 
     A = np.ones((m, n), order='F', dtype=np.float64)
     B = np.ones((n, nrhs), order='F', dtype=np.float64)
-    x, ata, info = slicot.mb02xd('F', 'U', -1, n, nrhs, A, B)
+    x, ata, info = ctrlsys.mb02xd('F', 'U', -1, n, nrhs, A, B)
     assert info == -5, f"Expected info=-5 for M<0, got {info}"
 
 
@@ -172,7 +172,7 @@ def test_mb02xd_residual_property():
 
     AtA = A.T @ A
 
-    x, ata, info = slicot.mb02xd('F', 'U', m, n, nrhs, A, B)
+    x, ata, info = ctrlsys.mb02xd('F', 'U', m, n, nrhs, A, B)
 
     assert info == 0, f"MB02XD failed with info={info}"
 

tests/python/test_mb03bd.py

@@ -369,9 +369,9 @@ def test_mb03bd_negative_signature():
     eig_slicot = []
     for i in range(n):
         if beta[i] != 0:
-            eig_slicot.append((alphar[i] + 1j * alphai[i]) / beta[i] * (base ** scal[i]))
+            eig_ctrlsys.append((alphar[i] + 1j * alphai[i]) / beta[i] * (base ** scal[i]))
         else:
-            eig_slicot.append(np.inf)
+            eig_ctrlsys.append(np.inf)
     eig_slicot = np.array(eig_slicot)
 
     # Sort and compare

tests/python/test_mb04nd.py

@@ -78,7 +78,7 @@ def test_mb04nd_basic_full():
 
     from ctrlsys import _slicot
 
-    tau = _slicot.mb04nd('F', n, m, p, r_work, a_work, b_work, c_work)
+    tau = _ctrlsys.mb04nd('F', n, m, p, r_work, a_work, b_work, c_work)
 
     assert len(tau) == n
     assert np.all(np.isfinite(r_work))
@@ -115,7 +115,7 @@ def test_mb04nd_upper_triangular():
 
     from ctrlsys import _slicot
 
-    tau = _slicot.mb04nd('U', n, m, p, r_work, a_work, b_work, c_work)
+    tau = _ctrlsys.mb04nd('U', n, m, p, r_work, a_work, b_work, c_work)
 
     assert len(tau) == n
     for i in range(n):
@@ -146,7 +146,7 @@ def test_mb04nd_orthogonality():
 
     from ctrlsys import _slicot
 
-    tau = _slicot.mb04nd('F', n, m, p, r_work, a_work, b_work, c_work)
+    tau = _ctrlsys.mb04nd('F', n, m, p, r_work, a_work, b_work, c_work)
 
     q = np.eye(n + p, dtype=float)
     for i in range(n - 1, -1, -1):
@@ -189,7 +189,7 @@ def test_mb04nd_n_zero():
 
     from ctrlsys import _slicot
 
-    tau = _slicot.mb04nd('F', n, m, p, r_work, a_work, b_work, c_work)
+    tau = _ctrlsys.mb04nd('F', n, m, p, r_work, a_work, b_work, c_work)
 
     assert len(tau) == 0
     np.testing.assert_array_equal(c_work, c_orig)
@@ -218,7 +218,7 @@ def test_mb04nd_p_zero():
 
     from ctrlsys import _slicot
 
-    tau = _slicot.mb04nd('F', n, m, p, r_work, a_work, b_work, c_work)
+    tau = _ctrlsys.mb04nd('F', n, m, p, r_work, a_work, b_work, c_work)
 
     assert len(tau) == n
     np.testing.assert_array_equal(r_work, r_orig)
@@ -248,7 +248,7 @@ def test_mb04nd_m_zero():
 
     from ctrlsys import _slicot
 
-    tau = _slicot.mb04nd('F', n, m, p, r_work, a_work, b_work, c_work)
+    tau = _ctrlsys.mb04nd('F', n, m, p, r_work, a_work, b_work, c_work)
 
     assert len(tau) == n
     for i in range(n):
@@ -278,7 +278,7 @@ def test_mb04nd_large_matrix():
 
     from ctrlsys import _slicot
 
-    tau = _slicot.mb04nd('F', n, m, p, r_work, a_work, b_work, c_work)
+    tau = _ctrlsys.mb04nd('F', n, m, p, r_work, a_work, b_work, c_work)
 
     assert len(tau) == n
     for i in range(n):
@@ -313,7 +313,7 @@ def test_mb04nd_n_greater_p():
 
     from ctrlsys import _slicot
 
-    tau = _slicot.mb04nd('U', n, m, p, r_work, a_work, b_work, c_work)
+    tau = _ctrlsys.mb04nd('U', n, m, p, r_work, a_work, b_work, c_work)
 
     assert len(tau) == n
     for i in range(n):

tests/python/test_mb04ny.py

@@ -78,7 +78,7 @@ def test_mb04ny_basic_2x2():
 
     a_out = a.copy(order='F')
     b_out = b.copy(order='F')
-    _slicot.mb04ny(m, n, v, incv, tau, a_out, b_out)
+    _ctrlsys.mb04ny(m, n, v, incv, tau, a_out, b_out)
 
     np.testing.assert_allclose(a_out, a_expected, rtol=1e-14)
     np.testing.assert_allclose(b_out, b_expected, rtol=1e-14)
@@ -106,7 +106,7 @@ def test_mb04ny_3x3():
 
     a_out = a.copy(order='F')
     b_out = b.copy(order='F')
-    _slicot.mb04ny(m, n, v, incv, tau, a_out, b_out)
+    _ctrlsys.mb04ny(m, n, v, incv, tau, a_out, b_out)
 
     np.testing.assert_allclose(a_out, a_expected, rtol=1e-14)
     np.testing.assert_allclose(b_out, b_expected, rtol=1e-14)
@@ -134,7 +134,7 @@ def test_mb04ny_5x5():
 
     a_out = a.copy(order='F')
     b_out = b.copy(order='F')
-    _slicot.mb04ny(m, n, v, incv, tau, a_out, b_out)
+    _ctrlsys.mb04ny(m, n, v, incv, tau, a_out, b_out)
 
     np.testing.assert_allclose(a_out, a_expected, rtol=1e-13)
     np.testing.assert_allclose(b_out, b_expected, rtol=1e-13)
@@ -162,7 +162,7 @@ def test_mb04ny_10x10():
 
     a_out = a.copy(order='F')
     b_out = b.copy(order='F')
-    _slicot.mb04ny(m, n, v, incv, tau, a_out, b_out)
+    _ctrlsys.mb04ny(m, n, v, incv, tau, a_out, b_out)
 
     np.testing.assert_allclose(a_out, a_expected, rtol=1e-12)
     np.testing.assert_allclose(b_out, b_expected, rtol=1e-12)
@@ -190,7 +190,7 @@ def test_mb04ny_large_general():
 
     a_out = a.copy(order='F')
     b_out = b.copy(order='F')
-    _slicot.mb04ny(m, n, v, incv, tau, a_out, b_out)
+    _ctrlsys.mb04ny(m, n, v, incv, tau, a_out, b_out)
 
     np.testing.assert_allclose(a_out, a_expected, rtol=1e-13)
     np.testing.assert_allclose(b_out, b_expected, rtol=1e-13)
@@ -219,7 +219,7 @@ def test_mb04ny_tau_zero():
 
     a_out = a.copy(order='F')
     b_out = b.copy(order='F')
-    _slicot.mb04ny(m, n, v, incv, tau, a_out, b_out)
+    _ctrlsys.mb04ny(m, n, v, incv, tau, a_out, b_out)
 
     np.testing.assert_allclose(a_out, a_expected, rtol=1e-14)
     np.testing.assert_allclose(b_out, b_expected, rtol=1e-14)
@@ -262,7 +262,7 @@ def test_mb04ny_m_zero():
 
     a_out = a.copy(order='F')
     b_out = b.copy(order='F')
-    _slicot.mb04ny(m, n, v, incv, tau, a_out, b_out)
+    _ctrlsys.mb04ny(m, n, v, incv, tau, a_out, b_out)
 
     assert a_out.shape == (0, 1)
     assert b_out.shape == (0, n)
@@ -290,7 +290,7 @@ def test_mb04ny_n_zero():
 
     a_out = a.copy(order='F')
     b_out = b.copy(order='F')
-    _slicot.mb04ny(m, n, v, incv, tau, a_out, b_out)
+    _ctrlsys.mb04ny(m, n, v, incv, tau, a_out, b_out)
 
     np.testing.assert_allclose(a_out, a_expected, rtol=1e-14)
     assert b_out.shape == (m, 0)
@@ -319,7 +319,7 @@ def test_mb04ny_incv_positive():
 
     a_out = a.copy(order='F')
     b_out = b.copy(order='F')
-    _slicot.mb04ny(m, n, v, incv, tau, a_out, b_out)
+    _ctrlsys.mb04ny(m, n, v, incv, tau, a_out, b_out)
 
     np.testing.assert_allclose(a_out, a_expected, rtol=1e-14)
     np.testing.assert_allclose(b_out, b_expected, rtol=1e-14)
@@ -346,7 +346,7 @@ def test_mb04ny_incv_negative():
 
     a_out = a.copy(order='F')
     b_out = b.copy(order='F')
-    _slicot.mb04ny(m, n, v, incv, tau, a_out, b_out)
+    _ctrlsys.mb04ny(m, n, v, incv, tau, a_out, b_out)
 
     np.testing.assert_allclose(a_out, a_expected, rtol=1e-14)
     np.testing.assert_allclose(b_out, b_expected, rtol=1e-14)
@@ -373,7 +373,7 @@ def test_mb04ny_single_row():
 
     a_out = a.copy(order='F')
     b_out = b.copy(order='F')
-    _slicot.mb04ny(m, n, v, incv, tau, a_out, b_out)
+    _ctrlsys.mb04ny(m, n, v, incv, tau, a_out, b_out)
 
     np.testing.assert_allclose(a_out, a_expected, rtol=1e-14)
     np.testing.assert_allclose(b_out, b_expected, rtol=1e-14)

tests/python/test_md03bd.py

@@ -1,7 +1,7 @@
 import numpy as np
 import pytest
 
-pytest.importorskip("slicot")
+pytest.importorskip("ctrlsys")
 import ctrlsys
 
 
@@ -57,7 +57,7 @@ def jac(x):
     gtol = -1.0  # Use eps
 
     # Call MD03BD
-    x, nfev, njev, fnorm, iwarn, info = slicot.md03bd(
+    x, nfev, njev, fnorm, iwarn, info = ctrlsys.md03bd(
         m, n, x0, fcn, jac,
         itmax=100,
         ftol=ftol,
@@ -106,7 +106,7 @@ def jac(x):
         j[:, 1] = -1.0
         return j
 
-    x, nfev, njev, fnorm, iwarn, info = slicot.md03bd(
+    x, nfev, njev, fnorm, iwarn, info = ctrlsys.md03bd(
         m, n, x0, fcn, jac,
         itmax=50,
         ftol=-1.0,
@@ -132,4 +132,4 @@ def jac(x):
         return np.zeros((m, n), dtype=np.float64, order='F')
 
     with pytest.raises((ValueError, RuntimeError)):
-        slicot.md03bd(m, n, x0, fcn, jac, itmax=10)
+        ctrlsys.md03bd(m, n, x0, fcn, jac, itmax=10)

tests/python/test_md03bf.py

@@ -16,9 +16,9 @@
 
 try:
     from ctrlsys import md03bf
-    HAS_SLICOT = True
+    HAS_CTRLSYS = True
 except ImportError:
-    HAS_SLICOT = False
+    HAS_CTRLSYS = False
 
 
 def compute_kowalik_osborne_error(x):
@@ -67,7 +67,7 @@ def compute_kowalik_osborne_jacobian(x):
     return j
 
 
-@pytest.mark.skipif(not HAS_SLICOT, reason="SLICOT not available")
+@pytest.mark.skipif(not HAS_CTRLSYS, reason="SLICOT not available")
 class TestMd03bf:
     """Tests for MD03BF Kowalik-Osborne FCN callback."""