Add missing handling of neighbor list (-n flag)

Author: j-fu

Project.toml

@@ -1,7 +1,7 @@
 name = "Triangulate"
 uuid = "f7e6ffb2-c36d-4f8f-a77e-16e897189344"
 authors = ["Juergen Fuhrmann <[email protected]>"]
-version = "2.1.2"
+version = "2.1.3"
 
 [deps]
 DocStringExtensions = "ffbed154-4ef7-542d-bbb7-c09d3a79fcae"

src/triangulateio.jl

@@ -337,6 +337,10 @@ function TriangulateIO(ctio::CTriangulateIO)
         tio.trianglearealist=convert(Array{Cdouble,1}, Base.unsafe_wrap(Array, ctio.trianglearealist, (Int(ctio.numberoftriangles)), own=true))
     end
 
+    if ctio.numberoftriangles>0 && ctio.neighborlist!=C_NULL
+        tio.neighborlist=convert(Array{Cint,2}, Base.unsafe_wrap(Array, ctio.neighborlist, (3,Int(ctio.numberoftriangles)), own=true))
+    end
+
     if ctio.numberofsegments>0  && ctio.segmentlist!=C_NULL && ctio.segmentmarkerlist!=C_NULL
         tio.segmentlist=convert(Array{Cint,2}, Base.unsafe_wrap(Array, ctio.segmentlist, (2,Int(ctio.numberofsegments)), own=true))
         tio.segmentmarkerlist=convert(Array{Cint,1}, Base.unsafe_wrap(Array, ctio.segmentmarkerlist, (Int(ctio.numberofsegments)), own=true))

test/test_neighbors.jl

@@ -0,0 +1,29 @@
+# thanks @tduretz
+module test_neighbors
+
+# Test correct handling of second order points
+
+using Triangulate
+
+function test()
+    triin=Triangulate.TriangulateIO()
+    triin.pointlist=Cdouble[0.0 0.0 ; 1.0 0.0 ; 1.0  1.0; 0.0 1.0]'
+    triin.segmentlist=Cint[1 2 ; 2 3 ; 3 4 ; 4 1 ]'
+    triin.segmentmarkerlist=Cint[1, 2, 3, 4]
+
+
+    (triout, vorout)=triangulate("pqQn", triin)  
+    test1 = triout.pointlist==[0.0  1.0  1.0  0.0;
+                               0.0  0.0  1.0  1.0]
+    
+    test1 = test1 && triout.trianglelist==[4  2;
+                                           1  3;
+                                           2  4]
+
+    test2=triout.neighborlist==[-1 -1;
+                                 2  1;
+                                -1 -1]
+    test1&&test2
+end
+
+end