Abstract
Matrix eigenvalue problems are frequently encountered in engineering and science. In this paper, we propose an extended shift-invert residual Arnoldi method for solving some smaller eigenvalues and the corresponding eigenvectors of a large and symmetric matrix, based on the idea of perfect Krylov subspace method proposed by Bai and Miao. We prove that the projected subspace of the proposed method is the same as the perfect Krylov subspace in exact computation. Based on this result, we also prove the convergence of the proposed method. Numerical experiments show that this method is superior to the perfect Krylov subspace method and the shift-invert residual Arnoldi method.
Similar content being viewed by others
References
Axelsson O (1985) A survey of preconditioned iterative methods for linear systems of algebraic equations. BIT Numer Math 25:165–187
Bai Z-Z, Miao C-Q (2017) On local quadratic convergence of inexact simplified Jacobi–Davidson method. Linear Algebra Appl 520:215–241
Bai Z-Z, Miao C-Q (2017) On local quadratic convergence of inexact simplified Jacobi–Davidson method for interior eigenpairs of Hermitian eigenproblems. Appl Math Lett 72:23–28
Bai Z-Z, Miao C-Q (2019) Computing eigenpairs of Hermitian matrices in perfect Krylov subspaces. Numer Algorithm 82:1251–1277
Belkin M, Niyogi P (2003) Laplacian eigenmaps for dimensionality reduction and data representation. Neural Comput 15:1373–1396
Grimes RG, Lewis JG, Simon HD (1994) A shifted block Lanczos algorithm for solving sparse symmetric generalized eigenproblems. SIAM J Matrix Anal Appl 15:228–272
Huang W-Q, Lin W-W, Lu HH-S, Yau S-T (2019) iSIRA: Integrated shift-invert residual Arnoldi method for graph Laplacian matrices from big data. J Comput Appl Math 346:518–531
Jia Z-X, Li C (2014) Inner iterations in the shift-invert residual Arnoldi method and the Jacobi-Davidson method. Sci China Math 57:1733–1752
Kunegis J (2013) Konect: the koblenz network collection. In: Proceedings of the 22nd International Conference on World Wide Web ACM, pp 1343–1350
Lanczos C (1950) An iteration method for the solution of the eigenvalue problem of linear differential and integral operators. J Res Nat Bur Stand 45:255–282
Lee C-R (2007) Residual Arnoldi method, theory, package and experiments. Ph.D. thesis, Department of Computer Science, University of Maryland at College Park
Lee C-R, Stewart GW (2007) Analysis of the residual Arnoldi method. Technical report, Department of Computer Science, University of Maryland at College Park
Meerbergen K, Spence A (1997) Implicitly restarted Arnoldi with purification for the shift-invert transformation. Math Comput 66:667–689
Molitierno JJ (2012) Applications of combinatorial matrix theory to Laplacian matrices of graphs. CRC Press, Boca Raton
Morgan RB, Scott DS (1986) Generalizations of Davidson’s method for computing eigenvalues of sparse symmetric matrices. SIAM J Sci Stat Comput 7:817–825
Pang J, Cheung G, Ortega A, Au OC (2015) Optimal graph Laplacian regularization for natural image denoising. In: 2015 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), pp 2294–2298
Parlett BN (1998) The symmetric eigenvalue problem. SIAM, Philadelphia
Ruhe A (1984) Rational Krylov sequence methods for eigenvalue computation. Linear Algebra Appl 58:391–405
Saad Y (1980) On the rates of convergence of the Lanczos and the block-Lanczos methods. SLAM J Numer Anal 17:687–706
Saad Y (2000) Numerical methods for large eigenvalue problems, 2nd edn. SIAM, Philadelphia
Sleijpen GLG, Van der Vorst HA (2000) A Jacobi–Davidson iteration method for linear eigenvalue problems. SIAM Rev 42:267–293
Stewart GW (2001) Matrix algorithms, vol. II, eigensystems. SIAM, Philadelphia
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Jinyun Yuan.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Yue, SF., Zhang, JJ. An extended shift-invert residual Arnoldi method. Comp. Appl. Math. 40, 63 (2021). https://doi.org/10.1007/s40314-021-01444-3
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s40314-021-01444-3