A pivoting strategy for symmetric tridiagonal matrices

The LBLT factorization of Bunch for solving linear systems involving a symmetric indefinite tridiagonal matrix T is a stable, efficient method. It computes a unit lower triangular matrix L and a block 1 x 1 and 2 x 2 matrix B such that T = LBLT. Choosing the pivot size requires knowing a priori the largest element a of T in magnitude. In some applications, it is required to factor T as it is formed without necessarily knowing sigma. In this paper, we present a modification of the Bunch algorithm that can satisfy this requirement. We demonstrate that this modification exhibits the same bound on the growth factor as the Bunch algorithm and is likewise normwise backward stable. Copyright (c) 2005 John Wiley & Sons, Ltd.