OpSparse: A Highly Optimized Framework for Sparse General Matrix Multiplication on GPUs

Abstract

Sparse general matrix multiplication (SpGEMM) is an important and expensive computation primitive in many real-world applications. Due to SpGEMM&#x2019;s inherent irregularity and the vast diversity of its input matrices, developing high-performance SpGEMM implementation on modern processors such as GPUs is challenging. The state-of-the-art SpGEMM libraries (i.e., <i>nsparse</i> and <i>spECK</i>) adopt several algorithms to tackle the challenges of global load balance, local load balance, and allocation of the result matrix. While these libraries focus on the high-level algorithm design for SpGEMM, they neglect several low-level architecture-specific optimizations, which causes inefficient implementations in their libraries. In this paper, we classify their inefficient implementations into several categories. Based on our observations, we propose a highly optimized SpGEMM library called <i>OpSparse</i>. The optimizations in <i>OpSparse</i> include 1) optimizing the binning method by improving the utilization of the shared memory, 2) optimizing the hashing method by reducing the access to the hash table, 3) improving the trade-off between hash collision rate and hardware utilization in the hashing method by setting appropriate binning ranges, 4) reducing the overheads of global memory utilization by minimizing the global memory usage of the metadata, and 5) improving the execution parallelism by overlapping global memory allocation with kernel execution. Performance evaluations with 26 commonly used matrices on an Nvidia Tesla V100 GPU show that <i>OpSparse</i> achieves on average 7.35&#x00D7; (up to 27.8&#x00D7;), 1.43&#x00D7; (up to 1.81&#x00D7;), and 1.52&#x00D7; (up to 2.04&#x00D7;) speedups over three state-of-the-art SpGEMM libraries: <i>cuSPARSE</i>, <i>nsparse</i>, and <i>spECK</i>, respectively.