11 May 2015 7:28:42.540 AM POLYGON_GRID_TEST: FORTRAN90 version Test the POLYGON_GRID library. POLYGON_GRID_COUNT_TEST: POLYGON_GRID_COUNT counts NG, the number of points in a grid defined with N+1 points on each side of a polygon of NV vertices. Polygonal vertex count NV = 3 N NG 0 1 1 4 2 10 3 19 4 31 5 46 Polygonal vertex count NV = 4 N NG 0 1 1 5 2 13 3 25 4 41 5 61 Polygonal vertex count NV = 5 N NG 0 1 1 6 2 16 3 31 4 51 5 76 POLYGON_GRID_POINTS_TEST01: POLYGON_GRID_POINTS returns grid points for a polygon of NV vertices, with N+1 points on a side For this test, the polygon is a triangle. Polygon vertices: Row 1 2 Col 1: 0.00000 0.00000 2: 1.00000 0.00000 3: 0.500000 0.866025 N = 5 Number of grid points will be NG = 46 The grid point array: Row 1 2 Col 1: 0.500000 0.288675 2: 0.400000 0.230940 3: 0.500000 0.173205 4: 0.600000 0.115470 5: 0.700000 0.577350E-01 6: 0.800000 0.00000 7: 0.300000 0.173205 8: 0.400000 0.115470 9: 0.500000 0.577350E-01 10: 0.600000 0.00000 11: 0.200000 0.115470 12: 0.300000 0.577350E-01 13: 0.400000 0.00000 14: 0.100000 0.577350E-01 15: 0.200000 0.00000 16: 0.00000 0.00000 17: 0.600000 0.230940 18: 0.600000 0.346410 19: 0.600000 0.461880 20: 0.600000 0.577350 21: 0.600000 0.692820 22: 0.700000 0.173205 23: 0.700000 0.288675 24: 0.700000 0.404145 25: 0.700000 0.519615 26: 0.800000 0.115470 27: 0.800000 0.230940 28: 0.800000 0.346410 29: 0.900000 0.577350E-01 30: 0.900000 0.173205 31: 1.00000 0.00000 32: 0.500000 0.404145 33: 0.400000 0.346410 34: 0.300000 0.288675 35: 0.200000 0.230940 36: 0.100000 0.173205 37: 0.500000 0.519615 38: 0.400000 0.461880 39: 0.300000 0.404145 40: 0.200000 0.346410 41: 0.500000 0.635085 42: 0.400000 0.577350 43: 0.300000 0.519615 44: 0.500000 0.750555 45: 0.400000 0.692820 46: 0.500000 0.866025 Created vertex file "triangle_vertex.txt". Created grid file "triangle_grid.txt". Created command file "triangle_commands.txt". Data written to the file "triangle.xy" POLYGON_GRID_POINTS_TEST02: POLYGON_GRID_POINTS returns grid points for a polygon of NV vertices, with N+1 points on a side For this test, the polygon is a convex quadrilateral with sides of varying length. Polygon vertices: Row 1 2 Col 1: 1.00000 1.00000 2: 2.00000 0.00000 3: 4.00000 3.00000 4: 0.00000 5.00000 N = 7 Number of grid points will be NG = 113 The grid point array: Row 1 2 Col 1: 1.75000 2.25000 2: 1.64286 2.07143 3: 1.67857 1.75000 4: 1.71429 1.42857 5: 1.75000 1.10714 6: 1.78571 0.785714 7: 1.82143 0.464286 8: 1.85714 0.142857 9: 1.53571 1.89286 10: 1.57143 1.57143 11: 1.60714 1.25000 12: 1.64286 0.928571 13: 1.67857 0.607143 14: 1.71429 0.285714 15: 1.42857 1.71429 16: 1.46429 1.39286 17: 1.50000 1.07143 18: 1.53571 0.750000 19: 1.57143 0.428571 20: 1.32143 1.53571 21: 1.35714 1.21429 22: 1.39286 0.892857 23: 1.42857 0.571429 24: 1.21429 1.35714 25: 1.25000 1.03571 26: 1.28571 0.714286 27: 1.10714 1.17857 28: 1.14286 0.857143 29: 1.00000 1.00000 30: 1.78571 1.92857 31: 2.10714 2.03571 32: 2.42857 2.14286 33: 2.75000 2.25000 34: 3.07143 2.35714 35: 3.39286 2.46429 36: 3.71429 2.57143 37: 1.82143 1.60714 38: 2.14286 1.71429 39: 2.46429 1.82143 40: 2.78571 1.92857 41: 3.10714 2.03571 42: 3.42857 2.14286 43: 1.85714 1.28571 44: 2.17857 1.39286 45: 2.50000 1.50000 46: 2.82143 1.60714 47: 3.14286 1.71429 48: 1.89286 0.964286 49: 2.21429 1.07143 50: 2.53571 1.17857 51: 2.85714 1.28571 52: 1.92857 0.642857 53: 2.25000 0.750000 54: 2.57143 0.857143 55: 1.96429 0.321429 56: 2.28571 0.428571 57: 2.00000 0.00000 58: 2.07143 2.35714 59: 1.82143 2.75000 60: 1.57143 3.14286 61: 1.32143 3.53571 62: 1.07143 3.92857 63: 0.821429 4.32143 64: 0.571429 4.71429 65: 2.39286 2.46429 66: 2.14286 2.85714 67: 1.89286 3.25000 68: 1.64286 3.64286 69: 1.39286 4.03571 70: 1.14286 4.42857 71: 2.71429 2.57143 72: 2.46429 2.96429 73: 2.21429 3.35714 74: 1.96429 3.75000 75: 1.71429 4.14286 76: 3.03571 2.67857 77: 2.78571 3.07143 78: 2.53571 3.46429 79: 2.28571 3.85714 80: 3.35714 2.78571 81: 3.10714 3.17857 82: 2.85714 3.57143 83: 3.67857 2.89286 84: 3.42857 3.28571 85: 4.00000 3.00000 86: 1.50000 2.64286 87: 1.39286 2.46429 88: 1.28571 2.28571 89: 1.17857 2.10714 90: 1.07143 1.92857 91: 0.964286 1.75000 92: 0.857143 1.57143 93: 1.25000 3.03571 94: 1.14286 2.85714 95: 1.03571 2.67857 96: 0.928571 2.50000 97: 0.821429 2.32143 98: 0.714286 2.14286 99: 1.00000 3.42857 100: 0.892857 3.25000 101: 0.785714 3.07143 102: 0.678571 2.89286 103: 0.571429 2.71429 104: 0.750000 3.82143 105: 0.642857 3.64286 106: 0.535714 3.46429 107: 0.428571 3.28571 108: 0.500000 4.21429 109: 0.392857 4.03571 110: 0.285714 3.85714 111: 0.250000 4.60714 112: 0.142857 4.42857 113: 0.00000 5.00000 Created vertex file "quad_vertex.txt". Created grid file "quad_grid.txt". Created command file "quad_commands.txt". Data written to the file "quad.xy" POLYGON_GRID_POINTS_TEST03: POLYGON_GRID_POINTS returns grid points for a polygon of NV vertices, with N+1 points on a side For this test, the polygon is nonconvex and six sided. Two degenerate triangles are created, and some grid points are generated several times. Polygon vertices: Row 1 2 Col 1: 0.00000 0.00000 2: 2.00000 0.00000 3: 2.00000 1.00000 4: 1.00000 1.00000 5: 1.00000 2.00000 6: 0.00000 2.00000 N = 5 Number of grid points will be NG = 91 The grid point array: Row 1 2 Col 1: 1.00000 1.00000 2: 0.800000 0.800000 3: 1.00000 0.600000 4: 1.20000 0.400000 5: 1.40000 0.200000 6: 1.60000 0.00000 7: 0.600000 0.600000 8: 0.800000 0.400000 9: 1.00000 0.200000 10: 1.20000 0.00000 11: 0.400000 0.400000 12: 0.600000 0.200000 13: 0.800000 0.00000 14: 0.200000 0.200000 15: 0.400000 0.00000 16: 0.00000 0.00000 17: 1.20000 0.800000 18: 1.40000 0.800000 19: 1.60000 0.800000 20: 1.80000 0.800000 21: 2.00000 0.800000 22: 1.40000 0.600000 23: 1.60000 0.600000 24: 1.80000 0.600000 25: 2.00000 0.600000 26: 1.60000 0.400000 27: 1.80000 0.400000 28: 2.00000 0.400000 29: 1.80000 0.200000 30: 2.00000 0.200000 31: 2.00000 0.00000 32: 1.20000 1.00000 33: 1.20000 1.00000 34: 1.20000 1.00000 35: 1.20000 1.00000 36: 1.20000 1.00000 37: 1.40000 1.00000 38: 1.40000 1.00000 39: 1.40000 1.00000 40: 1.40000 1.00000 41: 1.60000 1.00000 42: 1.60000 1.00000 43: 1.60000 1.00000 44: 1.80000 1.00000 45: 1.80000 1.00000 46: 2.00000 1.00000 47: 1.00000 1.00000 48: 1.00000 1.20000 49: 1.00000 1.40000 50: 1.00000 1.60000 51: 1.00000 1.80000 52: 1.00000 1.00000 53: 1.00000 1.20000 54: 1.00000 1.40000 55: 1.00000 1.60000 56: 1.00000 1.00000 57: 1.00000 1.20000 58: 1.00000 1.40000 59: 1.00000 1.00000 60: 1.00000 1.20000 61: 1.00000 1.00000 62: 1.00000 1.20000 63: 0.800000 1.40000 64: 0.600000 1.60000 65: 0.400000 1.80000 66: 0.200000 2.00000 67: 1.00000 1.40000 68: 0.800000 1.60000 69: 0.600000 1.80000 70: 0.400000 2.00000 71: 1.00000 1.60000 72: 0.800000 1.80000 73: 0.600000 2.00000 74: 1.00000 1.80000 75: 0.800000 2.00000 76: 1.00000 2.00000 77: 0.800000 1.20000 78: 0.600000 1.00000 79: 0.400000 0.800000 80: 0.200000 0.600000 81: 0.00000 0.400000 82: 0.600000 1.40000 83: 0.400000 1.20000 84: 0.200000 1.00000 85: 0.00000 0.800000 86: 0.400000 1.60000 87: 0.200000 1.40000 88: 0.00000 1.20000 89: 0.200000 1.80000 90: 0.00000 1.60000 91: 0.00000 2.00000 Created vertex file "ell_vertex.txt". Created grid file "ell_grid.txt". Created command file "ell_commands.txt". Data written to the file "ell.xy" POLYGON_GRID_TEST: Normal end of execution. 11 May 2015 7:28:42.638 AM