Typo from paper: should be + not - in front of integral for RHS of edge system. Also fixed a few typos. Almost working now - still fails for higher order deformed, but appears to be working in most cube cases.