A guide for implementing tridiagonal solvers on GPUs

Abstract

Tridiagonal solvers are important building blocks for many applications on GPUs. Although a wide range of algorithms and optimizations have been proposed for tridiagonal solvers, there are no comprehensive guidelines for building a high-performance tridiagonal solver for GPUs. In this chapter, we review and summarize several existing optimization strategies, and systematically analyze the purpose for each optimization. Finally, a case study, called SPIKE-CR, is given to demonstrate how to apply the guidelines to build a high-performance GPU tridiagonal solver.