In the last few years many numerical techniques for computing eigenvalues of structured rank matrices have been proposed. Most of them are based on QR iterations since, in the symmetric case, the rank structure is preserved and high accuracy is guaranteed. In the unsymmetric case, however, the QR algorithm destroys the rank structure, which is instead preserved if LR iterations are used. We consider a wide class of quasiseparable matrices which can be represented in terms of the same parameters involved in their Neville factorization. This class, if assumptions are made to prevent possible breakdowns, is closed under LR steps. Moreover, we propose an implicit shifted LR method with a linear cost per step, which resembles the qd method for tridiagonal matrices. We show that for totally nonnegative quasiseparable matrices the algorithm is stable and breakdowns cannot occur, if the Laguerre shift, or other shift strategy preserving nonnegativity, is used. Computational evidence shows that good accuracy is obtained also when applied to symmetric positive definite matrices.

Qd-type methods for quasiseparable matrices

BEVILACQUA, ROBERTO;DEL CORSO, GIANNA MARIA
2011

Abstract

In the last few years many numerical techniques for computing eigenvalues of structured rank matrices have been proposed. Most of them are based on QR iterations since, in the symmetric case, the rank structure is preserved and high accuracy is guaranteed. In the unsymmetric case, however, the QR algorithm destroys the rank structure, which is instead preserved if LR iterations are used. We consider a wide class of quasiseparable matrices which can be represented in terms of the same parameters involved in their Neville factorization. This class, if assumptions are made to prevent possible breakdowns, is closed under LR steps. Moreover, we propose an implicit shifted LR method with a linear cost per step, which resembles the qd method for tridiagonal matrices. We show that for totally nonnegative quasiseparable matrices the algorithm is stable and breakdowns cannot occur, if the Laguerre shift, or other shift strategy preserving nonnegativity, is used. Computational evidence shows that good accuracy is obtained also when applied to symmetric positive definite matrices.
Bevilacqua, Roberto; Bozzo, Enrico; DEL CORSO, GIANNA MARIA
File in questo prodotto:
File Dimensione Formato  
79448_p3.pdf

accesso aperto

Tipologia: Versione finale editoriale
Licenza: Tutti i diritti riservati (All rights reserved)
Dimensione 338.33 kB
Formato Adobe PDF
338.33 kB Adobe PDF Visualizza/Apri

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/11568/200409
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 12
  • ???jsp.display-item.citation.isi??? 12
social impact