Globally convergent Jacobi methods for positive definite matrix pairs.

The paper derives and investigates the Jacobi methods for the generalized eigenvalue problem Ax = λBx, where A is a symmetric and B is a symmetric positive definite matrix. The methods first “normalize” B to have the unit diagonal and then maintain that property during the iterative process. The global convergence is proved for all such methods. That result is obtained for the large class of generalized serial strategies from Hari and Begović Kovač (Trans. Numer. Anal. (ETNA) 47, 107-147, 2017). Preliminary numerical tests confirm a high relative accuracy of some of those methods, provided that both matrices are positive definite and the spectral condition numbers of ΔAAΔAand ΔBBΔBare small, for some nonsingular diagonal matrices ΔAand ΔB. [ABSTRACT FROM AUTHOR]