Function tetrahedralize

python/pymesh/tetrahedralize.py:15–186

Create a tetrahedral mesh from input triangle mesh. Arguments: mesh (:class:`Mesh`): Input triangular mesh. cell_size (``float``): Max radius of the circumscribed sphere of the output tet. radius_edge_ratio (``float``): Max radius of the circumscribed sphere to the

(mesh,
        cell_size,
        radius_edge_ratio=2.0,
        facet_distance=-1.0,
        feature_angle=120,
        engine="auto",
        with_timing=False)

from time import time

def tetrahedralize(mesh,
        cell_size,
        radius_edge_ratio=2.0,
        facet_distance=-1.0,
        feature_angle=120,
        engine="auto",
        with_timing=False):
    """ Create a tetrahedral mesh from input triangle mesh.

    Arguments:
        mesh (:class:`Mesh`): Input triangular mesh.
        cell_size (``float``): Max radius of the circumscribed sphere of the output tet.
        radius_edge_ratio (``float``): Max radius of the circumscribed sphere to the
            shortest edge length of each tet.
        facet_distance (``float``): Upper bound on the distance from the
            circumcenter of a facet to the center of its "Delaunay ball", where
            a Delaunay ball is defined as the smallest circumscribed sphere
            with center on the surface of the domain.
        feature_angle (``float``): Angle threshold (in degrees) for feature extraction.
        engine (``string``): The tetrahedralization engine to use.  Valid options are:

            * ``auto``: default to tetgen
            * ``cgal``: `CGAL 3D mesh generation`_, using Polyhedron domain with
              auto feature extraction.
            * ``cgal_no_features``: `CGAL 3D mesh generation`_,
              using Polyhedron domain without feature extraction.
            * ``cgal_implicit``: `CGAL 3D mesh generation`_,
              using implicit domain with winding number as oracle.
            * ``tetgen``: `TetGen`_ from Hang Si.
            * ``quartet``: `Quartet`_ from Robert Bridson and Crawford Doran
            * ``delpsc``: `DelPSC`_ from Tamal K Dey , Joshua A. Levine, Andrew
              Slatton
            * ``vegafem``: Tet mesher provided by `VegaFEM`_ library.
            * ``mmg``: Implicit domain meshing from `MMG3D`_.
            * ``tetwild``: `TetWild`_ engine based on our Siggraph paper.
        with_timing (``boolean``): whether to output timing info.

    Returns:
        Tetrahedral mesh (and running time if `with_timing` is True).

    .. _`CGAL 3D mesh generation`: https://doc.cgal.org/latest/Mesh_3/index.html
    .. _`TetGen`: http://wias-berlin.de/software/tetgen/
    .. _`Quartet`: https://github.com/crawforddoran/quartet
    .. _`DelPSC`: http://web.cse.ohio-state.edu/~dey.8/delpsc.html
    .. _`VegaFEM`: http://run.usc.edu/vega/
    .. _`MMG3D`: https://www.mmgtools.org/
    .. _`TetWild`: https://github.com/Yixin-Hu/TetWild
    """
    logger = logging.getLogger(__name__)
    bbox_min, bbox_max = mesh.bbox
    bbox_diagonal = numpy.linalg.norm(bbox_max - bbox_min)
    if cell_size <= 0.0:
        cell_size = bbox_diagonal / 20.0 * math.sqrt(6) / 4.0
    logger.info("Cell size: {}".format(cell_size))
    if facet_distance <= 0.0:
        facet_distance = bbox_diagonal / 1000.0
    logger.info("Facet distance: {}".format(facet_distance))
    if radius_edge_ratio <= 0.0: