Hostname: page-component-8448b6f56d-xtgtn Total loading time: 0 Render date: 2024-04-16T23:58:35.667Z Has data issue: false hasContentIssue false

Using randomization to make recursive matrix algorithms practical

Published online by Cambridge University Press:  01 November 1999

DINH LÊ
Affiliation:
Computer Science Department, University of California, Los Angeles, CA 90095-1596, USA (e-mail: stott@cs.ucla.edu)
D. STOTT PARKER
Affiliation:
Computer Science Department, University of California, Los Angeles, CA 90095-1596, USA (e-mail: stott@cs.ucla.edu)
Rights & Permissions [Opens in a new window]

Abstract

Core share and HTML view are not available for this content. However, as you have access to this content, a full PDF is available via the ‘Save PDF’ action button.

Recursive block decomposition algorithms (also known as quadtree algorithms when the blocks are all square) have been proposed to solve well-known problems such as matrix addition, multiplication, inversion, determinant computation, block LDU decomposition and Cholesky and QR factorization. Until now, such algorithms have been seen as impractical, since they require leading submatrices of the input matrix to be invertible (which is rarely guaranteed). We show how to randomize an input matrix to guarantee that submatrices meet these requirements, and to make recursive block decomposition methods practical on well-conditioned input matrices. The resulting algorithms are elegant, and we show the recursive programs can perform well for both dense and sparse matrices, although with randomization dense computations seem most practical. By ‘homogenizing’ the input, randomization provides a way to avoid degeneracy in numerical problems that permits simple recursive quadtree algorithms to solve these problems.

Type
Research Article
Copyright
© 1999 Cambridge University Press
Submit a response

Discussions

No Discussions have been published for this article.