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.
Merge branch 'master' into higherorderspaces

Conflicts:
	source/helmholtz.py