Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Ordering sparse matrices for cache-based systems

Conference ·
OSTI ID:787125
The Conjugate Gradient (CG) algorithm is the oldest and best-known Krylov subspace method used to solve sparse linear systems. Most of the coating-point operations within each CG iteration is spent performing sparse matrix-vector multiplication (SPMV). We examine how various ordering and partitioning strategies affect the performance of CG and SPMV when different programming paradigms are used on current commercial cache-based computers. However, a multithreaded implementation on the cacheless Cray MTA demonstrates high efficiency and scalability without any special ordering or partitioning.
Research Organization:
Lawrence Berkeley National Lab., CA (US)
Sponsoring Organization:
USDOE Director, Office of Science. Office of Advanced Scientific Computing Research. Mathematical, Information, and Computational Sciences Division (US)
DOE Contract Number:
AC03-76SF00098
OSTI ID:
787125
Report Number(s):
LBNL--47805
Country of Publication:
United States
Language:
English