https://hgpu.org/?p=6505
Efficient Two-Level Preconditionined Conjugate Gradient Method on the GPU