Trying grad(div) with gausseidel