From the thesis of Dafna Talmor, Well-Spaced Points for Numerical Methods (1997):
It is reasonable to conjecture that every well-spaced point set can be tetrahedralized to form a good aspect ratio mesh. However, the Delaunay tetrahedralization fails to oblige us in that respect. Currently, it is unknown if such a tetrahedralization algorithm exists, or if there exists a counter-example: a well-spaced point set with the property that none of its tetrahedralizations is of bounded aspect ratio.
This is probably hard to understand for most of you due to lack of familiarity with the subject. Let me try to provide a translation:
1. There are two kinds of point sets, namely those that are well-spaced and those that are not. Intuitively well-spacedness means that the points are distributed evenly over space. No two points are too close together and there are no gaps. In a crate of oranges the points formed by the centres of all oranges are well-spaced.
2. Delaunay is a criterion for completely triangulating a point set. Applying Delaunay we get a triangulation of any collection of points. When we apply Delaunay to a two-dimensional point set (points in a plane) then, in the case of well-spaced points, we get a triangulation that consists only of ‘voluminous’ triangles, that means triangles that are not arbitrarily flat. In three dimensions, however, a Delaunay triangulation of well-spaced points can contain highly flat tetrahedra (three-dimensional triangles; a pyramid with a triangle instead of a square as bottom face).
3. The Delaunay triangulation is not the only triangulation that is possible. In fact, the possible number of triangulations for a point set is usually very, very large. In the Delaunay triangulation of well-spaced points in 3D only relatively few tetrahedra are flat and with some manipulation they can often be removed by slightly changing the connectivity of the points. If that does not work then points can be moved around a little or even points can be added to remove the flat tetrahedron.
4. What Talmor remarks, and what I recently came across myself, is that is unknown if it is always possible to remove the flat tetrahedra from a Delaunay triangulation of well-spaced points by changing the connectivity. The huge number of possible triangulations makes it unfeasible to try them all. There are two possibilities:
- For some well-spaced point sets there is just no triangulation possible that does not contain a flat tetrahedron.
- Every well-spaced point set has at least one triangulation without flat tetrahedra. Talmor says that this option is ‘reasonable to conjecture‘ and I concur. However, if such a triangulation is always possible, that leaves us hanging with one big and obvious question:
Is there are algorithm that produces a triangulation without flat tetrahedra, without trying all possible combinations of triangulations?
That sure is a vexing issue.
I imagine that at this point some of you still do not understand the gist of the question. If I find time, I will add some pictures to illustrate some of the concepts.
My interest in this problem is not to solve the issue all together. I’m just interested in creating good triangulations, without flat tetrahedra that is. What I would like to have is an algorithm that for a well-spaced point set produces the triangulation in which the flattest tetrahedron, is more voluminous than the flattest tetrahedra of all other triangulations. Otherwise fomulated: if we measure the quality of a triangulation by the flattest tetrahedron in contains, then we want to find the triangulation with the highest quality that is possible for a particular point set.
It is obvious that Talmor’s problem and my problem relate: assuming that a triangulation without flat tetrahedra always exists, then the optimal triangulation would surely offer an example of such a triangulation.
update: I recently discovered that in the paper Sliver Exudation (2000) it is proven that a triangulation without slivers does indeed exist, if the points have a ‘nice enough’ distribution to begin with. I’m not sure whether their sliver exudation algorithm can actually find the (approximate) optimal triangulation. At least it can find a triangulation without slivers, though the theoretical guarantueed quality is very small.