Abstract
A new pivotal strategy for reducing storage and arithmetic operations in computing the inverse of a matrix is outlined. This algorithm uses recursive deletion and partition to generate an efficient structure. The optimal selection of rows and columns to be deleted (spikes) is greatly facilitated through the use of certain exclusion criteria. The algorithm produces significantly fewer fill-ins and requires significantly fewer operations than the alternative P3 and P4 algorithms for the product form of the inverse (PFI). It also produces satisfactory structures for problems for which P3 and P4 would generate zero pivots. Three variants of HP algorithm specifically tailored for the elimination form of the inverse (EFI) and for greater speed are also presented. The algorithm has been implemented on a CDC 6400 computer. Performance comparison on six sample problems is given.
Similar content being viewed by others
References
J.S. Duff and J.K. Reid, “A comparison of sparsity orderings for obtaining a pivotal sequence in Gaussian elimination”,Journal of the Institute of Mathematics and its Applications 14 (1974) 281–291.
E. Hellerman and D. Rarick, “Reinversion with the preassigned pivot sequence”,Mathematical Programming 1 (1971) 195–216.
E. Hellerman and D. Rarick, “The partitioned reassigned pivot procedure (p4)”, in: D.J. Rose and R.A. Willoughby, eds.,Sparse matrices and their applications (Plenum Press, New York, 1972).
H.Y. Hsieh, “Fill-in comparisons between Gauss—Jordan and Gaussian eliminations”,IEEE Transactions on Circuits and Systems CAS-21 (2) (1974) 230–233.
R.S.H. Mah, “Recent developments in process design”, in: W.R. Spillers, ed.,Basic questions of design theory (North-Holland, Amsterdam, 1974).
H.M. Markowitz, “The elimination form of the inverse and its application to linear programming”,Management Science 3 (3) (1957) 225–269.
D.J. Rose and J.R. Bunch, “The role of partitioning in the numerical solution of sparse systems”, in: D.J. Rose and R.A. Willoughby, eds.,Sparse matrices and their applications (Plenum Press, New York 1972).
R.W.H. Sargent and A.W. Westerberg, “SPEED-UP in chemical engineering design”,Transactions of the Institution of Chemical Engineers, 42 (1964) 190–197.
D.V. Steward, “On an approach to techniques for the analysis of the structures of large systems of equations”,SIAM Review 4 (4) (1962) 321–342.
D.V. Steward, “Partitioning and tearing systems of equations”,SIAM Journal on Numerical Analysis 2 (1965) 345–365.
R. Tarjan, “Enumeration of the elementary circuits of a directed graph”,SIAM Journal on Computing 2 (3) (1973) 211–216.
H. Weinblatt, “A new search algorithm for finding the simple cycles of a finite directed graph”,Journal of the Association for Computing Machinery 19 (1972) 43–56.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Lin, T.D., Mah, R.S.H. Hierarchical partition—a new optimal pivoting algorithm. Mathematical Programming 12, 260–278 (1977). https://doi.org/10.1007/BF01593792
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF01593792