Task:
Given the boundary of a 3-D solid, and a mesh-size function, generate
a triangulation by distributing a set of nodes inside and on the boundary,
and connecting them into tetrahedra.
Goals:
Robustness.
Good scalability.
Support for fragmentation simulations (decohesion
surfaces, crack growth).