Experiment with different block ordering for velocity solves. Here, currently only for cg serial, in numbering the dofs to solve for in petsc we number the components contiguously, instead of having all dofs for velocity component 1, followed by those of component 2, etc.
This is an experiment to see if we can improve multigrid solves on coupled velocity solves using bsor for smoothing and coarsening strategies based on blocks.

Currently fails:
Stokes_mms_cg_p1bp1.xml: F
Stokes_Subduction_VK_Case_1b.xml: FFFFFFFFF
channel_wind_drag_rotated.xml: F
The Stokes ones are probably expected, as we haven't improved our preconditioning yet. Not sure why that last one should fail (should be unaffected).