6 September 2009 2:33:23.681 PM
TRIANGULATION_DELAUNAY_DISCREPANCY
FORTRAN90 version
Read a node dataset of NODE_NUM points in 2 dimensions.
Read an associated triangulation dataset of
TRIANGLE_NUM triangles using 3 or 6 nodes.
Determine the Delaunay discrepancy, that is, the amount
by which the minimum angle in the triangulation could be
changed by a single adjustment of a pair of triangles.
If this discrepancy is negative,
then the triangulation is not a Delaunay triangulation.
If this discrepancy is 0 or essentially so,
then the triangulation is a Delaunay triangulation.
Read the header of "ten3_nodes.txt".
Spatial dimension DIM_NUM = 2
Number of nodes NODE_NUM = 10
Read the data in "ten3_nodes.txt".
Coordinates of first 5 nodes:
Row 1 2
Col
1 0.00000 4.00000
2 1.00000 13.0000
3 4.00000 7.00000
4 5.00000 2.00000
5 6.00000 15.0000
Read the header of "ten3_triangles.txt".
Triangle order TRIANGLE_ORDER = 3
Number of triangles TRIANGLE_NUM = 10
Read the data in "ten3_triangles.txt".
First 5 triangles:
Row 1 2 3
Col
1 1 4 7
2 1 7 3
3 3 7 9
4 1 3 2
5 3 10 6
First 5 triangle neighbors:
Row 1 2 3
Col
1 -1 2 -1
2 3 4 1
3 -1 7 2
4 6 -1 2
5 8 6 7
Delaunay discrepancy = -21.8014
Minimum angle (degrees) = 9.97771
occurred in triangle 5
Maximum angle (degrees) = 156.161
occurred in triangle 5
TRIANGULATION_DELAUNAY_DISCREPANCY:
Normal end of execution.
6 September 2009 2:33:23.687 PM