Error: Unable to orthonormalize vector basis. Reason: Vector space has linear dependency. Where: This error was encountered inside Vector Space Basis.cpp

Because projection requires a linear solver and legacy DOLFIN defaults to a direct solver.

See Where to find 'project'-function in dolfinx? - #2 by nate or something like project(*args, solver_type="cg").