Simultaneous computation of the row and column rank profilesReport as inadecuate

Simultaneous computation of the row and column rank profiles - Download this document for free, or read online. Document in PDF available to download.

* Corresponding author 1 CASYS - Calculs Algébriques et Systèmes Dynamiques LJK - Laboratoire Jean Kuntzmann 2 MOAIS - PrograMming and scheduling design fOr Applications in Interactive Simulation Inria Grenoble - Rhône-Alpes, LIG - Laboratoire d-Informatique de Grenoble

Abstract : Gaussian elimination with full pivoting generates a PLUQ matrix decomposition. Depending on the strategy used in the search for pivots, the permutation matrices can reveal some information about the row or the column rank profiles of the matrix. We propose a new pivoting strategy that makes it possible to recover at the same time both row and column rank profiles of the input matrix and of any of its leading sub-matrices. We propose a rank-sensitive and quad-recursive algorithm that computes the latter PLUQ triangular decomposition of an m × n matrix of rank r in Omnr^{ω−2} field operations, with ω the exponent of matrix multiplication. Compared to the LEU decomposition by Malashonock, sharing a similar recursive structure, its time complexity is rank sensitive and has a lower leading constant. Over a word size finite field, this algorithm also improveLs the practical efficiency of previously known implementations.

Keywords : In place Rank profile PLUQ triangularization Rank sensitive In place.

Author: Jean-Guillaume Dumas - Clément Pernet - Ziad Sultan -



Related documents