getEdges.py

edges = {}
F = meshes[0]["Faces"].copy()
F.sort(1)
nf = F.shape[0]
for i in range(nf):
    edges.setdefault(tuple(F[i,:2]), []).append(i)
    edges.setdefault(tuple(F[i,1:3]), []).append(i)
    edges.setdefault(tuple(F[i,[0,2]]), []).append(i)

neighbours = {}
for e in edges.values():
    neighbours.setdefault(e[0], []).append(e[1])
    neighbours.setdefault(e[1], []).append(e[0])

neighbours = array(list(neighbours.values()))
neighbours.sort(1)