An E cient Algorithm for the Real Root and Symmetric Tridiagonal Eigenvalue Problems
نویسنده
چکیده
Given a univariate complex polynomial f(x) of degree n with rational coe cients expressed as a ratio of two integers < 2, the root problem is to nd all the roots of f(x) up to speci ed precision 2 . In this paper we assume the arithmetic model for computation. We give an improved algorithm for nding a well-isolated splitting interval and for fast root proximity veri cation. Using these results, we give an algorithm for the real root problem: where all the roots of the polynomial are real. Our real root algorithm has time cost of O(n log n(logn+ log b)); where b = m + . Our arithmetic time cost is thus O(n log n) even in the case of high precision b n. This is within a small polylog factor of optimality, thus (perhaps surprisingly) upper bounding the arithmetic complexity of the real root problem to nearly the same as basic arithmetic operations on polynomials. The symmetric tridiagonal problem is: given an n n symmetric tridiagonal matrix, with 3n nonzero rational entries each expressed as a ratio of two integers < 2, to nd all the eigenvalues up to speci ed precision 2 . Using known e cient reductions from the symmetric tridiagonal eigenvalue problem to the real root problem, we also get an O(n log n(logn+log b)) arithmetic time bound for the symmetric tridiagonal eigenvalue problem.
منابع مشابه
An O(n Log 3 N) Algorithm for the Real Root and Symmetric Tridiagonal Eigenvalue Problems an O(n Log 3 N) Algorithm for the Real Root and Symmetric Tridiagonal Eigenvalue Problems
Given a univariate complex polynomial f(x) of degree n with rational coeecients expressed as a ratio of two integers < 2 m , the root problem is to nd all the roots of f(x) up to speciied precision 2 ?. In this paper we assume the arithmetic model for computation. We give an algorithm for the real root problem: where all the roots of the polynomial are real. Our real root algorithm has time cos...
متن کاملOn the nonnegative inverse eigenvalue problem of traditional matrices
In this paper, at first for a given set of real or complex numbers $sigma$ with nonnegative summation, we introduce some special conditions that with them there is no nonnegative tridiagonal matrix in which $sigma$ is its spectrum. In continue we present some conditions for existence such nonnegative tridiagonal matrices.
متن کاملQuasi-Laguerre Iteration in Solving Symmetric Tridiagonal Eigenvalue Problems
In this article, the quasi-Laguerre iteration is established in the spirit of Laguerre's iteration for solving polynomial f with all real zeros. The new algorithm, which maintains the monotonicity and global convergence of the Laguerre iteration, no longer needs to evaluate f". The ultimate convergence rate is + 1. When applied to approximate the eigenvalues of a symmetric tridiagonal matrix, t...
متن کاملSplitting of Expanded Tridiagonal Matrices
The article addresses a regular splitting of tridiagonal matrices. The given tridiagonal matrix A is rst expanded to an equivalent matrix e A and then split as e A = B R for which B is block-diagonal and every eigenvalue of B R is zero, i.e., (M N) = 0. The optimal splitting technique is applicable to various algorithms that incorporate one-dimensional solves or their approximations. Examples c...
متن کاملAn algorithm for symmetric generalized inverse eigenvalue problems
Using QR-like decomposition with column pivoting and least squares techniques, we propose a new and ecient algorithm for solving symmetric generalized inverse eigenvalue problems, and give its locally quadratic convergence analysis. We also present some numerical experiments which illustrate the behaviour of our algorithm. Ó 1999 Elsevier Science Inc. All rights reserved. AMS classi®cation: 65...
متن کامل