A Relaxed Algorithm for Online Matrix Inversion

Arne Storjohann,Shiyun Yang
DOI: https://doi.org/10.1145/2755996.2756672
2015-06-24
Abstract:We consider a variation of the well known problem of computing the unique solution to a nonsingular system Ax=b of n linear equations over a field K. The variation assumes that A has generic rank profile and requires as output not only the single solution vector A-1b Ε Kn x 1, but rather the solution to all leading principle subsystems. Most importantly, the rows of the augmented system [A||b] are given one at a time from first to last, and as soon as the next row is given the solution to the next leading principal subsystem should be produced. We call this problem ONLINESYSTEM. The obvious iterative algorithm for ONLINESYSTEM has a cost in terms of field operations that is cubic in the dimension of A. In this paper we introduce a relaxed representation for the inverse and show how to obtain an algorithm for ONLINESYSTEM that allows us to incorporate matrix multiplication. As an application we show how to introduce fast matrix multiplication into the inherently iterative algorithm for row rank profile computation presented previously by the authors.
What problem does this paper attempt to address?