Rank-profile revealing Gaussian elimination and the CUP matrix decomposition

Claude-Pierre Jeannerod 1, 2, * Clément Pernet 3, * Arne Storjohann 4, *
* Auteur correspondant
1 ARENAIRE - Computer arithmetic
Inria Grenoble - Rhône-Alpes, LIP - Laboratoire de l'Informatique du Parallélisme
3 MOAIS - PrograMming and scheduling design fOr Applications in Interactive Simulation
Inria Grenoble - Rhône-Alpes, LIG - Laboratoire d'Informatique de Grenoble
Abstract : Transforming a matrix over a field to echelon form, or decomposing the matrix as a product of structured matrices that reveal the rank profile, is a fundamental building block of computational exact linear algebra. This paper surveys the well known variations of such decompositions and transformations that have been proposed in the literature. We present an algorithm to compute the CUP decomposition of a matrix, adapted from the LSP algorithm of Ibarra, Moran and Hui (1982), and show reductions from the other most common Gaussian elimination based matrix transformations and decompositions to the CUP decomposition. We discuss the advantages of the CUP algorithm over other existing algorithms by studying time and space complexities: the asymptotic time complexity is rank sensitive, and comparing the constants of the leading terms, the algorithms for computing matrix invariants based on the CUP decomposition are always at least as good except in one case. We also show that the CUP algorithm, as well as the computation of other invariants such as transformation to reduced column echelon form using the CUP algorithm, all work in place, allowing for example to compute the inverse of a matrix on the same storage as the input matrix.
Type de document :
[Research Report] 2011, pp.32
Liste complète des métadonnées

Contributeur : Claude-Pierre Jeannerod <>
Soumis le : vendredi 30 décembre 2011 - 14:28:24
Dernière modification le : jeudi 11 octobre 2018 - 08:48:03

Lien texte intégral


  • HAL Id : hal-00655543, version 1
  • ARXIV : 1112.5717


Claude-Pierre Jeannerod, Clément Pernet, Arne Storjohann. Rank-profile revealing Gaussian elimination and the CUP matrix decomposition. [Research Report] 2011, pp.32. 〈hal-00655543〉



Consultations de la notice