Do not set the quadrature order to zero
Do not set the quadrature order to zero at the beginning of the computeH1HalfNormDifferenceSquared method. Instead, use the order requested by the caller.
I don't know how the line that sets the order to zero ever got into the code. It makes no sense at all.