You are here: MIMS > EPrints
MIMS EPrints

2006.151: A parallel algorithm for computing the polar decomposition

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

Full text available as:

PDF - Archive staff only - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
721 Kb

DOI: 10.1016/0167-8191(94)90073-6


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 2000 > 15 Linear and multilinear algebra; matrix theory
MSC 2000 > 65 Numerical analysis
MIMS number:2006.151
Deposited By:Miss Louise Stait
Deposited On:28 June 2006

Download Statistics: last 4 weeks
Repository Staff Only: edit this item