Ruppert's Refinement Algorithm.

CDT/include/CDTUtils.h

@@ -44,6 +44,7 @@ typedef char couldnt_parse_cxx_standard[-1]; ///< Error: couldn't parse standard
 #include <cassert>
 #include <cmath>
 #include <limits>
+#include <queue>
 #include <vector>
 
 #ifdef CDT_USE_STRONG_TYPING
@@ -243,6 +244,7 @@ inline Edge edge_make(VertInd iV1, VertInd iV2)
 }
 
 typedef std::vector<Edge> EdgeVec;                ///< Vector of edges
+typedef std::queue<Edge> EdgeQue;                 ///< Queue of edges
 typedef unordered_set<Edge> EdgeUSet;             ///< Hash table of edges
 typedef unordered_set<TriInd> TriIndUSet;         ///< Hash table of triangles
 typedef unordered_map<TriInd, TriInd> TriIndUMap; ///< Triangle hash map

CDT/include/Triangulation.h

@@ -83,6 +83,22 @@ struct CDT_EXPORT IntersectingConstraintEdges
     };
 };
 
+/**
+ * Enum of strategies for triangles refinement
+ */
+struct CDT_EXPORT RefineTriangles
+{
+    /**
+     * The Enum itself
+     * @note needed to pre c++11 compilers that don't support 'class enum'
+     */
+    enum Enum
+    {
+        ByAngle, ///< constraint minimum triangles angle
+        ByArea,  ///< constraint maximum triangles area
+    };
+};
+
 /**
  * Type used for storing layer depths for triangles
  * @note LayerDepth should support 60K+ layers, which could be to much or
@@ -110,8 +126,14 @@ class CDT_EXPORT Triangulation
 {
 public:
     typedef std::vector<V2d<T> > V2dVec; ///< Vertices vector
+    typedef std::vector<bool> boolVec;   ///< Steiner Vertices flag
     V2dVec vertices;                     ///< triangulation's vertices
-    TriangleVec triangles;               ///< triangulation's triangles
+    boolVec isSteinerVertex; ///< triangulation's vertices Steiner point flag
+    IndexSizeType
+        maxSteinerPoints; ///< triangulation's maximum number of Steiner
+    IndexSizeType numOfSteinerPoints; ///< triangulation's maximum number of
+                                      ///< Steiner points to be added
+    TriangleVec triangles;            ///< triangulation's triangles
     EdgeUSet fixedEdges; ///< triangulation's constraints (fixed edges)
 
     /** Stores count of overlapping boundaries for a fixed edge. If no entry is
@@ -275,6 +297,16 @@ class CDT_EXPORT Triangulation
      * @tparam edges edges to conform to
      */
     void conformToEdges(const std::vector<Edge>& edges);
+    /**
+     * Traingles refinement by removing bad triangles
+     * @note bad triangles don't fulfill constraints defined by the user
+     * @param refinement_constrain refinement strategy that is used to identify
+     * bad triangles
+     * @param threshold threshold value for refinement
+     */
+    void refineTriangles(
+        RefineTriangles::Enum refinementConstrain = RefineTriangles::ByAngle,
+        T threshold = 20 / 180.0 * M_PI);
     /**
      * Erase triangles adjacent to super triangle
      *
@@ -360,7 +392,7 @@ class CDT_EXPORT Triangulation
 private:
     /*____ Detail __*/
     void addSuperTriangle(const Box2d<T>& box);
-    void addNewVertex(const V2d<T>& pos, TriInd iT);
+    void addNewVertex(const V2d<T>& pos, TriInd iT, bool isSteiner = false);
     void insertVertex(VertInd iVert);
     void insertVertex(VertInd iVert, VertInd walkStart);
     void ensureDelaunayByEdgeFlips(
@@ -488,6 +520,32 @@ class CDT_EXPORT Triangulation
         VertInd iV2,
         VertInd iV3,
         VertInd iV4) const;
+    TriInd edgeTriangle(Edge edge) const;
+    /// Checks if edge e is encroached by vertex v
+    bool isEncroached(const V2d<T>& v, Edge e) const;
+    bool isBadTriangle(
+        const Triangle& tri,
+        RefineTriangles::Enum refinement,
+        T threshold) const;
+    V2d<T> circumcenter(const Triangle& tri) const;
+    /// Search in all fixed edges to find encroached edges, each fixed edge is
+    /// checked against its opposite vertices
+    /// Returns queue of encroached edges
+    EdgeQue detectEncroachedEdges();
+    /// Search in all fixed edges to find encroached edges, each fixed edge is
+    /// checked against its opposite vertices and vertex v
+    /// Returns queue of encroached edges
+    EdgeQue detectEncroachedEdges(const V2d<T>& v);
+    /// Recursively split encroached edges
+    /// returns vector of badly shaped triangles and number of splits
+    std::pair<TriIndVec, IndexSizeType> resolveEncroachedEdges(
+        EdgeQue encroachedEdges,
+        const V2d<T>& v = {},
+        bool validV = false,
+        bool fillBadTriangles = false,
+        RefineTriangles::Enum refinementConstrain = {},
+        T badTriangleThreshold = {});
+    VertInd splitEncroachedEdge(Edge e, TriInd iT, TriInd iTopo);
     void changeNeighbor(TriInd iT, TriInd oldNeighbor, TriInd newNeighbor);
     void changeNeighbor(
         TriInd iT,
@@ -744,7 +802,6 @@ void Triangulation<T, TNearPointLocator>::conformToEdges(
     }
     eraseDummies();
 }
-
 } // namespace CDT
 
 #ifndef CDT_USE_AS_COMPILED_LIBRARY