def radii(self, sys): # TODO: this recomputes the same control_gram multiple times over return [2 * np.sum(hankel(sub)) for sub in decompose_states(sys)]
def test_hankel(sys): assert np.allclose(hankel(sys), balreal(sys)[1])