Multigrid Methods on GPUs: High Performance GPGPU Solvers for Partial Differential Equations