Verified Error Bounds for Matrix Decompositions
Siegfried M. Rump,Takeshi Ogita
DOI: https://doi.org/10.1137/24m165096x
IF: 1.908
2024-11-13
SIAM Journal on Matrix Analysis and Applications
Abstract:SIAM Journal on Matrix Analysis and Applications, Volume 45, Issue 4, Page 2155-2183, December 2024. In this note we consider common matrix factorizations such as LU decomposition of a square and rectangular matrix, Cholesky and QR decomposition, singular value decomposition for square and rectangular matrices, and eigen-, Schur, and Takagi decomposition. We first note that well-conditioned factors tend to be sensitive to perturbations of the input matrix, while ill-conditioned factors tend to be insensitive. It seems that this behavior has not been recognized in numerical analysis. We develop a formula for the relation between the condition number of the factor and its sensitivity with respect to input perturbations, and give reasons for that. Our main focus is to describe verification methods for the factors of the mentioned decompositions. That means proving existence of the factorization together with rigorous entrywise error bounds for the factors. Our goal is to develop algorithms requiring [math] operations for an [math] matrix with [math] and [math]. Moreover, bounds of high quality are aimed for, often not far from maximal accuracy. A main tool to achieve that is accurate dot products based on error-free transformations. Since preconditioning based on approximate inverses is used, our methods are restricted to full matrices.
mathematics, applied