triangle_shewchuk, a C code which generates meshes, Delaunay triangulations and Voronoi diagrams for 2D pointsets, by Jonathan Shewchuk.
TRIANGLE generates exact Delaunay triangulations, constrained Delaunay triangulations, and quality conforming Delaunay triangulations. The latter can be generated with no small angles, and are thus suitable for finite element analysis.
TRIANGLE produces its own documentation. Complete instructions are printed by invoking the program with the `-h' switch:
triangle -h
The instructions are long; you'll probably want to pipe the output to more or redirect it to a file. The programs gives a short list of command line options if it is invoked without arguments.
You have to give TRIANGLE some input to work with. Usually, you start with a Planar Straight Line Graph (PSLG) which defines a set of points, and some nonintersecting segments that you want to be sure are included in the triangulation. This information can easily be stored in a POLY, which is one of the expected input formats.
Normally, TRIANGLE will triangulate the entire region defined by the convex hull of the points. You may define a portion of the convex hull to be a "hole", which is not to be triangulated. This is done by adding some hole information to your POLY file.
Try out TRIANGLE on the sample file, A.poly:
triangle -p A.poly
TRIANGLE will read the Planar Straight Line Graph defined by A.poly, and write its constrained Delaunay triangulation to A.1.node, A.1.ele and A.1.poly.
For contrast, try running
triangle -pq A.polyNow, click on the same "ele" button. A new triangulation will be loaded; this one having no angles smaller than 20 degrees.
To see a Voronoi diagram, try this:
cp A.poly A.node triangle -v A
triangle is available in a C version.
fem_to_triangle, a C code which reads FEM files defining a 2D mesh of triangles, namely a file of node coordinates and a file of elements defined by node indices, and creates a corresponding pair of node and element files for use by Jonathan Shewchuk's triangle() program.
POLY, a data directory which contains a description and examples of the POLY file format.
SHOWME, a C code which can display the POLY files uses as input to Jonathan Shewchuk's triangle() program, and the output files that define meshes and other objects.
TRIANGLE_BENCH, a script which times the execution of Jonathan Shewchuk's triangle() program on a sequence of sets of randomly generated nodes in the unit square.
TRIANGLE_FILES, a data directory of examples of files used by the Jonathan Shewchuk's triangle() and showme() programs.
TRIANGULATE, a C code which reads data defining a polygon and produces a PostScript file describing a triangulation of the polygon.
TRIANGULATION, a C code which performs various operations on order 3 (linear) or order 6 (quadratic) triangulations.
TRIANGULATION_TRIANGLE_NEIGHBORS, a C code which reads data defining a triangulation, determines the neighboring triangles of each triangle, and writes that information to a file.
Jonathan Shewchuk,
Computer Science Division,
University of California at Berkeley,
Berkeley, California, 94720-1776
jrs@cs.berkeley.edu.
A is a planar straight line graph of the capital letter A. We use it as input to get a constrained Delaunay triangulation.
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.
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.
triangle -p double_hex.polyrequests that we triangulate the current points:
triangle -pqY -a0.0015 double_hex.1.polyrequests 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.
triangle -p double_hex2.polyrequests that we triangulate the current points:
triangle -pqY -a0.0060 double_hex2.1.polyrequests 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:
SQUARE_CIRCLE_HOLE is a planar straight line graph of a square region with an off center circular hole, and 826 points computed by a CVT calculation, prepared by Hua Fei.