Skip to content

Fix test sg03ad #82

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
bnavigator wants to merge 5 commits into from
Closed

Fix test sg03ad #82

Changes from 1 commit
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
fd830ee
sg03ad: arange instead of range
Nov 8, 2019
d173984
remove unused module, clean pep8 warning
Nov 8, 2019
df3b944
fix sg03ad test for continuos time and add discrete time test
bnavigator Nov 15, 2019
caecd0b
docstrings for sg03ad tests
bnavigator Nov 15, 2019
7f393f2
Merge branch 'master' into fix_test_sg03ad_static
bnavigator Apr 3, 2020
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
fix sg03ad test for continuos time and add discrete time test
  • Loading branch information
@bnavigator
bnavigator committed Nov 15, 2019
commit df3b944e4ae8a5f2241f8cc354e4423d8d5ffc42
23 changes: 19 additions & 4 deletions slycot/tests/test_sg03ad.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,20 +16,35 @@

class test_sg03ad(unittest.TestCase):

def test_sg03ad_a(self):
def test_sg03ad_1c(self):
# Example 1
n = 100
Xref = np.ones((n, n))
U = np.tril(Xref)
for t in range(0, 50, 10):
A = 2.0**(-t) - np.eye(n) + np.diag(np.arange(1, n+1)) + U.T
E = np.eye(n) + 2.0**(-t)*U
A = (2**(-t) - 1) * np.eye(n) + np.diag(np.arange(1., n+1.)) + U.T
E = np.eye(n) + 2**(-t) * U
Y = A.T.dot(Xref).dot(E) + E.T.dot(Xref).dot(A)
Q = np.zeros((n, n))
Z = np.zeros((n, n))
A, E, Q, Z, X, scale, sep, ferr, alphar, alphai, beta = \
synthesis.sg03ad('C', 'B', 'N', 'N', 'L', n, A, E, Q, Z, Y)
assert_almost_equal(Xref, X)
assert_almost_equal(X, Xref)

def test_sg03ad_1d(self):
# Example 1
n = 100
Xref = np.ones((n, n))
U = np.tril(Xref)
for t in range(0, 50, 10):
A = 2**(-t) * np.eye(n) + np.diag(np.arange(1., n+1.)) + U.T
E = np.eye(n) + 2**(-t) * U
Y = A.T.dot(Xref).dot(A) - E.T.dot(Xref).dot(E)
Q = np.zeros((n, n))
Z = np.zeros((n, n))
A, E, Q, Z, X, scale, sep, ferr, alphar, alphai, beta = \
synthesis.sg03ad('D', 'B', 'N', 'N', 'L', n, A, E, Q, Z, Y)
assert_almost_equal(X, Xref)

def test_sg03ad_3(self):
n = 3
Expand Down