Performance Optimization for Sparse A(T)Ax in Parallel on Multicore Cpu

Yuan Tao,Yangdong Deng,Shuai Mu,Zhenzhong Zhang,Mingfa Zhu,Limin Xiao,Li Ruan
DOI: https://doi.org/10.1587/transinf.e97.d.315
2014-01-01
IEICE Transactions on Information and Systems
Abstract:The sparse matrix operation, y <- y+A(t)Ax, where A is a sparse matrix and x and y are dense vectors, is a widely used computing pattern in High Performance Computing (HPC) applications. The pattern poses challenge to efficient solutions because both a matrix and its transposed version are involved. An efficient sparse matrix format, Compressed Sparse Blocks (CSB), has been proposed to provide nearly the same performance for both Ax and A(t)x. We develop a multithreaded implementation for the CSB format and apply it to solve y <- y+A(t)Ax. Experiments show that our technique outperforms the Compressed Sparse Row (CSR) based solution in POSKI by up to 2.5 fold on over 70% of benchmarking matrices.
What problem does this paper attempt to address?