GJK error in interpolation between non matching second ordered 3D meshes

[chenyongxin]

I met the same GJK error when do non-matching interpolation between two 2nd order 3D meshes. See 3D ALE-FSI on a local refined 2nd order mesh: seek advices on remeshing for the my use case. With DOLFINx v0.9.0, only to replace

@Lingyue_Shen

def nmm_interpolate(f_out, f_in):
        interpolation_data=create_nonmatching_meshes_interpolation_data(
                                                                        f_out.function_space.mesh,
                                                                        f_out.function_space.element,
                                                                        f_in.function_space.mesh,
                                                                        padding=1e-11)
        f_out.interpolate(f_in,nmm_interpolation_data=interpolation_data)       

by

def nmm_interpolate(f_out, f_in):
        index_map = f_out.function_space.mesh.topology.index_map(3)    
        ncells = index_map.size_local + index_map.num_ghosts
        cells = np.arange(ncells, dtype=np.int32)
        interpolation_data = fem.create_interpolation_data(f_out.function_space, f_in.function_space, 
                                                           cells, padding=1e-11)
        f_out.interpolate_nonmatching(f_in, cells, interpolation_data=interpolation_data)

the error can be reproduced. Any suggestion on this? Or maybe at least how can we know which mesh point cause this error?