https://hgpu.org/?p=11809
Optimizing Krylov Subspace Solvers on Graphics Processing Units