Now the gauss seidel solver returns the errors and we plot them separately for pressure projection and implicit diffusion