Julia Delaunay 三角剖分无限期挂起

Question

我正在尝试掌握 Julia 中的 Delaunay 三角剖分。我正在尝试以下测试代码。

using VoronoiDelaunay
a = Point2D[Point(rand(),rand()) for i in 1:6]
tess = DelaunayTessellation()
push!(tess,a)

在将最终命令输入 REPL 时，代码似乎无限期挂起。我已经试了好几次了，放了好久运行，但一直没有完成。有谁知道我做错了什么吗？

Answer 1

看起来 VoronoiDelaunay package currently limits input points 到 [1 + eps, 2 - 2*eps] 的范围。通过此更改，我能够将您的代码设为运行：

a = Point2D[Point(1 + rand(),1 + rand()) for i in 1:6]
tess = DelaunayTessellation()
push!(tess,a)

跟随the example in their README，你还可以使用内置的max_coord和min_coord，这应该比我第一次尝试更健壮：

width = max_coord - min_coord
a = Point2D[Point(min_coord + rand() * width, min_coord + rand() * width) for i in 1:6]

注意我假设你的意思是 DelaunayTesselation 因为 DelaunayTriangulation 似乎不是包提供的方法。

我的结果：

Julia Delaunay triangulation hangs indefinitely