@pimms.value<EOL><INDENT>def contained_faces(surface, label):<DEDENT>
msk = np.where(label >= <NUM_LIT:0.5>)[<NUM_LIT:0>]<EOL>fs = np.where(np.sum(np.isin(surface.tess.indexed_faces, msk), axis=<NUM_LIT:0>) == <NUM_LIT:3>)[<NUM_LIT:0>]<EOL>return pimms.imm_array(surface.tess.faces[:,fs])<EOL>
path.contained_faces is a matrix of faces in path.surface that are completely contained
  by the given path; any face that intersects the path will not be included. Faces are
  returned in a (3 x n) matrix of vertex labels.