Draft: Implement subIndex for 3D edges and corners
Additionally, it fixes a bug for codim==0 subIndex of an edge: If the an 3D edge is passed as an argument requesting a 0 codim sub-index, the previous implementation would try to get the index of the 0-th vertex, rather than of the edge itself.
Close #155