The solution of large dense generalized eigenvalue problems on the Cray X-MP/24 with SSD
Journal of Computational Physics
The WY representation for products of householder matrices
SIAM Journal on Scientific and Statistical Computing - Papers from the Second Conference on Parallel Processing for Scientific Computin
A fully parallel algorithm for the symmetric eigenvalue problem
SIAM Journal on Scientific and Statistical Computing
A multiprocessor algorithm for the symmetric tridiagonal eigenvalue problem
SIAM Journal on Scientific and Statistical Computing
The symmetric eigenvalue problem
The symmetric eigenvalue problem
Basic Linear Algebra Subprograms for Fortran Usage
ACM Transactions on Mathematical Software (TOMS)
Banded Eigenvalue Solvers on Vector Machines
ACM Transactions on Mathematical Software (TOMS)
Proceedings of the fourth workshop on I/O in parallel and distributed systems: part of the federated computing research conference
A framework for symmetric band reduction
ACM Transactions on Mathematical Software (TOMS)
Proceedings of 2011 International Conference for High Performance Computing, Networking, Storage and Analysis
Toward a scalable multi-GPU eigensolver via compute-intensive kernels and efficient communication
Proceedings of the 27th international ACM conference on International conference on supercomputing
An improved parallel singular value algorithm and its implementation for multicore hardware
SC '13 Proceedings of the International Conference on High Performance Computing, Networking, Storage and Analysis
Hi-index | 0.00 |
This paper describes a new implementation of algorithms for solving large, dense symmetric eigen-problems AX = BX&Lgr;, where the matrices A and B are too large to fit in the central memory of the computer. Here A is assumed to be symmetric, and B symmetric positive definite. A combination of block Cholesky and block Householder transformations are used to reduce the problem to a symmetric banded eigenproblem whose eigenvalues can be computed in central memory. Inverse iteration is applied to the banded matrix to compute selected eigenvectors, which are then transformed back to eigenvectors of the original problem. This method is especially suitable for the solution of large eigenproblems arising in quantum physics, using a vector supercomputer with fast secondary storage device such as the Cray X-MP with SSD. Some numerical results demonstrate the efficiency of the new implementation.