JavaScript is disabled for your browser. Some features of this site may not work without it.

## The Sparse Null Space Basis Problem

#####
**Author**

Coleman, Thomas F.; Pothen, Alex

#####
**Abstract**

The sparse null space basis problem is the following: $A t \times n$ matrix $A (t less than n)$ is given. Find a matrix $N$, with the fewest nonzero entries in it, whose columns span the null space of $A$. This problem arises in the design of practical algorithms for large-scale numerical optimization problems. Surprisingly, this problem can be formulated as a combinatorial optimization problem under a non-degeneracy assumption on $A$. THe theory of matchings in bipartite graphs - marriage theorems - can then be used to obtain the nonzero positions in $N$. Numerically stable matrix factorizations are used in the next phase to compute $N$. We use conformal decompositions to characterize the columns of a sparsest null basis. Matroid theory is used to prove that a greedy algorithm constructs a sparsest null basis. We prove that finding a sparsest null basis is NP-hard by showing that associated matroidal and graph-theoretic problems are NP-complete. We propose two approximation algorithms to construct sparse null bases. Both of them make use of the Dulmage-Mendelsohn decomposition of rectangular matrices. One algorithm is a sparsity exploiting variant of the variable-reduction technique. The second is a locally greedy algorithm that constructs a null basis with an upper triangular submatrix. These results are extended to computing sparse orthogonal null bases. We discuss how this "two-phase" approach can construct sparser null bases than a purely numerical approach; it is also potentially faster than the latter. Finally, we classify all known methods for constructing null bases, and show some unexpected equivalences between some of them.

#####
**Date Issued**

1984-07#####
**Publisher**

Cornell University

#####
**Subject**

computer science; technical report

#####
**Previously Published As**

http://techreports.library.cornell.edu:8081/Dienst/UI/1.0/Display/cul.cs/TR84-598

#####
**Type**

technical report