On iterative refinement for the spectral decomposition of symmetric matrices
Résumé
A method is studied for computation of the spectral decomposition of a symmetric matrix with high precision. The approach used is based upon careful calculation of a solution of the Riccati equations arising in the refinement process. Detailed rounding error analysis and estimates of convergence rates are carried out.