@MISC{_oncomputing, author = {}, title = {On Computing an Eigenvector of a Tridiagonal Matrix}, year = {} }

Share

OpenURL

Abstract

We consider the solution of the homogeneous equation (J; I)x = 0 where J is a tridiagonal matrix, is a known eigenvalue, and x is the unknown eigenvector corresponding to. Since the system is under-determined, x could be obtained by setting xk = 1 and solving for the rest of the elements of x. This method is not entirely new and it can be traced back to the times of Cauchy (1829). In 1958, Wilkinson demonstrated that, in nite-precision arithmetic, the computed x is highly sensitive to the choice of k � the traditional practice of setting k =1ork = n can lead to disastrous results. We develop algorithms to nd optimal k which require a LDU and a UDL factorisation of J; I and are based on the theory developed by Fernando for general matrices. We have also discovered new formulae (valid also for more general Hessenberg matrices) for the determinant ofJ; I,whichgivebetter numerical results when the shifted matrix is nearly singular. These formulae could be used to compute eigenvalues (or to improve the accuracy of known estimates) based on standard zero nders such as Newton and Laguerre methods. The accuracy of the computed eigenvalues is crucial for obtaining small residuals for the computed eigenvectors. The algorithms for solving eigenproblems are embarrassingly parallel and hence suitable for modern architectures. 1