Nakatsukasa, Yuji and Higham, Nicholas J. (2012) Backward stability of iterations for computing the polar decomposition. SIAM Journal on Matrix Analysis and Applications, 33 (2). pp. 460479. ISSN 08954798
This is the latest version of this item.
PDF
SML000460(1).pdf Download (362kB) 
Abstract
Among the many iterations available for computing the polar decomposition the most practically useful are the scaled Newton iteration and the recently proposed dynamically weighted Halley iteration. Effective ways to scale these and other iterations are known, but their numerical stability is much less well understood. In this work we show that a general iteration $X_{k+1} = f(X_k)$ for computing the unitary polar factor is backward stable under two conditions. The first condition requires that the iteration is implemented in a mixed backwardforward stable manner and the second requires that the mapping $f$ does not significantly decrease the size of any singular value relative to the largest singular value. Using this result we show that the dynamically weighted Halley iteration is backward stable when it is implemented using Householder QR factorization with column pivoting and either row pivoting or row sorting. We also prove the backward stability of the scaled Newton iteration under the assumption that matrix inverses are computed in a mixed backwardforward stable fashion; our proof is much shorter than a previous one of Kielbasinski and Zietak. We also use our analysis to explain the instability of the inverse Newton iteration and to show that the NewtonSchulz iteration is only conditionally stable. This work shows that by carefully blending perturbation analysis with rounding error analysis it is possible to produce a general result that can prove the backward stability or predict or explain the instability (as the case may be) of a wide range of practically interesting iterations for the polar decomposition.
Item Type:  Article 

Uncontrolled Keywords:  polar decomposition, Newton iteration, inverse Newton iteration, NewtonSchulz iteration, dynamically weighted Halley iteration, QR factorization, row pivoting, row sorting, column pivoting, backward error analysis, rounding error analysis, numerical stability 
Subjects:  MSC 2010, the AMS's Mathematics Subject Classification > 15 Linear and multilinear algebra; matrix theory MSC 2010, the AMS's Mathematics Subject Classification > 65 Numerical analysis 
Depositing User:  Yuji Nakatsukasa 
Date Deposited:  12 Jun 2012 
Last Modified:  20 Oct 2017 14:13 
URI:  http://eprints.maths.manchester.ac.uk/id/eprint/1839 
Available Versions of this Item

Backward stability of iterations for computing the polar decomposition. (deposited 02 Dec 2011)

Backward stability of iterations for computing the polar decomposition. (deposited 22 Feb 2012)
 Backward stability of iterations for computing the polar decomposition. (deposited 12 Jun 2012) [Currently Displayed]

Backward stability of iterations for computing the polar decomposition. (deposited 22 Feb 2012)
Actions (login required)
View Item 