## 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

## Abstract

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