Moved the reduction functions to separate script to allow re-using in diffusion and pressure
Calculate max divergence at GPU
Minor cleanup, plus added function for CUDA error checking
Don't do unnecessary divisions at GPU. Can easily speed up a kernel with ~25%