A Subspace Error Estimate for Linear Systems

Abstract

This paper proposes a new method for estimating the error in the solution of linear systems. A condition number is defined for a linear function of the solution components. This definition of the condition number is quite versatile. It reduces to the component condition number proposed by Chandrasekaran and Ipsen [SIAM J. Matrix Anal. Appl., 16 (1995), pp. 93--112] and to Skeel's definition of condition number [J. ACM, 26 (1979), pp. 494--526] in some special cases, and it can be used to estimate the error in a subspace. The estimate is based on the adjoint equation in combination with small sample statistical theory. It can be implemented simply and is inexpensive to compute. Numerical examples are presented which illustrate the power and effectiveness of this error estimate.

Keywords

This publication has 13 references indexed in Scilit:

Accuracy and Stability of Numerical Algorithms
Published by Society for Industrial & Applied Mathematics (SIAM) ,2002
On the Complexity of Computing Error Bounds
Foundations of Computational Mathematics, 2001
Statistical Condition Estimation for Linear Systems
SIAM Journal on Scientific Computing, 1998
On the Sensitivity of Solution Components in Linear Systems of Equations
SIAM Journal on Matrix Analysis and Applications, 1995
Probabilistic Bounds on the Extremal Eigenvalues and Condition Number by the Lanczos Algorithm
SIAM Journal on Matrix Analysis and Applications, 1994
Small-Sample Statistical Condition Estimates for General Matrix Functions
SIAM Journal on Scientific Computing, 1994
Estimating the Largest Eigenvalue by the Power and Lanczos Algorithms with a Random Start
SIAM Journal on Matrix Analysis and Applications, 1992
Estimating Extremal Eigenvalues and Condition Numbers of Matrices
SIAM Journal on Numerical Analysis, 1983
An Estimate for the Condition Number of a Matrix
SIAM Journal on Numerical Analysis, 1979
Numerical Linear Algebra
Canadian Mathematical Bulletin, 1966