https://hgpu.org/?p=5983
Implementing a Preconditioned Iterative Linear Solver Using Massively Parallel Graphics Processing Units