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

Laura Grigori 1 Sebastien 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.
Complete list of metadatas

Cited literature [31 references]  Display  Hide  Download

https://hal.inria.fr/hal-01313856
Contributor : Laura Grigori <>
Submitted on : Tuesday, May 10, 2016 - 2:45:06 PM
Last modification on : Friday, September 20, 2019 - 4:34:04 PM

File

paper.pdf
Files produced by the author(s)

Identifiers

  • HAL Id : hal-01313856, version 1

Citation

Laura Grigori, Sebastien 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⟩

Share

Metrics

Record views

675

Files downloads

797