Approximations of expressions of the form If:=trace(W^Tf(A)W), where A ∈R^{m×m} is a large symmetric matrix, W∈R^{m×k} with k<<m, and f is a function, can be computed without evaluating f(A) by applying a few steps of the global block Lanczos method to A with initial block-vector W. This yields a partial global Lanczos decomposition of A. We show that for suitable functions f upper and lower bounds for If can be determined by exploiting the connection between the global block Lanczos method and Gauss-type quadrature rules. Our approach generalizes techniques advocated by Golub and Meurant for the standard Lanczos method (with block size one) to the global block Lanczos method. We describe applications to the computation of upper and lower bounds of the trace of f(A) and consider, in particular, the computation of upper and lower bounds for the Estrada index, which arises in network analysis. We also discuss an application to machine learning.

Bounding matrix functionals via partial global block Lanczos decomposition

RODRIGUEZ, GIUSEPPE;
2015

Abstract

Approximations of expressions of the form If:=trace(W^Tf(A)W), where A ∈R^{m×m} is a large symmetric matrix, W∈R^{m×k} with k<
Gauss quadrature; Global block Lanczos algorithm; Trace computation; Network analysis; Machine learning
File in questo prodotto:
File Dimensione Formato  
glanczos15.pdf

Solo gestori archivio

Tipologia: versione editoriale
Dimensione 469.35 kB
Formato Adobe PDF
469.35 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

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: http://hdl.handle.net/11584/58518
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 18
  • ???jsp.display-item.citation.isi??? 16
social impact