Packed storage extension for ScaLAPACK
- ORNL
The authors describe a new extension to ScaLAPACK for computing with symmetric (Hermitian) matrices stored in a packed form. The new code is built upon the ScaLAPACK routines for full dense storage for a high degree of software reuse. The original ScaLAPACK stores a symmetric matrix as a full matrix but accesses only the lower or upper triangular part. The new code enables more efficient use of memory by storing only the lower or upper triangular part of a symmetric (Hermitian) matrix. The packed storage scheme distributes the matrix by block column panels. Within each panel, the matrix is stored as a regular ScaLAPACK matrix. This storage arrangement simplifies the subroutine interface and code reuse. Routines PxPPTRF/PxPPTRS implement the Cholesky factorization and solution for symmetric (Hermitian) linear systems in packed storage. Routines PxSPEV/PxSPEVX (PxHPEV/PxHPEVX) implement the computation of eigenvalues and eigenvectors for symmetric (Hermitian) matrices in packed storage. Routines PxSPGVX (PxHPGVX) implement the expert driver for the generalized eigenvalue problem for symmetric (Hermitian) matrices in packed storage. Performance results on the Intel Paragon suggest that the packed storage scheme incurs only a small time overhead over the full storage scheme.
- Research Organization:
- Oak Ridge National Lab., TN (US)
- Sponsoring Organization:
- USDOE Office of Science (US)
- DOE Contract Number:
- AC05-00OR22725
- OSTI ID:
- 754353
- Report Number(s):
- ORNL/TM-13545
- Country of Publication:
- United States
- Language:
- English
Similar Records
The design and implementation of the parallel out-of-core ScaLAPACK LU, QR and Cholesky factorization routines
Toward an efficient parallel eigensolver for dense symmetric matrices
The performance of finding eigenvalues and eigenvectors of dense symmetric matrices on distributed memory computers
Technical Report
·
Mon Mar 31 23:00:00 EST 1997
·
OSTI ID:296722
Toward an efficient parallel eigensolver for dense symmetric matrices
Journal Article
·
Thu Dec 31 23:00:00 EST 1998
· SIAM Journal on Scientific Computing
·
OSTI ID:355657
The performance of finding eigenvalues and eigenvectors of dense symmetric matrices on distributed memory computers
Conference
·
Thu Nov 30 23:00:00 EST 1995
·
OSTI ID:125572