Fast Computation of the Rank Profile Matrix and the Generalized Bruhat Decomposition

by   Jean-Guillaume Dumas, et al.

The row (resp. column) rank profile of a matrix describes the stair-case shape of its row (resp. column) echelon form.We here propose a new matrix invariant, the rank profile matrix, summarizing all information on the row andcolumn rank profiles of all the leading sub-matrices.We show that this normal form exists and is unique over any ring, provided that the notion of McCoy's rank is used, in the presence of zero divisors.We then explore the conditions for a Gaussian elimination algorithm to compute all or part of this invariant,through the corresponding PLUQ decomposition. This enlarges the set of known Elimination variants that compute row or column rank profiles.As a consequence a new Crout base case variant significantly improves the practical efficiency of previously known implementations over a finite field.With matrices of very small rank, we also generalize the techniques ofStorjohann and Yang to the computation of the rank profile matrix, achieving an (r^ω+mn)^1+o(1) time complexity for an m × n matrix of rank r, where ω is the exponent of matrix multiplication. Finally, by give connections to the Bruhat decomposition, and severalof its variants and generalizations. Thus, our algorithmicimprovements for the PLUQ factorization, and their implementations,directly apply to these decompositions.In particular, we show how a PLUQ decomposition revealing the rankprofile matrix also reveals both a row and a column echelon form ofthe input matrix or of any of its leading sub-matrices, by a simplepost-processing made of row and column permutations.


page 1

page 2

page 3

page 4


Rank-Sensitive Computation of the Rank Profile of a Polynomial Matrix

Consider a matrix 𝐅∈𝕂[x]^m × n of univariate polynomials over a field 𝕂....

Elimination-based certificates for triangular equivalence and rank profiles

In this paper, we give novel certificates for triangular equivalence and...

Certificates for triangular equivalence and rank profiles

In this paper, we give novel certificates for triangular equivalence and...

Symmetric indefinite triangular factorization revealing the rank profile matrix

We present a novel recursive algorithm for reducing a symmetric matrix t...

Numerical considerations and a new implementation for ICS

Invariant Coordinate Selection (ICS) is a multivariate data transformati...

Matrix factoring by fraction-free reduction

We consider exact matrix decomposition by Gauss-Bareiss reduction. We in...

Perturbations of CUR Decompositions

The CUR decomposition is a factorization of a low-rank matrix obtained b...

Please sign up or login with your details

Forgot password? Click here to reset