pick 7ace4bc New function for LTI system norm computationpick 5103448 * Updated documentation of function norm * Added control/tests/sysnorm_test.pypick 4fb252e Update sysnorm.pypick 108817c Do not track changes in VS Code setup.pick 6683eb3 Lowered tolerances in tests.pick 32d38bf Added: * Use of warnings package. * Use routine statesp.linfnorm when slycot installed. * New routine internal _h2norm_slycot when slycot is installed.pick 793f0d6 Added: * type check when calling ct.norm * metod argument in ct.norm (slycot or scipy)pick fe22e33 Use function _slycot_or_scipy from ct.mateqn# Rebase 85328d9..fe22e33 onto 85328d9 (20 commands)## Commands:# p, pick <commit> = use commit# r, reword <commit> = use commit, but edit the commit message# e, edit <commit> = use commit, but stop for amending# s, squash <commit> = use commit, but meld into previous commit# f, fixup [-C | -c] <commit> = like "squash" but keep only the previous#                    commit's log message, unless -C is used, in which case#                    keep only this commit's message; -c is same as -C but#                    opens the editor# x, exec <command> = run command (the rest of the line) using shell# b, break = stop here (continue rebase later with 'git rebase --continue')# d, drop <commit> = remove commit# l, label <label> = label current HEAD with a name# t, reset <label> = reset HEAD to a label# m, merge [-C <commit> | -c <commit>] <label> [# <oneline>]#         create a merge commit using the original merge commit's#         message (or the oneline, if no original merge commit was#         specified); use -c <commit> to reword the commit message# u, update-ref <ref> = track a placeholder for the <ref> to be updated#                       to this position in the new commits. The <ref> is#                       updated at the end of the rebase## These lines can be re-ordered; they are executed from top to bottom.## If you remove a line here THAT COMMIT WILL BE LOST.## However, if you remove everything, the rebase will be aborted.#

ben@skylab:~/src/python-control> git statusOn branch system-normsYour branch is up to date with 'henriks76/system-norms'.nothing to commit, working tree cleanben@skylab:~/src/python-control> git diff system-norms2

diff --git a/control/tests/sysnorm_test.py b/control/tests/sysnorm_test.pyindex 917e98d..74e8de0 100644--- a/control/tests/sysnorm_test.py+++ b/control/tests/sysnorm_test.py@@ -13,35 +13,35 @@ def test_norm_1st_order_stable_system():     """First-order stable continuous-time system"""     s = ct.tf('s')     G1 = 1/(s+1)-    assert np.allclose(ct.norm(G1, p='inf', tol=1e-9), 1.0) # Comparison to norm computed in MATLAB-    assert np.allclose(ct.norm(G1, p=2), 0.707106781186547) # Comparison to norm computed in MATLAB+    assert np.allclose(ct.norm(G1, p='inf', tol=1e-9, print_warning=False), 1.0) # Comparison to norm computed in MATLAB+    assert np.allclose(ct.norm(G1, p=2, print_warning=False), 0.707106781186547) # Comparison to norm computed in MATLAB          Gd1 = ct.sample_system(G1, 0.1)-    assert np.allclose(ct.norm(Gd1, p='inf', tol=1e-9), 1.0) # Comparison to norm computed in MATLAB-    assert np.allclose(ct.norm(Gd1, p=2), 0.223513699524858) # Comparison to norm computed in MATLAB+    assert np.allclose(ct.norm(Gd1, p='inf', tol=1e-9, print_warning=False), 1.0) # Comparison to norm computed in MATLAB+    assert np.allclose(ct.norm(Gd1, p=2, print_warning=False), 0.223513699524858) # Comparison to norm computed in MATLAB   def test_norm_1st_order_unstable_system():     """First-order unstable continuous-time system"""     s = ct.tf('s')     G2 = 1/(1-s)-    assert np.allclose(ct.norm(G2, p='inf', tol=1e-9), 1.0) # Comparison to norm computed in MATLAB-    assert ct.norm(G2, p=2) == float('inf') # Comparison to norm computed in MATLAB+    assert np.allclose(ct.norm(G2, p='inf', tol=1e-9, print_warning=False), 1.0) # Comparison to norm computed in MATLAB+    assert ct.norm(G2, p=2, print_warning=False) == float('inf') # Comparison to norm computed in MATLAB          Gd2 = ct.sample_system(G2, 0.1)-    assert np.allclose(ct.norm(Gd2, p='inf', tol=1e-9), 1.0) # Comparison to norm computed in MATLAB-    assert ct.norm(Gd2, p=2) == float('inf') # Comparison to norm computed in MATLAB+    assert np.allclose(ct.norm(Gd2, p='inf', tol=1e-9, print_warning=False), 1.0) # Comparison to norm computed in MATLAB+    assert ct.norm(Gd2, p=2, print_warning=False) == float('inf') # Comparison to norm computed in MATLAB  def test_norm_2nd_order_system_imag_poles():     """Second-order continuous-time system with poles on imaginary axis"""     s = ct.tf('s')     G3 = 1/(s**2+1)-    assert ct.norm(G3, p='inf') == float('inf') # Comparison to norm computed in MATLAB-    assert ct.norm(G3, p=2) == float('inf') # Comparison to norm computed in MATLAB+    assert ct.norm(G3, p='inf', print_warning=False) == float('inf') # Comparison to norm computed in MATLAB+    assert ct.norm(G3, p=2, print_warning=False) == float('inf') # Comparison to norm computed in MATLAB          Gd3 = ct.sample_system(G3, 0.1)-    assert ct.norm(Gd3, p='inf') == float('inf') # Comparison to norm computed in MATLAB-    assert ct.norm(Gd3, p=2) == float('inf') # Comparison to norm computed in MATLAB+    assert ct.norm(Gd3, p='inf', print_warning=False) == float('inf') # Comparison to norm computed in MATLAB+    assert ct.norm(Gd3, p=2, print_warning=False) == float('inf') # Comparison to norm computed in MATLAB  def test_norm_3rd_order_mimo_system():     """Third-order stable MIMO continuous-time system"""@@ -55,9 +55,9 @@ def test_norm_3rd_order_mimo_system():                   [-0.863652821988714,  -1.214117043615409,  -0.006849328103348]])     D = np.zeros((2,2))     G4 = ct.ss(A,B,C,D) # Random system generated in MATLAB-    assert np.allclose(ct.norm(G4, p='inf', tol=1e-9), 4.276759162964244) # Comparison to norm computed in MATLAB-    assert np.allclose(ct.norm(G4, p=2), 2.237461821810309) # Comparison to norm computed in MATLAB+    assert np.allclose(ct.norm(G4, p='inf', tol=1e-9, print_warning=False), 4.276759162964244) # Comparison to norm computed in MATLAB+    assert np.allclose(ct.norm(G4, p=2, print_warning=False), 2.237461821810309) # Comparison to norm computed in MATLAB          Gd4 = ct.sample_system(G4, 0.1)-    assert np.allclose(ct.norm(Gd4, p='inf', tol=1e-9), 4.276759162964228) # Comparison to norm computed in MATLAB-    assert np.allclose(ct.norm(Gd4, p=2), 0.707434962289554) # Comparison to norm computed in MATLAB+    assert np.allclose(ct.norm(Gd4, p='inf', tol=1e-9, print_warning=False), 4.276759162964228) # Comparison to norm computed in MATLAB+    assert np.allclose(ct.norm(Gd4, p=2, print_warning=False), 0.707434962289554) # Comparison to norm computed in MATLAB