Abstract:We consider generalizations of the Sylvester matrix equation, consisting of the sum of a Sylvester operator and a linear operator $\Pi$ with a particular structure. More precisely, the commutator of the matrix coefficients of the operator $\Pi$ and the Sylvester operator coefficients are assumed to be matrices with low rank. We show (under certain additional conditions) low-rank approximability of this problem, i.e., the solution to this matrix equation can be approximated with a low-rank matrix. Projection methods have successfully been used to solve other matrix equations with low-rank approximability. We propose a new projection method for this class of matrix equations. The choice of subspace is a crucial ingredient for any projection method for matrix equations. Our method is based on an adaption and extension of the extended Krylov subspace method for Sylvester equations. A constructive choice of the starting vector/block is derived from the low-rank commutators. We illustrate the effectiveness of our method by solving large-scale matrix equations arising from applications in control theory and the discretization of PDEs. The advantages of our approach in comparison to other methods are also illustrated.

From Sylvester's Determinant Identity to Cramer's Rule

An efficient numerical algorithm for solving linear systems with cyclic tridiagonal coefficient matrices

Divided Differences and a Unified Approach to Evaluating Determinants Involving Generalized Factorials

Yet Another Proof of Sylvester's Determinant Identity

Non-commutative Sylvester's determinantal identity

A global variant of the COCR method for the complex symmetric Sylvester matrix equation AX + XB = C

Symbolic algorithm for solving cyclic penta-diagonal linear systems

Computing Hermitian determinantal representations of hyperbolic curves

A note on Dodgson's determinant condensation algorithm

A new projection method for solving large Sylvester equations

Computing the Determinant of a Dense Matrix over Z

A combinatorial proof of Cramer's Rule

The relaxed gradient based iterative algorithm for solving the generalized coupled complex conjugate and transpose Sylvester matrix equations

On the Determinant Evaluation of Quasi Penta-Diagonal Matrices and Quasi Penta-Diagonal Toeplitz Matrices

A Parallelizable Method for Computing the Singular Subspaces of a General Matrix

A general alternating direction implicit iteration method for solving complex symmetric linear systems

Krylov methods for low-rank commuting generalized Sylvester equations

An Explicit Identity to solve sums of Powers of Complex Functions

A Global Transpose-Free Method with Quasi-Minimal Residual Strategy for the Sylvester Equations

A NEW VERSION OF THE SMITH METHOD FOR SOLVING SYLVESTER EQUATION AND DISCRETE-TIME SYLVESTER EQUATION

A Hybrid Parallel Tridiagonal Solver on Multi-core Architectures