08-Jan-2022 10:40:27 triangulation_delaunay_discrepancy_test(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Test triangulation_delaunay_discrepancy(). 08-Jan-2022 10:40:27 TRIANGULATION_DELAUNAY_DISCREPANCY MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Read a node dataset of NODE_NUM points in 2 dimensions. Read an associated triangulation file 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 points NODE_NUM = 10 Read the data in "ten3_nodes.txt". First 5 nodes: Row: 1 2 Col 1 0.000000 4.000000 2 1.000000 13.000000 3 4.000000 7.000000 4 5.000000 2.000000 5 6.000000 15.000000 Read the header of ""ten3_elements.txt". 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 MESH_BASE_ONE: The element indexing appears to be 1-based! No conversion is necessary. 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 Discrepancy (degrees) = -21.801409 Minimum angle (degrees) = 9.977713 occurred in triangle 5 Maximum angle (degrees) = 156.161260 occurred in triangle 5 TRIANGULATION_DELAUNAY_DISCREPANCY: Normal end of execution. 08-Jan-2022 10:40:27 triangulation_delaunay_discrepancy_test(): Normal end of execution. 08-Jan-2022 10:40:27