Add project utility function
Add convergence tests of projection functionality

We check that an expression projected into CG1-3 and DG0-2 spaces
converges at the right order.