Department of Applied Mathematics, University of Twente
Doi
Abstract
When applying Gaussian elimination to a sparse matrix, it is desirable to avoid turning zeros into non-zeros to preserve the sparsity. The class of perfect elimination bipartite graphs is closely related to square matrices that Gaussian elimination can be applied to without turning any zero into a non-zero. Existing literature on the recognition of this class and finding suitable pivots mainly focusses on time complexity. For n×n matrices with m non-zero elements, the currently best known algorithm has a time complexity of O(n3/logn). However, when viewed from a practical perspective, the space complexity also deserves attention: it may not be worthwhile to look for a suitable set of pivots for a sparse matrix if this requires Ω(n2) space. We present two new algorithms for the recognition of sparse instances: one with a O(nm) time complexity in Θ(n2) space and one with a O(m2) time complexity in Θ(m) space. Furthermore, if we allow only pivots on the diagonal, our second algorithm can easily be adapted to run in time O(nm)