Fixed bug in residual norm calculation

The residual r_\phi was calculated incorrectly from the quantity M_\phi r_\phi, need to do mass matrix solve with RHS M_\phi r_\phi instead.