Check for new way of computing Dirichlet boundary conditions (replaced VectorTools::project_boundary_values_curl_conforming with my own version, which is included at the top in MyVectorTools). So far it works for a cube on higher degrees, but the convergence is off when the face is non-rectangular, although better than the library function