cdt using cgal

Author: tomvanmele

docs/examples/triangulation.py

@@ -0,0 +1,26 @@
+import numpy as np
+
+from compas_plotters import Plotter
+from compas.geometry import Polygon, Translation
+from compas.datastructures import Mesh
+from compas_cgal.triangulation import constrained_delaunay_triangulation
+
+
+boundary = Polygon.from_sides_and_radius_xy(4, 10)
+boundary = np.array(boundary, dtype=np.float64)
+
+hole = Polygon.from_sides_and_radius_xy(8, 2)
+
+hole1 = hole.transformed(Translation.from_vector([4, 0, 0]))
+hole2 = hole.transformed(Translation.from_vector([-4, 0, 0]))
+hole3 = hole.transformed(Translation.from_vector([0, 4, 0]))
+hole4 = hole.transformed(Translation.from_vector([0, -4, 0]))
+
+V, F = constrained_delaunay_triangulation(boundary, [hole1, hole2, hole3, hole4], [], maxlength=1.0)
+
+mesh = Mesh.from_vertices_and_faces(V, F)
+
+plotter = Plotter(figsize=(8, 8))
+plotter.add(mesh, sizepolicy='absolute', vertexsize=5)
+plotter.zoom_extents()
+plotter.show()

docs/examples/triangulation.rst

@@ -0,0 +1,11 @@
+********************************************************************************
+Triangulation
+********************************************************************************
+
+.. figure:: /_images/cgal_triangulation.png
+    :figclass: figure
+    :class: figure-img img-fluid
+
+
+.. literalinclude:: triangulation.py
+    :language: python

scripts/cdt.py

@@ -55,11 +55,11 @@ def get_library_dirs():
         sorted([
             'src/compas_cgal.cpp',
             'src/compas.cpp',
-            # 'src/meshing.cpp',
-            # 'src/booleans.cpp',
-            # 'src/slicer.cpp',
-            # 'src/intersections.cpp',
-            # 'src/measure.cpp',
+            'src/meshing.cpp',
+            'src/booleans.cpp',
+            'src/slicer.cpp',
+            'src/intersections.cpp',
+            'src/measure.cpp',
             'src/triangulations.cpp',
         ]),
         include_dirs=[

setup.py

@@ -3,11 +3,11 @@
 
 namespace py = pybind11;
 
-// void init_meshing(py::module&);
-// void init_booleans(py::module&);
-// void init_slicer(py::module&);
-// void init_intersections(py::module&);
-// void init_measure(py::module&);
+void init_meshing(py::module&);
+void init_booleans(py::module&);
+void init_slicer(py::module&);
+void init_intersections(py::module&);
+void init_measure(py::module&);
 void init_triangulations(py::module&);
 
 PYBIND11_MODULE(_cgal, m) {
@@ -17,10 +17,10 @@ PYBIND11_MODULE(_cgal, m) {
     	.def_readonly("vertices", &compas::Result::vertices)
     	.def_readonly("faces", &compas::Result::faces);
 
-    // init_meshing(m);
-    // init_booleans(m);
-    // init_slicer(m);
-    // init_intersections(m);
-    // init_measure(m);
+    init_meshing(m);
+    init_booleans(m);
+    init_slicer(m);
+    init_intersections(m);
+    init_measure(m);
     init_triangulations(m);
 }

src/compas_cgal.cpp

@@ -64,6 +64,15 @@
 
     slicer.slice_mesh
 
+Slicer
+------
+
+.. autosummary::
+    :toctree: generated/
+    :nosignatures:
+
+    triangulation.constrained_delaunay_triangulation
+
 """
 from __future__ import print_function
 

src/compas_cgal/__init__.py

@@ -0,0 +1,43 @@
+import numpy as np
+from compas_cgal._cgal import triangulations
+
+
+def constrained_delaunay_triangulation(boundary, holes=None, curves=None, maxlength=None):
+    """Construct a constrained Delaunay triangulation.
+
+    Parameters
+    ----------
+    boundary : :class:`compas.geometry.Polygon`
+        The boundary of the triangulation.
+    holes : list of :class:`compas.geometry.Polygon`, optional
+        Internal boundary polygons.
+    curves : list of :class:`compas.geometry.Polyline`, optional
+        Internal constraint curves.
+    maxlength : float, optional
+        The maximum length of the triangle edges.
+
+    Returns
+    -------
+    list
+        The vertices and faces of the triangulation.
+
+    Notes
+    -----
+    ...
+
+    """
+    boundary = np.array(boundary, dtype=np.float64)
+
+    if holes:
+        holes = [np.array(hole, dtype=np.float64) for hole in holes]
+    else:
+        holes = []
+
+    if curves:
+        curves = [np.array(curve, dtype=np.float64) for curve in curves]
+    else:
+        curves = []
+
+    maxlength = maxlength or 0.0
+
+    return triangulations.constrained_delaunay_triangulation(boundary, holes, curves, maxlength=maxlength)

src/compas_cgal/triangulation.py

@@ -0,0 +1,13 @@
+from typing import List
+from typing import Optional
+
+from compas.geometry import Polygon
+from compas.geometry import Polyline
+
+from compas_cgal.types import VerticesFacesNumpy
+
+
+def constrained_delaunay_triangulation(boundary: Polygon,
+                                       holes: Optional[List[Polygon]] = None,
+                                       curves: Optional[List[Polyline]] = None,
+                                       maxlength: Optional[float] = None) -> VerticesFacesNumpy: ...