Spectrum-Adapted Polynomial Approximation for Matrix Functions

Abstract

We propose and investigate two new methods to approximate f( A) b for large, sparse, Hermitian matrices A. The main idea behind both methods is to first estimate the spectral density of A, and then find polynomials of a fixed order that better approximate the function f on areas of the spectrum with a higher density of eigenvalues. Compared to state-of-the-art methods such as the Lanczos method and truncated Chebyshev expansion, the proposed methods tend to provide more accurate approximations of f( A) b at lower polynomial orders, and for matrices A with a large number of distinct interior eigenvalues and a small spectral width.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…