paper

On the growth factor upper bound for Aasen's algorithm

arXiv:1805.08994 · doi:10.1016/j.aml.2018.08.017

Abstract

Aasen's algorithm factorizes a symmetric indefinite matrix as , where is a permutation matrix, is unit lower triangular with its first column being the first column of the identity matrix, and is tridiagonal. In this note, we provide a growth factor upper bound for Aasen's algorithm which is much smaller than that given by Higham. We also show that the upper bound we have given is not tight when the matrix dimension is greater than or equal to .