Implement Delaunay propagation

Signed-off-by: pacien <pacien.trangirard@pacien.net>

2 files changed, 21 insertions, 0 deletions

diff --git a/include/morpher/quadrilateral.h b/include/morpher/quadrilateral.h
index c8cf3a1..8691263 100644
--- a/include/morpher/quadrilateral.h
+++ b/include/morpher/quadrilateral.h
@@ -36,4 +36,15 @@ void quadrilateral_flip_diagonal(TriangleMap *t1, TriangleMap *t2);
  */
 bool quadrilateral_is_delaunay(TriangleMap *t1, TriangleMap *t2);
 
+/**
+ * Function: quadrilateral_propagate_delaunay
+ * Ensures that the quadrilateral spawned by the given triangles fulfills the Delaunay criterion,
+ * flipping the diagonal if necessary and propagating the changes to the neighbouring triangles.
+ *
+ * Parameters:
+ *   *start    - the starting triangle
+ *   *neighbor - a neighboring triangle
+ */
+void quadrilateral_propagate_delaunay(TriangleMap *start, TriangleMap *neighbor);
+
 #endif
diff --git a/include/morpher/trianglemap.h b/include/morpher/trianglemap.h
index 05e7b87..0bb6a1d 100644
--- a/include/morpher/trianglemap.h
+++ b/include/morpher/trianglemap.h
@@ -106,6 +106,16 @@ void trianglemap_set_neighbors(TriangleMap *t, TriangleMap *n1, TriangleMap *n2,
 void trianglemap_replace_neighbor(TriangleMap *t, TriangleMap *old, TriangleMap *new);
 
 /**
+ * Function: trianglemap_foreach_neighbor
+ * Executes the given function for each existing (non-NULL) neighbour of the supplied triangle.
+ *
+ * Parameters:
+ *   *t - the base triangle
+ *   *f - the function
+ */
+void trianglemap_foreach_neighbor(TriangleMap *t, void (*f)(TriangleMap *current, TriangleMap *neighbor));
+
+/**
  * Function: trianglemap_split
  * Splits a triangle into three sub-triangles at the supplied center vertex, updating the surrounding triangles.
  * The first triangle resulting from the split is returned, with the two others chained as linear neighbours.