We can initialize based on a non-closed vector
Compute adjoint dirichlet constraints

And a bugfix for a case when some qois have constraints and others
don't
Remove overzealous asserts, print all constraints