Jacobi method for symmetric $4\times4$ matrices converges for every cyclic pivot strategy

Erna Begovic,Vjeran Hari
DOI: https://doi.org/10.1007/s11075-017-0396-8
2017-08-20
Abstract:The paper studies the global convergence of the Jacobi method for symmetric matrices of size $4$. We prove global convergence for all $720$ cyclic pivot strategies. Precisely, we show that inequality $S(A^{[t+3]})\leq\gamma S(A^{[t]})$, $t\geq1$, holds with the constant $\gamma<1$ that depends neither on the matrix $A$ nor on the pivot strategy. Here $A^{[t]}$ stands for the matrix obtained from $A$ after $t$ full cycles of the Jacobi method and $S(A)$ is the off-diagonal norm of $A$. We show why three consecutive cycles have to be considered. The result has a direct application on the $J$-Jacobi method.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the global convergence of the Jacobi method on symmetric 4×4 matrices. Specifically, the authors studied the Jacobi method under all possible cyclic pivot strategies and proved that these strategies can ensure the global convergence of the Jacobi method. ### Core problems of the paper 1. **Global convergence**: For any initial matrix \(A\), can the Jacobi method converge to a diagonal matrix after a certain number of iterations? 2. **Influence of cyclic pivot strategies**: What is the influence of different cyclic pivot strategies on the convergence of the Jacobi method? In particular, are there certain strategies that make the method unable to converge? ### Main contributions - The authors proved that for all 720 cyclic pivot strategies, the Jacobi method can converge globally. - Specifically, they proved that for any symmetric 4×4 matrix \(A\) and any cyclic pivot strategy, there exists a constant \(\gamma < 1\) such that the inequality \(S(A^{[t + 3]}) \leq \gamma S(A^{[t]})\) holds for all \(t\geq1\). Here \(S(A)\) represents the off - diagonal norm of matrix \(A\), which is defined as: \[ S(X)=\sqrt{\frac{2}{n(n - 1)}}\|X-\text{diag}(X)\|_F=\sqrt{\sum_{i = 1}^{n - 1}\sum_{j = i + 1}^n x_{ij}^2} \] where \(\|X\|_F\) is the Frobenius norm of matrix \(X\). ### Applications - This result can be directly applied to the J - Jacobi method, indicating that for 4×4 matrices, any cyclic J - Jacobi method can converge globally. - At the same time, this result also applies to block Jacobi methods, where each block is a 2×2 matrix, thereby further expanding the application range of the Jacobi method. ### Summary This paper solves the problem of the global convergence of the Jacobi method on symmetric 4×4 matrices through strict mathematical proofs and verifies the effectiveness of all 720 cyclic pivot strategies. This not only provides an important basis for theoretical research but also provides a reliable guarantee for numerical calculations in practical applications.