Abstract
In this paper a two year project whose principal goal is the development and implementation of efficient parallel algorithms for the computation of subsets of eigenvalues and associated eigenvectors of large, usually sparse, matrices is described. The eigenvalues may be the numerically largest, the numerically smallest, or those closest to a specified value. The main aims are to investigate and assess the comparative usefulness of new and existing algorithms for the solution of the partial eigenproblem, and to identify suitable parallel architectures on which they may be efficiently implemented. The project is currently in its first phase in which a detailed investigation of a number of variations of the basic Lanczos algorithm has been undertaken. A□standard version of the algorithm has been identified and implemented on two parallel machines and possibilities for enhancement of its performance based on the efficient monitoring of its computational progress are being studied. In this context a new version of the algorithm for the computation of the largest eigenvalue is proposed and evaluated.
Preview
Unable to display preview. Download preview PDF.
References
Kim, S.K., and A.T. Chronopoulos, A class of Lanczos-like algorithms implemented on parallel computers, Parallel Computing 17 (1991) 763–778.
Morgan, R. B., and D. S. Scott, Preconditioning the Lanczos Algorithm for Sparse Symmetric Eigenvalue Problems, SIAM J. Sci. Comput, Vol 14, No 3 (1993) 585–593.
Khelifi, M., Lanczos maximal algorithm for unsymmetric eigenvalue problems, Applied Numerical Mathematics 7 (1991) 179–193.
Kuczynski, J., and H. Wozniakowski, Estimating the Largest Eigenvalue by the Power and Lanczos Algorithms with a Random Start, Siam J. Matrix Anal Appl., Vol 13, No. 4 (1992) 1094–1122.
Jones, M.T., and M.L. Patrick, The Lanczos Algorithm for the Generalised Symmetric Eigenproblem on Shared-Memory Architectures, Preprint MCS-P182-0990 Mathematics and Computer Science Division, Argonne National Laboratory, December 1990.
Crouzeix, M., Philippe, B., and M. Sadkane, The Davidson method, Tech. Rep., Report TR/PA/90/45, CERFACS, Toulouse, 1990.
Morgan, R.B., and D.S. Scott, Generalizations of Davidson's method for computing eigenvalues of sparse symmetric matrices, SIAM J. Sci. Comput., Vol.7 (1986) 817–825.
Morgan, R. B., Generalizations of Davidson's Method for Computing Eigenvalues of Large Nonsymmetric Matrices, Journal of Computational Physics, 101 (1992) 287–291.
Sadkane, M., Block-Arnoldi and Davidson methods for unsymmetric large eigenvalue problems, Tech. Rep., Report TR/PA/92/45, CERFACS, Toulouse, 1990.
Arnoldi, W. E., The principle of minimized iterations in the solution of the matrix eigenvalue problem, Quart. Appl. Math., 9 (1951), 17–29.
Clint, M., and A. Jennings, The evaluation of eigenvalues and eigenvectors of real symmetric matrices by simultaneous iteration, Comput. J., (1970) 76–80.
Clint, M., and A. Jennings, A Simultaneous Iteration Method for the Unsymmetric Eigenproblem, J. Inst. Maths. Applics. 8, (1971) 111–121
Utku, S., and Y. Chang, Simultaneous Iteration Algorithm for General Eigenvalue Problems on Parallel Processors, in Proceedings of the 1986 International Conference on Parallel Processing, 19–22 August 1986, Hwang, K., Jacobs, S.M., and E.E. Swartzlander (Eds), IEEE Computer Society.
Stuart, E. J., and J. S. Weston, An Algorithm for the Parallel Computation of Subsets of Eigenvalues and Associated Eigenvectors of Large Symmetric Matrices using an Array Processor, in Proceedings Euromicro Workshop on Parallel and Distributed Processing, 27–29 January, 1993, Milligan, P., and A. Nunez (Eds.), IEEE Computer Society Press (1992), 211–217.
O'Leary, D. P., and P. Whitman, Parallel QR factorization by Householder and modified Gram-Schmidt algorithms, Parallel Computing, 16 (1991) 99–112.
Clint, M., and L. C. Waring, The orthogonalisation of small sets of very long vectors on massively parallel computers, (submitted to PARCO 1993) (1993).
Waring, L. C., and M. Clint, Parallel Gram-Schmidt orthogonalisation on a network of transputers, Parallel Computing 17 (1991) 1043–1050.
Weston, J. S., and M. Clint, Two algorithms for the parallel computation of eigenvalues and eigenvectors of large symmetric matrices using the ICL DAP, Parallel Computing 13 (1990) 281–288.
Weston, J. S., Clint, M., and C. W. Bleakney, The parallel computation of eigenvalues and eigenvectors of large Hermitian matrices using the AMT DAP 510, Concurrency: Practice and Experience, Vol 3(3) (1991) 179–185.
Clint, M., Weston, J. S., and C. W. Bleakney, A comparison of two Fortran dialects for expressing parallel solutions for a problem in linear algebra, Parallel Computing 18 (1992) 1325–1333.
Weston, J. S., The computation of eigensystems using an array processor, Ph D Thesis, The Queen's University of Belfast, 1988.
Clint, M., Weston, J. S., and C. W. Bleakney, Algorithms, languages and machines: An evaluation of a range of AMT DAP processors for the eigensolution of tridiagonal Hermitian matrices, Parallel Computing and Transputer Applications, Valero, M., Onate, E., Jane, M., Larriba J. L., and B. Suarez (Eds), IOS Press/CIMNE, Barcelona (1992) 1061–1069.
Parlett, B. N., and D. S. Scott, The Lanczos Algorithm with Selective Orthogonalisation, Mathematics of Computation, Vol 33, No 145 (1979) 217–238.
Basermann, A., and P. Weidner, A parallel algorithm for determining all eigenvalues of large real symmetric tridiagonal matrices, Parallel Computing 18 (1992) 1129–1141.
Dowell, M., and P. Jarratt, The “pegasus” method for computing the root of an equation, BIT 12 (1972) 503–508.
Lo, S., Philippe, B., and A. Sameh, A Multiprocessor Algorithm for the Symmetric Tridiagonal Eigenvalue Problem, SIAM J. Sci. Stat. Comput., 8 (1987) 155–165.
Krishnakumar, A.S., and M. Morf, Eigenvalues of a Symmetric Tridiagonal Matrix: A Divide-and-Conquer Approach', Numer. Math., 48 (1986) 349–368.
Author information
Authors and Affiliations
Editor information
Rights and permissions
Copyright information
© 1994 Springer-Verlag Berlin Heidelberg
About this paper
Cite this paper
Szularz, M., Weston, J., Murphy, K., Clint, M. (1994). Parallel algorithms for the partial eigensolution of large sparse matrices on novel architecture computers. In: Dongarra, J., Waśniewski, J. (eds) Parallel Scientific Computing. PARA 1994. Lecture Notes in Computer Science, vol 879. Springer, Berlin, Heidelberg. https://doi.org/10.1007/BFb0030174
Download citation
DOI: https://doi.org/10.1007/BFb0030174
Published:
Publisher Name: Springer, Berlin, Heidelberg
Print ISBN: 978-3-540-58712-5
Online ISBN: 978-3-540-49050-0
eBook Packages: Springer Book Archive