An Efficient Parallel Version of the Householder-QL Matrix Diagonalisation Algorithm
An Efficient Parallel Version of the Householder-QL Matrix Diagonalisation Algorithm
In this paper we report an effective parallelisation of the Householder routine for the reduction of a real symmetric matrix to tri-diagonal form and the QL algorithm for the diagonalisation of the resulting matrix. The Householder algorithm scales like $\alpha N^3/P+\beta N^2 \log_2(P)$ and the QL algorithm like $\gamma N^2 + \delta N^3/P$ as the number of processors $P$ is increased for fixed problem size. The constant parameters $\alpha$, $\beta$, $\gamma$ and $\delta$ are obtained empirically. When the eigenvalues only are required the Householder method scales as above while the QL algorithm remains sequential. The code is implemented in c in conjunction with the Message Passing Interface (MPI) libraries and verified on a sixteen node IBM SP2 and for real matrices that occur in the simulation of properties of crystalline materials
311-319
Reeve, JS
dd909010-7d44-44ea-83fe-a09e4d492618
Heath, M
aee4c988-75c5-4fdb-b21f-b571bb277384
1999
Reeve, JS
dd909010-7d44-44ea-83fe-a09e4d492618
Heath, M
aee4c988-75c5-4fdb-b21f-b571bb277384
Reeve, JS and Heath, M
(1999)
An Efficient Parallel Version of the Householder-QL Matrix Diagonalisation Algorithm.
Parallel Computing, 25 (3), .
Abstract
In this paper we report an effective parallelisation of the Householder routine for the reduction of a real symmetric matrix to tri-diagonal form and the QL algorithm for the diagonalisation of the resulting matrix. The Householder algorithm scales like $\alpha N^3/P+\beta N^2 \log_2(P)$ and the QL algorithm like $\gamma N^2 + \delta N^3/P$ as the number of processors $P$ is increased for fixed problem size. The constant parameters $\alpha$, $\beta$, $\gamma$ and $\delta$ are obtained empirically. When the eigenvalues only are required the Householder method scales as above while the QL algorithm remains sequential. The code is implemented in c in conjunction with the Message Passing Interface (MPI) libraries and verified on a sixteen node IBM SP2 and for real matrices that occur in the simulation of properties of crystalline materials
More information
Published date: 1999
Organisations:
EEE
Identifiers
Local EPrints ID: 251943
URI: http://eprints.soton.ac.uk/id/eprint/251943
ISSN: 0167-8191
PURE UUID: a3249294-1382-42d0-9686-ffc9a6ac3db7
Catalogue record
Date deposited: 14 Apr 2000
Last modified: 14 Mar 2024 05:14
Export record
Contributors
Author:
JS Reeve
Author:
M Heath
Download statistics
Downloads from ePrints over the past year. Other digital versions may also be available to download e.g. from the publisher's website.
View more statistics