Skip to Main content Skip to Navigation
New interface
Reports (Research report)

Low rank approximation of a sparse matrix based on LU factorization with column and row tournament pivoting

Laura Grigori 1 Sébastien Cayrols 1 James W. Demmel 2 
1 ALPINES - Algorithms and parallel tools for integrated numerical simulations
LJLL - Laboratoire Jacques-Louis Lions, Inria Paris-Rocquencourt, INSMI - Institut National des Sciences Mathématiques et de leurs Interactions
Abstract : In this paper we present an algorithm for computing a low rank approximation of a sparse matrix based on a truncated LU factorization with column and row permutations. We present various approaches for determining the column and row permutations that show a trade-off between speed versus deterministic/probabilistic accuracy. We show that if the permutations are chosen by using tournament pivoting based on QR factorization, then the obtained truncated LU factorization with column/row tournament pivoting, LU\_CRTP, satisfies bounds on the singular values which have similarities with the ones obtained by a communication avoiding rank revealing QR factorization. Experiments on challenging matrices show that LU_CRTP provides a good low rank approximation of the input matrix and it is less expensive than the rank revealing QR factorization in terms of computational and memory usage costs, while also minimizing the communication cost. We also compare the computational complexity of our algorithm with randomized algorithms and show that for sparse matrices and high enough but still modest accuracies, our approach is faster.
Document type :
Reports (Research report)
Complete list of metadata

Cited literature [31 references]  Display  Hide  Download
Contributor : Laura Grigori Connect in order to contact the contributor
Submitted on : Tuesday, May 10, 2016 - 2:45:06 PM
Last modification on : Tuesday, November 29, 2022 - 12:12:15 PM


Files produced by the author(s)


  • HAL Id : hal-01313856, version 1


Laura Grigori, Sébastien Cayrols, James W. Demmel. Low rank approximation of a sparse matrix based on LU factorization with column and row tournament pivoting. [Research Report] RR-8910, inria. 2016, pp.35. ⟨hal-01313856⟩



Record views


Files downloads