@pimms.value<EOL><INDENT>def intersected_edge_indices(edge_data, surface):<DEDENT>
idx = surface.tess.index<EOL>return pimms.imm_array([idx[(u,v)] for (u,v) in zip(edge_data[<NUM_LIT:0>], edge_data[<NUM_LIT:1>])])<EOL>
path.intersected_edge_indices is equivalent to to path.intersected_edges except that it
stores a list of edge indices instead of a matrix of edges.