In this paper we revisit the problem of performing a QR-step on an unreduced Hessenberg matrix H when we know an "exact" eigenvalue ?0 of H. Under exact arithmetic, this eigenvalue will appear on diagonal of the transformed Hessenberg matrix H~ and will be decoupled from the remaining part of the Hessenberg matrix, thus resulting in a deflation. But it is well known that in finite precision arithmetic the so-called perfect shift can get blurred and that the eigenvalue ?0 can then not be deflated and/or is perturbed significantly. In this paper, we develop a new strategy for computing such a QR step so that the deflation is almost always successful. We also show how to extend this technique to double QR-steps with complex conjugate shifts.

The QR Steps with Perfect Shifts

Nicola Mastronardi;
2018

Abstract

In this paper we revisit the problem of performing a QR-step on an unreduced Hessenberg matrix H when we know an "exact" eigenvalue ?0 of H. Under exact arithmetic, this eigenvalue will appear on diagonal of the transformed Hessenberg matrix H~ and will be decoupled from the remaining part of the Hessenberg matrix, thus resulting in a deflation. But it is well known that in finite precision arithmetic the so-called perfect shift can get blurred and that the eigenvalue ?0 can then not be deflated and/or is perturbed significantly. In this paper, we develop a new strategy for computing such a QR step so that the deflation is almost always successful. We also show how to extend this technique to double QR-steps with complex conjugate shifts.
2018
Istituto Applicazioni del Calcolo ''Mauro Picone''
Hessenberg form
QR step
eigenvalue problem
perfect shift
File in questo prodotto:
Non ci sono file associati a questo prodotto.

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/20.500.14243/370858
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 5
  • ???jsp.display-item.citation.isi??? 5
social impact