Abstract
In this paper, we consider how to factor symmetric totally nonpositive matrices and their inverses by taking advantage of the symmetric property. It is well-known that the Bunch-Kaufman algorithm is the most commonly used pivoting strategy which can, however, produce arbitrarily large entries in the lower triangular factor for such matrices as illustrated by our example. Therefore, it is interesting to show that when the Bunch-Parlett algorithm is simplified for these matrices, it only requires O(n 2) comparisons with the growth factor being nicely bounded by 4. These facts, together with a nicely bounded lower triangular factor and a pleasantly small relative backward error, show that the Bunch-Parlett algorithm is more preferable than the Bunch-Kaufman algorithm when dealing with these matrices.
Similar content being viewed by others
Notes
In fact, \(a^{(2)}_{ij}\) should be \(a^{(2)}_{i-1,j-1}\) or \(a^{(2)}_{ij}\). For the sake of the convenience, we perform a minor abuse of notation throughout the paper.
References
Alonso, P., Gasca, M., Peña, J.M.: Backward error analysis of Neville elimination. Appl. Numer. Math. 23, 193–204 (1997)
Ando, T.: Totally positive matrices. Linear Algebra Appl. 90, 165–219 (1987)
Ashcraft, C., Grimes, R.G., Lewis, J.G.: Accurate symmetric indefinite linear equation solvers. SIAM J. Matrix Anal. Appl. 20, 513–561 (1998)
Bunch, J.R.: Analysis of the diagonal pivot method. SIAM J. Numer. Anal. 8, 656–680 (1971)
Cortés, V., Peña, J.M.: A stable test for strict sign regularity. Math. Comput. 77, 2155–2171 (2008)
de Boor, C., Pinkus, A.: Backward error analysis for totally positive linear systems. Numer. Math. 27, 485–490 (1977)
Demmel, J., Koev, P.: The accurate and efficient solution of a totally positive generalized Vandermonde linear system. SIAM J. Matrix Anal. Appl. 27, 142–152 (2005)
Dopico, F.M., Koev, P.: Accurate symmetric rank revealing and eigendecompositions of symmetric structured matrices. SIAM J. Matrix Anal. Appl. 28, 1126–1156 (2006)
Fallat, S.M., van den Driessche, P.: On matrices with all minors negative. Electron. J. Linear Algebra 7, 92–99 (2000)
Fang, H.R., O’Leary, D.P.: Stable factorizations of symmetric tridiagonal and triadic matrices. SIAM J. Matrix Anal. Appl. 28, 576–595 (2006)
George, A., Ikramov, K.D.: Gaussian elimination is stable for the inverse of a diagonally dominant matrix. Math. Comput. 73, 653–657 (2003)
Higham, N.J.: Bounding the error in Gaussian elimination for tridiagonal systems. SIAM J. Matrix Anal. Appl. 11, 521–530 (1990)
Higham, N.J.: Stability of the diagonal pivoting method with partial pivoting. SIAM J. Matrix Anal. Appl. 18, 52–65 (1997)
Higham, N.J.: Stability of block LDL T factorization of a symmetric tridiagonal matrix. Linear Algebra Appl. 287, 181–189 (1999)
Higham, N.J.: Accuracy and Stability of Numerical Algorithms, 2nd edn. SIAM, Philadelphia (2002)
Huang, R.: Sign structure preserves for m-banded factorizations of sign regular matrices. Linear Algebra Appl. 436, 1990–2000 (2012)
Huang, R., Liu, J.Z.: On Schur complements of sign regular matrices of order k. Linear Algebra Appl. 433, 143–148 (2010)
Karlin, S.: Total Positivity. Stanford University Press, Stanford (1968)
Peláez, M.J., Moro, J.: Accurate factorization and eigenvalue algorithms for symmetric DSTU and TSC matrices. SIAM J. Matrix Anal. Appl. 28, 1173–1198 (2006)
Peña, J.M.: Backward stability of a pivoting strategy for sign regular linear systems. BIT Numer. Math. 37, 910–924 (1997)
Peña, J.M.: Shape Preserving Representations in Computer-Aided Geometric Design. Nova Science, Commack (1999)
Peña, J.M.: On nonsingular sign regular matrices. Linear Algebra Appl. 359, 91–100 (2003)
Slapničar, I.: Componentwise analysis of direct factorization of real symmetric and Hermitian matrices. Linear Algebra Appl. 272, 227–275 (1998)
Acknowledgements
The author is very grateful to the referees and Prof. D. Kressner for their very careful reading of the paper and valuable suggestions.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Daniel Kressner.
This work was supported by the National Natural Science Foundation for Youths of China (Grant No. 11001233), the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20094301120002), the Hunan Provincial Natural Science Foundation of China (Grant No. 11JJ4004).
Rights and permissions
About this article
Cite this article
Huang, R. Factoring symmetric totally nonpositive matrices and inverses with a diagonal pivoting method. Bit Numer Math 53, 443–458 (2013). https://doi.org/10.1007/s10543-012-0404-2
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10543-012-0404-2