Switch to cube mesh - this show expected convergence (exact at p=2) with the new function