Stable double LR algorithm and its error analysis (Q1333873)

A normative matrix \(A\) is a real tridiagonal matrix with unit super- diagonal, sub-diagonal elements \(b_ i\), \(i=2,\dots,n\), and diagonal elements \(a_ i\), \(i=1,\dots,n\). The symmetric tridiagonal matrix \(S\) with the same diagonal and off-diagonal elements \(\beta_ i\), \(i=2,\dots,n\), has the same eigenvalues as \(A\) if \(b_ i= \beta^ 2_ i\), \(i=2,\dots,n\). The author shows that the double LR transformation of \(A\) with origin shift is related to the QR transformation of \(S\) with the same origin shift. The double LR has operation count only 4/7 of that for the QR method. A stable version of the double LR algorithm is then given, followed by the result that a practical version of this is numerically stable.