|Efficiency of General Krylov Methods on GPUs – An Experimental Study
|Year of Publication
|Anzt, H., J. Dongarra, M. Kreutzer, G. Wellein, and M. Kohler
|The Sixth International Workshop on Accelerators and Hybrid Exascale Systems (AsHES)
|algorithmic bombardment, BiCGSTAB, CGS, gpu, IDR(s), Krylov solver, QMR
|This paper compares different Krylov methods based on short recurrences with respect to their efficiency when implemented on GPUs. The comparison includes BiCGSTAB, CGS, QMR, and IDR using different shadow space dimensions. These methods are known for their good convergence characteristics. For a large set of test matrices taken from the University of Florida Matrix Collection, we evaluate the methods’ performance against different target metrics: convergence, number of sparse matrix-vector multiplications, and execution time. We also analyze whether the methods are “orthogonal” in terms of problem suitability. We propose best practices for choosing methods in a “black box” scenario, where no information about the optimal solver is available.