Given a large square matrix A and a sufficiently regular function f so that f(A) is well defined, we are interested in the approximation of the leading singular values and corresponding left and right singular vectors of f(A), and in particular in the approximation of ‖f(A)‖, where ‖⋅‖ is the matrix norm induced by the Euclidean vector norm. Since neither f(A) nor f(A)v can be computed exactly, we introduce a new inexact Golub–Kahan–Lanczos bidiagonalization procedure, where the inexactness is related to the inaccuracy of the operations f(A)v, f(A)⁎v. Particular outer and inner stopping criteria are devised so as to cope with the lack of a true residual. Numerical experiments with the new algorithm on typical application problems are reported.

Gaaf, S.W., Simoncini, V. (2017). Approximating the leading singular triplets of a large matrix function. APPLIED NUMERICAL MATHEMATICS, 113, 26-43 [10.1016/j.apnum.2016.10.015].

Approximating the leading singular triplets of a large matrix function

SIMONCINI, VALERIA
2017

Abstract

Given a large square matrix A and a sufficiently regular function f so that f(A) is well defined, we are interested in the approximation of the leading singular values and corresponding left and right singular vectors of f(A), and in particular in the approximation of ‖f(A)‖, where ‖⋅‖ is the matrix norm induced by the Euclidean vector norm. Since neither f(A) nor f(A)v can be computed exactly, we introduce a new inexact Golub–Kahan–Lanczos bidiagonalization procedure, where the inexactness is related to the inaccuracy of the operations f(A)v, f(A)⁎v. Particular outer and inner stopping criteria are devised so as to cope with the lack of a true residual. Numerical experiments with the new algorithm on typical application problems are reported.
2017
Gaaf, S.W., Simoncini, V. (2017). Approximating the leading singular triplets of a large matrix function. APPLIED NUMERICAL MATHEMATICS, 113, 26-43 [10.1016/j.apnum.2016.10.015].
Gaaf, Sarah W.; Simoncini, Valeria
File in questo prodotto:
Eventuali allegati, non sono esposti

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11585/586013
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 5
  • ???jsp.display-item.citation.isi??? 5
social impact