8481 articles  [english version]

hal-00723564, version 1

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

Amal Khabou () 1, James W. Demmel () 2, Laura Grigori () 3, Ming Gu () 2

N° RR-7867 (2012)

  • 1 :  Laboratoire de Recherche en Informatique (LRI)
  • http://www.lri.fr/
    CNRS : UMR8623 – Université Paris XI - Paris Sud LRI - Bâtiment 490 Université Paris-Sud 91405 Orsay Cedex France
  • 2 :  Berkeley University of California (UC BERKELEY)
  • http://berkeley.edu/
    UC Berkeley University of California, Berkeley Department of Mathematics Berkeley, CA 94720, USA États-Unis
  • 3 :  GRAND-LARGE (INRIA Saclay - Ile de France)
  • http://grand-large.lri.fr/index.php/Main_Page
    INRIA – CNRS : UMR8623 – Université Paris XI - Paris Sud PCRI - Université de Paris-Sud, Bâtiment 650, 91405 ORSAY Cedex France

Références bibliographiques

  • Type de publication : Rapports
  • Domaine : Informatique/Analyse numérique
  • Titre : LU factorization with panel rank revealing pivoting and its communication avoiding version
  • Résumé : 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.
  • Langue du document : Anglais
  • Type de rapport : Rapport de recherche
  • Date de publication : 25/01/2012
  • Référence interne : RR-7867

Liste des fichiers attachés à ce document :

 
  • hal-00723564, version 1
  • oai:hal.inria.fr:hal-00723564
  • Contributeur : 
  • Soumis le : Vendredi 10 Août 2012, 18:08:51
  • Dernière modification le : Dimanche 12 Août 2012, 21:19:36