Isira: Integrated Shift-Invert Residual Arnoldi Method for Graph Laplacian Matrices from Big Data
Wei-Qiang Huang,Wen-Wei Lin,Henry Horng-Shing Lu,Shing-Tung Yau
DOI: https://doi.org/10.1016/j.cam.2018.07.031
IF: 2.872
2018-01-01
Journal of Computational and Applied Mathematics
Abstract:The eigenvalue problem of a graph Laplacian matrix L arising from a simple, connected and undirected graph has been given more attention due to its extensive applications, such as spectral clustering, community detection, complex network, image processing and so on. The associated matrix L is symmetric, positive semi-definite, and is usually large and sparse. It is often of interest for finding some smallest positive eigenvalues and corresponding eigenvectors.However, the singularity of L makes the classical eigensolvers inefficient since we need to factorize L for the purpose of solving large and sparse linear systems exactly. The next difficulty is that it is usually prohibitive to factorize L generated by real network problems from big data such as social media transactional databases, and sensor systems because there is in general not only local connections.In this paper, we propose a trimming to cure the singularity of L according to its special property: zero row/column sum. This remedy technique leads us to solve a positive definite linear system reduced in one dimension and then recover the result to get a suitable solution of the original system involved in an eigensolver. Besides, we apply a deflating approach to exclude the influence of converged eigenvalues. We show how to apply the idea of trimming to the graph Laplacian eigenvalue problem together with a deflated term and a target shift. Accordingly, based on the inexact residual Arnoldi (Lee, 2007; Lee and Stewart, 2007) method, we propose an integrated eigensolver for this kind of L combining with the implicit remedy of the singularity, an effective deflation for convergent eigenvalues and the shift–invert enhancement.Numerical experiments reveal that the integrated eigensolver outperforms the classical Arnoldi/Lanczos method for computing some smallest positive eigeninformation especially when the LU factorization is not available.