Applying Schur Complements for Handling General Up dates
of a Large, Sparse, Unsymmetric Matrix
Jacek Gondzio
Systems Research Institute, Polish Academy of Sciences,
Newelska 6, 01-447 Warsaw, Poland
e-mail: [email protected]
Technical Rep ort ZTSW-2-G244/93
September 1993, revised April 30, 1995
Abstract
We describ e a set of pro cedures for computing and up dating an inverse representation of
a large and sparse unsymmetric matrix A. The representation is built of two matrices: an
easily invertible, large and sparse matrix A and a dense Schur complement matrix S . An
0
ecient heuristic is given that nds this representation for any matrix A and keeps the size
of S as small as p ossible.
Equations with A are replaced with a sequence of equations that involve matrices A and
0
S . The former take full advantage of the sparsity of A ; the latter b ene t from applying
0
dense linear algebra techniques to a dense representation of the inverse of S .
We showhow to manage ve general up dates of A: row or column replacement, row and
column addition or deletion and a rank one correction. They all maintain an easily invertible
form of A and reduce to some low rank up date of matrix S .
0
An exp erimental implementation of the approach is describ ed and the preliminary com-
putational results of applying it to handling working basis up dates that arise in linear pro-
gramming are given.
Key words: large sparse unsymmetric matrix, inverse representation, Schur com-
plement, general up dates.
Supp orted by the Committe for Scienti c ResearchofPoland, grant No PB 8 S505 015 05. 1
1 Intro duction
In this pap er we describ e a metho d for the inverse representation of a large and sparse unsym-
metric matrix A. The representation is exp ected to b e used to solve equations
T
Ax = b and A y = d; 1
n
where x; b; y and d are vectors from R .
We present a metho d for up dating this representation after the following mo di cations of A:
{row replacement,
{ column replacement,
{row and column addition,
{row and column deletion,
{ rank one change.
LU decomp osition cf. Du et al. [5 ], Chapter 8 is probably the most widely used metho d
for representing an inverse of a sparse unsymmetric matrix A. Simplicity and eciency of its
up date after column replacement is one of its main advantages. The up dates by the metho d of
Bartels and Golub [1] preserve sparsity of LU factors and ensure accuracy of the representation
that is sucient for most practical applications. There exist many ecient implementations of
this metho d [4 , 29 , 32 ], including an extension of Gill et al. [13 ] for di erent up dates of a general
rectangular A.
The Bartels-Golub up dates of sparse LU factors have one basic disadvantage: they are
achieved by a purely sequential algorithm. Hence, relatively little can be done to sp ecialize
them for new computer architectures.
We feel motivated to lo ok for a paral lelisable metho d for the inverse representation of a
sparse unsymmetric matrix A. Additionally, in order to cover a wider eld of applications we
would like to b e able to up date this representation inexp ensively and in a stable way after simple
mo di cations of A.
The metho d presented in this pap er satis es these requirements.
We prop ose to represent the inverse of A with two matrices: an easily invertible fundamental
basis A and a Schur complement S hop efuly, of small size. The fundamental basis is an easily
0
invertible matrix close to A sub ject to the row and column p ermutations. The \closeness"
is measured with the number of columns of A that have to be replaced with appropriate unit
vectors in order to get a nonsingular matrix p ermutable to A . The numb er of columns that do
0
not t the required structure determines the size of the Schur complement. It is in our interest
to keep this numb er as small as p ossible.
Many di erent forms of A might be chosen. The most suitable one would clearly dep end
0
on the particular structure of A. The rst suggestion is to let A be triangular. Although it
0
seems very restrictive, this simple form already covers a large class of matrices, as e.g., linear
programming bases arising in network or sto chastic optimization [20 ]. The structure inherent
to these problems implies that all bases are \nearly triangular". 2
As shown in [17 ], the general linear programming bases can always be p ermuted to a form
close to a blo ck triangular matrix with at most 2x2 irreducible blo cks. This particular form of
A ensures that equations with it can b e solved very eciently since the sparsityof A as well
0 0
as that of the right hand side vectors can b e fully exploited.
It was also shown in the ab ovementioned pap er that the Schur complement inverse repre-