Add `permute_closure_dofs` function

[mscroggs]

Signature: permute_closure_dofs(span<int>, u8)

This should apply a permutation to all dofs including swapping positions of those on vertices.