Fix pathological case with CDT cavity doubling back on itself (#98)

JuliaGeometry · Nov 23, 2023 · 7762f1b · 7762f1b · DanielVandH · Nov 23, 2023
2 parents 7f0ed94 + 22a5128
commit 7762f1b
Show file tree

Hide file tree

Showing 3 changed files with 4 additions and 3 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "DelaunayTriangulation"
 uuid = "927a84f5-c5f4-47a5-9785-b46e178433df"
 authors = ["Daniel VandenHeuvel <[email protected]>"]
-version = "0.8.10"
+version = "0.8.11"
 
 [deps]
 DataStructures = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"

diff --git a/src/constrained_triangulation/segment_insertion.jl b/src/constrained_triangulation/segment_insertion.jl
@@ -172,6 +172,7 @@ end
 
 function add_point_cavity_cdt!(tri::Triangulation, marked_vertices, u, v, w)
     x = get_adjacent(tri, w, v)
+    (u == v || v == w || u == w) && return tri
     if !edge_exists(x)
         insert_flag = true
     else

diff --git a/test/triangulation/constrained.jl b/test/triangulation/constrained.jl
@@ -5,10 +5,10 @@ using StableRNGs
 include("../helper_functions.jl")
 
 @testset "Test random constrained Delaunay triangulations" begin
-    for i in 1:250
+    for i in 1:25000
         rng = StableRNG(i)
         points, edges, mat_edges = get_random_vertices_and_constrained_edges(40, 200, 20, rng)
-        tri = triangulate(points; edges, rng)
+        tri = triangulate(points; edges, rng=StableRNG(i))
         @test validate_triangulation(tri)
         empty!(get_all_constrained_edges(tri))
         @test !validate_triangulation(tri)