When restoring edge constraints, edges would be flipped between
adjacent triangles to try achieving a Delaunay conforming mesh.
However, after each flip, triangle IDs would not be updated, leading
to invalid further processing and producing self intersections.