{"about":{"site":"https://codewithpapers.app","non_affiliation":"Code with Papers and Syntology are not affiliated with, endorsed by, or sponsored by Papers with Code, Meta, or the pwc-archive mirror.","licence":"CC BY-SA 4.0","licence_url":"https://creativecommons.org/licenses/by-sa/4.0/legalcode","attribution":"https://codewithpapers.app/attribution","modified":"archive material modified by Syntology; see the attribution page"},"url":"/paper/lanczos-with-compression-for-symmetric-matrix","title":"Lanczos with compression for symmetric matrix Lyapunov equations","arxiv_id":"2505.22498","date":"2025-05-28","proceeding":null,"authors":["Angelo A. Casulli","Francesco Hrobat","Daniel Kressner"],"abstract":"This work considers large-scale Lyapunov matrix equations of the form $AX + XA = \\boldsymbol{c}\\boldsymbol{c}^T$, where $A$ is a symmetric positive definite matrix and $\\boldsymbol{c}$ is a vector. Motivated by the need to solve such equations in a wide range of applications, various numerical methods have been developed to compute low-rank approximations of the solution matrix $X$. In this work, we focus on the Lanczos method, which has the distinct advantage of requiring only matrix-vector products with $A$, making it broadly applicable. However, the Lanczos method may suffer from slow convergence when $A$ is ill-conditioned, leading to excessive memory requirements for storing the Krylov subspace basis generated by the algorithm. To address this issue, we propose a novel compression strategy for the Krylov subspace basis that significantly reduces memory usage without hindering convergence. This is supported by both numerical experiments and a convergence analysis. Our analysis also accounts for the loss of orthogonality due to round-off errors in the Lanczos process.","url_abs":"https://arxiv.org/abs/2505.22498v1","url_pdf":"https://arxiv.org/pdf/2505.22498v1.pdf","source":{"archive":"pwc-archive (Hugging Face), CC BY-SA 4.0","snapshot":"2025-07-28","licence_url":"https://creativecommons.org/licenses/by-sa/4.0/legalcode","row_kind":"links_only","authors_date_abstract":"arXiv metadata, CC0 1.0 (https://info.arxiv.org/help/license), from the Kaggle arXiv metadata snapshot of 2026-09-12"},"code_links":[{"paper_slug":"lanczos-with-compression-for-symmetric-matrix","repo_url":"https://github.com/fhrobat/lyap-compress","is_official":1,"mentioned_in_paper":0,"mentioned_in_github":0,"framework":"none","reach":null}],"tasks":[],"methods":[],"datasets_introduced":[],"methods_introduced":[],"results":[],"syntology":{"syntology_url":null,"atlas_url":null,"mcp":null,"developers":"https://syntology.ai/developers"},"arxiv_metadata":null,"syntology_extracted_results":null}