Forward stable eigenvalue decomposition of rank-one modifications of diagonal matrices

From MaRDI portal




Abstract: We present a new algorithm for solving an eigenvalue problem for a real symmetric matrix which is a rank-one modification of a diagonal matrix. The algorithm computes each eigenvalue and all components of the corresponding eigenvector with high relative accuracy in O(n) operations. The algorithm is based on a shift-and-invert approach. Only a single element of the inverse of the shifted matrix eventually needs to be computed with double the working precision. Each eigenvalue and the corresponding eigenvector can be computed separately, which makes the algorithm adaptable for parallel computing. Our results extend to the complex Hermitian case. The algorithm is similar to the algorithm for solving the eigenvalue problem for real symmetric arrowhead matrices from: N. Jakovv{c}evi'{c}~Stor, I. Slapniv{c}ar and J. L. Barlow, {Accurate eigenvalue decomposition of real symmetric arrowhead matrices and applications}, Lin. Alg. Appl., 464 (2015).





Describes a project that uses

Uses Software





This page was built for publication: Forward stable eigenvalue decomposition of rank-one modifications of diagonal matrices

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q890573)