SLICOT `C11 translation with Python wrapper` progress update #1192

jamestjsp · 2026-01-17T16:37:07Z

jamestjsp
Jan 17, 2026

Now ~150 routines remains to translate (Thanks to Claude Code), I will keep the repo private until 100% complete.
https://github.com/jamestjsp/slicot.c

SLICOT C11 Translation Status

Summary

Metric	Count
Translated routines	455
Test files	417
Pending translation	168

Well-Tested Translated Routines

State-Space Analysis & Transformation (AB family)

Routine	Description
AB01MD	Staircase form for state-space realization
AB01ND	Controllability staircase form
AB01OD	Orthogonal transformation to controllable form
AB04MD	Bilinear transformation (continuous ↔ discrete)
AB05MD	Series connection of systems
AB05ND	Parallel connection of systems
AB07MD	System inversion
AB08MD	System zeros via generalized eigenvalue problem
AB08ND	Right coprime factorization with pole assignment
AB09AD	Balanced truncation model reduction
AB09BD	Square-root balanced truncation
AB09CD	Frequency-weighted balanced truncation
AB09HD	Hankel-norm approximation
AB09JD	Singular perturbation approximation
AB09MD	Model reduction via modal truncation
AB13AD	Hankel singular values
AB13BD	H-infinity norm computation
AB13CD	Balanced transformation computation
AB13DD	L-infinity norm of transfer function
AB13ED	Controllability and observability Gramians
AB13FD	L2 gain computation
AB13MD	Stability margin computation

Controller Design (SB family)

Routine	Description
SB01BD	Pole placement via state feedback
SB02MD	Algebraic Riccati equation (continuous)
SB02ND	Algebraic Riccati equation (discrete)
SB02OD	Optimal Riccati solution
SB02RD	Riccati equation via Schur method
SB02SD	LQR/LQG optimal control
SB03MD	Lyapunov equation solver (continuous)
SB03OD	Lyapunov equation solver (discrete)
SB04MD	Sylvester equation solver
SB04ND	Discrete Sylvester equation
SB04QD	Generalized Sylvester equation
SB08DD	Spectral factorization
SB10AD	H-infinity controller (continuous)
SB10DD	H-infinity controller (discrete)
SB10HD	H2 optimal controller
SB10FD	Loop-shaping design
SB10PD	State-space H-infinity synthesis

Descriptor Systems (SG/TG family)

Routine	Description
SG02AD	Generalized Riccati equation
SG03AD	Generalized Lyapunov equation (continuous)
SG03BD	Generalized Lyapunov equation (discrete)
TG01AD	Descriptor system balancing
TG01BD	Descriptor to state-space conversion
TG01FD	Finite/infinite pole separation

System Identification (IB family)

Routine	Description
IB01AD	Subspace identification (MOESP/N4SID)
IB01BD	System matrices from subspace identification
IB01CD	Covariance sequence estimation
IB01PD	Kalman filter gain computation
IB01QD	Initial state estimation
IB01RD	Regularized subspace identification
IB03AD	Box-Jenkins model identification
IB03BD	ARMAX model identification

Transfer Functions (TB/TC/TD/TF family)

Routine	Description
TB01ID	State-space to transfer function
TB01KD	Minimal realization
TB01LD	Observability/controllability reduction
TB01PD	Pole-zero cancellation
TB01WD	Frequency response computation
TB03AD	Polynomial to state-space conversion
TB04AD	State-space to polynomial matrix
TB05AD	Transfer function evaluation at complex frequency
TC01OD	Polynomial matrix operations
TC04AD	Polynomial matrix to state-space
TD04AD	Right matrix fraction description
TF01MD	Transfer function matrix evaluation
TF01RD	Stable/unstable additive decomposition

Matrix Computations (MB family)

Routine	Description
MB02CD	QR factorization with column pivoting
MB02DD	Cholesky factorization
MB02ED	RQ factorization
MB02GD	Triangular matrix inversion
MB02JD	Band matrix operations
MB02MD	Least squares solution
MB02SD	Singular value decomposition
MB02TD	Block diagonal system solution
MB02VD	Linear system solution
MB03AD	QR algorithm for Hessenberg matrices
MB03BD	Real Schur form computation
MB03CD	Eigenvalue ordering in Schur form
MB03FD	Eigenvalue computation
MB03JD	Symplectic QR algorithm
MB03LD	Periodic Schur decomposition
MB03RD	Block diagonal form transformation
MB03TD	Hamiltonian Schur decomposition
MB03VD	Periodic eigenvalue computation
MB03XD	Skew-Hamiltonian/Hamiltonian pencil
MB04AD	Hessenberg-triangular reduction
MB04BD	Hessenberg form computation
MB04CD	Symplectic URV decomposition
MB04DD	Orthogonal balancing
MB04ED	Product Hessenberg form
MB04FD	Symplectic butterfly form
MB05MD	Matrix exponential (Padé approximation)
MB05ND	Matrix exponential and integral

Filtering (FB family)

Routine	Description
FB01QD	Square-root covariance filter
FB01RD	Information filter
FB01SD	Square-root information filter
FB01TD	Time-varying Kalman filter
FB01VD	One-step Kalman predictor

Polynomial Operations (MC family)

Routine	Description
MC01PD	Polynomial evaluation
MC01TD	Polynomial stability test

Notable Routines Pending Translation

High-Priority Controller Design

Routine	Description	Complexity
SB01DD	Partial pole placement	Medium
SB01MD	Deadbeat controller design	Medium
SB02PD	Riccati equation (Newton method)	High
SB03TD	Stein equation solver	Medium
SB03UD	Discrete Lyapunov (direct method)	Medium
SB16AD	Controller reduction (balanced truncation)	High
SB16BD	Controller reduction (frequency-weighted)	High
SB16CD	Coprime factorization controller reduction	High

Spectral Factorization & Coprime Methods

Routine	Description	Complexity
SB08CD	Left spectral factorization	High
SB08ED	J-spectral factorization	High
SB08GD	Normalized coprime factorization	High
SB08HD	Robust stabilization	High
SB08MD	Right coprime factorization	Medium
SB08ND	Left coprime factorization	Medium

Model Reduction

Routine	Description	Complexity
AB09FD	Frequency-weighted balanced reduction	High
AB09GD	Coprime factor model reduction	High
AB09ID	Inner-outer factorization model reduction	High
AB13HD	Impulse response energy	Medium
AB13ID	Peak gain computation	High

Polynomial & Matrix Fraction Operations

Routine	Description	Complexity
MC01MD	Polynomial GCD computation	Medium
MC01ND	Polynomial LCM computation	Medium
MC01OD	Polynomial division	Medium
MC01QD	Polynomial multiplication	Low
MC01RD	Polynomial root finding	High
MC01VD	Polynomial interpolation	Medium
MC03MD	Polynomial matrix determinant	High
MC03ND	Polynomial matrix inverse	High

Descriptor System Transformations

Routine	Description	Complexity
TG01CD	Descriptor coordinate transformation	Medium
TG01DD	Controllability reduction (descriptor)	Medium
TG01ED	Observability reduction (descriptor)	Medium
TG01GD	Irreducibility test	Medium
TG01HD	Descriptor to Weierstrass form	High
TG01ID	Impulse-free descriptor form	Medium
TG01JD	Proper/improper decomposition	High
TG01KD	Infinite eigenvalue removal	Medium
TG01LD	Descriptor system poles/zeros	High
TG01MD	Minimal realization (descriptor)	High
TG01ND	System normalization	Medium
TG01PD	Left null-space computation	Medium
TG01QD	Right null-space computation	Medium
TG01WD	Descriptor balancing	Medium

Transfer Function Conversion

Routine	Description	Complexity
TB04BD	State-space to transfer function (symbolic)	High
TB04CD	Transfer function matrix conversion	Medium
TD03AD	Left MFD conversion	Medium
TD05AD	Polynomial matrix factorization	High
TF01ND	Stable part extraction	Medium
TF01OD	Anti-stable part extraction	Medium
TF01PD	All-pass factorization	High
TF01QD	Inner-outer factorization	High

Complex/Periodic Matrix Operations

Routine	Description	Complexity
MB03WA	Complex periodic Schur form	High
MB03ZA	Complex Hamiltonian Schur form	High
MB03ZD	Complex skew-Hamiltonian decomposition	High
MB04WD	Symplectic QR factorization (blocked)	High
MB04XY	Periodic QR iteration	High
MB05OD	Matrix exponential (scaling & squaring)	Medium
MB05OY	Matrix exponential derivative	Medium

Translation Statistics by Family

Family	Translated	Pending	Coverage
AB (Analysis)	44	5	90%
AG (Generalized)	4	1	80%
BB/BD (Benchmarks)	6	0	100%
FB/FD (Filtering)	6	0	100%
IB (Identification)	17	0	100%
MA (Matrix utility)	38	0	100%
MB (Matrix ops)	257	27	90%
MC (Polynomials)	3	17	15%
NF (Nonlinear)	20	1	95%
SB (Synthesis)	72	45	62%
SG (Generalized)	14	6	70%
TB/TC/TD/TF (Transfer)	30	18	63%
TG (Descriptor)	4	26	13%
UD/UE (Utility)	0	7	0%

Getting Started

# Install
pip install slicot

# Example: Solve Riccati equation
import numpy as np
from slicot import sb02od

A = np.array([[0, 1], [-2, -3]], dtype=np.float64, order='F')
B = np.array([[0], [1]], dtype=np.float64, order='F')
Q = np.eye(2, dtype=np.float64, order='F')
R = np.array([[1]], dtype=np.float64, order='F')

X, rcond, wr, wi, S, U, info = sb02od('D', 'B', 'N', 'U', 'S', 'N', 2, 1, A, B, Q, R, L=None)
print(f"Riccati solution X:\n{X}")

bavcol · 2026-01-19T10:32:00Z

bavcol
Jan 19, 2026

So with this, python-control can get rid of the fortran dependency and the routines are becoming availabe with permissive open source licence, right?

2 replies

jamestjsp Jan 19, 2026
Author

Yes, it inherits the same license of the https://github.com/SLICOT/SLICOT-Reference which is BSD-3-Clause license

bnavigator Jan 31, 2026
Collaborator

The new license is definitely a plus. But the installation problems with slycot come more from the BLAS/LAPACK dependency than from the Fortran language.

jamestjsp · 2026-01-24T14:41:14Z

jamestjsp
Jan 24, 2026
Author

It is done and published a python package. Windows users can enjoy pip install.
https://pypi.org/project/slicot/
https://github.com/jamestjsp/slicot

OpenBLAS is bundled in the wheels (~10MB each). Users just run: pip install slicot
Only if installing from source (pip install slicot --no-binary :all:) would they need BLAS dev libraries.

0 replies

bnavigator · 2026-01-31T16:09:11Z

bnavigator
Jan 31, 2026
Collaborator

Out of curiosity, would it be possible to benchmark this and compare it to Fortran Slycot?

0 replies

jamestjsp · 2026-01-31T21:03:29Z

jamestjsp
Jan 31, 2026
Author

@bnavigator Benchmarks available at jamestjsp/ctrlsys#8

Summary (Apple M1 Pro, same BLAS backend):

C11 is 7.4x faster on average (63.4% of 71 tests)
BD01AD (Descriptor): up to 49.5x speedup
BB02AD (Discrete Riccati): 6.4x speedup
Fortran wins on some Lyapunov routines (~18% faster)

Benchmark code: https://github.com/jamestjsp/slicot/tree/main/benchmarks

0 replies

jamestjsp · 2026-01-31T21:05:20Z

jamestjsp
Jan 31, 2026
Author

@bnavigator The BLAS/LAPACK dependency is bundled in the wheels via scipy-openblas32 (~10MB per wheel). Users just pip install slicot - no system BLAS needed.

Wheel builds: https://github.com/jamestjsp/slicot/blob/main/pyproject.toml#L45-L68

Platform coverage:

Linux: manylinux_2_28, musllinux_1_2 (x86_64, aarch64)
macOS: arm64
Windows: AMD64

0 replies

jamestjsp · 2026-01-31T21:11:42Z

jamestjsp
Jan 31, 2026
Author

Correction on platform coverage:

Linux: manylinux_2_28, musllinux_1_2 (x86_64, aarch64)
macOS: arm64
Windows: AMD64 only (no ARM64 yet)

0 replies

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

SLICOT `C11 translation with Python wrapper` progress update #1192

Uh oh!

{{title}}

Uh oh!

Replies: 6 comments 2 replies

Uh oh!

{{title}}

Uh oh!

Uh oh!

{{title}}

Uh oh!

Uh oh!

{{title}}

Uh oh!

Uh oh!

{{title}}

Uh oh!

Uh oh!

{{title}}

Uh oh!

Uh oh!

{{title}}

Uh oh!

Uh oh!

{{title}}

Uh oh!

Uh oh!

{{editor}}'s edit

{{editor}}'s edit

Uh oh!

Uh oh!

{{title}}

Uh oh!

Select a reply

Uh oh!

SLICOT C11 translation with Python wrapper progress update #1192

Uh oh!

jamestjsp Jan 17, 2026

SLICOT C11 Translation Status

Summary

Well-Tested Translated Routines

State-Space Analysis & Transformation (AB family)

Controller Design (SB family)

Descriptor Systems (SG/TG family)

System Identification (IB family)

Transfer Functions (TB/TC/TD/TF family)

Matrix Computations (MB family)

Filtering (FB family)

Polynomial Operations (MC family)

Notable Routines Pending Translation

High-Priority Controller Design

Spectral Factorization & Coprime Methods

Model Reduction

Polynomial & Matrix Fraction Operations

Descriptor System Transformations

Transfer Function Conversion

Complex/Periodic Matrix Operations

Translation Statistics by Family

Getting Started

Replies: 6 comments · 2 replies

Uh oh!

bavcol Jan 19, 2026

Uh oh!

jamestjsp Jan 19, 2026 Author

Uh oh!

bnavigator Jan 31, 2026 Collaborator

Uh oh!

jamestjsp Jan 24, 2026 Author

Uh oh!

bnavigator Jan 31, 2026 Collaborator

Uh oh!

jamestjsp Jan 31, 2026 Author

Uh oh!

Uh oh!

jamestjsp Jan 31, 2026 Author

Uh oh!

jamestjsp Jan 31, 2026 Author

SLICOT `C11 translation with Python wrapper` progress update #1192

jamestjsp
Jan 17, 2026

Replies: 6 comments 2 replies

bavcol
Jan 19, 2026

jamestjsp Jan 19, 2026
Author

bnavigator Jan 31, 2026
Collaborator

jamestjsp
Jan 24, 2026
Author

bnavigator
Jan 31, 2026
Collaborator

jamestjsp
Jan 31, 2026
Author

jamestjsp
Jan 31, 2026
Author

jamestjsp
Jan 31, 2026
Author