A fully portable high performance minimal storage hybrid format Cholesky algorithm

  • Authors:
  • Bjarne S. Andersen;John A. Gunnels;Fred G. Gustavson;John K. Reid;Jerzy Waśniewski

  • Affiliations:
  • Danish Meteorological Institute, Copenhagen, Denmark;IBM T.J. Watson Research Center, Yorktown Heights, NY;IBM T.J. Watson Research Center, Yorktown Heights, NY;Atlas Centre, Rutherford Appleton Laboratory, Didcot Oxon, UK;Technical University of Denmark, Lyngby, Denmark

  • Venue:
  • ACM Transactions on Mathematical Software (TOMS)
  • Year:
  • 2005

Quantified Score

Hi-index 0.00

Visualization

Abstract

We consider the efficient implementation of the Cholesky solution of symmetric positive-definite dense linear systems of equations using packed storage. We take the same starting point as that of LINPACK and LAPACK, with the upper (or lower) triangular part of the matrix stored by columns. Following LINPACK and LAPACK, we overwrite the given matrix by its Cholesky factor. We consider the use of a hybrid format in which blocks of the matrices are held contiguously and compare this to the present LAPACK code. Code based on this format has the storage advantages of the present code but substantially outperforms it. Furthermore, it compares favorably to using conventional full format (LAPACK) and using the recursive format of Andersen et al. [2001].