TetGen example: Crystal


Return to TetGen examples page.
Return to TetGen home page.
The input surface mesh, 11706 nodes, 26399 triangles.


The constrained Delaunay tetrahedralization, 22271 nodes, 94569 tetrahedra, and 52676 boundary faces.

The quality Delaunay tetrahedralization (Statistics).

Runtime output and mesh quality statistics

si@a302n3: ~/meshmodels/crystal > tetgen -pqV crystal03_cut.d_cut.mesh 
Opening crystal03_cut.d_cut.mesh.
Constructing Delaunay tetrahedrization.
  Total flips: 209675, where T23 123795, T32 85854, T22 0, T44 26
Delaunay seconds:  2.21
Creating surface mesh.
  Constructing mapping from indices to points.
  Unifying segments.
  Constructing mapping from points to subfaces.
  Merging coplanar facets.
Perturbing vertices.
  Constructing mapping from points to tetrahedra.
  0 break points are inserted.
Delaunizing segments.
  Marking segments have acute corners.
  Constructing mapping from points to segments.
  Constructing mapping from points to tetrahedra.
  10738 protect points are inserted.
Constraining facets.
  Constructing mapping from points to tetrahedra.
  The biggest cavity: 132 faces, 68 vertices
  Enlarged 0 times
Segment and facet seconds:  5.54
Removing unwanted tetrahedra.
  Marking concavities for elimination.
  Marking neighbors of marked tetrahedra.
  Deleting marked tetrahedra.
Hole seconds:  0.07
Removing degenerate tets.
  0 degeneracies are removed.
Repair seconds:  0.06
Adding Steiner points to enforce quality.
  Marking sharp subfaces.
  30088 subfaces have been marked.
  Marking skinny subfaces.
  16812 subfaces have been marked.
  Initializing protecting spheres.
  Constructing mapping from points to tetrahedra.
  Splitting encroached subsegments.
  49904 points are inserted.
  Splitting encroached subfaces.
  35239 points are inserted.
  Splitting bad tetrahedra.
  52329 points are inserted.
  Removing slivers.
  3181 slivers are removed.
  87 sliver(s) survived.
  17 points are inserted.
  Total 137489 points are inserted.
Quality seconds:  98.17
Jettisoning redundants points.

Writing crystal03_cut.d_cut.1.node.
Writing crystal03_cut.d_cut.1.ele.
Writing crystal03_cut.d_cut.1.face.

Output seconds:  1.86
Total running seconds:  107.91

Statistics:

  Input points: 11706
  Input facets: 26399
  Input holes: 0
  Input regions: 0

  Mesh points: 159931
  Mesh tetrahedra: 709938
  Mesh faces: 1507632
  Mesh subfaces: 261850
  Mesh subsegments: 118591

Mesh quality statistics:

  Smallest volume:        3.546e-29   |  Largest volume:       1.0562e-05
  Shortest edge:          4.732e-10   |  Longest edge:           0.064647
  Smallest dihedral:     7.0331e-05   |  Largest dihedral:       179.9430

  Radius-edge ratio histogram:
         < 0.707     :      5897      |    1.6 - 1.8        :     53539
   0.707 - 1         :    182032      |    1.8 - 2          :     28981
       1 - 1.1       :     81138      |      2 - 2.5        :     20617
     1.1 - 1.2       :     76580      |    2.5 - 3          :     11043
     1.2 - 1.4       :    129627      |      3 - 10         :     25776
     1.4 - 1.6       :     89128      |     10 -            :      5580
  (A tetrahedron's radius-edge ratio is its radius of circumsphere divided
    by its shortest edge length)


  Aspect ratio histogram:
  1.1547 - 1.5       :         0      |     15 - 25         :     15819
     1.5 - 2         :         0      |     25 - 50         :      5901
       2 - 2.5       :         9      |     50 - 100        :      1696
     2.5 - 3         :     13621      |    100 - 300        :       799
       3 - 4         :    155255      |    300 - 1000       :       154
       4 - 6         :    350612      |   1000 - 10000      :        23
       6 - 10        :    135883      |  10000 - 100000     :        15
      10 - 15        :     30130      | 100000 -            :        21
  (A tetrahedron's aspect ratio is its longest edge length divided by the
    diameter of its inscribed sphere)

  Dihedral Angle histogram:
       0 - 10 degrees:     28391      |     90 - 100 degrees:    131439
      10 - 20 degrees:     77962      |    100 - 110 degrees:    159962
      20 - 30 degrees:    168328      |    110 - 120 degrees:    136148
      30 - 40 degrees:    196923      |    120 - 130 degrees:     88043
      40 - 50 degrees:    152600      |    130 - 140 degrees:     54006
      50 - 60 degrees:     76832      |    140 - 150 degrees:     33984
      60 - 70 degrees:      8901      |    150 - 160 degrees:     21704
      70 - 80 degrees:      4324      |    160 - 170 degrees:     12771
      80 - 90 degrees:     66164      |    170 - 180 degrees:      1394

Return to TetGen examples page.
Return to TetGen home page.
Hang Si