Commit makes solver use the new div and divT operators that were added
to the P2P1TaylorHoodStokesOperator. These replace the individual apply()
calls on div_x and div_y etc.