// Find edge between points v0 and v1. label findEdge(const primitiveMesh& mesh, const label v0, const label v1) { const labelList& pEdges = mesh.pointEdges()[v0]; forAll(pEdges, pEdgeI) { label edgeI = pEdges[pEdgeI]; const edge& e = mesh.edges()[edgeI]; if (e.otherVertex(v0) == v1) { return edgeI; } }