Abstract:Owing to the growing size of the eigenvalue problem and the growing number of eigenvalues desired, solution methods of iterative nature are becoming more popular than ever, which however suffer from low efficiency and lack of proper convergence criteria. In this paper, three efficient iterative eigenvalue algorithms are considered, i.e., subspace iteration method, iterative Ritz vector method and iterative Lanczos method based on the cell sparse fast solver and loop-unrolling. They are examined under the mode error criterion, i.e., the ratio of the out-of-balance nodal forces and the maximum elastic nodal point forces. Averagely speaking, the iterative Ritz vector method is the most efficient one among the three. Based on the mode error convergence criteria, the eigenvalue solvers are shown to be more stable than those based on eigenvalues only. Compared with ANSYS's subspace iteration and block Lanczos approaches, the subspace iteration presented here appears to be more efficient, while the Lanczos approach has roughly equal efficiency. The methods proposed are robust and efficient. Large size tests show that the improvement in terms of CPU time and storage is tremendous. Also reported is an aggressive shifting technique for the subspace iteration method, based on the mode error convergence criteria. A backward technique is introduced when the shift is not located in the right region. The efficiency of such a technique was demonstrated in the numerical tests.

Root-Finding with Eigen-Solving

A Parallel Algorithm For Finding All The Roots Of A Polynomial

Survey of Real Root Finding of Univariate Polynomial Equation in CAGD/CG

Advances in Solution of Classical Generalized Eigenvalue Problem

A novel Newton-Raphson style root finding algorithm

Finding roots of complex analytic functions via generalized colleague matrices

Highly Efficient Parallel Computing Schemes

Fast and backward stable computation of roots of polynomials, Part II: backward error analysis; companion matrix and companion pencil

Randomized methods for computing joint eigenvalues, with applications to multiparameter eigenvalue problems and root finding

Solving a Class of Nonlinear Eigenvalue Problems by Newton's Method

Fast and backward stable computation of the eigenvalues and eigenvectors of matrix polynomials

A robust iterative family for multiple roots of nonlinear equations: Enhancing accuracy and handling critical points

Root finding techniques that work

Root Radii and Subdivision for Polynomial Root-Finding

Simple and Nearly Optimal Polynomial Root-finding by Means of Root Radii Approximation

The Cell Index Concept in Fast Sparse Static and Eigenvalue Solvers for Finite Element Analysis

Fast Sparse Eigen Solution In Finite Element Analysis

Solving Rational Eigenvalue Problems Via Linearization

Solving large-scale nonlinear eigenvalue problems by rational interpolation approach and resolvent sampling based Rayleigh-Ritz method

Polynomial Root Isolation by Means of Root Radii Approximation

Real root refinements for univariate polynomial equations