13 September 2009 10:10:06.395 AM 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_elements.txt". Triangle order TRIANGLE_ORDER = 3 Number of triangles TRIANGLE_NUM = 10 Read the data in "ten3_elements.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. 13 September 2009 10:10:06.401 AM