pytetwild 0.1.dev0__cp312-cp312-win_amd64.whl → 0.1.1__cp312-cp312-win_amd64.whl

pytetwild/__init__.py

@@ -2,16 +2,15 @@
 
 
 # start delvewheel patch
-def _delvewheel_patch_1_5_2():
+def _delvewheel_patch_1_11_2():
     import os
-    libs_dir = os.path.abspath(os.path.join(os.path.dirname(__file__), os.pardir, 'pytetwild.libs'))
-    if os.path.isdir(libs_dir):
+    if os.path.isdir(libs_dir := os.path.abspath(os.path.join(os.path.dirname(__file__), os.pardir, 'pytetwild.libs'))):
         os.add_dll_directory(libs_dir)
 
 
-_delvewheel_patch_1_5_2()
-del _delvewheel_patch_1_5_2
+_delvewheel_patch_1_11_2()
+del _delvewheel_patch_1_11_2
 # end delvewheel patch
 
 from ._version import __version__  # noqa: F401
-from .pytetwild import tetrahedralize, tetrahedralize_pv  # noqa: F401
+from .pytetwild import tetrahedralize, tetrahedralize_pv, tetrahedralize_csg  # noqa: F401

pytetwild/_version.py

@@ -1,4 +1,5 @@
 """Contains the pytetwild version."""
+
 from importlib import metadata
 
 __version__ = metadata.version("pytetwild")

pytetwild/py.typed

@@ -0,0 +1 @@
+partial

pytetwild/pytetwild.py

@@ -1,14 +1,33 @@
 """Wrapper for fTetWild."""
+
 import warnings
 import numpy as np
-from pytetwild import PyfTetWildWrapper
+from pytetwild import PyfTetWildWrapper  # type: ignore
 from typing import Tuple, TYPE_CHECKING
+import numpy.typing as npt
 
 if TYPE_CHECKING:
     import pyvista as pv
 
 
-def tetrahedralize_pv(mesh: "pv.PolyData") -> "pv.UnstructuredGrid":
+def _check_edge_length(edge_length_fac: float) -> None:
+    """Check edge length.
+
+    Parameters
+    ----------
+    edge_length_fac : float, default: 0.05
+        Tetrahedral edge length as a function of bounding box diagonal. The
+        default ideal edge length is bb/20 (bounding box divided by 20).
+
+
+    """
+    if edge_length_fac > 1.0 or edge_length_fac < 0.0 + 1e-16:
+        raise ValueError("Edge length factor must be between 1e-16 and 1.0")
+
+
+def tetrahedralize_pv(
+    mesh: "pv.PolyData", edge_length_fac: float = 0.05, optimize: bool = True
+) -> "pv.UnstructuredGrid":
     """
     Convert a PyVista surface mesh to a PyVista unstructured grid.
 
@@ -16,19 +35,44 @@ def tetrahedralize_pv(mesh: "pv.PolyData") -> "pv.UnstructuredGrid":
     ----------
     mesh : pv.PolyData
         The input surface mesh.
+    edge_length_fac : float, default: 0.05
+        Tetrahedral edge length as a function of bounding box diagonal. The
+        default ideal edge length is bb/20 (bounding box divided by 20).
+    optimize : bool
+        Improve the minimum scaled Jacobean for each cell. This leads to higher
+        cell quality at the expense of computation time.
 
     Returns
     -------
     pv.UnstructuredGrid
-        The converted unstructured grid.
+        The converted unstructured grid containing only tetrahedra.
+
+    Examples
+    --------
+    >>> import pyvista as pv
+    >>> import pytetwild
+    >>> surface_mesh = pv.Sphere()
+    >>> tetrahedral_mesh = pytetwild.tetrahedralize_pv(surface_mesh)
+    >>> tetrahedral_mesh
+    UnstructuredGrid (0x7ff568593d00)
+      N Cells:    7247
+      N Points:   1658
+      X Bounds:   -4.993e-01, 4.993e-01
+      Y Bounds:   -4.965e-01, 4.965e-01
+      Z Bounds:   -5.000e-01, 4.999e-01
+      N Arrays:   0
+
     """
     try:
         import pyvista as pv
     except:
         raise ModuleNotFoundError(
-            "Install PyVista to use this feature with:\n\n" "pip install pytetwild[all]"
+            "Install PyVista to use this feature with:\n\npip install pytetwild[all]"
         )
 
+    if not isinstance(mesh, pv.PolyData):
+        raise TypeError(f"`mesh` must be a pyvista.PolyData, got {type(mesh)}")
+
     if not mesh.is_all_triangles:
         warnings.warn(
             "Input mesh is not all triangles. Either call `.triangulate()`"
@@ -36,13 +80,14 @@ def tetrahedralize_pv(mesh: "pv.PolyData") -> "pv.UnstructuredGrid":
         )
         mesh = mesh.triangulate()
 
-    vertices = np.array(mesh.points, dtype=np.float64)
-    faces = np.array(mesh.faces.reshape((-1, 4))[:, 1:4], dtype=np.int32)
+    _check_edge_length(edge_length_fac)
+    vertices = mesh.points.astype(np.float64, copy=False)
+    faces = mesh.faces.reshape((-1, 4))[:, 1:4].astype(np.float32, copy=False)
 
     (
         tetrahedral_mesh_vertices,
         tetrahedral_mesh_tetrahedra,
-    ) = PyfTetWildWrapper.tetrahedralize_mesh(vertices, faces)
+    ) = PyfTetWildWrapper.tetrahedralize_mesh(vertices, faces, optimize, edge_length_fac)
 
     cells = np.hstack(
         [
@@ -52,10 +97,15 @@ def tetrahedralize_pv(mesh: "pv.PolyData") -> "pv.UnstructuredGrid":
     )
     cell_types = np.full(tetrahedral_mesh_tetrahedra.shape[0], 10, dtype=np.uint8)
 
-    return pv.UnstructuredGrid(cells, cell_types, tetrahedral_mesh_vertices)
+    return pv.UnstructuredGrid(cells, cell_types, tetrahedral_mesh_vertices).clean()
 
 
-def tetrahedralize(vertices: np.ndarray, faces: np.ndarray) -> Tuple[np.ndarray, np.ndarray]:
+def tetrahedralize(
+    vertices: npt.NDArray[np.float64],
+    faces: npt.NDArray[np.int32],
+    optimize: bool = True,
+    edge_length_fac: float = 0.05,
+) -> Tuple[npt.NDArray[np.float64], npt.NDArray[np.int32]]:
     """
     Convert mesh vertices and faces to a tetrahedral mesh.
 
@@ -65,10 +115,84 @@ def tetrahedralize(vertices: np.ndarray, faces: np.ndarray) -> Tuple[np.ndarray,
         The vertices of the mesh.
     faces : np.ndarray
         The faces of the mesh.
+    optimize : bool
+        Improve the minimum scaled Jacobean for each cell. This leads to higher
+        cell quality at the expense of computation time.
+    edge_length_fac : float, default: 0.05
+        Tetrahedral edge length as a function of bounding box diagonal. The
+        default ideal edge length is bb/20 (bounding box divided by 20).
 
     Returns
     -------
     Tuple[np.ndarray, np.ndarray]
         A tuple containing the vertices and tetrahedra of the tetrahedral mesh.
     """
-    return PyfTetWildWrapper.tetrahedralize_mesh(vertices, faces)
+    _check_edge_length(edge_length_fac)
+    if not isinstance(vertices, np.ndarray):
+        raise TypeError("`vertices` must be a numpy array")
+    if not isinstance(faces, np.ndarray):
+        raise TypeError("`faces` must be a numpy array")
+    vertices = vertices.astype(np.float64, copy=False)
+    faces = faces.astype(np.int32, copy=False)
+    return PyfTetWildWrapper.tetrahedralize_mesh(vertices, faces, optimize, edge_length_fac)
+
+
+def tetrahedralize_csg(
+    csg_file: str,
+    epsilon: float = 1e-3,
+    edge_length_r: float = 0.05,
+    stop_energy: float = 10.0,
+    coarsen: bool = True,
+    num_threads: int = 0,
+    loglevel: int = 3,
+) -> "pv.UnstructuredGrid":
+    """
+    Generate a tetrahedral mesh based on a the CSG tree specified in the csf_file.
+
+    Parameters
+    ----------
+    csg_file : str
+        Path to the input json file.
+    epsilon : float, default 1e-3
+        Envelop size, specifying the maximum distance of the output surface from the input surface,
+        relative to the bounding box size.
+    edge_length_r : float, default: 0.05
+        Tetrahedral edge length as a function of bounding box diagonal. The
+        default ideal edge length is bb/20 (bounding box divided by 20).
+    stop_energy : float, default: 10.0
+        The mesh optimization stops when the  conformal AMIPS energy reaches 'stop_energy'.
+    coarsen : bool, default: true
+        Coarsen the output as much as possible, while maintaining the mesh quality.
+    num_threads : int, default: 0
+        Set number of threads used (0 means all available cores).
+    loglevel : int, default: 6
+        Set log level (0 = most verbose, 6 = minimal output).
+
+    Returns
+    -------
+    pv.UnstructuredGrid
+        The converted unstructured grid containing only tetrahedra,
+        with a cell attribute 'marker' indicating which of the input surfaces the cell belongs to.
+    """
+    try:
+        import pyvista as pv
+    except:
+        raise ModuleNotFoundError(
+            "Install PyVista to use this feature with:\n\npip install pytetwild[all]"
+        )
+    (tetrahedral_mesh_vertices, tetrahedral_mesh_tetrahedra, tetrahedral_marker) = (
+        PyfTetWildWrapper.tetrahedralize_csg(
+            csg_file, epsilon, edge_length_r, stop_energy, coarsen, num_threads, loglevel
+        )
+    )
+    cells = np.hstack(
+        [
+            np.full((tetrahedral_mesh_tetrahedra.shape[0], 1), 4, dtype=np.int32),
+            tetrahedral_mesh_tetrahedra,
+        ]
+    )
+    cell_types = np.full(tetrahedral_mesh_tetrahedra.shape[0], 10, dtype=np.uint8)
+
+    grid = pv.UnstructuredGrid(cells, cell_types, tetrahedral_mesh_vertices).clean()
+    grid["marker"] = tetrahedral_marker
+    return grid

pytetwild-0.1.1.dist-info/DELVEWHEEL

@@ -0,0 +1,2 @@
+Version: 1.11.2
+Arguments: ['C:\\Users\\runneradmin\\AppData\\Local\\Temp\\cibw-run-aufwxr3e\\cp312-win_amd64\\build\\venv\\Scripts\\delvewheel', 'repair', '-w', 'C:\\Users\\runneradmin\\AppData\\Local\\Temp\\cibw-run-aufwxr3e\\cp312-win_amd64\\repaired_wheel', 'C:\\Users\\runneradmin\\AppData\\Local\\Temp\\cibw-run-aufwxr3e\\cp312-win_amd64\\built_wheel\\pytetwild-0.1.1-cp312-cp312-win_amd64.whl', '--add-path', '.']

{pytetwild-0.1.dev0.dist-info → pytetwild-0.1.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
-Metadata-Version: 2.1
+Metadata-Version: 2.2
 Name: pytetwild
-Version: 0.1.dev0
+Version: 0.1.1
 Summary: Python wrapper of fTetWild
 Author-Email: Alex Kaszynski <akascap@gmail.com>
 Classifier: Development Status :: 3 - Alpha
@@ -10,22 +10,20 @@ Classifier: License :: OSI Approved :: Mozilla Public License 2.0 (MPL 2.0)
 Classifier: Operating System :: Microsoft :: Windows
 Classifier: Operating System :: POSIX
 Classifier: Operating System :: MacOS
-Classifier: Programming Language :: Python :: 3.8
-Classifier: Programming Language :: Python :: 3.9
 Classifier: Programming Language :: Python :: 3.10
 Classifier: Programming Language :: Python :: 3.11
 Classifier: Programming Language :: Python :: 3.12
+Classifier: Programming Language :: Python :: 3.13
 Requires-Dist: numpy
+Provides-Extra: all
 Requires-Dist: pyvista; extra == "all"
+Provides-Extra: dev
 Requires-Dist: pytest; extra == "dev"
 Requires-Dist: pre-commit; extra == "dev"
 Requires-Dist: pyvista; extra == "dev"
 Requires-Dist: scipy; extra == "dev"
 Requires-Dist: meshio; extra == "dev"
-Provides-Extra: all
-Provides-Extra: dev
 Description-Content-Type: text/x-rst
-License-File: LICENSE
 
 pytetwild
 #########
@@ -62,28 +60,72 @@ more details.
 Usage
 *****
 
-To tetrahedralize a surface mesh:
+To tetrahedralize a surface mesh from `PyVista <https://docs.pyvista.org>`_:
+
+.. code:: py
+
+   import pyvista as pv
+   import pytetwild
+
+   # Load or create a PyVista PolyData surface mesh
+   # Here, we'll create a simple sphere mesh as an example
+   surface_mesh = pv.Icosphere(nsub=2)
+
+   # Convert the surface mesh to a tetrahedral mesh. For this example let's
+   # use a coarse mesh
+   tetrahedral_mesh = pytetwild.tetrahedralize_pv(surface_mesh, edge_length_fac=1)
+
+   # Visualize the tetrahedral mesh in an "exploded" view
+   tetrahedral_mesh.explode(1).plot(show_edges=True)
+
+.. image:: https://github.com/pyvista/pytetwild/raw/main/exploded-sphere.png
+
+You can also work with raw arrays. Here's a simple cube that we turn into tetrahedra.
+
+.. code:: py
+
+   import numpy as np
 
-.. code:: pycon
+   # Define vertices of the cube
+   vertices = np.array([
+       [0, 0, 0],  # Vertex 0
+       [1, 0, 0],  # Vertex 1
+       [1, 1, 0],  # Vertex 2
+       [0, 1, 0],  # Vertex 3
+       [0, 0, 1],  # Vertex 4
+       [1, 0, 1],  # Vertex 5
+       [1, 1, 1],  # Vertex 6
+       [0, 1, 1]   # Vertex 7
+   ])
 
-   >>> import pytetwild
-   >>> mesh = pytetwild.tetrahedralize("input_mesh.ply")
-   >>> mesh.vertices
-   array([[x1, y1, z1],
-          [x2, y2, z2],
-          ...], dtype=float32)
-   >>> mesh.tetrahedra
-   array([[i1, j1, k1, l1],
-          [i2, j2, k2, l2],
-          ...], dtype=int32)
+   # Define faces using vertex indices
+   # Each face is a rectangle (also accepts triangles)
+   faces = np.array([
+       [0, 1, 2, 3],  # Front face
+       [1, 5, 6, 2],  # Right face
+       [5, 4, 7, 6],  # Back face
+       [4, 0, 3, 7],  # Left face
+       [4, 5, 1, 0],  # Bottom face
+       [3, 2, 6, 7]   # Top face
+   ])
+   v_out, tetra = pytetwild.tetrahedralize(vertices, faces, optimize=False)
 
-You can also load a mesh, perform tetrahedralization, and export the tetrahedral mesh:
 
-.. code:: pycon
+Usage - Options
+---------------
+We've surfaced a handful of parameters to each of our interfaces
+``tetrahedralize`` and ``tetrahedralize_pv``. Here are the optional parameters.
 
-    >>> import pytetwild
-    >>> pytetwild.tetrahedralize("input_mesh.ply", "output_mesh.ply")
+.. code::
 
+    Additional Parameters
+    ---------------------
+    edge_length_fac : float, default: 0.05
+        Tetrahedral edge length as a function of bounding box diagonal. The
+        default ideal edge length is bb/20 (bounding box divided by 20).
+    optimize : bool
+        Improve the minimum scaled Jacobean for each cell. This leads to higher
+        cell quality at the expense of computation time.
 
 
 License and Acknowledgments

pytetwild-0.1.1.dist-info/RECORD

@@ -0,0 +1,13 @@
+pytetwild/py.typed,sha256=mDShSrm8qg9qjacQc2F-rI8ATllqP6EdgHuEYxuCXZ0,7
+pytetwild/PyfTetWildWrapper.pyd,sha256=ol9Knvqr5Kf4MHq22-X3yRVJUgKYYswgHV0x8bdjqAk,5940736
+pytetwild/pytetwild.py,sha256=JZNJzIsoC9fOLqb4xDJ7oREJarEYkdvCgbIgRCvp52Q,6977
+pytetwild/_version.py,sha256=A4_1DWpVQi-Rp5-lH3Itz0MNxg-VoQh0v0s0Wub4PQ4,120
+pytetwild/__init__.py,sha256=l0Dwm58loQqSSrYqdxUx8lRHxeEBYyxr9gBvmvo-QIA,498
+pytetwild-0.1.1.dist-info/DELVEWHEEL,sha256=c7PLQUz98MRU1_k7pfUBO60FuE1AfgWCkxr52fi2ayU,421
+pytetwild-0.1.1.dist-info/METADATA,sha256=4NA0MlyTgECJengiMK3QHzM89oMFU_IPMCscS8Mg4u0,4508
+pytetwild-0.1.1.dist-info/RECORD,,
+pytetwild-0.1.1.dist-info/WHEEL,sha256=chqeLhPBtPdrOoreR34YMcofSk3yWDQhkrsDJ2n48LU,106
+pytetwild-0.1.1.dist-info/licenses/LICENSE,sha256=k0P1sbYupRWlxoESvbYq2UzoMW6MXA84LltqR-4mEo4,17096
+pytetwild.libs/concrt140-63a42d2def1023b22584817146cee33d.dll,sha256=tWTNrcqIbCkALGnMNEoKRuaSbXscOIw0F5gWYZzfmDY,304128
+pytetwild.libs/mpir-68915c4420c2a122a0aba56cd715bd26.dll,sha256=96pOmyi8-b4gAyXakkz7w0arxMoDHxmlAU33OZXMKAY,612352
+pytetwild.libs/msvcp140-8f141b4454fa78db34bc1f28c571b4da.dll,sha256=jxQbRFT6eNs0vB8oxXG02g4AzSxD960OKC8xMDaCaq4,557648

{pytetwild-0.1.dev0.dist-info → pytetwild-0.1.1.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: scikit-build-core 0.8.2
+Generator: scikit-build-core 0.11.6
 Root-Is-Purelib: false
 Tag: cp312-cp312-win_amd64
 

pytetwild-0.1.dev0.dist-info/DELVEWHEEL

@@ -1,2 +0,0 @@
-Version: 1.5.2
-Arguments: ['C:\\Users\\runneradmin\\AppData\\Local\\Temp\\cibw-run-a87uf2fd\\cp312-win_amd64\\build\\venv\\Scripts\\delvewheel', 'repair', '-w', 'C:\\Users\\runneradmin\\AppData\\Local\\Temp\\cibw-run-a87uf2fd\\cp312-win_amd64\\repaired_wheel', 'C:\\Users\\runneradmin\\AppData\\Local\\Temp\\cibw-run-a87uf2fd\\cp312-win_amd64\\built_wheel\\pytetwild-0.1.dev0-cp312-cp312-win_amd64.whl', '--add-path', '.']

pytetwild-0.1.dev0.dist-info/RECORD

@@ -1,13 +0,0 @@
-pytetwild/PyfTetWildWrapper.pyd,sha256=YPKOtIzkZ0Xb03svtxo6Sy7QqAn6hoQ8sk2NUaVcW00,4992000
-pytetwild/pytetwild.py,sha256=fA0lXIqIf97TONhpX2vME7_3vS_9kVh8ianOVIC44EU,2169
-pytetwild/_version.py,sha256=hbg9UkagWDGV6Zg23kEpex_otuUYHUzMCHXYTnFnHmE,118
-pytetwild/__init__.py,sha256=01OqAVF8WQmKyJjKPOTwgJCDAqYmP-KwbCF9VhN--xo,486
-pytetwild-0.1.dev0.dist-info/DELVEWHEEL,sha256=DIwy-gDi9myKYwjvZ8kcJOz4oKS6TGfwFn6Vk9_6g6E,423
-pytetwild-0.1.dev0.dist-info/entry_points.txt,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-pytetwild-0.1.dev0.dist-info/METADATA,sha256=tAozjdbNnB40sjP1A_7R-maND-a9aNqKpNibMwI0k3c,3038
-pytetwild-0.1.dev0.dist-info/RECORD,,
-pytetwild-0.1.dev0.dist-info/WHEEL,sha256=I4fjUu5CKXtMlETsB-F4Q7DROW-oHB8rNWT-iDmQHok,105
-pytetwild-0.1.dev0.dist-info/licenses/LICENSE,sha256=k0P1sbYupRWlxoESvbYq2UzoMW6MXA84LltqR-4mEo4,17096
-pytetwild.libs/concrt140-ddbfd8b75976783430ced289c6b989c3.dll,sha256=BZPGv1uz0YBgBOgTWBBL0elt7C3clmR348_2Zlqh47Q,302080
-pytetwild.libs/mpir-0f58ffaa159ebf7efb9359a55c2d7f23.dll,sha256=qd6vEIL0JPtRZ39EzelGE9rApfJ6duoMgV8jT8UHGns,612352
-pytetwild.libs/msvcp140-3499a3a8427eec4915749ee9a8991ddd.dll,sha256=oDDcLf0uyiipN1ySmJrfTa8WH5iNteFrnhBnjrDf9Mc,573008