Abstract:The solution of a large-scale Sylvester matrix equation plays an important role in control and large scientific computations. In this paper, we are interested in the large Sylvester matrix equation with large dimensionA and small dimension B, and a popular approach is to use the global Krylov subspace method. In this paper, we propose three new algorithms for this problem. We first consider the global GMRES algorithm with weighting strategy, which can be viewed as a precondition method. We present three new schemes to update the weighting matrix during iterations. Due to the growth of memory requirements and computational cost, it is necessary to restart the algorithm effectively. The deflation strategy is efficient for the solution of large linear systems and large eigenvalue problems; to the best of our knowledge, little work is done on applying deflation to the (weighted) global GMRES algorithm for large Sylvester matrix equations. We then consider how to combine the weighting strategy with deflated restarting, and propose a weighted global GMRES algorithm with deflation for solving large Sylvester matrix equations. In particular, we are interested in the global GMRES algorithm with deflation, which can be viewed as a special case when the weighted matrix is chosen as the identity. Theoretical analysis is given to show rationality of the new algorithms. Numerical experiments illustrate the numerical behavior of the proposed algorithms.

A new projection method for solving large Sylvester equations

A Projection Method Based on Extended Krylov Subspaces for Solving Sylvester Equations

Krylov Subspace Methods for the Generalized Sylvester Equation.

Krylov methods for low-rank commuting generalized Sylvester equations

Weighted and deflated global GMRES algorithms for solving large Sylvester matrix equations

Matrix Sign Function Methods for Solving Projected Generalized Continuous-Time Sylvester Equations.

Low-Rank Methods for Solving Discrete-Time Projected Lyapunov Equations

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

A Non-Krylov subspace Method for Solving Large and Sparse Linear System of Equations

A sketch-and-project method for solving the matrix equation AXB = C

Accurate solutions of M -matrix Sylvester equations

Restarted Generalized Krylov Subspace Methods for Solving Large-Scale Polynomial Eigenvalue Problems

Sketched and Truncated Polynomial Krylov Subspace Methods: Matrix Sylvester Equations

A Global Harmonic Arnoldi Method for Large Non-Hermitian Eigenproblems with an Application to Multiple Eigenvalue Problems

A General Approach to Sylvester-Polynomial-Conjugate Matrix Equations

A derivative-free memoryless BFGS hyperplane projection method for solving large-scale nonlinear monotone equations

Developing iterative algorithms to solve Sylvester tensor equations

Lifting Sylvester equations: singular value decay for non-normal coefficients

On Analysis Of Projection Methods For Rational Function Approximation To The Matrix Exponential

Preconditioning techniques for generalized Sylvester matrix equations

A class of multiplicative splitting iterations for solving the continuous Sylvester equation