Pivotal interchanges are commonly used in the solution of large and sparse systems of linear algebraic equations by Gaussian elimination (in order to preserve the sparsity of the matrix and to prevent the appearance of large roundoff errors during the computations). The Markowitz strategy (see [H. M. Markowitz, The elimination form of inverse and its applications to linear programming, Management Sci., 3 (1957), pp. 255–269]) is often used to determine the pivotal sequence. An efficient implementation of this strategy is given by Curtis and Reid (see [A. R. Curtis and J. K. Reid, Fortran subroutines for the solution of sparse sets of linear equations, A.E.R.E., Report R.6844, HMSO, London, 1971]) and improved by Duff (see [I. S. Duff, MA28—a set of Fortran subroutines for sparse unsymmetric matrices, A.E.R.E., Report R.8730, HMSO, London, 1977]). In this paper it is shown how the classical Markowitz idea can be generalized. Consider the following parameters: u —the stability factor and $p(s)$—the number o...
Read full abstract