Systems Optimization Laboratory
Stanford, CA 94305-4121 USA
LUMOD: Updating a dense square factorization L*C = U
M. A. Saunders.
- CONTENTS: Fortran software for updating a dense square factorization
L*C = U when rows and columns of C are added, deleted or
(Suitable as basis package for dense simplex method,
or for updating sparse factorizations via the
L is square, stored by rows in a 1-D array.
It is a product of stabilized elementary transformations.
U is upper triangular, stored by rows in a 1-D array.
The dimension of C, L and U may change. If maxn
is the largest C allowed for, the total storage is
maxn^2 for L and maxn(maxn+1)/2 for U.
- REFERENCES: Stabilized elementary transformations are described in
J. H. Wilkinson, The Algebraic Eigenvalue Problem,
Clarendon Press, Oxford (1965).
- RELEASE: f77 files are well tested.