Low-rank updates of matrix functions II: Rational Krylov methods

08/26/2020

∙

This work develops novel rational Krylov methods for updating a large-scale matrix function f(A) when A is subject to low-rank modifications. It extends our previous work in this context on polynomial Krylov methods, for which we present a simplified convergence analysis. For the rational case, our convergence analysis is based on an exactness result that is connected to work by Bernstein and Van Loan on rank-one updates of rational matrix functions. We demonstrate the usefulness of the derived error bounds for guiding the choice of poles in the rational Krylov method for the exponential function and Markov functions. Low-rank updates of the matrix sign function require additional attention; we develop and analyze a combination of our methods with a squaring trick for this purpose. A curious connection between such updates and existing rational Krylov subspace methods for Sylvester matrix equations is pointed out.

READ FULL TEXT

Low-rank updates of matrix functions II: Rational Krylov methods

Divide and conquer methods for functions of matrices with banded or hierarchical low-rank structure

Low-memory Krylov subspace methods for optimal rational matrix function approximation

Rational Krylov for Stieltjes matrix functions: convergence and pole selection

An effcient block rational Krylov solver for Sylvester equations with adaptive pole selection

Computing low-rank approximations of the Fréchet derivative of a matrix function using Krylov subspace methods

Interpolatory projection technique for Riccati-based feedback stabilization of index-1 descriptor systems

A review of the Separation Theorem of Chebyshev-Markov-Stieltjes for polynomial and some rational Krylov subspaces

Low-rank updates of matrix functions II: Rational Krylov methods

Related Research

Divide and conquer methods for functions of matrices with banded or hierarchical low-rank structure

Low-memory Krylov subspace methods for optimal rational matrix function approximation

Rational Krylov for Stieltjes matrix functions: convergence and pole selection

An effcient block rational Krylov solver for Sylvester equations with adaptive pole selection

Computing low-rank approximations of the Fréchet derivative of a matrix function using Krylov subspace methods

Interpolatory projection technique for Riccati-based feedback stabilization of index-1 descriptor systems

A review of the Separation Theorem of Chebyshev-Markov-Stieltjes for polynomial and some rational Krylov subspaces