Finished writing constraints code. Created new testcase for BDFM1
Helmholtz equation, it runs and correctly enforces the continuity
constraints and elementwise constraints (i.e. normal components are linear
rather than quadratic). Next, a convergence test.