remove super-slow mesh.tets.reserve(mesh.tets.size() + new_tets.size());