From f78257eeb19b6585daa6470ba124516fedec26cc Mon Sep 17 00:00:00 2001 From: "Art J. R. Pelling" Date: Wed, 5 Jan 2022 15:46:30 +0100 Subject: [PATCH] increase ldwork for dico=='D' --- slycot/synthesis.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/slycot/synthesis.py b/slycot/synthesis.py index c232c5ce..f2117bb3 100644 --- a/slycot/synthesis.py +++ b/slycot/synthesis.py @@ -788,7 +788,7 @@ def sb03md57(A, U=None, C=None, if C is None: C = _np.zeros((n, n)) if ldwork is None: - ldwork = max(2*n*n, 3*n) + ldwork = max(2*n*n, 3*n) if dico == 'C' else 2*n*n + 2*n if dico != 'C' and dico != 'D': raise SlycotParameterError('dico must be either D or C', -1) out = _wrapper.sb03md(dico, n, A, U, C,