Wei-Qiang Huang

2papers

2 Papers

NAOct 23, 2017
iSIRA: Integrated Shift-Invert Residual Arnoldi Method for Graph Laplacian Matrices from Big Data

Wei-Qiang Huang, Wen-Wei Lin, Henry Horng-Shing Lu et al.

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 graph Laplacian matrix is symmetric, positive semi-definite, and is usually large and sparse. Computing some smallest positive eigenvalues and corresponding eigenvectors is often of interest. 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 time consuming or even unavailable to factorize a large and sparse matrix arising from 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 an eignsolver based on the inexact residual Arnoldi method together with an implicit remedy of the singularity and an effective deflation for convergent eigenvalues. Numerical experiments reveal that the integrated eigensolver outperforms the classical Arnoldi/Lanczos method for computing some smallest positive eigeninformation provided the LU factorization is not available.

NAOct 24, 2013
A Novel Skew-Hamiltonian Isotropic Lanczos Algorithm for Spectral Conformal Parameterizations

Wei-Qiang Huang, Xianfeng David Gu, Wen-Wei Lin et al.

Numerous methods for computing conformal mesh paramterizations has been developed due to the vast applications in the field of geometry processing. Spectral conformal parameterization (SCP) is one of these methods to computing a quality conformal parameterization based on the spectral technique. SCP focus on a generalized eigenvalue problem (GEP) $L_{C}\mathbf{f} = λB\mathbf{f}$ whose eigenvector(s) associated with the smallest positive eigenvalue(s) will provide the parameterization result. This paper devotes to study a novel eigensolver for this GEP. Based on the structures of matrix pair $(L_{C},B)$, we show that this GEP can be transformed into a small-scaled compressed deating standard eigenvalue problem with a symmetric positive definite skew-Hamiltonian operator. We then propose a skew-Hamiltonian isotropic Lanczos algorithm (SHILA) to solve the reducing problem. Numerical experiments show that our compressed deating skill remove the inuence of the kernel of $L_{C}$ and transform the original problem to a more robust system. The novel SHILA method can effective avoid the disturbance of duplicate eigenvalues. As a result, our numerical eigensolver can accurately and efficiently compute the conformal parameterization based on the spectral model of SCP.