A large unsymmetric linear system problem is transformed into the problem of computing the eigenvector of a large symmetric nonnegative definite matrix associated with the eigenvalue zero, i.e., the computation of the...A large unsymmetric linear system problem is transformed into the problem of computing the eigenvector of a large symmetric nonnegative definite matrix associated with the eigenvalue zero, i.e., the computation of the elgenvector of the cross-product matrix of an augmented matrix associated with the eigenvalue zero. The standard Lanczos method and an improved refined Lanczos method are proposed that compute approximate eigenvectors and return approximate solutions of the linear system. An implicitly restarted Lanczos algorithm and its refined version are developed. Theoretical analysis and numerical experiments show the refined method is better than the standard one. If the large matrix has small eigenvalues, the two new algorithms are much faster than the unpreconditioned restarted GMRES.展开更多
In this paper, we develop an implicitly restarted block Arnoldi algorithm in a vector-wise fashion. The vector-wise construction greatly simplifies both the detection of necessary deflation and the actual deflation it...In this paper, we develop an implicitly restarted block Arnoldi algorithm in a vector-wise fashion. The vector-wise construction greatly simplifies both the detection of necessary deflation and the actual deflation itself, so it is preferable to the block-wise construction. The numerical experiment shows that our algorithm is effective.展开更多
The singular value decomposition problem is mathematically equivalent to the eigenproblem of an argumented matrix. Golub et al. give a bidiagonalization Lanczos method for computing a number of largest or smallest sin...The singular value decomposition problem is mathematically equivalent to the eigenproblem of an argumented matrix. Golub et al. give a bidiagonalization Lanczos method for computing a number of largest or smallest singular values and corresponding singular vertors, but the method may encounter some convergence problems. In this paper we analyse the convergence of the method and show why it may fail to converge. To correct this possible nonconvergence, we propose a refined bidiagonalization Lanczos method and apply the implicitly restarting technique to it, and we then present an implicitly restarted bidiagonalization Lanczos algorithm(IRBL) and an implicitly restarted refined bidiagonalization Lanczos algorithm (IRRBL). A new implicitly restarting scheme and a reliable and efficient algorithm for computing refined shifts are developed for this special structure eigenproblem.Theoretical analysis and numerical experiments show that IRRBL performs much better than IRBL.展开更多
文摘A large unsymmetric linear system problem is transformed into the problem of computing the eigenvector of a large symmetric nonnegative definite matrix associated with the eigenvalue zero, i.e., the computation of the elgenvector of the cross-product matrix of an augmented matrix associated with the eigenvalue zero. The standard Lanczos method and an improved refined Lanczos method are proposed that compute approximate eigenvectors and return approximate solutions of the linear system. An implicitly restarted Lanczos algorithm and its refined version are developed. Theoretical analysis and numerical experiments show the refined method is better than the standard one. If the large matrix has small eigenvalues, the two new algorithms are much faster than the unpreconditioned restarted GMRES.
基金This work is supported by National Natural Science Foundation of China No. 10531080.
文摘In this paper, we develop an implicitly restarted block Arnoldi algorithm in a vector-wise fashion. The vector-wise construction greatly simplifies both the detection of necessary deflation and the actual deflation itself, so it is preferable to the block-wise construction. The numerical experiment shows that our algorithm is effective.
文摘The singular value decomposition problem is mathematically equivalent to the eigenproblem of an argumented matrix. Golub et al. give a bidiagonalization Lanczos method for computing a number of largest or smallest singular values and corresponding singular vertors, but the method may encounter some convergence problems. In this paper we analyse the convergence of the method and show why it may fail to converge. To correct this possible nonconvergence, we propose a refined bidiagonalization Lanczos method and apply the implicitly restarting technique to it, and we then present an implicitly restarted bidiagonalization Lanczos algorithm(IRBL) and an implicitly restarted refined bidiagonalization Lanczos algorithm (IRRBL). A new implicitly restarting scheme and a reliable and efficient algorithm for computing refined shifts are developed for this special structure eigenproblem.Theoretical analysis and numerical experiments show that IRRBL performs much better than IRBL.
