Compress‐and‐restart block Krylov subspace methods for Sylvester matrix equations. Issue 1 (13th October 2020)

• Record Type: Journal Article
• Authors: Kressner, Daniel; Lund, Kathryn; Massei, Stefano; Palitta, Davide
• Is Part Of: Numerical linear algebra with applications. Volume 28:Issue 1(2021)
• Page Start: n/a
• Page End: n/a
• Publication Date: 2020-10-13
• Subjects: linear matrix equations -- block Krylov subspace methods -- low‐rank compression -- restarts; Algebras, Linear -- Periodicals; 512.5
• DOI: 10.1002/nla.2339
• Languages: English
• ISSNs: 1070-5325

Abstract

Summary: Block Krylov subspace methods (KSMs) comprise building blocks in many state‐of‐the‐art solvers for large‐scale matrix equations as they arise, for example, from the discretization of partial differential equations. While extended and rational block Krylov subspace methods provide a major reduction in iteration counts over polynomial block KSMs, they also require reliable solvers for the coefficient matrices, and these solvers are often iterative methods themselves. It is not hard to devise scenarios in which the available memory, and consequently the dimension of the Krylov subspace, is limited. In such scenarios for linear systems and eigenvalue problems, restarting is a well‐explored technique for mitigating memory constraints. In this work, such restarting techniques are applied to polynomial KSMs for matrix equations with a compression step to control the growing rank of the residual. An error analysis is also performed, leading to heuristics for dynamically adjusting the basis size in each restart cycle. A panel of numerical experiments demonstrates the effectiveness of the new method with respect to extended block KSMs.