Stability of the Diagonal Pivoting Method with Partial Pivoting
- 1 January 1997
- journal article
- research article
- Published by Society for Industrial & Applied Mathematics (SIAM) in SIAM Journal on Matrix Analysis and Applications
- Vol. 18 (1) , 52-65
- https://doi.org/10.1137/s0895479895290371
Abstract
LAPACK and LINPACK both solve symmetric indefinite linear systems using the diagonal pivoting method with the partial pivoting strategy of Bunch and Kaufman [Math. Comp., 31 (1977), pp. 163--179]. No proof of the stability of this method has appeared in the literature. It is tempting to argue that the diagonal pivoting method is stable for a given pivoting strategy if the growth factor is small. We show that this argument is false in general and give a sufficient condition for stability. This condition is not satisfied by the partial pivoting strategy because the multipliers are unbounded. Nevertheless, using a more specific approach we are able to prove the stability of partial pivoting, thereby filling a gap in the body of theory supporting LAPACK and LINPACK.Keywords
This publication has 11 references indexed in Scilit:
- Accurate Symmetric Indefinite Linear Equation SolversSIAM Journal on Matrix Analysis and Applications, 1998
- Stability of block LU factorizationNumerical Linear Algebra with Applications, 1995
- The Factorization of Sparse Symmetric Indefinite MatricesIMA Journal of Numerical Analysis, 1991
- The Strong Stability of Algorithms for Solving Symmetric Linear SystemsSIAM Journal on Matrix Analysis and Applications, 1989
- Direct Solution of Sets of Linear Equations whose Matrix is Sparse, Symmetric and IndefiniteIMA Journal of Applied Mathematics, 1979
- Some stable methods for calculating inertia and solving symmetric linear systemsMathematics of Computation, 1977
- Decomposition of a symmetric matrixNumerische Mathematik, 1976
- Analysis of the Diagonal Pivoting MethodSIAM Journal on Numerical Analysis, 1971
- Direct Methods for Solving Symmetric Indefinite Systems of Linear EquationsSIAM Journal on Numerical Analysis, 1971
- Compatibility of approximate solution of linear equations with given error bounds for coefficients and right-hand sidesNumerische Mathematik, 1964