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.

Isira: Integrated Shift-Invert Residual Arnoldi Method for Graph Laplacian Matrices from Big Data

A Refined Shift-And-Invert Arnoldi Algorithm for Large Unsymmetric Generalized Eigenproblems

An Inexact Shift-And-Invert Arnoldi Algorithm for Large Non-Hermitian Generalised Toeplitz Eigenproblems

Riemannian Linearized Proximal Algorithms for Nonnegative Inverse Eigenvalue Problem.

Fast Calculation of Laplace-Beltrami Eigenproblems Via Subdivision Linear Subspace

A Refined Arnoldi Type Method for Large Scale Eigenvalue Problems

Scalable Algorithms for Laplacian Pseudo-inverse Computation

An Invert-Free Arnoldi Method for Computing Interior Eigenpairs of Large Matrices

The shift-and-invert Arnoldi method for singular matrix pencils

Inverse eigenvalue problem for Laplacian matrices of a graph

Generalized Block Lanczos Methods for Large Unsymmetric Eigenproblems

Absolute-value based preconditioner for complex-shifted Laplacian systems

The Shift-Invert Residual Arnoldi Method and the Jacobi-Davidson Method: Theory and Algorithms ∗

Sparsified Cholesky and Multigrid Solvers for Connection Laplacians

Parallel Eigenvalue Calculation Based On Multiple Shift-Invert Lanczos And Contour Integral Based Spectral Projection Method

Lean Algebraic Multigrid (LAMG): Fast Graph Laplacian Linear Solver

Eigenvalue-free Iterative Shrinkage-Thresholding Algorithm for Solving the Linear Inverse Problems

A Quantum Algorithm for Solving Eigenproblem of the Laplacian Matrix of a Fully Connected Weighted Graph

The High Relative Accuracy of Computations with Laplacian Matrices

A Globally Convergent Inexact Newton-Like Cayley Transform Method for Inverse Eigenvalue Problems

A Refined Iterative Algorithm Based on the Block Arnoldi Process for Large Unsymmetric Eigenproblems