ADD boundaries wrapped method.

Author: petrasvestartas

CMakeLists.txt

@@ -161,14 +161,14 @@ function(add_nanobind_module module_name source_file)
 endfunction()
 
 # Add modules
-add_nanobind_module(nanobind src/nanobind.cpp)
-add_nanobind_module(types_std src/types_std.cpp)
-add_nanobind_module(boundaries src/boundaries.cpp)
-add_nanobind_module(curvature src/curvature.cpp)
-add_nanobind_module(geodistance src/geodistance.cpp)
-add_nanobind_module(intersections src/intersections.cpp)
-add_nanobind_module(isolines src/isolines.cpp)
-add_nanobind_module(massmatrix src/massmatrix.cpp)
-add_nanobind_module(meshing src/meshing.cpp)
-add_nanobind_module(parametrisation src/parametrisation.cpp)
-add_nanobind_module(planarize src/planarize.cpp)
+add_nanobind_module(_nanobind src/nanobind.cpp)
+add_nanobind_module(_types_std src/types_std.cpp)
+add_nanobind_module(_boundaries src/boundaries.cpp)
+add_nanobind_module(_curvature src/curvature.cpp)
+add_nanobind_module(_geodistance src/geodistance.cpp)
+add_nanobind_module(_intersections src/intersections.cpp)
+add_nanobind_module(_isolines src/isolines.cpp)
+add_nanobind_module(_massmatrix src/massmatrix.cpp)
+add_nanobind_module(_meshing src/meshing.cpp)
+add_nanobind_module(_parametrisation src/parametrisation.cpp)
+add_nanobind_module(_planarize src/planarize.cpp)

docs/examples/boundaries.py

@@ -43,6 +43,8 @@
 )
 
 for vertices in boundaries:
+    vertices = list(vertices)
+    vertices.append(vertices[0])
     points = mesh.vertices_attributes("xyz", keys=vertices)
     polyline = Polyline(points)
     viewer.scene.add(polyline, linecolor=Color.red(), linewidth=3)

src/boundaries.cpp

@@ -2,22 +2,16 @@
 
 using namespace Eigen;
 
-std::vector<std::vector<int>>
-trimesh_boundaries(
-    compas::RowMatrixXi F)
-{
+std::vector<std::vector<int>> trimesh_boundaries(const MatrixXi& F) {
     std::vector<std::vector<int>> L;
-
-	igl::boundary_loop(F, L);
-
+    igl::boundary_loop(F, L);
     return L;
 }
 
-NB_MODULE(compas_libigl_boundaries, m) {
+NB_MODULE(_boundaries, m) {
     m.def(
         "trimesh_boundaries",
         &trimesh_boundaries,
         "Compute list of ordered boundary loops for a manifold mesh.",
-        "F"_a
-    );
+        "F"_a);
 }

src/compas_libigl/__init__.py

@@ -1,6 +1,7 @@
 import os
 import compas
-from .nanobind import add, __doc__
+from ._nanobind import add, __doc__
+from .boundaries import trimesh_boundaries
 
 # from .boundaries import trimesh_boundaries
 # from .curvature import trimesh_gaussian_curvature, trimesh_principal_curvature
@@ -83,7 +84,7 @@ def get_armadillo():
 
 
 __all_plugins__ = [
-    "compas_libigl.nanobind"
+    "compas_libigl._nanobindcompas_libigl._boundaries",
     # "compas_libigl.curvature",
     # "compas_libigl.geodistance",
     # "compas_libigl.intersections",
@@ -101,7 +102,7 @@ def get_armadillo():
     "TEMP",
     add,
     __doc__,
-
+    trimesh_boundaries,
     # "get",
     # "get_beetle",
     # "get_armadillo",

src/compas_libigl/boundaries.py

@@ -0,0 +1,41 @@
+import numpy as np
+
+from compas_libigl import _boundaries
+from compas_libigl._types_std import VectorVectorInt
+
+
+def trimesh_boundaries(M):
+    """Compute all (ordered) boundary loops of a manifold triangle mesh.
+
+    Parameters
+    ----------
+    M : (list, list)
+        A mesh represented by a list of vertices and a list of faces.
+
+    Returns
+    -------
+    array
+        The ordered boundary loops of the triangle mesh.
+        Each loop is a sequence of vertex indices.
+
+    Notes
+    -----
+    The input mesh should be manifold.
+
+    Examples
+    --------
+    >>> import compas
+    >>> import compas_libigl
+    >>> from compas.datastructures import Mesh
+    >>> mesh = Mesh.from_off(compas.get("tubemesh.off"))
+    >>> mesh.quads_to_triangles()
+    >>> M = mesh.to_vertices_and_faces()
+    >>> boundaries = compas_libigl.trimesh_boundaries(M)
+    >>> len(boundaries) == 1
+    True
+
+    """
+    V, F = M
+    F = np.asarray(F, dtype=np.int32)
+    result: VectorVectorInt = _boundaries.trimesh_boundaries(F)
+    return [list(loop) for loop in result]

src/curvature.cpp

@@ -2,7 +2,7 @@
 
 void function(const Eigen::MatrixXd& V, const Eigen::MatrixXi& F) {}
 
-NB_MODULE(curvature, m) {
+NB_MODULE(_curvature, m) {
 
     m.def(
         "function_name",

src/geodistance.cpp

@@ -2,7 +2,7 @@
 
 void function(const Eigen::MatrixXd& V, const Eigen::MatrixXi& F) {}
 
-NB_MODULE(geodistance, m) {
+NB_MODULE(_geodistance, m) {
 
     m.def(
         "function_name",

src/intersections.cpp

@@ -2,7 +2,7 @@
 
 void function(const Eigen::MatrixXd& V, const Eigen::MatrixXi& F) {}
 
-NB_MODULE(intersections, m) {
+NB_MODULE(_intersections, m) {
 
     m.def(
         "function_name",

src/isolines.cpp

@@ -2,7 +2,7 @@
 
 void function(const Eigen::MatrixXd& V, const Eigen::MatrixXi& F) {}
 
-NB_MODULE(isolines, m) {
+NB_MODULE(_isolines, m) {
 
     m.def(
         "function_name",

src/massmatrix.cpp

@@ -2,7 +2,7 @@
 
 void function(const Eigen::MatrixXd& V, const Eigen::MatrixXi& F) {}
 
-NB_MODULE(massmatrix, m) {
+NB_MODULE(_massmatrix, m) {
 
     m.def(
         "function_name",