When I compare the L2 error produced with my Julia code versus mach's on the steady isentropic vortex, there is a small difference. I have isolated the issue to the metric/mapping Jacobian terms present in the LPS terms. The difference may be due to how mfem computes the mapping Jacobian, so I am not that worried about it...for now. However, I wanted to document the issue somewhere.