Yuquan Sun

Yuquan Sun, Fanghui Gong, Igor V. Ovchinnikov et al.

This paper proposes an efficient method for computing partial eigenvalues of large sparse matrices what can be called the inexact inverse power method (IIPM). It is similar to the inexact Rayleigh quotient method and inexact Jacobi-Davidson method that it uses only a low precision approximate solution for the inner iteration. But this method uses less memory than inexact Jacobi-Davidson method and has stronger convergence performance than inexact Rayleigh quotient method. We exemplify the advantages of IIPM by applying it to find the ground state in theory of stochastics. Here we need to solve hundreds of large-scale matrix. The computational results show that this approach is a particularly useful method.

Zhongxiao Jia, Yuquan Sun

We investigate the generalized second-order Arnoldi (GSOAR) method, a generalization of the SOAR method proposed by Bai and Su [{\em SIAM J. Matrix Anal. Appl.}, 26 (2005): 640--659.], and the Refined GSOAR (RGSOAR) method for the quadratic eigenvalue problem (QEP). The two methods use the GSOAR procedure to generate an orthonormal basis of a given generalized second-order Krylov subspace, and with such basis they project the QEP onto the subspace and compute the Ritz pairs and the refined Ritz pairs, respectively. We develop implicitly restarted GSOAR and RGSOAR algorithms, in which we propose certain exact and refined shifts for respective use within the two algorithms. Numerical experiments on real-world problems illustrate the efficiency of the restarted algorithms and the superiority of the restarted RGSOAR to the restarted GSOAR. The experiments also demonstrate that both IGSOAR and IRGSOAR generally perform much better than the implicitly restarted Arnoldi method applied to the corresponding linearization problems, in terms of the accuracy and the computational efficiency.

Liu Yang, Yuquan Sun, Fanghui Gong

In this paper, we first establish the convergence criteria of the residual iteration method for solving quadratic eigenvalue problem- s. We analyze the impact of shift point and the subspace expansion on the convergence of this method. In the process of expanding subspace, this method needs to solve a linear system at every step. For large scale problems in which the equations cannot be solved directly, we propose an inner and outer iteration version of the residual iteration method. The new method uses the iterative method to solve the equations and uses the approximate solution to expand the subspace. We analyze the relationship between inner and outer iterations and provide a quantita- tive criterion for the inner iteration which can ensure the convergence of the outer iteration. Finally, our numerical experiments provide proof of our analysis and demonstrate the effectiveness of the inexact residual iteration method.

FangHui Gong, Yuquan Sun

In this paper, a new shift strategy for the implicitly restarted generalized second-order Arnoldi (GSOAR) method is proposed. In implicitly restarted processes, we can get a $k$-step GSOAR decomposition from a $m$-step GSOAR decomposition by performing $p = m-k$ implicit shifted QR iterations. The problem of the implicitly restarted GSOAR is the mismatch between the number of shifts and the dimension of the subspace. There are $2p$ shifts for $p$ QR iterations. We use the shifts to filter out the unwanted information in the current subspace; when more shifts are used, one obtains a better updated subspace. But, if we use more than $p$ shifts, the structure of the GSOAR decomposition will be destroyed. We propose a novel method which can use all $2p$ candidates and preserve the special structure. The new method vastly enhances the overall efficiency of the algorithm. Numerical experiments illustrate the efficiency of every restart process.