English

Convergence of the Cyclic and Quasi-cyclic Block Jacobi Methods

Numerical Analysis 2017-06-27 v3

Abstract

The paper studies the global convergence of the block Jacobi me\-thod for symmetric matrices. Given a symmetric matrix AA of order nn, the method generates a sequence of matrices by the rule A(k+1)=UkTA(k)UkA^{(k+1)}=U_k^TA^{(k)}U_k, k0k\geq0, where UkU_k are orthogonal elementary block matrices. A class of generalized serial pivot strategies is introduced, significantly enlarging the known class of weak wavefront strategies, and appropriate global convergence proofs are obtained. The results are phrased in the stronger form: S(A)cS(A)S(A')\leq c S(A), where AA' is the matrix obtained from AA after one full cycle, c<1c<1 is a constant and S(A)S(A) is the off-norm of AA. Hence, using the theory of block Jacobi operators, one can apply the obtained results to prove convergence of block Jacobi methods for other eigenvalue problems, such as the generalized eigenvalue problem. As an example, the results are applied to the block JJ-Jacobi method. Finally, all results are extended to the corresponding quasi-cyclic strategies.

Keywords

Cite

@article{arxiv.1604.05825,
  title  = {Convergence of the Cyclic and Quasi-cyclic Block Jacobi Methods},
  author = {Vjeran Hari and Erna Begovic},
  journal= {arXiv preprint arXiv:1604.05825},
  year   = {2017}
}

Comments

41 pages