Skip to content

Bugfix ab13bd #200

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

Merged
bnavigator merged 4 commits into from
Aug 24, 2023
Merged

Bugfix ab13bd #200

Changes from 1 commit
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
75b1fdf
Fix datatype and value check of parameter (m,p)
KybernetikJo Jul 19, 2023
14472a6
Add shallow copy to ab13bd, fix reference call of state space paramet…
KybernetikJo Jul 19, 2023
f93ef47
Remove shallow copy of A,B,C,D; Add intent(in,out,copy) in pyf
KybernetikJo Jul 20, 2023
823fda9
Add unittest to ad13bd
KybernetikJo Jul 30, 2023
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
Add shallow copy to ab13bd, fix reference call of state space paramet… 
…ers (A,B,C,D)
  • Loading branch information
@KybernetikJo
KybernetikJo committed Aug 23, 2023
commit 14472a620857f4628c61db775fcde9c500a505fa
10 changes: 8 additions & 2 deletions slycot/analysis.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1378,13 +1378,19 @@ def ab13bd(dico, jobn, n, m, p, A, B, C, D, tol = 0.0):
denominator of `G` (see the SLICOT subroutine SB08DD).
"""

out = _wrapper.ab13bd(dico, jobn, n, m, p, A, B, C, D, tol)

hidden = ' (hidden by the wrapper)'
arg_list = ('dico', 'jobn', 'n', 'm', 'p',
'A', 'lda' + hidden, 'B', 'ldb' + hidden, 'C', 'ldc' + hidden,
'D', 'ldd' + hidden, 'nq' + hidden,'tol', 'dwork' + hidden,
'ldwork' + hidden, 'iwarn', 'info')

a = A.copy()
b = B.copy()
c = C.copy()
d = D.copy()

out = _wrapper.ab13bd(dico, jobn, n, m, p, a, b, c, d, tol)

raise_if_slycot_error(out[-2:], arg_list, ab13bd.__doc__, locals())
return out[0]

Expand Down