squeeze: return shape of the frequency response.

[bnavigator]

Something strange about the return shape of the frequency response. For a square 2x2 system:

sys(np.array(0.0), squeeze=False).shape
(2, 2, 1)

sys(0.0, squeeze=False).shape
(2, 2)

Originally posted with the wrong result mapping by @forgi86 in #593 (comment)

[bnavigator]

BTW, omega has to be complex.

>>> import numpy as np
>>> from control import rss
>>> sys = rss(2,2,2)
>>> sys(np.array(0.0j), squeeze=False).shape
(2, 2, 1)
>>> sys(0.0j, squeeze=False).shape
(2, 2)

Now I am even more confused and will stop for this Sunday evening.

[murrayrm]

Sorry, accidentally closed.

[forgi86]

@bnavigator as far as I can see looks like you can pass real/complex/imaginary argument. I would expect the output to be the Laplace/Z transform for ct/dt in general, but haven't really tested!