The parallel multisection method for solving algebraic eigenproblem has been presented in recent years with the development of the parallel computers, but all the research work is limited in standard eigenproblems of ...The parallel multisection method for solving algebraic eigenproblem has been presented in recent years with the development of the parallel computers, but all the research work is limited in standard eigenproblems of symmetric tridiagonal matrix. The multisection method for solving the generalized eigenproblem applied significantly in many science and engineering domains has not been studied. The parallel region preserving multisection method (PRM for short) for solving generalized eigenproblems of large sparse and real symmetric matrix is presented in this paper. This method not only retains the advantages of the conventional determinant search method (DS for short), but also overcomes its disadvantages such as leaking roots and disconvergence. We have tested the method on the YH 1 vector computer, and compared it with the parallel region preserving determinant search method the parallel region preserving bisection method (PRB for short). The numerical results show that PRM has a higher speed up, for instance, it attains the speed up of 7.7 when the scale of the problem is 2 114 and the eigenpair found is 3, and PRM is superior to PRB when the scale of the problem is large.展开更多
In(relativistic)electronic structure methods,the quaternion matrix eigenvalue problem and the linear response(Bethe-Salpeter)eigenvalue problem for excitation energies are two frequently encoun-tered structured eigenv...In(relativistic)electronic structure methods,the quaternion matrix eigenvalue problem and the linear response(Bethe-Salpeter)eigenvalue problem for excitation energies are two frequently encoun-tered structured eigenvalue problems.While the former problem was thoroughly studied,the later problem in its most general form,namely,the complex case without assuming the positive definiteness of the electronic Hessian,was not fully understood.In view of their very similar mathematical structures,we examined these two problems from a unified point of view.We showed that the identification of Lie group structures for their eigenvectors provides a framework to design diagonalization algorithms as well as numerical optimizations techniques on the corresponding manifolds.By using the same reduction algorithm for the quaternion matrix eigenvalue problem,we provided a necessary and sufficient condition to characterize the different scenarios,where the eigenvalues of the original linear response eigenvalue problem are real,purely imaginary,or complex.The result can be viewed as a natural generalization of the well-known condition for the real matrix case.展开更多
文摘The parallel multisection method for solving algebraic eigenproblem has been presented in recent years with the development of the parallel computers, but all the research work is limited in standard eigenproblems of symmetric tridiagonal matrix. The multisection method for solving the generalized eigenproblem applied significantly in many science and engineering domains has not been studied. The parallel region preserving multisection method (PRM for short) for solving generalized eigenproblems of large sparse and real symmetric matrix is presented in this paper. This method not only retains the advantages of the conventional determinant search method (DS for short), but also overcomes its disadvantages such as leaking roots and disconvergence. We have tested the method on the YH 1 vector computer, and compared it with the parallel region preserving determinant search method the parallel region preserving bisection method (PRB for short). The numerical results show that PRM has a higher speed up, for instance, it attains the speed up of 7.7 when the scale of the problem is 2 114 and the eigenpair found is 3, and PRM is superior to PRB when the scale of the problem is large.
基金supported by the National Natural Science Foundation of China (No.21973003)the Beijing Normal University Startup Package
文摘In(relativistic)electronic structure methods,the quaternion matrix eigenvalue problem and the linear response(Bethe-Salpeter)eigenvalue problem for excitation energies are two frequently encoun-tered structured eigenvalue problems.While the former problem was thoroughly studied,the later problem in its most general form,namely,the complex case without assuming the positive definiteness of the electronic Hessian,was not fully understood.In view of their very similar mathematical structures,we examined these two problems from a unified point of view.We showed that the identification of Lie group structures for their eigenvectors provides a framework to design diagonalization algorithms as well as numerical optimizations techniques on the corresponding manifolds.By using the same reduction algorithm for the quaternion matrix eigenvalue problem,we provided a necessary and sufficient condition to characterize the different scenarios,where the eigenvalues of the original linear response eigenvalue problem are real,purely imaginary,or complex.The result can be viewed as a natural generalization of the well-known condition for the real matrix case.