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])
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])
