LU factorization with panel rank revealing pivoting and its communication avoiding version - Inria - Institut national de recherche en sciences et technologies du numérique Access content directly
Reports (Research Report) Year : 2012

LU factorization with panel rank revealing pivoting and its communication avoiding version

Amal Khabou
  • Function : Author
  • PersonId : 918056
James W. Demmel
  • Function : Author
  • PersonId : 849003
Laura Grigori
Ming Gu
  • Function : Author
  • PersonId : 918057

Abstract

We present the LU decomposition with panel rank revealing pivoting (LU_PRRP), an LU factorization algorithm based on strong rank revealing QR panel factorization. LU_PRRP is more stable than Gaussian elimination with partial pivoting (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 and easily solves the Wilkinson matrix and the Foster matrix. We also present CALU_PRRP, a communication avoiding version of LU_PRRP that minimizes communication. 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. CALU_PRRP is more stable than CALU, the communication avoiding version of GEPP. CALU_PRRP is also more stable in practice and is resistant to pathological cases on which GEPP and CALU fail.
Fichier principal
Vignette du fichier
RR-7867.pdf (1 Mo) Télécharger le fichier
Origin : Files produced by the author(s)
Loading...

Dates and versions

hal-00723564 , version 1 (10-08-2012)

Identifiers

Cite

Amal Khabou, James W. Demmel, Laura Grigori, Ming Gu. LU factorization with panel rank revealing pivoting and its communication avoiding version. [Research Report] RR-7867, INRIA. 2012. ⟨hal-00723564⟩
209 View
204 Download

Altmetric

Share

Gmail Facebook X LinkedIn More