Computation of generalized matrix functions with rational Krylov methods – Igor Simunec (Scuola Normale Superiore)


Venue: Aula Magna, Dipartimento di Matematica. Generalized matrix functions [3] are an extension of the notion of standard matrix functions to rectangular matrices, defined using the singular value decomposition instead of an eigenvalue decomposition. In this talk, we consider the computation of the action of a generalized matrix function on a vector and we present a class of algorithms based on rational Krylov methods [2]. These algorithms incorporate as a special case previous methods based on the Golub-Kahan bidiagonalization [1]. By exploiting the quasiseparable structure of the projected matrices, we show that the basis vectors can be updated using a short recurrence, which can be seen as a generalization to the rational case of the Golub-Kahan bidiagonalization. We also prove error bounds that relate the error of these methods to uniform rational approximation on an interval containing the singular values of the matrix. The effectiveness of the algorithms and the accuracy of the bounds is illustrated with numerical experiments. This is joint work with Angelo A. Casulli (Scuola Normale Superiore). [1] F. Arrigo, M. Benzi, and C. Fenu, Computation of generalized matrix functions, SIAM J. Matrix Anal. Appl. 37 (2016), no. 3, 836–860. [2] A. A. Casulli, I. Simunec, Computation of generalized matrix functions with rational Krylov methods, arXiv:2107.12074 (2021). [3] J. B. Hawkins and A. Ben-Israel, On generalized matrix functions, Linear and Multilinear Algebra 1 (1973), no. 2, 163–171.

Torna in cima