TRIANGLE_FILES
Files Used by the TRIANGLE Program


TRIANGLE_FILES is a data directory which contains some examples of the groups of files used by the TRIANGLE and SHOWME programs to produce triangulations and Voronoi decompositions of a region.

Licensing:

The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.

Related Data and Programs:

POLY, a data directory which contains examples of POLY files, a 2D format used by triangle and showme for Planar Straight Line Graphs.

SHOWME, a C program which uses the X Window library to display triangulatons, by Jonathan Shewchuk.

TRIANGLE, a C program which computes a triangulation of a geometric region, by Jonathan Shewchuk.

TRIANGLE_DISPLAY, a MATLAB program which displays the nodes and elements of a triangulation on the MATLAB graphics screen, assuming the data has been stored in NODE and ELE files by the TRIANGLE program.

Reference:

  1. Jonathan Shewchuk,
    Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator,
    in Applied Computational Geometry: Towards Geometric Engineering,
    edited by Ming Lin, Dinesh Manocha,
    Lecture Notes in Computer Science, Volume 1148,
    Springer, 1996,
    ISBN: 354061785X,
    LC: QA448.D38.A635.
  2. Jonathan Shewchuk,
    Delaunay Refinement Algorithms for Triangular Mesh Generation,
    Computational Geometry, Theory and Applications,
    Volume 23, pages 21-74, May 2002.
  3. http://www-2.cs.cmu.edu/~quake/triangle.html, the TRIANGLE web site;

Example Files:

A is a planar straight line graph of the capital letter A. We use it as input to get a constrained Delaunay triangulation.

AIRFOIL_EXTERIOR defines the boundary of an airfoil, and a box. We want to mesh the region inside the box and exterior to the airfoil. We use the command

        triangle -a0.05 airfoil_exterior.poly
      

AIRFOIL_INTERIOR defines the boundary of an airfoil, whose interior we want to mesh. We use the command

        triangle -a0.05 airfoil_interior.poly
      

BOX is a planar straight line graph of a double box. We use it as input to get a constrained Delaunay triangulation.

The command "triangle -pc box" creates the first mesh. (We need the -pc because we didn't surround our initial area with segments.) The command "triangle -rpa0.2 box.1" refines the box.1 mesh with an area constraint of 0.2: The command "triangle -rpa0.05 box.2" refines the box.2 mesh with an area constraint of 0.05, 1/4 of the previous maximum area: The command "triangle -rpa0.0125 box.3" refines the box.3 mesh with an area constraint of 0.0125, 1/4 of the previous maximum area:

DIAMOND_02_00009 is another set of test data, for which we want the Voronoi diagram.

DOUBLE_HEX describes a unit square with two hexagonal holes. 72 points are listed on the outer boundary, and 12 on each of the holes. It is desired to create a nice looking mesh of about 500 nodes, and no additional nodes on the boundary segments.

Our first command
triangle -p double_hex.poly
requests that we triangulate the current points: Our second command
triangle -pqY -a0.0015 double_hex.1.poly
requests that we triangulate the current points, adding new nodes as necessary to make a nice mesh, with no triangle being larger than 0.0015 in area, and with no points added on boundary segments. We end up with 525 nodes and 956 elements:

DOUBLE_HEX2 describes a unit square with two hexagonal holes. 36 points are listed on the outer boundary, and 6 on each of the holes. It is desired to create a nice looking mesh of about 235 elements, and no additional nodes on the boundary segments.

Our first command
triangle -p double_hex2.poly
requests that we triangulate the current points: Our second command
triangle -pqY -a0.0060 double_hex2.1.poly
requests that we triangulate the current points, adding new nodes as necessary to make a nice mesh, with no triangle being larger than 0.0060 in area, and with no points added on boundary segments. We end up with 141 nodes and 236 elements:

DOUBLE_HEX3 describes a unit square with two hexagonal holes. 4 points are listed on the outer boundary, and 6 on each of the holes. We want triangle to triangulate this region.

Our command
triangle -pq -a0.0015 double_hex3.poly
requests that we triangulate the region, adding points as necessary so that no triangle has an area greater than 0.0015.

ELL is a triangulation of an L-shaped region, using a mesh of 21 nodes and 24 elements.

GREENLAND is a triangulation of Greenland, using a graded (varying-size) mesh of 33,343 nodes and 64,125 elements. The resulting image is essentially a red blob the shape of Greenland. However, using the ZOOM tool, it is possible to see interesting details of the mesh.

LA is a POLY file containing information representing soil layers. The data includes points that are bounded by line segments defining the different layers. The intent is that certain layers will be triangulated with smaller area requirements.

SPIRAL is a node file containing points that form a spiral.

The command "triangle spiral" produces a Delaunay triangulation of the points, in the following node and element files: