https://hgpu.org/?p=8682
A Parallel Preconditioned Bi-Conjugate Gradient Stabilized Solver for the Poisson Problem