[vtkusers] vtkDelaunay3D skips points?
Joerg Rings
Joerg.Rings at iup.uni-heidelberg.de
Thu Apr 22 08:32:51 EDT 2004
Dear VTK users,
I use VTK for a 3-D tet mesh generator. Here I first create a point cloud and
triangulate this cloud with vtkDelaunay3D.
My geometry is described with coarse premesh tetrahedra, then the points are
placed within these tetrahedra, the mesh for each tetrahedron is calculated
and then collected and merged into a single unstructued grid. This procedure
has some important advantages over triangulating the whole point cloud in one
step.
If two premesh tetrahedra touch with one facet, this results in having that
facet triangulated twice, but the triangulations have to match for collecting
and merging and should if the input points are the same.
But I discover this not to be true, and it seems like the Delaunay3D
triangulation skips some points that are in the resultung unstructured grid
data, but cannot be found as edge vertices of tetrahedra.
It's especially easy when looking at the points of a facet, as the point grid
is somewhat regular, and I can see points missing.
Now the questions are:
- Does Delaunay3D do such things, so that I can blame it on the filter, or are
some of the filter options helpful here? I've tinkered with them but without
any visible change.
- If Delaunay3D is unreliable here, is their a way to create the surface of
the premesh tets otherwise and initialize Delaunay3D with these triangles?
It's only important that the triangles match on the premesh tet surfaces, so
it wouldn't be that bad if the Delaunay3D isn't perfect inside.
Any suggestion is warmly welcome,
thanks in advance,
Joerg Rings
More information about the vtkusers
mailing list