// Find label of face. label findFace(const primitiveMesh& mesh, const face& f) { const labelList& pFaces = mesh.pointFaces()[f[0]]; forAll(pFaces, i) { label faceI = pFaces[i]; if (mesh.faces()[faceI] == f) { return faceI; } }