add mesh split endpoint

Author: Licini

src/booleans.cpp

@@ -1,5 +1,6 @@
 #include "booleans.h"
 #include <CGAL/Polygon_mesh_processing/corefinement.h>
+#include <CGAL/Polygon_mesh_processing/clip.h>
 
 
 namespace PMP = CGAL::Polygon_mesh_processing;
@@ -68,6 +69,28 @@ pmp_boolean_intersection(
 };
 
 
+std::tuple<compas::RowMatrixXd, compas::RowMatrixXi>
+pmp_split(
+    Eigen::Ref<const compas::RowMatrixXd> &VA,
+    Eigen::Ref<const compas::RowMatrixXi> &FA,
+    Eigen::Ref<const compas::RowMatrixXd> &VB,
+    Eigen::Ref<const compas::RowMatrixXi> &FB)
+{
+    Mesh A = compas::mesh_from_vertices_and_faces(VA, FA);
+    Mesh B = compas::mesh_from_vertices_and_faces(VB, FB);
+    // Mesh C;
+
+    // PMP::clip(A, B, CGAL::parameters::clip_volume(false), CGAL::parameters::clip_volume(false));
+    PMP::split(A, B);
+
+    // Result
+
+    std::tuple<compas::RowMatrixXd, compas::RowMatrixXi> R = compas::mesh_to_vertices_and_faces(A);
+
+    return R;
+};
+
+
 void init_booleans(pybind11::module &m)
 {
     pybind11::module submodule = m.def_submodule("booleans");
@@ -95,4 +118,12 @@ void init_booleans(pybind11::module &m)
         pybind11::arg("FA").noconvert(),
         pybind11::arg("VB").noconvert(),
         pybind11::arg("FB").noconvert());
+    
+    submodule.def(
+        "split",
+        &pmp_split,
+        pybind11::arg("VA").noconvert(),
+        pybind11::arg("FA").noconvert(),
+        pybind11::arg("VB").noconvert(),
+        pybind11::arg("FB").noconvert());
 };

src/booleans.h

@@ -28,4 +28,12 @@ pmp_boolean_intersection(
     Eigen::Ref<const compas::RowMatrixXi> &FB);
 
 
+std::tuple<compas::RowMatrixXd, compas::RowMatrixXi>
+pmp_clip(
+    Eigen::Ref<const compas::RowMatrixXd> &VA,
+    Eigen::Ref<const compas::RowMatrixXi> &FA,
+    Eigen::Ref<const compas::RowMatrixXd> &VB,
+    Eigen::Ref<const compas::RowMatrixXi> &FB);
+
+
 #endif /* COMPAS_BOOLEANS_H */

src/compas_cgal/booleans.py

@@ -41,6 +41,8 @@ def _boolean(A, B, operation):
         result = booleans.boolean_difference(VA, FA, VB, FB)
     elif operation == 'intersection':
         result = booleans.boolean_intersection(VA, FA, VB, FB)
+    elif operation == 'split':
+        result = booleans.split(VA, FA, VB, FB)
     else:
         raise NotImplementedError
 
@@ -159,3 +161,9 @@ def boolean_intersection(A, B):
 
     """
     return _boolean(A, B, 'intersection')
+
+@plugin(category='booleans', pluggable_name='split_mesh_mesh')
+def split(A, B):
+    """
+    """
+    return _boolean(A, B, 'split')

src/compas_cgal/booleans.pyi

@@ -4,7 +4,7 @@ from compas_cgal.types import VerticesFaces
 from compas_cgal.types import VerticesFacesNumpy
 
 
-def _boolean(A: VerticesFaces, B: VerticesFaces, operation: Literal['union', 'difference', 'intersection']) -> VerticesFacesNumpy:
+def _boolean(A: VerticesFaces, B: VerticesFaces, operation: Literal['union', 'difference', 'intersection', 'split']) -> VerticesFacesNumpy:
     pass
 
 
@@ -18,3 +18,6 @@ def boolean_difference(A: VerticesFaces, B: VerticesFaces) -> VerticesFacesNumpy
 
 def boolean_intersection(A: VerticesFaces, B: VerticesFaces) -> VerticesFacesNumpy:
     pass
+
+def split(A: VerticesFaces, B: VerticesFaces) -> VerticesFacesNumpy:
+    pass