`scipy.linalg._flapack:strsyl`

Signature

def strsyl(a, b, c, trana = N, tranb = N, isgn = 1, overwrite_c = 0)

Solve the Sylvester equation op(a) @ x + isgn * x @ op(b) = scale * c (LAPACK strsyl).

a : ndarray: Upper (quasi-)triangular matrix of shape (m, m), in Schur canonical form.
b : ndarray: Upper (quasi-)triangular matrix of shape (n, n), in Schur canonical form.
c : ndarray: Right-hand side of shape (m, n).
trana : str, optional: 'N', 'T' or 'C' for op(a). Default is 'N'.
tranb : str, optional: 'N', 'T' or 'C' for op(b). Default is 'N'.
isgn : int, optional: Sign of the x @ op(b) term: 1 or -1. Default is 1.
overwrite_c : int, optional: If nonzero, c may be overwritten in place. Default is 0.

x : ndarray: The solution, overwriting c.
scale : float: Scale factor, at most 1, chosen to avoid overflow in the solution.
info : int: 0 on success; if negative, the -info-th argument had an illegal value; 1 if a and b have common or very close eigenvalues, in which case perturbed values were used.