add c_grad_p to forcing function
use reformulated forcing function, remove c_grad_p from pressure gradient force
revised velocity prediction, velocity correction and forcing function
c_grad_p fixes, do not specify device, add c=0.1 neumann test
do not run tests with c=0.1