More work on Laplace3DSingleElement.  Some fixes for edge/vertex dof matching, but still not working.  Remains the case that we solve correctly for constant solutions, and don't for linears.  There is directionality in the bug: we get different L^2 error in phi if phi=x than if phi=y.