Approximating sparse matrices and their functions using matrix-vector products

Journal article

The computation of a matrix function 𝑓(𝐴) is an important task in scientific computing appearing in machine learning, network analysis and the solution of partial differential equations. In this work, we use only matrix-vector products 𝑥 ↦ 𝐴𝑥 to approximate functions of sparse matrices and matrices with similar structures such as sparse matrices 𝐴 themselves or matrices that have a similar decay property as matrix functions. We show that when 𝐴 is a sparse matrix with an unknown sparsity pattern, techniques from compressed sensing can be used under natural assumptions. Moreover, if 𝐴 is a banded matrix then certain deterministic matrix-vector products can efficiently recover the large entries of 𝑓(𝐴). We describe an algorithm for each of the two cases and give error analysis based on the decay bound for the entries of 𝑓(𝐴). We finish with numerical experiments showing the accuracy of our algorithms.

Authors: Park, T; Nakatsukasa, Y

• Publication status: Published
• Peer review status: Peer reviewed
• Publisher: Elsevier
• Journal: Applied and Computational Harmonic Analysis
• Volume: 83
• Article number: 101869
• Publication date: 2026-02-26
• Acceptance date: 2026-02-24
• DOI: 10.1016/j.acha.2026.101869
• EISSN: 1096-603X
• ISSN: 1063-5203
• Language: English
• Keywords: matrix function; banded matrix; sparse matrix; matrix-vector products; compressed sensing