... and Schulz (2011) used shared memory of GPU to store the nodal data of three successive 2D slices of nodes for each thread block. By successively moving in the third direction, SpMV for the entire mesh can be efficiently performed. Using an NbN − based matrix-free PCG solver, the authors achieved...