BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Ricardo S. Leite contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We provide a detailed description of the maps associated with spectral interlacing, for rank one perturbations and bordering of symmetric and Hermitian matrices. The arguments rely on standard techniques of nonlinear analysis.
The computation of eigenvalues of real symmetric tridiagonal matrices frequently proceeds by a sequence of QR steps with shifts. We introduce simple shift strategies, functions sigma satisfying natural conditions, taking each n x n matrix T to a real number sigma(T). The strategy specifies the shift to be applied by the QR step at T. Rayleigh and Wilkinson's are examples of simple shift strategies. We show that if sigma is continuous then there exist initial conditions for which deflation does not occur, i.e., subdiagonal entries do not tend to zero. In case of deflation, we consider the rate of convergence to zero of the (n, n-1) entry: for simple shift strategies this is always at least quadratic. If the function sigma is smooth in a suitable region and the spectrum of T does not include three consecutive eigenvalues in arithmetic progression then convergence is cubic. This implies cubic convergence to deflation of Wilkinson's shift for generic spectra. The study of the algorithm near deflation uses tubular coordinates, under which QR steps with shifts are given by a simple formula.
One of the most widely used methods for eigenvalue computation is the $QR$ iteration with Wilkinson's shift: here the shift $s$ is the eigenvalue of the bottom $2\times 2$ principal minor closest to the corner entry. It has been a long-standing conjecture that the rate of convergence of the algorithm is cubic. In contrast, we show that there exist matrices for which the rate of convergence is strictly quadratic. More precisely, let $T_X$ be the $3 \times 3$ matrix having only two nonzero entries $(T_X)_{12} = (T_X)_{21} = 1$ and let $T_L$ be the set of real, symmetric tridiagonal matrices with the same spectrum as $T_X$. There exists a neighborhood $U \subset T_L$ of $T_X$ which is invariant under Wilkinson's shift strategy with the following properties. For $T_0 \in U$, the sequence of iterates $(T_k)$ exhibits either strictly quadratic or strictly cubic convergence to zero of the entry $(T_k)_{23}$. In fact, quadratic convergence occurs exactly when $\lim T_k = T_X$. Let $X$ be the union of such quadratically convergent sequences $(T_k)$: the set $X$ has Hausdorff dimension 1 and is a union of disjoint arcs $X^σ$ meeting at $T_X$, where $σ$ ranges over a Cantor set.