TRIANGULATION_Q2L
6-Node Triangulation to 3-Node Triangulation

TRIANGULATION_Q2L is a FORTRAN90 program which reads information describing a triangulation of a set of points using 6-node ("quadratic") triangles, and creates a 3-node ("linear") triangulation.

The same nodes are used, but each 6-node triangle is broken up into four smaller 3-node triangles.

Thus, the program might be given the following 4 triangles:

        11-12-13-14-15
         |\    |\    |
         | \   | \   |
         6  7  8  9 10
         |   \ |   \ |
         |    \|    \|
         1--2--3--4--5
      

It would make a new triangulation involving 16 triangles:

        11-12-13-14-15
         |\ |\ |\ |\ |
         | \| \| \| \|
         6--7--8--9-10
         |\ |\ |\ |\ |
         | \| \| \| \|
         1--2--3--4--5
      

The input and output files use the simple TABLE format; comment lines begin with a "#" character. Otherwise, each line of the file contains one set of information, either the coordinates of a node (for a node file), or the indices of nodes that make up a triangle, (for a triangle file).

The input file prefix_elements.txt contains the triangle information for the 6-node triangulation. Each line contains the indices of six nodes that form a triangle, in counterclockwise order. The first three indices are the vertices, in counterclockwise order. The fourth index is the midside node between vertices 1 and 2, the fifth the midside between vertices 2 and 3, and the sixth the midside between vertices 3 and 1.

The output file prefix_q2l_elements.txt contains the triangle information for the 3-node triangulation. The vertices for each triangle are listed in counterclockwise order. There are 4 times as many triangles in this triangulation.

Usage:

triangulation_q2l prefix

• prefix_nodes.txt contains the node coordinates (not needed by this program),
• prefix_elements.txt contains the quadratic element definitions.
• prefix_q2l_elements.txt will contain the linear element definitions, based on the same set of nodes.

Licensing:

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

Languages:

TRIANGULATION_Q2L is available in a C++ version and a FORTRAN90 version and a MATLAB version.

Related Data and Programs:

MESH_TO_XML, a FORTRAN90 program which reads information defining a 1D, 2D or 3D mesh, namely a file of node coordinates and a file of elements defined by node indices, and creates a corresponding XML file for input to DOLFIN or FENICS.

TABLE_DELAUNAY, a FORTRAN90 program which triangulates a set of nodes whose coordinates are stored in a file.

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

TRIANGULATION, a FORTRAN90 library which carries out various operations on order 3 ("linear") or order 6 ("quadratic") triangulations.

TRIANGULATION_BOUNDARY_NODES, a FORTRAN90 program which reads data defining a triangulation, determines which nodes lie on the boundary, and writes their coordinates to a file.

TRIANGULATION_CORNER, a FORTRAN90 program which patches triangulations so that no triangle has two sides on the boundary.

TRIANGULATION_DELAUNAY_DISCREPANCY, a FORTRAN90 program which measures the amount by which a triangulation fails the local Delaunay test;

TRIANGULATION_HISTOGRAM, a FORTRAN90 program which computes histograms of data over a triangulation.

TRIANGULATION_L2Q, a FORTRAN90 program which reads data defining a 3-node triangulation and generates midside nodes and writes out the corresponding 6-node triangulation.

TRIANGULATION_MASK, a FORTRAN90 program which takes an existing triangulation and deletes triangles and their corresponding nodes as requested by the user.

TRIANGULATION_NODE_TO_ELEMENT, a FORTRAN90 program which reads files describing a set of nodes, their triangulation, and the value of one or more quantities at each node, and outputs a file that averages the quantities for each element. This operation in effect creates an "order1" finite element model of the data.

TRIANGULATION_ORDER3, a data directory which contains a description and examples of order 3 triangulations.

TRIANGULATION_ORDER6, a data directory which contains a description and examples of order 6 triangulations.

TRIANGULATION_ORIENT, a FORTRAN90 program which reads data defining a triangulation, makes sure that every triangle has positive orientation, and if not, writes a corrected triangle file.

TRIANGULATION_PLOT, a FORTRAN90 program which reads data defining a triangulation and creates a PostScript image of the nodes and triangles.

TRIANGULATION_QUAD, a FORTRAN90 program which estimates the integral of a function over a triangulated region.

TRIANGULATION_QUALITY, a FORTRAN90 program which reads data defining a triangulation and computes a number of quality measures.

TRIANGULATION_RCM, a FORTRAN90 program which reads data defining a triangulation, determines an ordering of the nodes that will reduce the bandwidth of the adjacency matrix, and writes the new triangulation information to a file.

TRIANGULATION_REFINE, a FORTRAN90 program which reads data defining a triangulation, replaces each triangle by four congruent smaller ones, and writes the new triangulation information to a file.

TRIANGULATION_TRIANGLE_NEIGHBORS, a FORTRAN90 program which reads data defining a triangulation, determines the neighboring triangles of each triangle, and writes that information to a file.

Reference:

1. Marc deBerg, Marc Krevald, Mark Overmars, Otfried Schwarzkopf,
   Computational Geometry,
   Springer, 2000.
2. Joseph ORourke,
   Computational Geometry,
   Cambridge University Press,
   Second Edition, 1998.

Source Code:

• triangulation_q2l.f90, the source code.

Examples and Tests:

EXAMPLE is a set of demonstration data.

• example_nodes.txt, the nodes.
• example_elements.txt, the triangles for the 6-node triangulation.
• example.png, a PNG image of the 6-node triangulation, created by TRIANGULATION_PLOT.
• example_q2l_elements.txt, the triangles for the 3-node triangulation.
• example_q2l.png, a PNG image of the 3-node triangulation, created by TRIANGULATION_PLOT.

List of Routines:

• MAIN is the main program for TRIANGULATION_Q2L.
• CH_CAP capitalizes a single character.
• CH_EQI is a case insensitive comparison of two characters for equality.
• CH_TO_DIGIT returns the integer value of a base 10 digit.
• FILE_COLUMN_COUNT counts the number of columns in the first line of a file.
• FILE_NAME_EXT_GET determines the "extension" of a file name.
• FILE_NAME_EXT_SWAP replaces the current "extension" of a file name.
• FILE_ROW_COUNT counts the number of row records in a file.
• GET_UNIT returns a free FORTRAN unit number.
• I4_MODP returns the nonnegative remainder of integer division.
• I4_WRAP forces an integer to lie between given limits by wrapping.
• I4MAT_TRANSPOSE_PRINT prints an I4MAT, transposed.
• I4MAT_TRANSPOSE_PRINT_SOME prints some of the transpose of an I4MAT.
• ITABLE_DATA_READ reads data from an integer table file.
• ITABLE_DATA_WRITE writes data to an integer table file.
• ITABLE_HEADER_READ reads the header from an integer table file.
• ITABLE_HEADER_WRITE writes the header to an integer table file.
• ITABLE_WRITE writes an integer table file.
• S_BLANK_DELETE removes blanks from a string, left justifying the remainder.
• S_INDEX_LAST finds the LAST occurrence of a given substring.
• S_TO_I4 reads an I4 from a string.
• S_TO_I4VEC reads an I4VEC from a string.
• S_WORD_COUNT counts the number of "words" in a string.
• TIMESTAMP prints the current YMDHMS date as a time stamp.
• TIMESTRING writes the current YMDHMS date into a string.
• TRIANGULATION_ORDER6_TO_ORDER3 linearizes a quadratic triangulation.

You can go up one level to the FORTRAN90 source codes.