Would it make sense to add an epsilon here? So you can find mesh nodes that are not exactly on the surface but within the defined distance epsilon.