Khabou, Amal and Demmel, James W. and Grigori, Laura and Gu, Ming (2013) LU FACTORIZATION WITH PANEL RANK REVEALING PIVOTING AND ITS COMMUNICATION AVOIDING VERSION. [MIMS Preprint]
![[thumbnail of luprrp_caluprrp.pdf]](https://eprints.maths.manchester.ac.uk/style/images/fileicons/application_pdf.png) |
PDF
luprrp_caluprrp.pdf Download (513kB) |
Abstract
We present the block LU factorization with panel rank revealing pivoting (block LU_PRRP), a decomposition algorithm based on strong rank revealing QR panel factorization. Block LU_PRRP is more stable than Gaussian elimination with partial pivoting (GEPP), with a theoretical upper bound of the growth factor of $(1+ \tau b)^{(n/ b)-1}$, where $b$ is the size of the panel used during the block factorization, $\tau$ is a parameter of the strong rank revealing QR factorization, and $n$ is the number of columns of the matrix. For example, if the size of the panel is $b = 64$, and $\tau = 2$, then $(1+2b)^{(n/b)-1} = (1.079)^{n-64} \ll 2^{n-1}$, where $2^{n-1}$ is the upper bound of the growth factor of GEPP. Our extensive numerical experiments show that the new factorization scheme is as numerically stable as GEPP in practice, but it is more resistant to pathological cases. The block LU_PRRP factorization does only $O(n^2 b)$ additional floating point operations compared to GEPP. We also present block CALU_PRRP, a communication avoiding version of block LU_PRRP that minimizes communication. Block CALU_PRRP is based on tournament pivoting, with the selection of the pivots at each step of the tournament being performed via strong rank revealing QR factorization. Block CALU_PRRP is more stable than CALU, the communication avoiding version of GEPP, with a theoretical upper bound of the growth factor of $(1+ \tau b)^{{n\over b}(H+1)-1}$, where $H$ is the height of the reduction tree used during tournament pivoting. The upper bound of the growth factor of CALU is $2^{n(H+1)-1}$. Block CALU_PRRP is also more stable in practice and is resistant to pathological cases on which GEPP and CALU fail.
| Item Type: | MIMS Preprint |
|---|---|
| Additional Information: | submitted to SIAM J. Matrix Anal. & Appl., in revision |
| Uncontrolled Keywords: | LU factorization, growth factor, numerical stability, communication avoiding, strong rank revealing QR factorization |
| Subjects: | MSC 2010, the AMS's Mathematics Subject Classification > 65 Numerical analysis |
| Depositing User: | Amal Khabou |
| Date Deposited: | 13 Mar 2013 |
| Last Modified: | 08 Nov 2017 18:18 |
| URI: | https://eprints.maths.manchester.ac.uk/id/eprint/1958 |
Available Versions of this Item
- LU FACTORIZATION WITH PANEL RANK REVEALING PIVOTING AND ITS COMMUNICATION AVOIDING VERSION. (deposited 13 Mar 2013) [Currently Displayed]
Actions (login required)
 |
View Item |