A comparison of limited-memory Krylov methods for Stieltjes functions of Hermitian matrices

06/10/2020
by   Stefan Güttel, et al.
0

Given a limited amount of memory and a target accuracy, we propose and compare several polynomial Krylov methods for the approximation of f(A)b, the action of a Stieltjes matrix function of a large Hermitian matrix on a vector. Using new error bounds and estimates, as well as existing results, we derive predictions of the practical performance of the methods, and rank them accordingly. As by-products, we derive new results on inexact Krylov iterations for matrix functions in order to allow for a fair comparison of rational Krylov methods with polynomial inner solves.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
02/05/2020

Limited-memory polynomial methods for large-scale matrix functions

Matrix functions are a central topic of linear algebra, and problems req...
research
06/10/2023

Sketched and truncated polynomial Krylov methods: Evaluation of matrix functions

Among randomized numerical linear algebra strategies, so-called sketchin...
research
07/02/2019

Modeling and Computation of Kubo Conductivity for 2D Incommensurate Bilayers

This paper presents a unified approach to the modeling and computation o...
research
07/26/2021

Computation of generalized matrix functions with rational Krylov methods

We present a class of algorithms based on rational Krylov methods to com...
research
01/31/2020

A study of defect-based error estimates for the Krylov approximation of φ-functions

Prior recent work, devoted to the study of polynomial Krylov techniques ...
research
08/06/2019

Rational Krylov for Stieltjes matrix functions: convergence and pole selection

Evaluating the action of a matrix function on a vector, that is x=f( M)v...
research
01/18/2021

Randomised preconditioning for the forcing formulation of weak constraint 4D-Var

There is growing awareness that errors in the model equations cannot be ...

Please sign up or login with your details

Forgot password? Click here to reset