Back to Search Start Over

A treatment of breakdowns and near breakdowns in a reduction of a matrix to upper [formula omitted]-Hessenberg form and related topics.

Authors :
Salam, Ahmed
Ben Kahla, Haithem
Source :
Journal of Computational & Applied Mathematics. Oct2019, Vol. 359, p16-29. 14p.
Publication Year :
2019

Abstract

The reduction of a matrix to an upper J -Hessenberg form is a crucial step in the S R -algorithm (which is a Q R -like algorithm), structure-preserving, for computing eigenvalues and vectors, for a class of structured matrices. This reduction may be handled via the algorithm JHESS or via the recent algorithm JHMSH and its variants. The main drawback of JHESS (or JHMSH) is that it may suffer from a fatal breakdown, causing a brutal stop of the computations and hence, the S R -algorithm does not run. JHESS may also encounter near-breakdowns, source of serious numerical instability. In this paper, we focus on these aspects. We first bring light on the necessary and sufficient condition for the existence of the S R -decomposition, which is intimately linked to J -Hessenberg reduction. Then we will derive a strategy for curing fatal breakdowns and also for treating near breakdowns. Hence, the J -Hessenberg form may be obtained. Numerical experiments are given, demonstrating the efficiency of our strategies to cure and treat breakdowns or near breakdowns. [ABSTRACT FROM AUTHOR]

Details

Language :
English
ISSN :
03770427
Volume :
359
Database :
Academic Search Index
Journal :
Journal of Computational & Applied Mathematics
Publication Type :
Academic Journal
Accession number :
136344237
Full Text :
https://doi.org/10.1016/j.cam.2019.03.011