Repository logo

Infoscience

  • English
  • French
Log In
Logo EPFL, École polytechnique fédérale de Lausanne

Infoscience

  • English
  • French
Log In
  1. Home
  2. Academic and Research Output
  3. Journal articles
  4. A Low-Memory Lanczos Method with Rational Krylov Compression for Matrix Functions
 
research article

A Low-Memory Lanczos Method with Rational Krylov Compression for Matrix Functions

Casulli, Angelo A.
•
Simunec, Igor  
May 2, 2025
SIAM Journal on Scientific Computing

In this work, we introduce a memory-efficient method for computing the action of a Hermitian matrix function on a vector. Our method consists of a rational Lanczos algorithm combined with a basis compression procedure based on rational Krylov subspaces that only involve small matrices. The cost of the compression procedure is negligible with respect to the cost of the Lanczos algorithm. This enables us to avoid storing the whole Krylov basis, leading to substantial reductions in memory requirements. This method is particularly effective when the rational Lanczos algorithm needs a significant number of iterations to converge and each iteration involves a low computational effort. This scenario often occurs when polynomial Lanczos as well as extended and shift-and-invert Lanczos are employed. Theoretical results prove that, for a wide variety of functions, the proposed algorithm differs from rational Lanczos by an error term that is usually negligible. The algorithm is compared with other low-memory Krylov methods from the literature on a variety of test problems, showing competitive performance.

  • Details
  • Metrics
Type
research article
DOI
10.1137/24m1644699
Author(s)
Casulli, Angelo A.

Gran Sasso Science Institute, 67100 L’Aquila, Italy.

Simunec, Igor  

École Polytechnique Fédérale de Lausanne

Date Issued

2025-05-02

Publisher

Society for Industrial & Applied Mathematics (SIAM)

Published in
SIAM Journal on Scientific Computing
Volume

47

Issue

3

Start page

A1358

End page

A1382

Editorial or Peer reviewed

REVIEWED

Written at

EPFL

EPFL units
HPNALGS  
FunderFunding(s)Grant NumberGrant URL

INdAM-GNCS

Italian Ministry of University and Research

20227PCCKZ

RelationRelated workURL/DOI

IsNewVersionOf

A low-memory Lanczos method with rational Krylov compression for matrix functions

https://doi.org/10.48550/arXiv.2403.04390

IsSupplementedBy

[Code] Ratkrylov-compress-matfun

https://github.com/casulli/ratkrylov-compress-matfun
Available on Infoscience
May 5, 2025
Use this identifier to reference this record
https://infoscience.epfl.ch/handle/20.500.14299/249782
Logo EPFL, École polytechnique fédérale de Lausanne
  • Contact
  • infoscience@epfl.ch

  • Follow us on Facebook
  • Follow us on Instagram
  • Follow us on LinkedIn
  • Follow us on X
  • Follow us on Youtube
AccessibilityLegal noticePrivacy policyCookie settingsEnd User AgreementGet helpFeedback

Infoscience is a service managed and provided by the Library and IT Services of EPFL. © EPFL, tous droits réservés