DOCS

Author: petrasvestartas

docs/api/compas_libigl.rst

@@ -20,5 +20,4 @@ Functions
     trimesh_massmatrix
     trimesh_harmonic
     trimesh_lscm
-    trimesh_remesh_along_isoline
     quadmesh_planarize

docs/examples/curvature.rst

@@ -6,5 +6,13 @@ Curvature
     :figclass: figure
     :class: figure-img img-fluid
 
-.. literalinclude:: curvature.py
+.. literalinclude:: curvature_gaussian.py
+    :language: python
+
+
+.. figure:: /_images/curvature.png
+    :figclass: figure
+    :class: figure-img img-fluid
+
+.. literalinclude:: curvature_principal.py
     :language: python

docs/examples/planarize.py

@@ -3,7 +3,6 @@
 from compas.colors import Color
 from compas.colors import ColorMap
 from compas.datastructures import Mesh
-from compas.datastructures import mesh_flatness
 from compas_viewer import Viewer
 
 # ==============================================================================
@@ -14,31 +13,26 @@
 MAXDEV = 0.005
 KMAX = 500
 
-mesh = Mesh.from_obj(compas.get("tubemesh.obj"))
+mesh_not_planarized = Mesh.from_obj(compas.get("tubemesh.obj"))
+mesh_not_planarized.name = "Not Planarized"
 
 # ==============================================================================
 # Planarize
 # ==============================================================================
 
-V, F = mesh.to_vertices_and_faces()
+V, F = mesh_not_planarized.to_vertices_and_faces()
 V2 = igl.quadmesh_planarize((V, F), KMAX, MAXDEV)
 
 # ==============================================================================
 # Visualize
 # ==============================================================================
 
-mesh = Mesh.from_vertices_and_faces(V2, F)
-dev = mesh_flatness(mesh, maxdev=TOL)
+mesh_planarized = Mesh.from_vertices_and_faces(V2, F)
+mesh_planarized.name = "Planarized"
 
 cmap = ColorMap.from_two_colors(Color.white(), Color.blue())
 
 viewer = Viewer(width=1600, height=900)
-# viewer.view.camera.position = [1, -6, 2]
-# viewer.view.camera.look_at([1, 1, 1])
-
-viewer.scene.add(
-    mesh,
-    facecolor={face: (cmap(dev[face]) if dev[face] <= 1.0 else Color.red()) for face in mesh.faces()},
-    show_points=False,
-)
-viewer.show()
+viewer.scene.add(mesh_not_planarized, show_faces=False)
+viewer.scene.add(mesh_planarized)
+viewer.show()

src/compas_libigl/__init__.py

@@ -8,8 +8,8 @@
 from .isolines import trimesh_isolines, groupsort_isolines
 from .massmatrix import trimesh_massmatrix
 from .parametrisation import trimesh_harmonic, trimesh_lscm
+from .planarize import quadmesh_planarize
 
-# from .planarize import quadmesh_planarize
 # from .meshing import trimesh_remesh_along_isoline
 
 
@@ -115,6 +115,6 @@ def get_armadillo():
     "trimesh_massmatrix",
     "trimesh_harmonic",
     "trimesh_lscm",
-    # "quadmesh_planarize",
+    "quadmesh_planarize",
     # "trimesh_remesh_along_isoline",
 ]

src/compas_libigl/planarize.py

@@ -0,0 +1,36 @@
+import numpy as np
+from compas.plugins import plugin
+
+from compas_libigl import _planarize
+
+
+@plugin(category="quadmesh")
+def quadmesh_planarize(M, kmax=500, maxdev=0.005):
+    """Planarize the faces of a quad mesh.
+
+    Parameters
+    ----------
+    M : tuple or :class:`compas.datastructures.Mesh`
+        A quad mesh represented by a list of vertices and a list of faces
+        or by a COMPAS mesh object.
+    kmax : int, optional
+        The maximum number of iterations.
+        Default is ``500``.
+    maxdev : float, optional
+        The maximum deviation from planar.
+        Default is ``0.005``.
+
+    Returns
+    -------
+    list
+        The coordinates of the new vertices.
+
+    Examples
+    --------
+    >>>
+
+    """
+    V, F = M
+    V = np.asarray(V, dtype=np.float64)
+    F = np.asarray(F, dtype=np.int32)
+    return _planarize.planarize_quads(V, F, kmax, maxdev)

src/planarize.cpp

@@ -1,13 +1,23 @@
 #include "planarize.hpp"
 
-void function(const compas::RowMatrixXd& V, const compas::RowMatrixXi& F) {}
+
+compas::RowMatrixXd planarize_quads(compas::RowMatrixXd V, compas::RowMatrixXi F, int maxiter = 100, double threshold = 0.005) {
+  compas::RowMatrixXd Vplanar;
+
+  igl::planarize_quad_mesh(V, F, maxiter, threshold, Vplanar);
+
+  return Vplanar;
+}
+
 
 NB_MODULE(_planarize, m) {
 
     m.def(
-        "function_name",
-        &function,
-        "Description.",
+        "planarize_quads",
+        &planarize_quads,
+        "Planarize a quad mesh.",
         "V"_a,
-        "F"_a);
+        "F"_a,
+        "maxiter"_a,
+        "threshold"_a);
 }