Sparse matrix solvers on the GPU: conjugate gradients and multigrid
Vector Models for Data-Parallel Computing