Michael Ray, Yves C. Martin
Proceedings of SPIE - The International Society for Optical Engineering
Inverse iteration is widely used to compute the eigenvectors of a matrix once accurate eigenvalues are known. We discuss various issues involved in any implementation of inverse iteration for real, symmetric matrices. Current implementations resort to reorthogonalization when eigenvalues agree to more than three digits relative to the norm. Such reorthogonalization can have unexpected consequences. Indeed, as we show in this paper, the implementations in EISPACK and LAPACK may fail. We illustrate with both theoretical and empirical failures.
Michael Ray, Yves C. Martin
Proceedings of SPIE - The International Society for Optical Engineering
Salvatore Certo, Anh Pham, et al.
Quantum Machine Intelligence
Moutaz Fakhry, Yuri Granik, et al.
SPIE Photomask Technology + EUV Lithography 2011
I.K. Pour, D.J. Krajnovich, et al.
SPIE Optical Materials for High Average Power Lasers 1992