A parallel algorithm for computing the polar decomposition

Higham, Nicholas J. and Papadimitriou, Pythagoras (1994) A parallel algorithm for computing the polar decomposition. Elsevier, Parallel Computing, 20. pp. 1161-1173. ISSN 0024-3795

[thumbnail of sdarticle[2].pdf] PDF
Restricted to Repository staff only

Download (738kB)


The polar decomposition A = UH of a rectangular matrix A, where U is unitary and H is Hermitian positive semidefinite, is an important tool in various applications, including aerospace computations, factor analysis and signal processing. We consider a pth order iteration for computing U that involves p independent matrix inversions per step and which is hence very amenable to parallel computation. We show that scaling the iterates speeds convergence of the iteration but makes the iteration only conditionally stable, with the backward error typically κappa_2(A) times bigger than the unit roundoff. In our implementation of the iteration on the Kendall Square Research KSR1 virtaul shared memory MIMD computer we take p to be the number of processors (p ≤ 16 in our experiments). Our code is found to be significantly faster than two existing techniques for computing the polar decomposition: one a Newton iteration, the other based on the singular value decomposition.

Item Type: Article
Uncontrolled Keywords: Polar decomposition; Singular value decomposition; Numerical stability; LA-PACK; Level 3 BLAS;Kendall Square Research KSR1 computer
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: Ms Lucy van Russelt
Date Deposited: 28 Jun 2006
Last Modified: 20 Oct 2017 14:12
URI: https://eprints.maths.manchester.ac.uk/id/eprint/329

Actions (login required)

View Item View Item