Summary
In this paper a Gauss-Jordan algorithm with column interchanges is presented and analysed. We show that, in contrast with Gaussian elimination, the Gauss-Jordan algorithm has essentially differing properties when using column interchanges instead of row interchanges for improving the numerical stability. For solutions obtained by Gauss-Jordan with column interchanges, a more satisfactory bound for the residual norm can be given. The analysis gives theoretical evidence that the algorithm yields numerical solutions as good as those obtained by Gaussian elimination and that, in most practical situations, the residuals are equally small. This is confirmed by numerical experiments. Moreover, timing experiments on a Cyber 205 vector computer show that the algorithm presented has good vectorisation properties.
Similar content being viewed by others
References
Businger, P.A.: Monitoring the numerical stability of Gaussian elimination. Numer. Math.16, 360–361 (1971)
Dongarra, J.J., Moler, C.F., Bunch, J.R., Stewart, G.W.: LINPACK User's guide. Philadelphia: SIAM 1979
Golub, G.H., Van Loan, C.F.: Matrix computations. Oxford: North Oxford Academic 1983
Hoffmann, W.: Solving linear systems on a vector computer. J. Comput. Appl. Math.18, 353–367 (1987)
Hoffmann, W., Lioen, W.M.: Chapter simultaneous linear equations. Report NM-R8614. In: NUMVEC FORTRAN Library Manual. Amsterdam: Center for Mathematics and Computer Science 1986
Peters, G., Wilkinson, J.H.: On the stability of Gauss-Jordan elimination with pivoting. Commun. ACM18, 20–24 (1975)
Stewart, G.W.: Introduction to matrix computations. New York London: Academic Press 1973
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Dekker, T.J., Hoffmann, W. Rehabilitation of the Gauss-Jordan algorithm. Numer. Math. 54, 591–599 (1989). https://doi.org/10.1007/BF01396364
Issue Date:
DOI: https://doi.org/10.1007/BF01396364