I have a symmetric real matrix $A$. Iteratively I am doing the following:
$0) A_1=A$
$1)$ Calculate the largest eigenvalue $\lambda_i$ and corresponding unit norm eigenvector $v_i$ (you may assume a single largest eigenvalue).
$2)$ Then $A_i$ is constructed from $A_{i-1}$ by setting the diagonal element $A_{i}(j,j)=\lambda_i (v_i(j))^2$.
Repeat until convergence i.e. $\lambda_{i+1}=\lambda_{i}$ for some $i$. The simulations $i$ did shows that this always converge, sometimes slower sometimes faster. Do anyone know why this procedure is converging, if so to what it is converging?
Thanks