Modified p_norm_of_rel_error function to include epsilon term to prevent issues when v_new[i] == 0
Add acceleration using Eddington factors for MB3 for homogeneous case.

Updated edge flux acceleration method to be compatible with MB3.
Move variable status update line into the update_Qhat method.
Move update_Qhat to before acceleration step.
Only check for Qhat convergence if not accelerating.