I noticed that if I have complex numbers in my A,B,C,D state space matrices, there's casting complex numbers to real in control.StateSpace(A,B,C,D) which does not appear to be right (?) if there's imaginary part. Is it intentional? Is there a way around it?

Specifically, I think this line is responsible: https://github.com/python-control/python-control/blob/main/control/statesp.py#L2170. Note dtype=float. I have matrices with dtype=complex128. Minimally reproducible example (matrices are completely made up here, of course the real state space I'm working with is more realistic but also much bigger):

from control import StateSpace
import numpy as np
A = np.array([[1+1j,1+1j],[1+1j,1+1j]])
B = np.array([[1+1j,1+1j],[1+1j,1+1j]])
C = np.array([[1+1j,1+1j],[1+1j,1+1j]])
D = np.array([[1+1j,1+1j],[1+1j,1+1j]])
ss = StateSpace(A,B,C,D)

results in
statesp.py:2170: ComplexWarning: Casting complex values to real discards the imaginary part arr = np.array(data, dtype=float)