22 December 2015 11:08:34 AM QUAD_MESH_PRB C++ version Test the QUAD_MESH library. ADJ_SIZE_Q4_MESH_TEST ADJ_SIZE_Q4_MEXH counts the node adjacencies. Number of adjacency entries is 186 Adjacency pointers: 0 0 3 1 4 9 2 10 15 3 16 21 4 22 25 5 26 31 6 32 40 7 41 49 8 50 58 9 59 66 10 67 70 11 71 76 12 77 84 13 85 92 14 93 101 15 102 110 16 111 116 17 117 122 18 123 130 19 131 138 20 139 147 21 148 155 22 156 159 23 160 163 24 164 169 25 170 175 26 176 181 27 182 185 AREA_Q4_MESH_TEST AREA_Q4_MESH computes the area of each element in a Q4 mesh. Element areas: 0 1 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1 Mesh = 17 AREA_QUAD_TEST AREA_QUAD computes the area of a quadrilateral. Area = 4.5 BOUNDARY_EDGE_COUNT_Q4_MESH_TEST BOUNDARY_EDGE_COUNT_Q4_MESH counts the boundary edges. Number of boundary edges = 22 Correct number = 22 BOUNDARY_EDGE_COUNT_EULER_Q4_MESH_TEST BOUNDARY_EDGE_COUNT_EULER_Q4_MESH counts the boundary edges using Euler's formula. Number of boundary edges = 22 Correct number = 22 EXAMPLE1_Q4_MESH_TEST 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 0 2 1 0 3 2 0 4 3 0 5 4 0 6 0 1 7 1 1 8 2 1 9 3 1 10 4 1 11 5 1 12 0 2 13 1 2 14 2 2 15 3 2 16 4 2 17 5 2 18 0 3 19 1 3 20 2 3 21 3 3 22 4 3 23 5 3 24 0 4 25 1 4 26 2 4 27 3 4 28 4 4 Elements: Row: 1 2 3 4 Col 1 0 1 6 5 2 1 2 7 6 3 2 3 8 7 4 3 4 9 8 5 5 6 12 11 6 6 7 13 12 7 7 8 14 13 8 8 9 15 14 9 9 10 16 15 10 11 12 18 17 11 13 14 20 19 12 14 15 21 20 13 15 16 22 21 14 17 18 24 23 15 18 19 25 24 16 19 20 26 25 17 20 21 27 26 Element neighbors Row: 1 2 3 4 Col 1 -1 1 4 -1 2 -1 2 5 0 3 -1 3 6 1 4 -1 -1 7 2 5 0 5 9 -1 6 1 6 -2 4 7 2 7 10 5 8 3 8 11 6 9 -1 -1 12 7 10 4 -2 13 -1 11 6 11 15 -2 12 7 12 16 10 13 8 -1 -1 11 14 9 14 -1 -1 15 -2 15 -1 13 16 10 16 -1 14 17 11 -1 -1 15 EXAMPLE2_Q4_MESH_TEST 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 0 2 -1 1.22465e-16 3 -0.92388 0.382683 4 -0.707107 0.707107 5 -0.382683 0.92388 6 6.12323e-17 1 7 0.382683 0.92388 8 0.707107 0.707107 9 0.92388 0.382683 10 1 0 11 -2 2.44929e-16 12 -1.84776 0.765367 13 -1.41421 1.41421 14 -0.765367 1.84776 15 1.22465e-16 2 16 0.765367 1.84776 17 1.41421 1.41421 18 1.84776 0.765367 19 2 0 20 -3 3.67394e-16 21 -2.77164 1.14805 22 -2.12132 2.12132 23 -1.14805 2.77164 24 1.83697e-16 3 25 1.14805 2.77164 26 2.12132 2.12132 27 2.77164 1.14805 28 3 0 29 -4 4.89859e-16 30 -3.69552 1.53073 31 -2.82843 2.82843 32 -1.53073 3.69552 33 2.44929e-16 4 34 1.53073 3.69552 35 2.82843 2.82843 36 3.69552 1.53073 37 4 0 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 1 8 -1 2 -1 2 9 0 3 -1 3 10 1 4 -1 4 11 2 5 -1 5 12 3 6 -1 6 13 4 7 -1 7 14 5 8 -1 -1 15 6 9 0 9 16 -1 10 1 10 17 8 11 2 11 18 9 12 3 12 19 10 13 4 13 20 11 14 5 14 21 12 15 6 15 22 13 16 7 -1 23 14 17 8 17 24 -1 18 9 18 25 16 19 10 19 26 17 20 11 20 27 18 21 12 21 28 19 22 13 22 29 20 23 14 23 30 21 24 15 -1 31 22 25 16 25 -1 -1 26 17 26 -1 24 27 18 27 -1 25 28 19 28 -1 26 29 20 29 -1 27 30 21 30 -1 28 31 22 31 -1 29 32 23 -1 -1 30 TEST08 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 1 4 -1 2 -1 2 5 0 3 -1 3 6 1 4 -1 -1 7 2 5 0 5 9 -1 6 1 6 -2 4 7 2 7 10 5 8 3 8 11 6 9 -1 -1 12 7 10 4 -2 13 -1 11 6 11 15 -2 12 7 12 16 10 13 8 -1 -1 11 14 9 14 -1 -1 15 -2 15 -1 13 16 10 16 -1 14 17 11 -1 -1 15 Element neighbors computed by NEIGHBOR_ELEMENTS_Q4_MESH: Row: 1 2 3 4 Col 1 -1 1 4 -1 2 -1 2 5 0 3 -1 3 6 1 4 -1 -1 7 2 5 0 5 9 -1 6 1 6 -1 4 7 2 7 10 5 8 3 8 11 6 9 -1 -1 12 7 10 4 -1 13 -1 11 6 11 15 -1 12 7 12 16 10 13 8 -1 -1 11 14 9 14 -1 -1 15 -1 15 -1 13 16 10 16 -1 14 17 11 -1 -1 15 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". TEST10 NODE_ORDER_4_MESH computes the order of the nodes in a Q4 mesh. NODE ORDER 0 1 1 2 2 2 3 2 4 1 5 2 6 4 7 4 8 4 9 3 10 1 11 2 12 3 13 3 14 4 15 4 16 2 17 2 18 3 19 3 20 4 21 3 22 1 23 1 24 2 25 2 26 2 27 1 TEST105 SAMPLE_Q4_MESH returns uniform sample points from a Q4 mesh. X Y Element 0 3.08923 0.488422 4 1 3.50752 1.05581 8 2 0.75155 0.888312 1 3 2.10709 1.99836 7 4 2.30748 3.00419 16 5 1.64914 3.99736 15 6 0.662859 1.88145 5 7 0.831905 3.46725 14 8 1.95497 3.571 15 9 3.36993 0.883239 4 10 0.785624 2.28403 10 11 3.65692 0.897305 4 12 1.75128 3.82321 15 13 1.13669 0.733518 2 14 2.12723 3.64481 16 15 2.86651 3.93311 16 16 0.521271 0.447549 1 17 0.390364 0.0446259 1 18 2.45459 1.09679 7 19 2.23378 2.83489 11 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_PRB Normal end of execution. 22 December 2015 11:08:34 AM