Monday 6-th May 2019
Bernhard Beckermann (Univ. Lille).

Title: Low-rank updates of matrix functions

Abstract: We consider the task of updating a matrix function f(A) when the square matrix A is subject to a low-rank modification. The approach proposed in this talk attains efficiency by projecting onto tensorized Krylov subspaces produced by matrix-vector multiplications with A and its adjoint. We prove the approximations obtained from m steps of the proposed methods are exact if f is a polynomial of degree at most m, and use this as a basis for proving a variety of convergence results, in particular for the matrix exponential and for Markov functions. We illustrate the performance of our method by considering various examples from network analysis, where our approach can be used to cheaply update centrality and communicability measures.

Joint work with M. Schweitzer and D. Kressner (EPFL).

Comments are closed.