In the past few decades, Delaunay-based mesh generation algorithms for 3D domains have also been implemented in various software packages and toolboxes, such as the functions delaunayTriangulaion, delaunay and distmesh/distmeshnd [1]} in MATLAB, QualMesh [2]} and the CGAL library [3]}, [4]} for C++, TetGen[5]} and iso2mesh [6]} built on CGAL as well as TetGen for MATLAB/Octave. Note that even though all the algorithms can be used to generate triangular meshes, the implementation details and usages differ. For example, the functions distmesh and distmeshnd take distance functions as representations of boundary geometry, while delaunayTriangulation generates Delaunay 2D/3D triangulation based on given point sets.
In practice, when the mathematical representations of boundary geometry or the vertex sets for Delaunay triangulations are not immediately available, one can employ TetGen instead, whose inputs are triangular surface mesh on domain boundaries, as shown in Figure REF .
