11 September 2021 2:59:21.119 PM QUAD_MESH_TEST FORTRAN90 version Test the QUAD_MESH library. ADJ_SIZE_Q4_MESH_TEST ADJ_SIZE_Q4_MESH counts the node adjacencies. Number of adjacency entries is 186 Adjacency pointers: 1 1 4 2 5 10 3 11 16 4 17 22 5 23 26 6 27 32 7 33 41 8 42 50 9 51 59 10 60 67 11 68 71 12 72 77 13 78 85 14 86 93 15 94 102 16 103 111 17 112 117 18 118 123 19 124 131 20 132 139 21 140 148 22 149 156 23 157 160 24 161 164 25 165 170 26 171 176 27 177 182 28 183 186 AREA_Q4_MESH_TEST AREA_Q4_MESH computes the area of each element in a Q4 mesh. Element areas: 1 0.19134172 2 0.19134172 3 0.19134172 4 0.19134172 5 0.19134172 6 0.19134172 7 0.19134172 8 0.19134172 9 0.57402515 10 0.57402515 11 0.57402515 12 0.57402515 13 0.57402515 14 0.57402515 15 0.57402515 16 0.57402515 17 0.95670858 18 0.95670858 19 0.95670858 20 0.95670858 21 0.95670858 22 0.95670858 23 0.95670858 24 0.95670858 25 1.3393920 26 1.3393920 27 1.3393920 28 1.3393920 29 1.3393920 30 1.3393920 31 1.3393920 32 1.3393920 Mesh = 24.491740 AREA_QUAD_TEST AREA_QUAD computes the area of a quadrilateral. Area = 4.5000000 BOUNDARY_EDGE_COUNT_Q4_MESH_TEST BOUNDARY_EDGE_COUNT_Q4_MESH counts the boundary edges by looking at the mesh and counting unpaired edges. Number of boundary edges = 22 Correct number = 22 TEST05 BOUNDARY_EDGE_COUNT_EULER_Q4_MESH counts the boundary edges using Euler's formula. Number of boundary edges = 22 Correct number = 22 TEST06 EXAMPLE1_Q4_MESH sets up example #1 Q4 mesh. Number of nodes = 28 Number of elements = 17 Number of holes = 1 Node coordinates: Row 1 2 Col 1 0.00000 0.00000 2 1.00000 0.00000 3 2.00000 0.00000 4 3.00000 0.00000 5 4.00000 0.00000 6 0.00000 1.00000 7 1.00000 1.00000 8 2.00000 1.00000 9 3.00000 1.00000 10 4.00000 1.00000 11 5.00000 1.00000 12 0.00000 2.00000 13 1.00000 2.00000 14 2.00000 2.00000 15 3.00000 2.00000 16 4.00000 2.00000 17 5.00000 2.00000 18 0.00000 3.00000 19 1.00000 3.00000 20 2.00000 3.00000 21 3.00000 3.00000 22 4.00000 3.00000 23 5.00000 3.00000 24 0.00000 4.00000 25 1.00000 4.00000 26 2.00000 4.00000 27 3.00000 4.00000 28 4.00000 4.00000 Elements: Row 1 2 3 4 Col 1 1 2 7 6 2 2 3 8 7 3 3 4 9 8 4 4 5 10 9 5 6 7 13 12 6 7 8 14 13 7 8 9 15 14 8 9 10 16 15 9 10 11 17 16 10 12 13 19 18 11 14 15 21 20 12 15 16 22 21 13 16 17 23 22 14 18 19 25 24 15 19 20 26 25 16 20 21 27 26 17 21 22 28 27 Element neighbors Row 1 2 3 4 Col 1 -1 2 5 -1 2 -1 3 6 1 3 -1 4 7 2 4 -1 -1 8 3 5 1 6 10 -1 6 2 7 -2 5 7 3 8 11 6 8 4 9 12 7 9 -1 -1 13 8 10 5 -2 14 -1 11 7 12 16 -2 12 8 13 17 11 13 9 -1 -1 12 14 10 15 -1 -1 15 -2 16 -1 14 16 11 17 -1 15 17 12 -1 -1 16 TEST062 EXAMPLE1_Q4_MESH sets up example #1 Q4 mesh. PLOT_Q4_MESH plots it. Number of nodes = 28 Number of elements = 17 Number of holes = 1 TEST07 EXAMPLE2_Q4_MESH sets up example #2 Q4 mesh. Number of nodes = 37 Number of elements = 32 Number of holes = 0 Node coordinates: Row 1 2 Col 1 0.00000 0.00000 2 -1.00000 0.122465E-15 3 -0.923880 0.382683 4 -0.707107 0.707107 5 -0.382683 0.923880 6 0.612323E-16 1.00000 7 0.382683 0.923880 8 0.707107 0.707107 9 0.923880 0.382683 10 1.00000 0.00000 11 -2.00000 0.244929E-15 12 -1.84776 0.765367 13 -1.41421 1.41421 14 -0.765367 1.84776 15 0.122465E-15 2.00000 16 0.765367 1.84776 17 1.41421 1.41421 18 1.84776 0.765367 19 2.00000 0.00000 20 -3.00000 0.367394E-15 21 -2.77164 1.14805 22 -2.12132 2.12132 23 -1.14805 2.77164 24 0.183697E-15 3.00000 25 1.14805 2.77164 26 2.12132 2.12132 27 2.77164 1.14805 28 3.00000 0.00000 29 -4.00000 0.489859E-15 30 -3.69552 1.53073 31 -2.82843 2.82843 32 -1.53073 3.69552 33 0.244929E-15 4.00000 34 1.53073 3.69552 35 2.82843 2.82843 36 3.69552 1.53073 37 4.00000 0.00000 Elements: Row 1 2 3 4 Col 1 1 1 3 2 2 1 1 4 3 3 1 1 5 4 4 1 1 6 5 5 1 1 7 6 6 1 1 8 7 7 1 1 9 8 8 1 1 10 9 9 2 3 12 11 10 3 4 13 12 11 4 5 14 13 12 5 6 15 14 13 6 7 16 15 14 7 8 17 16 15 8 9 18 17 16 9 10 19 18 17 11 12 21 20 18 12 13 22 21 19 13 14 23 22 20 14 15 24 23 21 15 16 25 24 22 16 17 26 25 23 17 18 27 26 24 18 19 28 27 25 20 21 30 29 26 21 22 31 30 27 22 23 32 31 28 23 24 33 32 29 24 25 34 33 30 25 26 35 34 31 26 27 36 35 32 27 28 37 36 Element neighbors Row 1 2 3 4 Col 1 -1 2 9 -1 2 -1 3 10 1 3 -1 4 11 2 4 -1 5 12 3 5 -1 6 13 4 6 -1 7 14 5 7 -1 8 15 6 8 -1 -1 16 7 9 1 10 17 -1 10 2 11 18 9 11 3 12 19 10 12 4 13 20 11 13 5 14 21 12 14 6 15 22 13 15 7 16 23 14 16 8 -1 24 15 17 9 18 25 -1 18 10 19 26 17 19 11 20 27 18 20 12 21 28 19 21 13 22 29 20 22 14 23 30 21 23 15 24 31 22 24 16 -1 32 23 25 17 26 -1 -1 26 18 27 -1 25 27 19 28 -1 26 28 20 29 -1 27 29 21 30 -1 28 30 22 31 -1 29 31 23 32 -1 30 32 24 -1 -1 31 TEST072 EXAMPLE2_Q4_MESH sets up example #2 Q4 mesh. PLOT_Q4_MESH plots it. Number of nodes = 37 Number of elements = 32 Number of holes = 0 NEIGHBOR_ELEMENTS_Q4_MESH_TEST NEIGHBOR_ELEMENTS_Q4_MESH determines the adjacency relationships between elements. Element neighbors as reported by EXAMPLE1_Q4_MESH: Row 1 2 3 4 Col 1 -1 2 5 -1 2 -1 3 6 1 3 -1 4 7 2 4 -1 -1 8 3 5 1 6 10 -1 6 2 7 -2 5 7 3 8 11 6 8 4 9 12 7 9 -1 -1 13 8 10 5 -2 14 -1 11 7 12 16 -2 12 8 13 17 11 13 9 -1 -1 12 14 10 15 -1 -1 15 -2 16 -1 14 16 11 17 -1 15 17 12 -1 -1 16 Element neighbors computed by NEIGHBOR_ELEMENTS_Q4_MESH: Row 1 2 3 4 Col 1 -1 2 5 -1 2 -1 3 6 1 3 -1 4 7 2 4 -1 -1 8 3 5 1 6 10 -1 6 2 7 -1 5 7 3 8 11 6 8 4 9 12 7 9 -1 -1 13 8 10 5 -1 14 -1 11 7 12 16 -1 12 8 13 17 11 13 9 -1 -1 12 14 10 15 -1 -1 15 -1 16 -1 14 16 11 17 -1 15 17 12 -1 -1 16 TEST09 Write Q4 Mesh Example #1 to files. Element neighbors written to "q4_mesh_ex2_element_neighbors.txt". Elements written to "q4_mesh_ex2_elements.txt". Node coordinates written to "q4_mesh_ex2_xy.txt". NODE_ORDER_Q4_MESH_TEST NODE_ORDER_Q4_MESH computes the order of the nodes in a Q4 mesh. NODE ORDER 1 1 2 2 3 2 4 2 5 1 6 2 7 4 8 4 9 4 10 3 11 1 12 2 13 3 14 3 15 4 16 4 17 2 18 2 19 3 20 3 21 4 22 3 23 1 24 1 25 2 26 2 27 2 28 1 SAMPLE_Q4_MESH_TEST SAMPLE_Q4_MESH returns uniform sample points from a Q4 mesh. X Y Element 1 1.5289 3.1749 15 2 3.3200 1.9424 8 3 0.3362 0.5043 1 4 0.9748 2.5192 10 5 4.6930 2.0288 13 6 2.7831 2.9488 11 7 4.7097 2.1586 13 8 2.2468 3.6211 16 9 2.5209 3.2711 16 10 0.7547 0.9733 1 11 1.5424 3.7752 15 12 2.3345 3.1840 16 13 0.1029 2.2004 10 14 3.1786 0.1119 4 15 0.0036 2.5329 10 16 2.2232 3.1481 16 17 0.0749 2.5827 10 18 3.9438 0.6197 4 19 1.4811 3.8957 15 20 0.5854 1.5024 5 SAMPLE_QUAD_TEST SAMPLE_QUAD computes N random points in a quadrilateral. Write them to a file. Point coordinates written to "sample_quad.txt". QUAD_MESH_TEST Normal end of execution. 11 September 2021 2:59:21.125 PM