It is astonishing that add_face works for the PolyMesh but not for the TriMesh. The function is_isolated(V) does not the right thing for a TriMesh