Adaptive precision in block‐Jacobi preconditioning for iterative sparse linear system solvers
Citations Over TimeTop 10% of 2018 papers
Abstract
Summary We propose an adaptive scheme to reduce communication overhead caused by data movement by selectively storing the diagonal blocks of a block‐Jacobi preconditioner in different precision formats (half, single, or double). This specialized preconditioner can then be combined with any Krylov subspace method for the solution of sparse linear systems to perform all arithmetic in double precision. We assess the effects of the adaptive precision preconditioner on the iteration count and data transfer cost of a preconditioned conjugate gradient solver. A preconditioned conjugate gradient method is, in general, a memory bandwidth‐bound algorithm, and therefore its execution time and energy consumption are largely dominated by the costs of accessing the problem's data in memory. Given this observation, we propose a model that quantifies the time and energy savings of our approach based on the assumption that these two costs depend linearly on the bit length of a floating point number. Furthermore, we use a number of test problems from the SuiteSparse matrix collection to estimate the potential benefits of the adaptive block‐Jacobi preconditioning scheme.
Related Papers
- → Communication Avoiding ILU0 Preconditioner(2015)30 cited
- → Fast and improved scaled HSS preconditioner for steady-state space-fractional diffusion equations(2020)2 cited
- → Application of the preconditioned conjugate-gradient algorithm to the edge FEM for electromagnetic boundary-value problems(2000)23 cited
- → The Construction of Parallel Preconditioner in Preconditioned Conjugate Gradient Method(2010)
- A New Adaptive Preconditioned Global CGS Algorithm(2011)