Coverage decreased (-0.08%) to 75.927% when pullinga68ad6c on roryyorke:rory-ss-static-pole intocdd3e73 on python-control:master.

Reran this on Travis today, everything passes [1].

[1]https://travis-ci.org/roryyorke/python-control/builds/156620754

I will review it within 12 hours.

In315a1ea you add the comment "TODO: use super here?" above

LTI.__init__(self,inputs=D.shape[1],outputs=D.shape[0],dt=dt)

in the definition of__init__ ofStateSpace in control/statesp.py. I think that usingsuper here will ease code maintenance because it provides a reference to the parent without naming the parent, which is one of the motivations provided in the Python documentation athttps://docs.python.org/3/library/functions.html#super

However, changing to usesuper is orthogonal to this pull request, so we can wait until later to do so.

@roryyorke Did you write "BugFix?" in the first commit to indicate that the change might be thought of as adding a new feature, instead of fixing a bug?

For posterity and because I closed#107 in favor of this PR, I want to reference the comment#107 (comment) from there, which raises concerns aboutnumpy.matrix([]) being a 1\times 0 matrix instead of 0\times 0. A solution proposed in that comment and in one of the commits of that PR (e640b17) is to create a wrapper aroundnumpy.matrix.

Regarding the behavior ofnumpy.matrix([]), I did not find any written justification. It might be worthwhile for someone else to also search.

I studied the constructor code formatrix objects, which shows thatnumpy.matrix([]) returns a value with shape(1,0) because

1. none of theisinstance calls match for[] (typelist);
2. [] is passed throughnumpy.array, which returns anndarray of shape(0,) [1];
3. np.array([]).ndim == 1, so thematrix object is givenshape = (1, np.array([]).shape[0]) [2].

Thus, I guess that one design motivation is to map 1-dimensional arrays to row matrices, without specially treating empty arrays.

[1] line 270 of numpy/matrixlib/defmatrix.py,https://github.com/numpy/numpy/blob/2498b74420074dffa3a062ab7fc69dfeb88a5050/numpy/matrixlib/defmatrix.py#L270
[2] lines 277-278 of numpy/matrixlib/defmatrix.py,https://github.com/numpy/numpy/blob/2498b74420074dffa3a062ab7fc69dfeb88a5050/numpy/matrixlib/defmatrix.py#L277-L278

My recommendations:

• Assertself.states == 0 intest_remove_useless_states of statesp_test.py because StateSpace._remove_useless_states contains the assignmentself.states = self.A.shape[0], and one consequence of changes from your pull request is thatself.states = 0 when all states are useless (instead ofself.states = 1 from before).

• Remove the new import oftf2ss in statesp_test.py. If you can quickly rebase commits so as to entirely avoid the noise (i.e., accidental addition, then removal), please do. Otherwise, OK to just add a correction commit.

Once these are addressed, this PR is ready to be merged.

Thanks for the thorough review! I'll put all my responses here.

TODO for super: I think we agree this change doesn't belong in this PR; I suppose the Right Thing would have been to open a new issue.

"BugFix?": indeed, whether this is a bugfix or enhancement is not very clear. I can edit the log (I guess? never done that with git) if you like.

On empty SS objects: I have implemented this, unfortunately didn't push it to the PR. It follows the proposed solution. See [1].

tf2ss: I don't remember why that's there; some sort of exploratory debugging, I guess. Will remove. I'll see if I can do the rebase you suggest.

I'll add the recommended assertion, thanks.

OK, so for now I will:

• rebase this onc498d3d
• remove the ? from the "BugFix?" log message
• remove the unnecessary tf2ss import
• add the extra assertion

The empty SS object extension of [1] is on the same theme as the rest of the PR, but would require additional review. I could either add it here, or open a new PR for it.

[1]roryyorke@ba75ae5

Changes Unknown when pullingc84debb on roryyorke:rory-ss-static-pole into ** on python-control:master**.

I squashed the "messagefix" commit in the rebase, otherwise changes as in previous comment.

I squashed the messagefix commit, otherwise rebase is as described in previous comment.

Thanks; the rebase is good. To better organize discussion, I decided to leave consideration of a wrapper ofnumpy.matrix([]) (as in#110 (comment) androryyorke@ba75ae5) to another issue or PR.