From: Dominik Still Date: Thu, 21 Mar 2024 08:47:20 +0000 (+0100) Subject: Improve refinement of simplex meshes X-Git-Tag: v9.6.0-rc1~462^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=40f683839226f7b73336e54bf821f26a1d14118f;p=dealii.git Improve refinement of simplex meshes Improve the refinement strategy of Tetrahedrons. Formerly one way to refine a tetrahedron was hard-coded, now to cut the tetrahedron, the line between the 2 opposite midpoints of the partent's edges, which is the shortest, is chosen. This leads to more regular children. --- diff --git a/doc/news/changes/minor/20240320StillFehnKronbichler b/doc/news/changes/minor/20240320StillFehnKronbichler new file mode 100644 index 0000000000..ea5c2fb2cf --- /dev/null +++ b/doc/news/changes/minor/20240320StillFehnKronbichler @@ -0,0 +1,8 @@ +Improved: The isotropic refinement function +TriangulationImplementation::Implementation::execute_refinement_isotropic() +now supports 2 additional ways to split a tetrahedron into its children. +Cutting the parent between the 2 edge midpoints, which are the shortest +distance apart, leads to more regular refined cells and optimal convergence +under mesh refinement. +
+(Dominik Still, Niklas Fehn, Martin Kronbichler, 2024/03/20) diff --git a/source/grid/tria.cc b/source/grid/tria.cc index ce80da5ec1..bc38869e87 100644 --- a/source/grid/tria.cc +++ b/source/grid/tria.cc @@ -6704,8 +6704,16 @@ namespace internal const std::array line_indices = TriaAccessorImplementation::Implementation:: get_line_indices_of_cell(*hex); - // avoid a compiler warning by fixing the max number of - // loop iterations to 12 + + // For the tetrahedron the parent consists of the vertices + // 0,1,2,3, the new vertices 4-9 are defined as the + // midpoints fo the edges: 4 -> (0,1), 5 -> (1,2), 6 -> + // (2,0), 7 -> (0,3), 8 -> (1,3), 9 -> (2,3). + // Order is defined by the reference cell, see + // https://dealii.org/developer/doxygen/deal.II/group__simplex.html#simplex_reference_cells. + + // Avoid a compiler warning by fixing the max number of loop + // iterations to 12 const unsigned int n_lines = std::min(hex->n_lines(), 12u); for (unsigned int l = 0; l < n_lines; ++l) { @@ -6732,34 +6740,55 @@ namespace internal } } + unsigned int chosen_line_tetrahedron = 0; // set up new lines - { - static constexpr dealii::ndarray - new_line_vertices_hex = {{{{22, 26}}, - {{26, 23}}, - {{20, 26}}, - {{26, 21}}, - {{24, 26}}, - {{26, 25}}}}; - - static constexpr dealii::ndarray - new_line_vertices_tet = {{{{6, 8}}, - {{X, X}}, - {{X, X}}, - {{X, X}}, - {{X, X}}, - {{X, X}}}}; - - const auto &new_line_vertices = - (reference_cell_type == ReferenceCells::Hexahedron) ? - new_line_vertices_hex : - new_line_vertices_tet; + if (reference_cell_type == ReferenceCells::Hexahedron) + { + static constexpr dealii::ndarray + new_line_vertices = {{{{22, 26}}, + {{26, 23}}, + {{20, 26}}, + {{26, 21}}, + {{24, 26}}, + {{26, 25}}}}; + for (unsigned int i = 0; i < n_new_lines; ++i) + new_lines[i]->set_bounding_object_indices( + {vertex_indices[new_line_vertices[i][0]], + vertex_indices[new_line_vertices[i][1]]}); + } + else if (reference_cell_type == ReferenceCells::Tetrahedron) + { + // in the tetrahedron case, we have the three + // possibilities (6,8), (5,7), (4,9) -> pick the + // shortest line to guarantee the best possible aspect + // ratios + static constexpr dealii::ndarray + new_line_vertices = {{{{6, 8}}, {{5, 7}}, {{4, 9}}}}; + + const auto &vertices = triangulation.get_vertices(); + double min_distance = + std::numeric_limits::infinity(); + for (unsigned int i = 0; i < new_line_vertices.size(); + ++i) + { + const double current_distance = + vertices[vertex_indices[new_line_vertices[i][0]]] + .distance( + vertices + [vertex_indices[new_line_vertices[i][1]]]); + if (current_distance < min_distance) + { + chosen_line_tetrahedron = i; + min_distance = current_distance; + } + } - for (unsigned int i = 0; i < n_new_lines; ++i) - new_lines[i]->set_bounding_object_indices( - {vertex_indices[new_line_vertices[i][0]], - vertex_indices[new_line_vertices[i][1]]}); - } + new_lines[0]->set_bounding_object_indices( + {vertex_indices + [new_line_vertices[chosen_line_tetrahedron][0]], + vertex_indices + [new_line_vertices[chosen_line_tetrahedron][1]]}); + } // set up new quads { @@ -6799,6 +6828,17 @@ namespace internal } else if (reference_cell_type == ReferenceCells::Tetrahedron) { + // The order of the lines is defined by the ordering + // of the faces of the reference cell and the ordering + // of the lines within a face. + // Each face is split into 4 child triangels, the + // relevant lines are defined by the vertices of the + // center triangles: 0 -> (4,5), 1 -> (5,6), 2 -> (4,6), + // 3 -> (4,7), 4 -> (7,8), 5 -> (4,8), 6 -> (6,9), 7 -> + // (9,7), 8 -> (6,7), 9 -> (5,8), 10 -> (8,9), 11 -> + // (5,9), Line 12 is determined by + // chosen_line_tetrahedron i.e. (6,8), (5,7) or (4,9) + relevant_lines.resize(13); unsigned int k = 0; @@ -6828,7 +6868,6 @@ namespace internal } relevant_lines[k++] = new_lines[0]; - AssertDimension(k, 13); } else @@ -6854,19 +6893,61 @@ namespace internal {{3, 25, 26, 12}}, {{25, 7, 27, 13}}}}; - static constexpr dealii::ndarray - new_quad_lines_tet = {{{{2, 3, 8, X}}, - {{0, 9, 5, X}}, - {{1, 6, 11, X}}, - {{4, 10, 7, X}}, - {{2, 12, 5, X}}, - {{1, 9, 12, X}}, - {{4, 8, 12, X}}, - {{6, 12, 10, X}}, - {{X, X, X, X}}, - {{X, X, X, X}}, - {{X, X, X, X}}, - {{X, X, X, X}}}}; + // It is easierst to start at table cell_vertices_tet, + // there the vertices are listed which build up the + // 8 child tets. To build the child tets, 8 new faces are + // needed. The the vertices, which define the lines of these + // new faces are listed in table_tet. Now only the + // coresponding index of the lines and quads have to be + // listed in new_quad_lines_tet and cell_quads_tet. + + // The first 4 define the faces which cut off the + // parent tetrahedron at the edges. the numbers are the + // index of the relevant_lines defined above the last 4 + // faces cut apart the remaining octahedron, such that all + // of these contain line number 12. the ordering of the + // faces is arbitrary, the ordering within the faces has to + // follow the righthand convention for triangles + static constexpr dealii::ndarray + new_quad_lines_tet = {{// new line is (6,8) + {{{{2, 3, 8, X}}, + {{0, 9, 5, X}}, + {{1, 6, 11, X}}, + {{4, 10, 7, X}}, + {{2, 12, 5, X}}, + {{1, 9, 12, X}}, + {{4, 8, 12, X}}, + {{6, 12, 10, X}}, + {{X, X, X, X}}, + {{X, X, X, X}}, + {{X, X, X, X}}, + {{X, X, X, X}}}}, + // new line is (5,7) + {{{{2, 3, 8, X}}, + {{0, 9, 5, X}}, + {{1, 6, 11, X}}, + {{4, 10, 7, X}}, + {{0, 3, 12, X}}, + {{1, 12, 8, X}}, + {{4, 12, 9, X}}, + {{7, 11, 12, X}}, + {{X, X, X, X}}, + {{X, X, X, X}}, + {{X, X, X, X}}, + {{X, X, X, X}}}}, + // new line is (4,9) + {{{{2, 3, 8, X}}, + {{0, 9, 5, X}}, + {{1, 6, 11, X}}, + {{4, 10, 7, X}}, + {{0, 12, 11, X}}, + {{2, 6, 12, X}}, + {{3, 12, 7, X}}, + {{5, 10, 12, X}}, + {{X, X, X, X}}, + {{X, X, X, X}}, + {{X, X, X, X}}, + {{X, X, X, X}}}}}}; static constexpr dealii::ndarray table_hex = { @@ -6883,30 +6964,60 @@ namespace internal {{{{20, 18}}, {{26, 23}}, {{20, 26}}, {{18, 23}}}}, {{{{26, 23}}, {{21, 19}}, {{26, 21}}, {{23, 19}}}}}}; - static constexpr dealii::ndarray + // The table defines the vertices of the lines above + // see relevant_lines for mapping between line indices and + // vertex numbering + static constexpr dealii::ndarray table_tet = { - {{{{{6, 4}}, {{4, 7}}, {{7, 6}}, {{X, X}}}}, - {{{{4, 5}}, {{5, 8}}, {{8, 4}}, {{X, X}}}}, - {{{{5, 6}}, {{6, 9}}, {{9, 5}}, {{X, X}}}}, - {{{{7, 8}}, {{8, 9}}, {{9, 7}}, {{X, X}}}}, - {{{{4, 6}}, {{6, 8}}, {{8, 4}}, {{X, X}}}}, - {{{{6, 5}}, {{5, 8}}, {{8, 6}}, {{X, X}}}}, - {{{{8, 7}}, {{7, 6}}, {{6, 8}}, {{X, X}}}}, - {{{{9, 6}}, {{6, 8}}, {{8, 9}}, {{X, X}}}}, - {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}, - {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}, - {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}, - {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}}}; + {// new line is (6, 8) + {{{{{{6, 4}}, {{4, 7}}, {{7, 6}}, {{X, X}}}}, + {{{{4, 5}}, {{5, 8}}, {{8, 4}}, {{X, X}}}}, + {{{{5, 6}}, {{6, 9}}, {{9, 5}}, {{X, X}}}}, + {{{{7, 8}}, {{8, 9}}, {{9, 7}}, {{X, X}}}}, + {{{{4, 6}}, {{6, 8}}, {{8, 4}}, {{X, X}}}}, + {{{{6, 5}}, {{5, 8}}, {{8, 6}}, {{X, X}}}}, + {{{{8, 7}}, {{7, 6}}, {{6, 8}}, {{X, X}}}}, + {{{{9, 6}}, {{6, 8}}, {{8, 9}}, {{X, X}}}}, + {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}, + {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}, + {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}, + {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}}}, + // new line is (5, 7) + {{{{{{6, 4}}, {{4, 7}}, {{7, 6}}, {{X, X}}}}, + {{{{4, 5}}, {{5, 8}}, {{8, 4}}, {{X, X}}}}, + {{{{5, 6}}, {{6, 9}}, {{9, 5}}, {{X, X}}}}, + {{{{7, 8}}, {{8, 9}}, {{9, 7}}, {{X, X}}}}, + {{{{5, 4}}, {{4, 7}}, {{7, 5}}, {{X, X}}}}, + {{{{6, 5}}, {{5, 7}}, {{7, 6}}, {{X, X}}}}, + {{{{8, 7}}, {{7, 5}}, {{5, 8}}, {{X, X}}}}, + {{{{7, 9}}, {{9, 5}}, {{5, 7}}, {{X, X}}}}, + {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}, + {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}, + {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}, + {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}}}, + // new line is (4, 9) + {{{{{{6, 4}}, {{4, 7}}, {{7, 6}}, {{X, X}}}}, + {{{{4, 5}}, {{5, 8}}, {{8, 4}}, {{X, X}}}}, + {{{{5, 6}}, {{6, 9}}, {{9, 5}}, {{X, X}}}}, + {{{{7, 8}}, {{8, 9}}, {{9, 7}}, {{X, X}}}}, + {{{{5, 4}}, {{4, 9}}, {{9, 5}}, {{X, X}}}}, + {{{{4, 6}}, {{6, 9}}, {{9, 4}}, {{X, X}}}}, + {{{{7, 4}}, {{4, 9}}, {{9, 7}}, {{X, X}}}}, + {{{{4, 8}}, {{8, 9}}, {{9, 4}}, {{X, X}}}}, + {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}, + {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}, + {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}, + {{{{X, X}}, {{X, X}}, {{X, X}}, {{X, X}}}}}}}}; const auto &new_quad_lines = (reference_cell_type == ReferenceCells::Hexahedron) ? new_quad_lines_hex : - new_quad_lines_tet; + new_quad_lines_tet[chosen_line_tetrahedron]; const auto &table = (reference_cell_type == ReferenceCells::Hexahedron) ? table_hex : - table_tet; + table_tet[chosen_line_tetrahedron]; static constexpr dealii::ndarray representative_lines{ @@ -6997,6 +7108,12 @@ namespace internal } else if (reference_cell_type == ReferenceCells::Tetrahedron) { + // list of the indices of the surfaces which define the + // 8 new tets. the indices 0-7 are the new quads defined + // above (so 0-3 cut off the corners and 4-7 sperate the + // remaining octahedral), the indices between 8-11 are + // the children of the first face, from 12-15 of the + // second, etc. for (unsigned int i = 0; i < n_new_quads; ++i) quad_indices[i] = new_quads[i]->index(); @@ -7035,70 +7152,93 @@ namespace internal {{3, 19, 7, 27, 11, 35}} // }}; - static constexpr dealii::ndarray - cell_quads_tet{{{{8, 13, 16, 0, X, X}}, - {{9, 12, 1, 21, X, X}}, - {{10, 2, 17, 20, X, X}}, - {{3, 14, 18, 22, X, X}}, - {{11, 1, 4, 5, X, X}}, - {{15, 0, 4, 6, X, X}}, - {{19, 7, 6, 3, X, X}}, - {{23, 5, 2, 7, X, X}}}}; - - static constexpr dealii::ndarray - cell_face_vertices_tet{{{{{{0, 4, 6, X}}, - {{4, 0, 7, X}}, - {{0, 6, 7, X}}, - {{6, 4, 7, X}}, - {{X, X, X, X}}, - {{X, X, X, X}}}}, - {{{{4, 1, 5, X}}, - {{1, 4, 8, X}}, - {{4, 5, 8, X}}, - {{5, 1, 8, X}}, - {{X, X, X, X}}, - {{X, X, X, X}}}}, - {{{{6, 5, 2, X}}, - {{5, 6, 9, X}}, - {{6, 2, 9, X}}, - {{2, 5, 9, X}}, - {{X, X, X, X}}, - {{X, X, X, X}}}}, - {{{{7, 8, 9, X}}, - {{8, 7, 3, X}}, - {{7, 9, 3, X}}, - {{9, 8, 3, X}}, - {{X, X, X, X}}, - {{X, X, X, X}}}}, - {{{{4, 5, 6, X}}, - {{5, 4, 8, X}}, - {{4, 6, 8, X}}, - {{6, 5, 8, X}}, - {{X, X, X, X}}, - {{X, X, X, X}}}}, - {{{{4, 7, 8, X}}, - {{7, 4, 6, X}}, - {{4, 8, 6, X}}, - {{8, 7, 6, X}}, - {{X, X, X, X}}, - {{X, X, X, X}}}}, - {{{{6, 9, 7, X}}, - {{9, 6, 8, X}}, - {{6, 7, 8, X}}, - {{7, 9, 8, X}}, - {{X, X, X, X}}, - {{X, X, X, X}}}}, - {{{{5, 8, 9, X}}, - {{8, 5, 6, X}}, - {{5, 9, 6, X}}, - {{9, 8, 6, X}}, - {{X, X, X, X}}, - {{X, X, X, X}}}}}}; + // indices of the faces which define the new tets + // the ordering of the tets is arbitrary + // the first 4 determine the tets cutting of the corners + // the last 4 are ordered after their appearance in the + // faces. + // the ordering within the faces is determined by + // convention for the tetrahedron unit cell, see + // cell_vertices_tet below + static constexpr dealii::ndarray + cell_quads_tet = {{// new line is (6, 8) + {{{{8, 13, 16, 0, X, X}}, + {{9, 12, 1, 21, X, X}}, + {{10, 2, 17, 20, X, X}}, + {{3, 14, 18, 22, X, X}}, + {{11, 1, 4, 5, X, X}}, + {{15, 0, 4, 6, X, X}}, + {{19, 7, 6, 3, X, X}}, + {{23, 5, 2, 7, X, X}}}}, + // new line is (5, 7) + {{ + {{8, 13, 16, 0, X, X}}, + {{9, 12, 1, 21, X, X}}, + {{10, 2, 17, 20, X, X}}, + {{3, 14, 18, 22, X, X}}, + {{11, 4, 0, 5, X, X}}, + {{15, 4, 1, 6, X, X}}, + {{19, 2, 5, 7, X, X}}, + {{23, 6, 7, 3, X, X}}, + }}, + // new line is (4, 9) + {{ + {{8, 13, 16, 0, X, X}}, + {{9, 12, 1, 21, X, X}}, + {{10, 2, 17, 20, X, X}}, + {{3, 14, 18, 22, X, X}}, + {{11, 4, 5, 2, X, X}}, + {{15, 6, 7, 3, X, X}}, + {{19, 5, 0, 6, X, X}}, + {{23, 1, 4, 7, X, X}}, + }}}}; + + // the 4 faces of the reference tet are defined by 3 + // vertices each and follow the order of the unit cell + static constexpr dealii::ndarray + reference_cell_face_vertices = { + {{{0, 1, 2}}, {{1, 0, 3}}, {{0, 2, 3}}, {{2, 1, 3}}}}; + + // the 8 child tets are each defined by 4 vertices + // the ordering of the tets has to be consistent with above + // the ordering within the tets is given by the reference + // tet i.e. looking at the fifth line the first 3 vertices + // are given by face 11, the last vertex is the remaining of + // the tet + static constexpr dealii::ndarray + cell_vertices_tet = {{// new line is (6,8) + {{{{0, 4, 6, 7}}, + {{4, 1, 5, 8}}, + {{6, 5, 2, 9}}, + {{7, 8, 9, 3}}, + {{4, 5, 6, 8}}, + {{4, 7, 8, 6}}, + {{6, 9, 7, 8}}, + {{5, 8, 9, 6}}}}, + // new line is (5,7) + {{{{0, 4, 6, 7}}, + {{4, 1, 5, 8}}, + {{6, 5, 2, 9}}, + {{7, 8, 9, 3}}, + {{4, 5, 6, 7}}, + {{4, 7, 8, 5}}, + {{6, 9, 7, 5}}, + {{5, 8, 9, 7}}}}, + // new line is (4,9) + {{{{0, 4, 6, 7}}, + {{4, 1, 5, 8}}, + {{6, 5, 2, 9}}, + {{7, 8, 9, 3}}, + {{4, 5, 6, 9}}, + {{4, 7, 8, 9}}, + {{6, 9, 7, 4}}, + {{5, 8, 9, 4}}}}}}; + const auto &cell_quads = (reference_cell_type == ReferenceCells::Hexahedron) ? cell_quads_hex : - cell_quads_tet; + cell_quads_tet[chosen_line_tetrahedron]; for (unsigned int c = 0; c < GeometryInfo::max_children_per_cell; @@ -7114,6 +7254,7 @@ namespace internal quad_indices[cell_quads[c][2]], quad_indices[cell_quads[c][3]]}); + // for tets, we need to go through the faces and // figure the orientation out the hard way for (const auto f : new_hex->face_indices()) @@ -7128,14 +7269,19 @@ namespace internal face->vertex_index(1), face->vertex_index(2)}}; + const auto &new_hex_vertices = + cell_vertices_tet[chosen_line_tetrahedron][c]; + const auto face_f_vertices = + reference_cell_face_vertices[f]; + const std::array vertices_1 = { { - vertex_indices[cell_face_vertices_tet[c][f] - [0]], - vertex_indices[cell_face_vertices_tet[c][f] - [1]], - vertex_indices[cell_face_vertices_tet[c][f] - [2]], + vertex_indices + [new_hex_vertices[face_f_vertices[0]]], + vertex_indices + [new_hex_vertices[face_f_vertices[1]]], + vertex_indices + [new_hex_vertices[face_f_vertices[2]]], }}; new_hex->set_combined_face_orientation( diff --git a/tests/simplex/convert_hypercube_to_simplex_mesh_01.output b/tests/simplex/convert_hypercube_to_simplex_mesh_01.output index 6497021e42..f25a2a3ef9 100644 --- a/tests/simplex/convert_hypercube_to_simplex_mesh_01.output +++ b/tests/simplex/convert_hypercube_to_simplex_mesh_01.output @@ -1102,50 +1102,50 @@ CELLS 1248 5472 4 33 1 39 38 4 98 39 12 99 4 35 38 99 10 -4 33 39 98 38 -4 33 35 38 98 -4 98 99 35 38 -4 39 38 99 98 +4 33 39 98 99 +4 33 35 38 99 +4 98 99 35 33 +4 39 38 99 33 4 2 42 91 45 4 42 3 50 51 4 91 50 11 94 4 45 51 94 12 -4 42 50 91 51 -4 42 45 51 91 -4 91 94 45 51 -4 50 51 94 91 +4 42 50 91 94 +4 42 45 51 94 +4 91 94 45 42 +4 50 51 94 42 4 7 74 92 76 4 74 6 68 69 4 92 68 11 95 4 76 69 95 13 -4 74 68 92 69 -4 74 76 69 92 -4 92 95 76 69 -4 68 69 95 92 +4 74 68 92 95 +4 74 76 69 95 +4 92 95 76 74 +4 68 69 95 74 4 5 61 101 63 4 61 4 59 58 4 101 59 13 104 4 63 58 104 10 -4 61 59 101 58 -4 61 63 58 101 -4 101 104 63 58 -4 59 58 104 101 +4 61 59 101 104 +4 61 63 58 104 +4 101 104 63 61 +4 59 58 104 61 4 0 34 81 98 4 34 2 44 45 4 81 44 8 84 4 98 45 84 12 -4 34 44 81 45 -4 34 98 45 81 -4 81 84 98 45 -4 44 45 84 81 +4 34 44 81 84 +4 34 98 45 84 +4 81 84 98 34 +4 44 45 84 34 4 4 56 59 57 4 56 6 69 67 4 59 69 13 102 4 57 67 102 8 -4 56 69 59 67 -4 56 57 67 59 -4 59 102 57 67 -4 69 67 102 59 +4 56 69 59 102 +4 56 57 67 102 +4 59 102 57 56 +4 69 67 102 56 4 5 101 73 62 4 101 13 76 103 4 73 76 7 75 @@ -1174,58 +1174,58 @@ CELLS 1248 5472 4 36 5 62 63 4 37 62 9 86 4 38 63 86 10 -4 36 62 37 63 -4 36 38 63 37 -4 37 86 38 63 -4 62 63 86 37 +4 36 62 37 86 +4 36 38 63 86 +4 37 86 38 36 +4 62 63 86 36 4 3 49 50 85 4 49 7 92 75 4 50 92 11 93 4 85 75 93 9 -4 49 92 50 75 -4 49 85 75 50 -4 50 93 85 75 -4 92 75 93 50 +4 49 92 50 93 +4 49 85 75 93 +4 50 93 85 49 +4 92 75 93 49 4 2 43 44 91 4 43 6 67 68 4 44 67 8 83 4 91 68 83 11 -4 43 67 44 68 -4 43 91 68 44 -4 44 83 91 68 -4 67 68 83 44 +4 43 67 44 83 +4 43 91 68 83 +4 44 83 91 43 +4 67 68 83 43 4 12 100 99 87 4 100 13 104 103 4 99 104 10 86 4 87 103 86 9 -4 100 104 99 103 -4 100 87 103 99 -4 99 86 87 103 -4 104 103 86 99 +4 100 104 99 86 +4 100 87 103 86 +4 99 86 87 100 +4 104 103 86 100 4 12 100 87 94 4 100 13 103 95 4 87 103 9 93 4 94 95 93 11 -4 100 103 87 95 -4 100 94 95 87 -4 87 93 94 95 -4 103 95 93 87 +4 100 103 87 93 +4 100 94 95 93 +4 87 93 94 100 +4 103 95 93 100 4 12 100 94 84 4 100 13 95 102 4 94 95 11 83 4 84 102 83 8 -4 100 95 94 102 -4 100 84 102 94 -4 94 83 84 102 -4 95 102 83 94 +4 100 95 94 83 +4 100 84 102 83 +4 94 83 84 100 +4 95 102 83 100 4 12 100 84 99 4 100 13 102 104 4 84 102 8 82 4 99 104 82 10 -4 100 102 84 104 -4 100 99 104 84 -4 84 82 99 104 -4 102 104 82 84 +4 100 102 84 82 +4 100 99 104 82 +4 84 82 99 100 +4 102 104 82 100 4 13 102 104 59 4 102 8 82 57 4 104 82 10 58 @@ -1254,10 +1254,10 @@ CELLS 1248 5472 4 95 11 83 68 4 102 83 8 67 4 69 68 67 6 -4 95 83 102 68 -4 95 69 68 102 -4 102 67 69 68 -4 83 68 67 102 +4 95 83 102 67 +4 95 69 68 67 +4 102 67 69 95 +4 83 68 67 95 4 10 99 86 38 4 99 12 87 39 4 86 87 9 37 @@ -1294,50 +1294,50 @@ CELLS 1248 5472 4 108 15 113 112 4 139 113 21 140 4 109 112 140 20 -4 108 113 139 112 -4 108 109 112 139 -4 139 140 109 112 -4 113 112 140 139 +4 108 113 139 140 +4 108 109 112 140 +4 139 140 109 108 +4 113 112 140 108 4 1 48 37 41 4 48 3 85 53 4 37 85 9 89 4 41 53 89 21 -4 48 85 37 53 -4 48 41 53 37 -4 37 89 41 53 -4 85 53 89 37 +4 48 85 37 89 +4 48 41 53 89 +4 37 89 41 48 +4 85 53 89 48 4 7 73 75 79 4 73 5 62 65 4 75 62 9 90 4 79 65 90 22 -4 73 62 75 65 -4 73 79 65 75 -4 75 90 79 65 -4 62 65 90 75 +4 73 62 75 90 +4 73 79 65 90 +4 75 90 79 73 +4 62 65 90 73 4 17 123 143 125 4 123 16 120 119 4 143 120 22 146 4 125 119 146 20 -4 123 120 143 119 -4 123 125 119 143 -4 143 146 125 119 -4 120 119 146 143 +4 123 120 143 146 +4 123 125 119 146 +4 143 146 125 123 +4 120 119 146 123 4 14 107 130 139 4 107 1 40 41 4 130 40 18 132 4 139 41 132 21 -4 107 40 130 41 -4 107 139 41 130 -4 130 132 139 41 -4 40 41 132 130 +4 107 40 130 132 +4 107 139 41 132 +4 130 132 139 107 +4 40 41 132 107 4 16 116 120 118 4 116 5 65 64 4 120 65 22 144 4 118 64 144 18 -4 116 65 120 64 -4 116 118 64 120 -4 120 144 118 64 -4 65 64 144 120 +4 116 65 120 144 +4 116 118 64 144 +4 120 144 118 116 +4 65 64 144 116 4 17 143 77 124 4 143 22 79 145 4 77 79 7 78 @@ -1366,66 +1366,66 @@ CELLS 1248 5472 4 110 17 124 125 4 111 124 19 135 4 112 125 135 20 -4 110 124 111 125 -4 110 112 125 111 -4 111 135 112 125 -4 124 125 135 111 +4 110 124 111 135 +4 110 112 125 135 +4 111 135 112 110 +4 124 125 135 110 4 3 49 85 133 4 49 7 75 78 4 85 75 9 88 4 133 78 88 19 -4 49 75 85 78 -4 49 133 78 85 -4 85 88 133 78 -4 75 78 88 85 +4 49 75 85 88 +4 49 133 78 88 +4 85 88 133 49 +4 75 78 88 49 4 1 36 40 37 4 36 5 64 62 4 40 64 18 129 4 37 62 129 9 -4 36 64 40 62 -4 36 37 62 40 -4 40 129 37 62 -4 64 62 129 40 +4 36 64 40 129 +4 36 37 62 129 +4 40 129 37 36 +4 64 62 129 36 4 21 141 140 136 4 141 22 146 145 4 140 146 20 135 4 136 145 135 19 -4 141 146 140 145 -4 141 136 145 140 -4 140 135 136 145 -4 146 145 135 140 +4 141 146 140 135 +4 141 136 145 135 +4 140 135 136 141 +4 146 145 135 141 4 21 141 136 89 4 141 22 145 90 4 136 145 19 88 4 89 90 88 9 -4 141 145 136 90 -4 141 89 90 136 -4 136 88 89 90 -4 145 90 88 136 +4 141 145 136 88 +4 141 89 90 88 +4 136 88 89 141 +4 145 90 88 141 4 21 141 89 132 4 141 22 90 144 4 89 90 9 129 4 132 144 129 18 -4 141 90 89 144 -4 141 132 144 89 -4 89 129 132 144 -4 90 144 129 89 +4 141 90 89 129 +4 141 132 144 129 +4 89 129 132 141 +4 90 144 129 141 4 21 141 132 140 4 141 22 144 146 4 132 144 18 131 4 140 146 131 20 -4 141 144 132 146 -4 141 140 146 132 -4 132 131 140 146 -4 144 146 131 132 +4 141 144 132 131 +4 141 140 146 131 +4 132 131 140 141 +4 144 146 131 141 4 22 144 146 120 4 144 18 131 118 4 146 131 20 119 4 120 118 119 16 -4 144 131 146 118 -4 144 120 118 146 -4 146 119 120 118 -4 131 118 119 146 +4 144 131 146 119 +4 144 120 118 119 +4 146 119 120 144 +4 131 118 119 144 4 22 146 145 143 4 146 20 135 125 4 145 135 19 124 @@ -1438,10 +1438,10 @@ CELLS 1248 5472 4 145 19 88 78 4 90 88 9 75 4 79 78 75 7 -4 145 88 90 78 -4 145 79 78 90 -4 90 75 79 78 -4 88 78 75 90 +4 145 88 90 75 +4 145 79 78 75 +4 90 75 79 145 +4 88 78 75 145 4 22 90 144 65 4 90 9 129 62 4 144 129 18 64 @@ -1454,82 +1454,82 @@ CELLS 1248 5472 4 140 21 136 113 4 135 136 19 111 4 112 113 111 15 -4 140 136 135 113 -4 140 112 113 135 -4 135 111 112 113 -4 136 113 111 135 +4 140 136 135 112 +4 140 112 113 136 +4 135 111 112 136 +4 136 113 111 112 4 19 136 88 133 4 136 21 89 53 4 88 89 9 85 4 133 53 85 3 -4 136 89 88 53 -4 136 133 53 88 -4 88 85 133 53 -4 89 53 85 88 +4 136 89 88 85 +4 136 133 53 85 +4 88 85 133 136 +4 89 53 85 136 4 9 89 129 37 4 89 21 132 41 4 129 132 18 40 4 37 41 40 1 -4 89 132 129 41 -4 89 37 41 129 -4 129 40 37 41 -4 132 41 40 129 +4 89 132 129 37 +4 89 37 41 132 +4 129 40 37 132 +4 132 41 40 37 4 18 132 131 130 4 132 21 140 139 4 131 140 20 109 4 130 139 109 14 -4 132 140 131 139 -4 132 130 139 131 -4 131 109 130 139 -4 140 139 109 131 +4 132 140 131 109 +4 132 130 139 109 +4 131 109 130 132 +4 140 139 109 132 4 23 148 162 150 4 148 2 47 46 4 162 47 27 163 4 150 46 163 26 -4 148 47 162 46 -4 148 150 46 162 -4 162 163 150 46 -4 47 46 163 162 +4 148 47 162 163 +4 148 150 46 163 +4 162 163 150 148 +4 47 46 163 148 4 15 52 111 115 4 52 3 133 54 4 111 133 19 137 4 115 54 137 27 -4 52 133 111 54 -4 52 115 54 111 -4 111 137 115 54 -4 133 54 137 111 +4 52 133 111 137 +4 52 115 54 137 +4 111 137 115 52 +4 133 54 137 52 4 7 77 78 80 4 77 17 124 127 4 78 124 19 138 4 80 127 138 28 -4 77 124 78 127 -4 77 80 127 78 -4 78 138 80 127 -4 124 127 138 78 +4 77 124 78 138 +4 77 80 127 138 +4 78 138 80 77 +4 124 127 138 77 4 6 70 165 71 4 70 24 155 154 4 165 155 28 170 4 71 154 170 26 -4 70 155 165 154 -4 70 71 154 165 -4 165 170 71 154 -4 155 154 170 165 +4 70 155 165 170 +4 70 71 154 170 +4 165 170 71 70 +4 155 154 170 70 4 23 149 159 162 4 149 15 114 115 4 159 114 25 161 4 162 115 161 27 -4 149 114 159 115 -4 149 162 115 159 -4 159 161 162 115 -4 114 115 161 159 +4 149 114 159 161 +4 149 162 115 161 +4 159 161 162 149 +4 114 115 161 149 4 24 151 155 153 4 151 17 127 126 4 155 127 28 169 4 153 126 169 25 -4 151 127 155 126 -4 151 153 126 155 -4 155 169 153 126 -4 127 126 169 155 +4 151 127 155 169 +4 151 153 126 169 +4 155 169 153 151 +4 127 126 169 151 4 6 165 74 68 4 165 28 80 166 4 74 80 7 92 @@ -1558,58 +1558,58 @@ CELLS 1248 5472 4 43 6 68 71 4 91 68 11 96 4 46 71 96 26 -4 43 68 91 71 -4 43 46 71 91 -4 91 96 46 71 -4 68 71 96 91 +4 43 68 91 96 +4 43 46 71 96 +4 91 96 46 43 +4 68 71 96 43 4 3 49 133 50 4 49 7 78 92 4 133 78 19 134 4 50 92 134 11 -4 49 78 133 92 -4 49 50 92 133 -4 133 134 50 92 -4 78 92 134 133 +4 49 78 133 134 +4 49 50 92 134 +4 133 134 50 49 +4 78 92 134 49 4 15 110 114 111 4 110 17 126 124 4 114 126 25 158 4 111 124 158 19 -4 110 126 114 124 -4 110 111 124 114 -4 114 158 111 124 -4 126 124 158 114 +4 110 126 114 158 +4 110 111 124 158 +4 114 158 111 110 +4 126 124 158 110 4 27 164 163 97 4 164 28 170 166 4 163 170 26 96 4 97 166 96 11 -4 164 170 163 166 -4 164 97 166 163 -4 163 96 97 166 -4 170 166 96 163 +4 164 170 163 96 +4 164 97 166 96 +4 163 96 97 164 +4 170 166 96 164 4 27 164 97 137 4 164 28 166 138 4 97 166 11 134 4 137 138 134 19 -4 164 166 97 138 -4 164 137 138 97 -4 97 134 137 138 -4 166 138 134 97 +4 164 166 97 134 +4 164 137 138 134 +4 97 134 137 164 +4 166 138 134 164 4 27 164 137 161 4 164 28 138 169 4 137 138 19 158 4 161 169 158 25 -4 164 138 137 169 -4 164 161 169 137 -4 137 158 161 169 -4 138 169 158 137 +4 164 138 137 158 +4 164 161 169 158 +4 137 158 161 164 +4 138 169 158 164 4 27 164 161 163 4 164 28 169 170 4 161 169 25 160 4 163 170 160 26 -4 164 169 161 170 -4 164 163 170 161 -4 161 160 163 170 -4 169 170 160 161 +4 164 169 161 160 +4 164 163 170 160 +4 161 160 163 164 +4 169 170 160 164 4 28 169 170 155 4 169 25 160 153 4 170 160 26 154 @@ -1622,10 +1622,10 @@ CELLS 1248 5472 4 170 26 96 71 4 166 96 11 68 4 165 71 68 6 -4 170 96 166 71 -4 170 165 71 166 -4 166 68 165 71 -4 96 71 68 166 +4 170 96 166 68 +4 170 165 71 68 +4 166 68 165 170 +4 96 71 68 170 4 28 166 138 80 4 166 11 134 92 4 138 134 19 78 @@ -1638,90 +1638,90 @@ CELLS 1248 5472 4 138 19 158 124 4 169 158 25 126 4 127 124 126 17 -4 138 158 169 124 -4 138 127 124 169 -4 169 126 127 124 -4 158 124 126 169 +4 138 158 169 126 +4 138 127 124 126 +4 169 126 127 138 +4 158 124 126 138 4 26 163 96 46 4 163 27 97 47 4 96 97 11 91 4 46 47 91 2 -4 163 97 96 47 -4 163 46 47 96 -4 96 91 46 47 -4 97 47 91 96 +4 163 97 96 91 +4 163 46 47 91 +4 96 91 46 163 +4 97 47 91 163 4 11 97 134 50 4 97 27 137 54 4 134 137 19 133 4 50 54 133 3 -4 97 137 134 54 -4 97 50 54 134 -4 134 133 50 54 -4 137 54 133 134 +4 97 137 134 50 +4 97 50 54 137 +4 134 133 50 137 +4 137 54 133 50 4 19 137 158 111 4 137 27 161 115 4 158 161 25 114 4 111 115 114 15 -4 137 161 158 115 -4 137 111 115 158 -4 158 114 111 115 -4 161 115 114 158 +4 137 161 158 114 +4 137 111 115 114 +4 158 114 111 137 +4 161 115 114 137 4 25 161 160 159 4 161 27 163 162 4 160 163 26 150 4 159 162 150 23 -4 161 163 160 162 -4 161 159 162 160 -4 160 150 159 162 -4 163 162 150 160 +4 161 163 160 159 +4 161 159 162 163 +4 160 150 159 163 +4 163 162 150 159 4 4 61 59 60 4 61 5 101 66 4 59 101 13 105 4 60 66 105 31 -4 61 101 59 66 -4 61 60 66 59 -4 59 105 60 66 -4 101 66 105 59 +4 61 101 59 105 +4 61 60 66 105 +4 59 105 60 61 +4 101 66 105 61 4 6 74 165 69 4 74 7 80 76 4 165 80 28 167 4 69 76 167 13 -4 74 80 165 76 -4 74 69 76 165 -4 165 167 69 76 -4 80 76 167 165 +4 74 80 165 167 +4 74 69 76 167 +4 165 167 69 74 +4 80 76 167 74 4 17 151 127 128 4 151 24 155 157 4 127 155 28 171 4 128 157 171 32 -4 151 155 127 157 -4 151 128 157 127 -4 127 171 128 157 -4 155 157 171 127 +4 151 155 127 171 +4 151 128 157 171 +4 127 171 128 151 +4 155 157 171 151 4 16 121 181 122 4 121 29 176 175 4 181 176 32 184 4 122 175 184 31 -4 121 176 181 175 -4 121 122 175 181 -4 181 184 122 175 -4 176 175 184 181 +4 121 176 181 184 +4 121 122 175 184 +4 181 184 122 121 +4 176 175 184 121 4 4 56 177 59 4 56 6 72 69 4 177 72 30 178 4 59 69 178 13 -4 56 72 177 69 -4 56 59 69 177 -4 177 178 59 69 -4 72 69 178 177 +4 56 72 177 178 +4 56 59 69 178 +4 177 178 59 56 +4 72 69 178 56 4 29 173 176 174 4 173 24 157 156 4 176 157 32 183 4 174 156 183 30 -4 173 157 176 156 -4 173 174 156 176 -4 176 183 174 156 -4 157 156 183 176 +4 173 157 176 183 +4 173 174 156 183 +4 176 183 174 173 +4 157 156 183 173 4 16 181 123 120 4 181 32 128 182 4 123 128 17 143 @@ -1750,90 +1750,90 @@ CELLS 1248 5472 4 116 16 120 122 4 65 120 22 147 4 66 122 147 31 -4 116 120 65 122 -4 116 66 122 65 -4 65 147 66 122 -4 120 122 147 65 +4 116 120 65 147 +4 116 66 122 147 +4 65 147 66 116 +4 120 122 147 116 4 7 77 80 79 4 77 17 127 143 4 80 127 28 168 4 79 143 168 22 -4 77 127 80 143 -4 77 79 143 80 -4 80 168 79 143 -4 127 143 168 80 +4 77 127 80 168 +4 77 79 143 168 +4 80 168 79 77 +4 127 143 168 77 4 6 70 72 165 4 70 24 156 155 4 72 156 30 179 4 165 155 179 28 -4 70 156 72 155 -4 70 165 155 72 -4 72 179 165 155 -4 156 155 179 72 +4 70 156 72 179 +4 70 165 155 179 +4 72 179 165 70 +4 156 155 179 70 4 13 106 105 142 4 106 32 184 182 4 105 184 31 147 4 142 182 147 22 -4 106 184 105 182 -4 106 142 182 105 -4 105 147 142 182 -4 184 182 147 105 +4 106 184 105 147 +4 106 142 182 147 +4 105 147 142 106 +4 184 182 147 106 4 13 106 142 167 4 106 32 182 171 4 142 182 22 168 4 167 171 168 28 -4 106 182 142 171 -4 106 167 171 142 -4 142 168 167 171 -4 182 171 168 142 +4 106 182 142 168 +4 106 167 171 168 +4 142 168 167 106 +4 182 171 168 106 4 13 106 167 178 4 106 32 171 183 4 167 171 28 179 4 178 183 179 30 -4 106 171 167 183 -4 106 178 183 167 -4 167 179 178 183 -4 171 183 179 167 +4 106 171 167 179 +4 106 178 183 179 +4 167 179 178 106 +4 171 183 179 106 4 13 106 178 105 4 106 32 183 184 4 178 183 30 180 4 105 184 180 31 -4 106 183 178 184 -4 106 105 184 178 -4 178 180 105 184 -4 183 184 180 178 +4 106 183 178 180 +4 106 105 184 180 +4 178 180 105 106 +4 183 184 180 106 4 32 183 184 176 4 183 30 180 174 4 184 180 31 175 4 176 174 175 29 -4 183 180 184 174 -4 183 176 174 184 -4 184 175 176 174 -4 180 174 175 184 +4 183 180 184 176 +4 183 176 174 180 +4 184 175 176 180 +4 180 174 175 176 4 32 184 182 181 4 184 31 147 122 4 182 147 22 120 4 181 122 120 16 -4 184 147 182 122 -4 184 181 122 182 -4 182 120 181 122 -4 147 122 120 182 +4 184 147 182 181 +4 184 181 122 147 +4 182 120 181 147 +4 147 122 120 181 4 32 182 171 128 4 182 22 168 143 4 171 168 28 127 4 128 143 127 17 -4 182 168 171 143 -4 182 128 143 171 -4 171 127 128 143 -4 168 143 127 171 +4 182 168 171 128 +4 182 128 143 168 +4 171 127 128 168 +4 168 143 127 128 4 32 171 183 157 4 171 28 179 155 4 183 179 30 156 4 157 155 156 24 -4 171 179 183 155 -4 171 157 155 183 -4 183 156 157 155 -4 179 155 156 183 +4 171 179 183 157 +4 171 157 155 179 +4 183 156 157 179 +4 179 155 156 157 4 31 105 147 66 4 105 13 142 101 4 147 142 22 65 diff --git a/tests/simplex/data_out_write_gnuplot_01.output b/tests/simplex/data_out_write_gnuplot_01.output index c56be9e692..01d3388410 100644 --- a/tests/simplex/data_out_write_gnuplot_01.output +++ b/tests/simplex/data_out_write_gnuplot_01.output @@ -2061,643 +2061,643 @@ # For a description of the GNUPLOT format see the GNUPLOT manual. # # -0.00000 0.00000 0.00000 -0.00854702 -0.250000 0.00000 0.00000 0.0507372 -0.00000 0.250000 0.00000 -0.00723101 -0.00000 0.00000 0.00000 -0.00854702 -0.00000 0.00000 0.250000 -0.00557879 -0.00000 0.250000 0.00000 -0.00723101 - -0.00000 0.00000 0.250000 -0.00557879 -0.250000 0.00000 0.00000 0.0507372 - -0.250000 0.00000 0.00000 0.0507372 -0.500000 0.00000 0.00000 0.244680 -0.250000 0.250000 0.00000 0.0504510 -0.250000 0.00000 0.00000 0.0507372 -0.250000 0.00000 0.250000 0.0480283 -0.250000 0.250000 0.00000 0.0504510 - -0.250000 0.00000 0.250000 0.0480283 -0.500000 0.00000 0.00000 0.244680 - -0.00000 0.250000 0.00000 -0.00723101 -0.250000 0.250000 0.00000 0.0504510 -0.00000 0.500000 0.00000 -0.00764870 -0.00000 0.250000 0.00000 -0.00723101 -0.00000 0.250000 0.250000 -0.0211579 -0.00000 0.500000 0.00000 -0.00764870 - -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.250000 0.00000 0.0504510 - -0.00000 0.00000 0.250000 -0.00557879 -0.250000 0.00000 0.250000 0.0480283 -0.00000 0.250000 0.250000 -0.0211579 -0.00000 0.00000 0.250000 -0.00557879 -0.00000 0.00000 0.500000 -0.00337312 -0.00000 0.250000 0.250000 -0.0211579 - -0.00000 0.00000 0.500000 -0.00337312 -0.250000 0.00000 0.250000 0.0480283 - -0.250000 0.00000 0.00000 0.0507372 -0.250000 0.250000 0.00000 0.0504510 -0.00000 0.250000 0.00000 -0.00723101 -0.250000 0.00000 0.00000 0.0507372 -0.250000 0.00000 0.250000 0.0480283 -0.00000 0.250000 0.00000 -0.00723101 - -0.250000 0.00000 0.250000 0.0480283 -0.250000 0.250000 0.00000 0.0504510 - -0.250000 0.00000 0.00000 0.0507372 -0.00000 0.00000 0.250000 -0.00557879 -0.250000 0.00000 0.250000 0.0480283 -0.250000 0.00000 0.00000 0.0507372 -0.00000 0.250000 0.00000 -0.00723101 -0.250000 0.00000 0.250000 0.0480283 - -0.00000 0.250000 0.00000 -0.00723101 -0.00000 0.00000 0.250000 -0.00557879 - -0.00000 0.250000 0.00000 -0.00723101 -0.00000 0.250000 0.250000 -0.0211579 -0.00000 0.00000 0.250000 -0.00557879 -0.00000 0.250000 0.00000 -0.00723101 -0.250000 0.00000 0.250000 0.0480283 -0.00000 0.00000 0.250000 -0.00557879 - -0.250000 0.00000 0.250000 0.0480283 -0.00000 0.250000 0.250000 -0.0211579 - -0.250000 0.250000 0.00000 0.0504510 -0.250000 0.00000 0.250000 0.0480283 -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.250000 0.00000 0.0504510 -0.00000 0.250000 0.00000 -0.00723101 -0.00000 0.250000 0.250000 -0.0211579 - -0.00000 0.250000 0.00000 -0.00723101 -0.250000 0.00000 0.250000 0.0480283 - -0.500000 0.00000 0.00000 0.244680 -0.750000 0.00000 0.00000 0.554111 -0.500000 0.250000 0.00000 0.243049 -0.500000 0.00000 0.00000 0.244680 -0.500000 0.00000 0.250000 0.227790 -0.500000 0.250000 0.00000 0.243049 - -0.500000 0.00000 0.250000 0.227790 -0.750000 0.00000 0.00000 0.554111 - -0.750000 0.00000 0.00000 0.554111 -1.00000 0.00000 0.00000 0.983987 -0.750000 0.250000 0.00000 0.550862 -0.750000 0.00000 0.00000 0.554111 -0.750000 0.00000 0.250000 0.552054 -0.750000 0.250000 0.00000 0.550862 - -0.750000 0.00000 0.250000 0.552054 -1.00000 0.00000 0.00000 0.983987 - -0.500000 0.250000 0.00000 0.243049 -0.750000 0.250000 0.00000 0.550862 -0.500000 0.500000 0.00000 0.238853 -0.500000 0.250000 0.00000 0.243049 -0.500000 0.250000 0.250000 0.240449 -0.500000 0.500000 0.00000 0.238853 - -0.500000 0.250000 0.250000 0.240449 -0.750000 0.250000 0.00000 0.550862 - -0.500000 0.00000 0.250000 0.227790 -0.750000 0.00000 0.250000 0.552054 -0.500000 0.250000 0.250000 0.240449 -0.500000 0.00000 0.250000 0.227790 -0.500000 0.00000 0.500000 0.241196 -0.500000 0.250000 0.250000 0.240449 - -0.500000 0.00000 0.500000 0.241196 -0.750000 0.00000 0.250000 0.552054 - -0.750000 0.00000 0.00000 0.554111 -0.750000 0.250000 0.00000 0.550862 -0.500000 0.250000 0.00000 0.243049 -0.750000 0.00000 0.00000 0.554111 -0.750000 0.00000 0.250000 0.552054 -0.500000 0.250000 0.00000 0.243049 - -0.750000 0.00000 0.250000 0.552054 -0.750000 0.250000 0.00000 0.550862 - -0.750000 0.00000 0.00000 0.554111 -0.500000 0.00000 0.250000 0.227790 -0.750000 0.00000 0.250000 0.552054 -0.750000 0.00000 0.00000 0.554111 -0.500000 0.250000 0.00000 0.243049 -0.750000 0.00000 0.250000 0.552054 - -0.500000 0.250000 0.00000 0.243049 -0.500000 0.00000 0.250000 0.227790 - -0.500000 0.250000 0.00000 0.243049 -0.500000 0.250000 0.250000 0.240449 -0.500000 0.00000 0.250000 0.227790 -0.500000 0.250000 0.00000 0.243049 -0.750000 0.00000 0.250000 0.552054 -0.500000 0.00000 0.250000 0.227790 - -0.750000 0.00000 0.250000 0.552054 -0.500000 0.250000 0.250000 0.240449 - -0.750000 0.250000 0.00000 0.550862 -0.750000 0.00000 0.250000 0.552054 -0.500000 0.250000 0.250000 0.240449 -0.750000 0.250000 0.00000 0.550862 -0.500000 0.250000 0.00000 0.243049 -0.500000 0.250000 0.250000 0.240449 - -0.500000 0.250000 0.00000 0.243049 -0.750000 0.00000 0.250000 0.552054 - -0.00000 0.500000 0.00000 -0.00764870 -0.250000 0.500000 0.00000 0.0516556 -0.00000 0.750000 0.00000 -0.0114345 -0.00000 0.500000 0.00000 -0.00764870 -0.00000 0.500000 0.250000 -0.00982944 -0.00000 0.750000 0.00000 -0.0114345 - -0.00000 0.500000 0.250000 -0.00982944 -0.250000 0.500000 0.00000 0.0516556 - -0.250000 0.500000 0.00000 0.0516556 -0.500000 0.500000 0.00000 0.238853 -0.250000 0.750000 0.00000 0.0516143 -0.250000 0.500000 0.00000 0.0516556 -0.250000 0.500000 0.250000 0.0517263 -0.250000 0.750000 0.00000 0.0516143 - -0.250000 0.500000 0.250000 0.0517263 -0.500000 0.500000 0.00000 0.238853 - -0.00000 0.750000 0.00000 -0.0114345 -0.250000 0.750000 0.00000 0.0516143 -0.00000 1.00000 0.00000 -0.00563512 -0.00000 0.750000 0.00000 -0.0114345 -0.00000 0.750000 0.250000 -0.00807627 -0.00000 1.00000 0.00000 -0.00563512 - -0.00000 0.750000 0.250000 -0.00807627 -0.250000 0.750000 0.00000 0.0516143 - -0.00000 0.500000 0.250000 -0.00982944 -0.250000 0.500000 0.250000 0.0517263 -0.00000 0.750000 0.250000 -0.00807627 -0.00000 0.500000 0.250000 -0.00982944 -0.00000 0.500000 0.500000 -0.00919997 -0.00000 0.750000 0.250000 -0.00807627 - -0.00000 0.500000 0.500000 -0.00919997 -0.250000 0.500000 0.250000 0.0517263 - -0.250000 0.500000 0.00000 0.0516556 -0.250000 0.750000 0.00000 0.0516143 -0.00000 0.750000 0.00000 -0.0114345 -0.250000 0.500000 0.00000 0.0516556 -0.250000 0.500000 0.250000 0.0517263 -0.00000 0.750000 0.00000 -0.0114345 - -0.250000 0.500000 0.250000 0.0517263 -0.250000 0.750000 0.00000 0.0516143 - -0.250000 0.500000 0.00000 0.0516556 -0.00000 0.500000 0.250000 -0.00982944 -0.250000 0.500000 0.250000 0.0517263 -0.250000 0.500000 0.00000 0.0516556 -0.00000 0.750000 0.00000 -0.0114345 -0.250000 0.500000 0.250000 0.0517263 - -0.00000 0.750000 0.00000 -0.0114345 -0.00000 0.500000 0.250000 -0.00982944 - -0.00000 0.750000 0.00000 -0.0114345 -0.00000 0.750000 0.250000 -0.00807627 -0.00000 0.500000 0.250000 -0.00982944 -0.00000 0.750000 0.00000 -0.0114345 -0.250000 0.500000 0.250000 0.0517263 -0.00000 0.500000 0.250000 -0.00982944 - -0.250000 0.500000 0.250000 0.0517263 -0.00000 0.750000 0.250000 -0.00807627 - -0.250000 0.750000 0.00000 0.0516143 -0.250000 0.500000 0.250000 0.0517263 -0.00000 0.750000 0.250000 -0.00807627 -0.250000 0.750000 0.00000 0.0516143 -0.00000 0.750000 0.00000 -0.0114345 -0.00000 0.750000 0.250000 -0.00807627 - -0.00000 0.750000 0.00000 -0.0114345 -0.250000 0.500000 0.250000 0.0517263 - -0.00000 0.00000 0.500000 -0.00337312 -0.250000 0.00000 0.500000 0.0499617 -0.00000 0.250000 0.500000 -0.00803541 -0.00000 0.00000 0.500000 -0.00337312 -0.00000 0.00000 0.750000 -0.00883340 -0.00000 0.250000 0.500000 -0.00803541 - -0.00000 0.00000 0.750000 -0.00883340 -0.250000 0.00000 0.500000 0.0499617 - -0.250000 0.00000 0.500000 0.0499617 -0.500000 0.00000 0.500000 0.241196 -0.250000 0.250000 0.500000 0.0524966 -0.250000 0.00000 0.500000 0.0499617 -0.250000 0.00000 0.750000 0.0485802 -0.250000 0.250000 0.500000 0.0524966 - -0.250000 0.00000 0.750000 0.0485802 -0.500000 0.00000 0.500000 0.241196 - -0.00000 0.250000 0.500000 -0.00803541 -0.250000 0.250000 0.500000 0.0524966 -0.00000 0.500000 0.500000 -0.00919997 -0.00000 0.250000 0.500000 -0.00803541 -0.00000 0.250000 0.750000 -0.00870722 -0.00000 0.500000 0.500000 -0.00919997 - -0.00000 0.250000 0.750000 -0.00870722 -0.250000 0.250000 0.500000 0.0524966 - -0.00000 0.00000 0.750000 -0.00883340 -0.250000 0.00000 0.750000 0.0485802 -0.00000 0.250000 0.750000 -0.00870722 -0.00000 0.00000 0.750000 -0.00883340 -0.00000 0.00000 1.00000 -0.00510315 -0.00000 0.250000 0.750000 -0.00870722 - -0.00000 0.00000 1.00000 -0.00510315 -0.250000 0.00000 0.750000 0.0485802 - -0.250000 0.00000 0.500000 0.0499617 -0.250000 0.250000 0.500000 0.0524966 -0.00000 0.250000 0.500000 -0.00803541 -0.250000 0.00000 0.500000 0.0499617 -0.250000 0.00000 0.750000 0.0485802 -0.00000 0.250000 0.500000 -0.00803541 - -0.250000 0.00000 0.750000 0.0485802 -0.250000 0.250000 0.500000 0.0524966 - -0.250000 0.00000 0.500000 0.0499617 -0.00000 0.00000 0.750000 -0.00883340 -0.250000 0.00000 0.750000 0.0485802 -0.250000 0.00000 0.500000 0.0499617 -0.00000 0.250000 0.500000 -0.00803541 -0.250000 0.00000 0.750000 0.0485802 - -0.00000 0.250000 0.500000 -0.00803541 -0.00000 0.00000 0.750000 -0.00883340 - -0.00000 0.250000 0.500000 -0.00803541 -0.00000 0.250000 0.750000 -0.00870722 -0.00000 0.00000 0.750000 -0.00883340 -0.00000 0.250000 0.500000 -0.00803541 -0.250000 0.00000 0.750000 0.0485802 -0.00000 0.00000 0.750000 -0.00883340 - -0.250000 0.00000 0.750000 0.0485802 -0.00000 0.250000 0.750000 -0.00870722 - -0.250000 0.250000 0.500000 0.0524966 -0.250000 0.00000 0.750000 0.0485802 -0.00000 0.250000 0.750000 -0.00870722 -0.250000 0.250000 0.500000 0.0524966 -0.00000 0.250000 0.500000 -0.00803541 -0.00000 0.250000 0.750000 -0.00870722 - -0.00000 0.250000 0.500000 -0.00803541 -0.250000 0.00000 0.750000 0.0485802 - -0.500000 0.00000 0.00000 0.244680 -0.500000 0.250000 0.00000 0.243049 -0.250000 0.250000 0.00000 0.0504510 -0.500000 0.00000 0.00000 0.244680 -0.500000 0.00000 0.250000 0.227790 -0.250000 0.250000 0.00000 0.0504510 - -0.500000 0.00000 0.250000 0.227790 -0.500000 0.250000 0.00000 0.243049 - -0.500000 0.250000 0.00000 0.243049 -0.500000 0.500000 0.00000 0.238853 -0.250000 0.500000 0.00000 0.0516556 -0.500000 0.250000 0.00000 0.243049 -0.500000 0.250000 0.250000 0.240449 -0.250000 0.500000 0.00000 0.0516556 - -0.500000 0.250000 0.250000 0.240449 -0.500000 0.500000 0.00000 0.238853 - -0.250000 0.250000 0.00000 0.0504510 -0.250000 0.500000 0.00000 0.0516556 -0.00000 0.500000 0.00000 -0.00764870 -0.250000 0.250000 0.00000 0.0504510 -0.250000 0.250000 0.250000 0.0523980 -0.00000 0.500000 0.00000 -0.00764870 - -0.250000 0.250000 0.250000 0.0523980 -0.250000 0.500000 0.00000 0.0516556 - -0.500000 0.00000 0.250000 0.227790 -0.500000 0.250000 0.250000 0.240449 -0.250000 0.250000 0.250000 0.0523980 -0.500000 0.00000 0.250000 0.227790 -0.500000 0.00000 0.500000 0.241196 -0.250000 0.250000 0.250000 0.0523980 - -0.500000 0.00000 0.500000 0.241196 -0.500000 0.250000 0.250000 0.240449 - -0.500000 0.250000 0.00000 0.243049 -0.250000 0.500000 0.00000 0.0516556 -0.250000 0.250000 0.00000 0.0504510 -0.500000 0.250000 0.00000 0.243049 -0.500000 0.250000 0.250000 0.240449 -0.250000 0.250000 0.00000 0.0504510 - -0.500000 0.250000 0.250000 0.240449 -0.250000 0.500000 0.00000 0.0516556 - -0.500000 0.250000 0.00000 0.243049 -0.500000 0.00000 0.250000 0.227790 -0.500000 0.250000 0.250000 0.240449 -0.500000 0.250000 0.00000 0.243049 -0.250000 0.250000 0.00000 0.0504510 -0.500000 0.250000 0.250000 0.240449 - -0.250000 0.250000 0.00000 0.0504510 -0.500000 0.00000 0.250000 0.227790 - -0.250000 0.250000 0.00000 0.0504510 -0.250000 0.250000 0.250000 0.0523980 -0.500000 0.00000 0.250000 0.227790 -0.250000 0.250000 0.00000 0.0504510 -0.500000 0.250000 0.250000 0.240449 -0.500000 0.00000 0.250000 0.227790 - -0.500000 0.250000 0.250000 0.240449 -0.250000 0.250000 0.250000 0.0523980 - -0.250000 0.500000 0.00000 0.0516556 -0.500000 0.250000 0.250000 0.240449 -0.250000 0.250000 0.250000 0.0523980 -0.250000 0.500000 0.00000 0.0516556 -0.250000 0.250000 0.00000 0.0504510 -0.250000 0.250000 0.250000 0.0523980 - -0.250000 0.250000 0.00000 0.0504510 -0.500000 0.250000 0.250000 0.240449 - -0.500000 0.00000 0.00000 0.244680 -0.250000 0.00000 0.250000 0.0480283 -0.500000 0.00000 0.250000 0.227790 -0.500000 0.00000 0.00000 0.244680 -0.250000 0.250000 0.00000 0.0504510 -0.500000 0.00000 0.250000 0.227790 - -0.250000 0.250000 0.00000 0.0504510 -0.250000 0.00000 0.250000 0.0480283 - -0.250000 0.00000 0.250000 0.0480283 -0.00000 0.00000 0.500000 -0.00337312 -0.250000 0.00000 0.500000 0.0499617 -0.250000 0.00000 0.250000 0.0480283 -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.00000 0.500000 0.0499617 - -0.00000 0.250000 0.250000 -0.0211579 -0.00000 0.00000 0.500000 -0.00337312 - -0.500000 0.00000 0.250000 0.227790 -0.250000 0.00000 0.500000 0.0499617 -0.500000 0.00000 0.500000 0.241196 -0.500000 0.00000 0.250000 0.227790 -0.250000 0.250000 0.250000 0.0523980 -0.500000 0.00000 0.500000 0.241196 - -0.250000 0.250000 0.250000 0.0523980 -0.250000 0.00000 0.500000 0.0499617 - -0.250000 0.250000 0.00000 0.0504510 -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.250000 0.250000 0.0523980 -0.250000 0.250000 0.00000 0.0504510 -0.00000 0.500000 0.00000 -0.00764870 -0.250000 0.250000 0.250000 0.0523980 - -0.00000 0.500000 0.00000 -0.00764870 -0.00000 0.250000 0.250000 -0.0211579 - -0.250000 0.00000 0.250000 0.0480283 -0.250000 0.00000 0.500000 0.0499617 -0.500000 0.00000 0.250000 0.227790 -0.250000 0.00000 0.250000 0.0480283 -0.00000 0.250000 0.250000 -0.0211579 -0.500000 0.00000 0.250000 0.227790 - -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.00000 0.500000 0.0499617 - -0.250000 0.00000 0.250000 0.0480283 -0.250000 0.250000 0.00000 0.0504510 -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.00000 0.250000 0.0480283 -0.500000 0.00000 0.250000 0.227790 -0.00000 0.250000 0.250000 -0.0211579 - -0.500000 0.00000 0.250000 0.227790 -0.250000 0.250000 0.00000 0.0504510 - -0.500000 0.00000 0.250000 0.227790 -0.250000 0.250000 0.250000 0.0523980 -0.250000 0.250000 0.00000 0.0504510 -0.500000 0.00000 0.250000 0.227790 -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.250000 0.00000 0.0504510 - -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.250000 0.250000 0.0523980 - -0.250000 0.00000 0.500000 0.0499617 -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.250000 0.250000 0.0523980 -0.250000 0.00000 0.500000 0.0499617 -0.500000 0.00000 0.250000 0.227790 -0.250000 0.250000 0.250000 0.0523980 - -0.500000 0.00000 0.250000 0.227790 -0.00000 0.250000 0.250000 -0.0211579 - -0.00000 0.500000 0.00000 -0.00764870 -0.00000 0.500000 0.250000 -0.00982944 -0.00000 0.250000 0.250000 -0.0211579 -0.00000 0.500000 0.00000 -0.00764870 -0.250000 0.250000 0.250000 0.0523980 -0.00000 0.250000 0.250000 -0.0211579 - -0.250000 0.250000 0.250000 0.0523980 -0.00000 0.500000 0.250000 -0.00982944 - -0.00000 0.500000 0.250000 -0.00982944 -0.00000 0.500000 0.500000 -0.00919997 -0.00000 0.250000 0.500000 -0.00803541 -0.00000 0.500000 0.250000 -0.00982944 -0.250000 0.250000 0.500000 0.0524966 -0.00000 0.250000 0.500000 -0.00803541 - -0.250000 0.250000 0.500000 0.0524966 -0.00000 0.500000 0.500000 -0.00919997 - -0.00000 0.250000 0.250000 -0.0211579 -0.00000 0.250000 0.500000 -0.00803541 -0.00000 0.00000 0.500000 -0.00337312 -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.00000 0.500000 0.0499617 -0.00000 0.00000 0.500000 -0.00337312 - -0.250000 0.00000 0.500000 0.0499617 -0.00000 0.250000 0.500000 -0.00803541 - -0.250000 0.250000 0.250000 0.0523980 -0.250000 0.250000 0.500000 0.0524966 -0.250000 0.00000 0.500000 0.0499617 -0.250000 0.250000 0.250000 0.0523980 -0.500000 0.00000 0.500000 0.241196 -0.250000 0.00000 0.500000 0.0499617 - -0.500000 0.00000 0.500000 0.241196 -0.250000 0.250000 0.500000 0.0524966 - -0.00000 0.500000 0.250000 -0.00982944 -0.00000 0.250000 0.500000 -0.00803541 -0.00000 0.250000 0.250000 -0.0211579 -0.00000 0.500000 0.250000 -0.00982944 -0.250000 0.250000 0.500000 0.0524966 -0.00000 0.250000 0.250000 -0.0211579 - -0.250000 0.250000 0.500000 0.0524966 -0.00000 0.250000 0.500000 -0.00803541 - -0.00000 0.500000 0.250000 -0.00982944 -0.250000 0.250000 0.250000 0.0523980 -0.250000 0.250000 0.500000 0.0524966 -0.00000 0.500000 0.250000 -0.00982944 -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.250000 0.500000 0.0524966 - -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.250000 0.250000 0.0523980 - -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.00000 0.500000 0.0499617 -0.250000 0.250000 0.250000 0.0523980 -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.250000 0.500000 0.0524966 -0.250000 0.250000 0.250000 0.0523980 - -0.250000 0.250000 0.500000 0.0524966 -0.250000 0.00000 0.500000 0.0499617 - -0.00000 0.250000 0.500000 -0.00803541 -0.250000 0.250000 0.500000 0.0524966 -0.250000 0.00000 0.500000 0.0499617 -0.00000 0.250000 0.500000 -0.00803541 -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.00000 0.500000 0.0499617 - -0.00000 0.250000 0.250000 -0.0211579 -0.250000 0.250000 0.500000 0.0524966 - -0.500000 0.500000 0.00000 0.238853 -0.500000 0.250000 0.250000 0.240449 -0.250000 0.500000 0.250000 0.0517263 -0.500000 0.500000 0.00000 0.238853 -0.250000 0.500000 0.00000 0.0516556 -0.250000 0.500000 0.250000 0.0517263 - -0.250000 0.500000 0.00000 0.0516556 -0.500000 0.250000 0.250000 0.240449 - -0.500000 0.250000 0.250000 0.240449 -0.500000 0.00000 0.500000 0.241196 -0.250000 0.250000 0.500000 0.0524966 -0.500000 0.250000 0.250000 0.240449 -0.250000 0.250000 0.250000 0.0523980 -0.250000 0.250000 0.500000 0.0524966 - -0.250000 0.250000 0.250000 0.0523980 -0.500000 0.00000 0.500000 0.241196 - -0.250000 0.500000 0.250000 0.0517263 -0.250000 0.250000 0.500000 0.0524966 -0.00000 0.500000 0.500000 -0.00919997 -0.250000 0.500000 0.250000 0.0517263 -0.00000 0.500000 0.250000 -0.00982944 -0.00000 0.500000 0.500000 -0.00919997 - -0.00000 0.500000 0.250000 -0.00982944 -0.250000 0.250000 0.500000 0.0524966 - -0.250000 0.500000 0.00000 0.0516556 -0.250000 0.250000 0.250000 0.0523980 -0.00000 0.500000 0.250000 -0.00982944 -0.250000 0.500000 0.00000 0.0516556 -0.00000 0.500000 0.00000 -0.00764870 -0.00000 0.500000 0.250000 -0.00982944 - -0.00000 0.500000 0.00000 -0.00764870 -0.250000 0.250000 0.250000 0.0523980 - -0.500000 0.250000 0.250000 0.240449 -0.250000 0.250000 0.500000 0.0524966 -0.250000 0.500000 0.250000 0.0517263 -0.500000 0.250000 0.250000 0.240449 -0.250000 0.250000 0.250000 0.0523980 -0.250000 0.500000 0.250000 0.0517263 - -0.250000 0.250000 0.250000 0.0523980 -0.250000 0.250000 0.500000 0.0524966 - -0.500000 0.250000 0.250000 0.240449 -0.250000 0.500000 0.00000 0.0516556 -0.250000 0.250000 0.250000 0.0523980 -0.500000 0.250000 0.250000 0.240449 -0.250000 0.500000 0.250000 0.0517263 -0.250000 0.250000 0.250000 0.0523980 - -0.250000 0.500000 0.250000 0.0517263 -0.250000 0.500000 0.00000 0.0516556 - -0.250000 0.500000 0.250000 0.0517263 -0.00000 0.500000 0.250000 -0.00982944 -0.250000 0.500000 0.00000 0.0516556 -0.250000 0.500000 0.250000 0.0517263 -0.250000 0.250000 0.250000 0.0523980 -0.250000 0.500000 0.00000 0.0516556 - -0.250000 0.250000 0.250000 0.0523980 -0.00000 0.500000 0.250000 -0.00982944 - -0.250000 0.250000 0.500000 0.0524966 -0.250000 0.250000 0.250000 0.0523980 -0.00000 0.500000 0.250000 -0.00982944 -0.250000 0.250000 0.500000 0.0524966 -0.250000 0.500000 0.250000 0.0517263 -0.00000 0.500000 0.250000 -0.00982944 - -0.250000 0.500000 0.250000 0.0517263 -0.250000 0.250000 0.250000 0.0523980 +0.00000 0.00000 0.00000 -0.00590900 +0.250000 0.00000 0.00000 0.0512603 +0.00000 0.250000 0.00000 -0.0105936 +0.00000 0.00000 0.00000 -0.00590900 +0.00000 0.00000 0.250000 -0.00801537 +0.00000 0.250000 0.00000 -0.0105936 + +0.00000 0.00000 0.250000 -0.00801537 +0.250000 0.00000 0.00000 0.0512603 + +0.250000 0.00000 0.00000 0.0512603 +0.500000 0.00000 0.00000 0.239570 +0.250000 0.250000 0.00000 0.0511687 +0.250000 0.00000 0.00000 0.0512603 +0.250000 0.00000 0.250000 0.0519658 +0.250000 0.250000 0.00000 0.0511687 + +0.250000 0.00000 0.250000 0.0519658 +0.500000 0.00000 0.00000 0.239570 + +0.00000 0.250000 0.00000 -0.0105936 +0.250000 0.250000 0.00000 0.0511687 +0.00000 0.500000 0.00000 -0.0103118 +0.00000 0.250000 0.00000 -0.0105936 +0.00000 0.250000 0.250000 -0.0117534 +0.00000 0.500000 0.00000 -0.0103118 + +0.00000 0.250000 0.250000 -0.0117534 +0.250000 0.250000 0.00000 0.0511687 + +0.00000 0.00000 0.250000 -0.00801537 +0.250000 0.00000 0.250000 0.0519658 +0.00000 0.250000 0.250000 -0.0117534 +0.00000 0.00000 0.250000 -0.00801537 +0.00000 0.00000 0.500000 -0.00809848 +0.00000 0.250000 0.250000 -0.0117534 + +0.00000 0.00000 0.500000 -0.00809848 +0.250000 0.00000 0.250000 0.0519658 + +0.250000 0.00000 0.00000 0.0512603 +0.250000 0.250000 0.00000 0.0511687 +0.00000 0.250000 0.00000 -0.0105936 +0.250000 0.00000 0.00000 0.0512603 +0.250000 0.00000 0.250000 0.0519658 +0.00000 0.250000 0.00000 -0.0105936 + +0.250000 0.00000 0.250000 0.0519658 +0.250000 0.250000 0.00000 0.0511687 + +0.250000 0.00000 0.00000 0.0512603 +0.00000 0.00000 0.250000 -0.00801537 +0.250000 0.00000 0.250000 0.0519658 +0.250000 0.00000 0.00000 0.0512603 +0.00000 0.250000 0.00000 -0.0105936 +0.250000 0.00000 0.250000 0.0519658 + +0.00000 0.250000 0.00000 -0.0105936 +0.00000 0.00000 0.250000 -0.00801537 + +0.00000 0.250000 0.00000 -0.0105936 +0.00000 0.250000 0.250000 -0.0117534 +0.00000 0.00000 0.250000 -0.00801537 +0.00000 0.250000 0.00000 -0.0105936 +0.250000 0.00000 0.250000 0.0519658 +0.00000 0.00000 0.250000 -0.00801537 + +0.250000 0.00000 0.250000 0.0519658 +0.00000 0.250000 0.250000 -0.0117534 + +0.250000 0.250000 0.00000 0.0511687 +0.250000 0.00000 0.250000 0.0519658 +0.00000 0.250000 0.250000 -0.0117534 +0.250000 0.250000 0.00000 0.0511687 +0.00000 0.250000 0.00000 -0.0105936 +0.00000 0.250000 0.250000 -0.0117534 + +0.00000 0.250000 0.00000 -0.0105936 +0.250000 0.00000 0.250000 0.0519658 + +0.500000 0.00000 0.00000 0.239570 +0.750000 0.00000 0.00000 0.552558 +0.500000 0.250000 0.00000 0.241213 +0.500000 0.00000 0.00000 0.239570 +0.500000 0.00000 0.250000 0.239642 +0.500000 0.250000 0.00000 0.241213 + +0.500000 0.00000 0.250000 0.239642 +0.750000 0.00000 0.00000 0.552558 + +0.750000 0.00000 0.00000 0.552558 +1.00000 0.00000 0.00000 0.984642 +0.750000 0.250000 0.00000 0.552908 +0.750000 0.00000 0.00000 0.552558 +0.750000 0.00000 0.250000 0.550251 +0.750000 0.250000 0.00000 0.552908 + +0.750000 0.00000 0.250000 0.550251 +1.00000 0.00000 0.00000 0.984642 + +0.500000 0.250000 0.00000 0.241213 +0.750000 0.250000 0.00000 0.552908 +0.500000 0.500000 0.00000 0.239596 +0.500000 0.250000 0.00000 0.241213 +0.500000 0.250000 0.250000 0.238578 +0.500000 0.500000 0.00000 0.239596 + +0.500000 0.250000 0.250000 0.238578 +0.750000 0.250000 0.00000 0.552908 + +0.500000 0.00000 0.250000 0.239642 +0.750000 0.00000 0.250000 0.550251 +0.500000 0.250000 0.250000 0.238578 +0.500000 0.00000 0.250000 0.239642 +0.500000 0.00000 0.500000 0.238105 +0.500000 0.250000 0.250000 0.238578 + +0.500000 0.00000 0.500000 0.238105 +0.750000 0.00000 0.250000 0.550251 + +0.750000 0.00000 0.00000 0.552558 +0.750000 0.250000 0.00000 0.552908 +0.500000 0.250000 0.00000 0.241213 +0.750000 0.00000 0.00000 0.552558 +0.750000 0.00000 0.250000 0.550251 +0.500000 0.250000 0.00000 0.241213 + +0.750000 0.00000 0.250000 0.550251 +0.750000 0.250000 0.00000 0.552908 + +0.750000 0.00000 0.00000 0.552558 +0.500000 0.00000 0.250000 0.239642 +0.750000 0.00000 0.250000 0.550251 +0.750000 0.00000 0.00000 0.552558 +0.500000 0.250000 0.00000 0.241213 +0.750000 0.00000 0.250000 0.550251 + +0.500000 0.250000 0.00000 0.241213 +0.500000 0.00000 0.250000 0.239642 + +0.500000 0.250000 0.00000 0.241213 +0.500000 0.250000 0.250000 0.238578 +0.500000 0.00000 0.250000 0.239642 +0.500000 0.250000 0.00000 0.241213 +0.750000 0.00000 0.250000 0.550251 +0.500000 0.00000 0.250000 0.239642 + +0.750000 0.00000 0.250000 0.550251 +0.500000 0.250000 0.250000 0.238578 + +0.750000 0.250000 0.00000 0.552908 +0.750000 0.00000 0.250000 0.550251 +0.500000 0.250000 0.250000 0.238578 +0.750000 0.250000 0.00000 0.552908 +0.500000 0.250000 0.00000 0.241213 +0.500000 0.250000 0.250000 0.238578 + +0.500000 0.250000 0.00000 0.241213 +0.750000 0.00000 0.250000 0.550251 + +0.00000 0.500000 0.00000 -0.0103118 +0.250000 0.500000 0.00000 0.0521580 +0.00000 0.750000 0.00000 -0.0110003 +0.00000 0.500000 0.00000 -0.0103118 +0.00000 0.500000 0.250000 -0.0109271 +0.00000 0.750000 0.00000 -0.0110003 + +0.00000 0.500000 0.250000 -0.0109271 +0.250000 0.500000 0.00000 0.0521580 + +0.250000 0.500000 0.00000 0.0521580 +0.500000 0.500000 0.00000 0.239596 +0.250000 0.750000 0.00000 0.0511859 +0.250000 0.500000 0.00000 0.0521580 +0.250000 0.500000 0.250000 0.0518325 +0.250000 0.750000 0.00000 0.0511859 + +0.250000 0.500000 0.250000 0.0518325 +0.500000 0.500000 0.00000 0.239596 + +0.00000 0.750000 0.00000 -0.0110003 +0.250000 0.750000 0.00000 0.0511859 +0.00000 1.00000 0.00000 -0.00565835 +0.00000 0.750000 0.00000 -0.0110003 +0.00000 0.750000 0.250000 -0.00803550 +0.00000 1.00000 0.00000 -0.00565835 + +0.00000 0.750000 0.250000 -0.00803550 +0.250000 0.750000 0.00000 0.0511859 + +0.00000 0.500000 0.250000 -0.0109271 +0.250000 0.500000 0.250000 0.0518325 +0.00000 0.750000 0.250000 -0.00803550 +0.00000 0.500000 0.250000 -0.0109271 +0.00000 0.500000 0.500000 -0.00807187 +0.00000 0.750000 0.250000 -0.00803550 + +0.00000 0.500000 0.500000 -0.00807187 +0.250000 0.500000 0.250000 0.0518325 + +0.250000 0.500000 0.00000 0.0521580 +0.250000 0.750000 0.00000 0.0511859 +0.00000 0.750000 0.00000 -0.0110003 +0.250000 0.500000 0.00000 0.0521580 +0.250000 0.500000 0.250000 0.0518325 +0.00000 0.750000 0.00000 -0.0110003 + +0.250000 0.500000 0.250000 0.0518325 +0.250000 0.750000 0.00000 0.0511859 + +0.250000 0.500000 0.00000 0.0521580 +0.00000 0.500000 0.250000 -0.0109271 +0.250000 0.500000 0.250000 0.0518325 +0.250000 0.500000 0.00000 0.0521580 +0.00000 0.750000 0.00000 -0.0110003 +0.250000 0.500000 0.250000 0.0518325 + +0.00000 0.750000 0.00000 -0.0110003 +0.00000 0.500000 0.250000 -0.0109271 + +0.00000 0.750000 0.00000 -0.0110003 +0.00000 0.750000 0.250000 -0.00803550 +0.00000 0.500000 0.250000 -0.0109271 +0.00000 0.750000 0.00000 -0.0110003 +0.250000 0.500000 0.250000 0.0518325 +0.00000 0.500000 0.250000 -0.0109271 + +0.250000 0.500000 0.250000 0.0518325 +0.00000 0.750000 0.250000 -0.00803550 + +0.250000 0.750000 0.00000 0.0511859 +0.250000 0.500000 0.250000 0.0518325 +0.00000 0.750000 0.250000 -0.00803550 +0.250000 0.750000 0.00000 0.0511859 +0.00000 0.750000 0.00000 -0.0110003 +0.00000 0.750000 0.250000 -0.00803550 + +0.00000 0.750000 0.00000 -0.0110003 +0.250000 0.500000 0.250000 0.0518325 + +0.00000 0.00000 0.500000 -0.00809848 +0.250000 0.00000 0.500000 0.0520679 +0.00000 0.250000 0.500000 -0.00930138 +0.00000 0.00000 0.500000 -0.00809848 +0.00000 0.00000 0.750000 -0.00846907 +0.00000 0.250000 0.500000 -0.00930138 + +0.00000 0.00000 0.750000 -0.00846907 +0.250000 0.00000 0.500000 0.0520679 + +0.250000 0.00000 0.500000 0.0520679 +0.500000 0.00000 0.500000 0.238105 +0.250000 0.250000 0.500000 0.0511676 +0.250000 0.00000 0.500000 0.0520679 +0.250000 0.00000 0.750000 0.0488913 +0.250000 0.250000 0.500000 0.0511676 + +0.250000 0.00000 0.750000 0.0488913 +0.500000 0.00000 0.500000 0.238105 + +0.00000 0.250000 0.500000 -0.00930138 +0.250000 0.250000 0.500000 0.0511676 +0.00000 0.500000 0.500000 -0.00807187 +0.00000 0.250000 0.500000 -0.00930138 +0.00000 0.250000 0.750000 -0.00817340 +0.00000 0.500000 0.500000 -0.00807187 + +0.00000 0.250000 0.750000 -0.00817340 +0.250000 0.250000 0.500000 0.0511676 + +0.00000 0.00000 0.750000 -0.00846907 +0.250000 0.00000 0.750000 0.0488913 +0.00000 0.250000 0.750000 -0.00817340 +0.00000 0.00000 0.750000 -0.00846907 +0.00000 0.00000 1.00000 -0.00570776 +0.00000 0.250000 0.750000 -0.00817340 + +0.00000 0.00000 1.00000 -0.00570776 +0.250000 0.00000 0.750000 0.0488913 + +0.250000 0.00000 0.500000 0.0520679 +0.250000 0.250000 0.500000 0.0511676 +0.00000 0.250000 0.500000 -0.00930138 +0.250000 0.00000 0.500000 0.0520679 +0.250000 0.00000 0.750000 0.0488913 +0.00000 0.250000 0.500000 -0.00930138 + +0.250000 0.00000 0.750000 0.0488913 +0.250000 0.250000 0.500000 0.0511676 + +0.250000 0.00000 0.500000 0.0520679 +0.00000 0.00000 0.750000 -0.00846907 +0.250000 0.00000 0.750000 0.0488913 +0.250000 0.00000 0.500000 0.0520679 +0.00000 0.250000 0.500000 -0.00930138 +0.250000 0.00000 0.750000 0.0488913 + +0.00000 0.250000 0.500000 -0.00930138 +0.00000 0.00000 0.750000 -0.00846907 + +0.00000 0.250000 0.500000 -0.00930138 +0.00000 0.250000 0.750000 -0.00817340 +0.00000 0.00000 0.750000 -0.00846907 +0.00000 0.250000 0.500000 -0.00930138 +0.250000 0.00000 0.750000 0.0488913 +0.00000 0.00000 0.750000 -0.00846907 + +0.250000 0.00000 0.750000 0.0488913 +0.00000 0.250000 0.750000 -0.00817340 + +0.250000 0.250000 0.500000 0.0511676 +0.250000 0.00000 0.750000 0.0488913 +0.00000 0.250000 0.750000 -0.00817340 +0.250000 0.250000 0.500000 0.0511676 +0.00000 0.250000 0.500000 -0.00930138 +0.00000 0.250000 0.750000 -0.00817340 + +0.00000 0.250000 0.500000 -0.00930138 +0.250000 0.00000 0.750000 0.0488913 + +0.500000 0.00000 0.00000 0.239570 +0.500000 0.250000 0.00000 0.241213 +0.250000 0.250000 0.00000 0.0511687 +0.500000 0.00000 0.00000 0.239570 +0.500000 0.00000 0.250000 0.239642 +0.250000 0.250000 0.00000 0.0511687 + +0.500000 0.00000 0.250000 0.239642 +0.500000 0.250000 0.00000 0.241213 + +0.500000 0.250000 0.00000 0.241213 +0.500000 0.500000 0.00000 0.239596 +0.250000 0.500000 0.00000 0.0521580 +0.500000 0.250000 0.00000 0.241213 +0.500000 0.250000 0.250000 0.238578 +0.250000 0.500000 0.00000 0.0521580 + +0.500000 0.250000 0.250000 0.238578 +0.500000 0.500000 0.00000 0.239596 + +0.250000 0.250000 0.00000 0.0511687 +0.250000 0.500000 0.00000 0.0521580 +0.00000 0.500000 0.00000 -0.0103118 +0.250000 0.250000 0.00000 0.0511687 +0.250000 0.250000 0.250000 0.0539741 +0.00000 0.500000 0.00000 -0.0103118 + +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.500000 0.00000 0.0521580 + +0.500000 0.00000 0.250000 0.239642 +0.500000 0.250000 0.250000 0.238578 +0.250000 0.250000 0.250000 0.0539741 +0.500000 0.00000 0.250000 0.239642 +0.500000 0.00000 0.500000 0.238105 +0.250000 0.250000 0.250000 0.0539741 + +0.500000 0.00000 0.500000 0.238105 +0.500000 0.250000 0.250000 0.238578 + +0.500000 0.250000 0.00000 0.241213 +0.250000 0.500000 0.00000 0.0521580 +0.250000 0.250000 0.00000 0.0511687 +0.500000 0.250000 0.00000 0.241213 +0.500000 0.250000 0.250000 0.238578 +0.250000 0.250000 0.00000 0.0511687 + +0.500000 0.250000 0.250000 0.238578 +0.250000 0.500000 0.00000 0.0521580 + +0.500000 0.250000 0.00000 0.241213 +0.500000 0.00000 0.250000 0.239642 +0.500000 0.250000 0.250000 0.238578 +0.500000 0.250000 0.00000 0.241213 +0.250000 0.250000 0.00000 0.0511687 +0.500000 0.250000 0.250000 0.238578 + +0.250000 0.250000 0.00000 0.0511687 +0.500000 0.00000 0.250000 0.239642 + +0.250000 0.250000 0.00000 0.0511687 +0.250000 0.250000 0.250000 0.0539741 +0.500000 0.00000 0.250000 0.239642 +0.250000 0.250000 0.00000 0.0511687 +0.500000 0.250000 0.250000 0.238578 +0.500000 0.00000 0.250000 0.239642 + +0.500000 0.250000 0.250000 0.238578 +0.250000 0.250000 0.250000 0.0539741 + +0.250000 0.500000 0.00000 0.0521580 +0.500000 0.250000 0.250000 0.238578 +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.500000 0.00000 0.0521580 +0.250000 0.250000 0.00000 0.0511687 +0.250000 0.250000 0.250000 0.0539741 + +0.250000 0.250000 0.00000 0.0511687 +0.500000 0.250000 0.250000 0.238578 + +0.500000 0.00000 0.00000 0.239570 +0.250000 0.00000 0.250000 0.0519658 +0.500000 0.00000 0.250000 0.239642 +0.500000 0.00000 0.00000 0.239570 +0.250000 0.250000 0.00000 0.0511687 +0.500000 0.00000 0.250000 0.239642 + +0.250000 0.250000 0.00000 0.0511687 +0.250000 0.00000 0.250000 0.0519658 + +0.250000 0.00000 0.250000 0.0519658 +0.00000 0.00000 0.500000 -0.00809848 +0.250000 0.00000 0.500000 0.0520679 +0.250000 0.00000 0.250000 0.0519658 +0.00000 0.250000 0.250000 -0.0117534 +0.250000 0.00000 0.500000 0.0520679 + +0.00000 0.250000 0.250000 -0.0117534 +0.00000 0.00000 0.500000 -0.00809848 + +0.500000 0.00000 0.250000 0.239642 +0.250000 0.00000 0.500000 0.0520679 +0.500000 0.00000 0.500000 0.238105 +0.500000 0.00000 0.250000 0.239642 +0.250000 0.250000 0.250000 0.0539741 +0.500000 0.00000 0.500000 0.238105 + +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.00000 0.500000 0.0520679 + +0.250000 0.250000 0.00000 0.0511687 +0.00000 0.250000 0.250000 -0.0117534 +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.250000 0.00000 0.0511687 +0.00000 0.500000 0.00000 -0.0103118 +0.250000 0.250000 0.250000 0.0539741 + +0.00000 0.500000 0.00000 -0.0103118 +0.00000 0.250000 0.250000 -0.0117534 + +0.250000 0.00000 0.250000 0.0519658 +0.250000 0.00000 0.500000 0.0520679 +0.500000 0.00000 0.250000 0.239642 +0.250000 0.00000 0.250000 0.0519658 +0.250000 0.250000 0.250000 0.0539741 +0.500000 0.00000 0.250000 0.239642 + +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.00000 0.500000 0.0520679 + +0.250000 0.00000 0.250000 0.0519658 +0.250000 0.250000 0.00000 0.0511687 +0.00000 0.250000 0.250000 -0.0117534 +0.250000 0.00000 0.250000 0.0519658 +0.250000 0.250000 0.250000 0.0539741 +0.00000 0.250000 0.250000 -0.0117534 + +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.250000 0.00000 0.0511687 + +0.500000 0.00000 0.250000 0.239642 +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.250000 0.00000 0.0511687 +0.500000 0.00000 0.250000 0.239642 +0.250000 0.00000 0.250000 0.0519658 +0.250000 0.250000 0.00000 0.0511687 + +0.250000 0.00000 0.250000 0.0519658 +0.250000 0.250000 0.250000 0.0539741 + +0.250000 0.00000 0.500000 0.0520679 +0.00000 0.250000 0.250000 -0.0117534 +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.00000 0.500000 0.0520679 +0.250000 0.00000 0.250000 0.0519658 +0.250000 0.250000 0.250000 0.0539741 + +0.250000 0.00000 0.250000 0.0519658 +0.00000 0.250000 0.250000 -0.0117534 + +0.00000 0.500000 0.00000 -0.0103118 +0.00000 0.500000 0.250000 -0.0109271 +0.00000 0.250000 0.250000 -0.0117534 +0.00000 0.500000 0.00000 -0.0103118 +0.250000 0.250000 0.250000 0.0539741 +0.00000 0.250000 0.250000 -0.0117534 + +0.250000 0.250000 0.250000 0.0539741 +0.00000 0.500000 0.250000 -0.0109271 + +0.00000 0.500000 0.250000 -0.0109271 +0.00000 0.500000 0.500000 -0.00807187 +0.00000 0.250000 0.500000 -0.00930138 +0.00000 0.500000 0.250000 -0.0109271 +0.250000 0.250000 0.500000 0.0511676 +0.00000 0.250000 0.500000 -0.00930138 + +0.250000 0.250000 0.500000 0.0511676 +0.00000 0.500000 0.500000 -0.00807187 + +0.00000 0.250000 0.250000 -0.0117534 +0.00000 0.250000 0.500000 -0.00930138 +0.00000 0.00000 0.500000 -0.00809848 +0.00000 0.250000 0.250000 -0.0117534 +0.250000 0.00000 0.500000 0.0520679 +0.00000 0.00000 0.500000 -0.00809848 + +0.250000 0.00000 0.500000 0.0520679 +0.00000 0.250000 0.500000 -0.00930138 + +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.250000 0.500000 0.0511676 +0.250000 0.00000 0.500000 0.0520679 +0.250000 0.250000 0.250000 0.0539741 +0.500000 0.00000 0.500000 0.238105 +0.250000 0.00000 0.500000 0.0520679 + +0.500000 0.00000 0.500000 0.238105 +0.250000 0.250000 0.500000 0.0511676 + +0.00000 0.500000 0.250000 -0.0109271 +0.00000 0.250000 0.500000 -0.00930138 +0.00000 0.250000 0.250000 -0.0117534 +0.00000 0.500000 0.250000 -0.0109271 +0.250000 0.250000 0.500000 0.0511676 +0.00000 0.250000 0.250000 -0.0117534 + +0.250000 0.250000 0.500000 0.0511676 +0.00000 0.250000 0.500000 -0.00930138 + +0.00000 0.500000 0.250000 -0.0109271 +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.250000 0.500000 0.0511676 +0.00000 0.500000 0.250000 -0.0109271 +0.00000 0.250000 0.250000 -0.0117534 +0.250000 0.250000 0.500000 0.0511676 + +0.00000 0.250000 0.250000 -0.0117534 +0.250000 0.250000 0.250000 0.0539741 + +0.00000 0.250000 0.250000 -0.0117534 +0.250000 0.00000 0.500000 0.0520679 +0.250000 0.250000 0.250000 0.0539741 +0.00000 0.250000 0.250000 -0.0117534 +0.250000 0.250000 0.500000 0.0511676 +0.250000 0.250000 0.250000 0.0539741 + +0.250000 0.250000 0.500000 0.0511676 +0.250000 0.00000 0.500000 0.0520679 + +0.00000 0.250000 0.500000 -0.00930138 +0.250000 0.250000 0.500000 0.0511676 +0.250000 0.00000 0.500000 0.0520679 +0.00000 0.250000 0.500000 -0.00930138 +0.00000 0.250000 0.250000 -0.0117534 +0.250000 0.00000 0.500000 0.0520679 + +0.00000 0.250000 0.250000 -0.0117534 +0.250000 0.250000 0.500000 0.0511676 + +0.500000 0.500000 0.00000 0.239596 +0.500000 0.250000 0.250000 0.238578 +0.250000 0.500000 0.250000 0.0518325 +0.500000 0.500000 0.00000 0.239596 +0.250000 0.500000 0.00000 0.0521580 +0.250000 0.500000 0.250000 0.0518325 + +0.250000 0.500000 0.00000 0.0521580 +0.500000 0.250000 0.250000 0.238578 + +0.500000 0.250000 0.250000 0.238578 +0.500000 0.00000 0.500000 0.238105 +0.250000 0.250000 0.500000 0.0511676 +0.500000 0.250000 0.250000 0.238578 +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.250000 0.500000 0.0511676 + +0.250000 0.250000 0.250000 0.0539741 +0.500000 0.00000 0.500000 0.238105 + +0.250000 0.500000 0.250000 0.0518325 +0.250000 0.250000 0.500000 0.0511676 +0.00000 0.500000 0.500000 -0.00807187 +0.250000 0.500000 0.250000 0.0518325 +0.00000 0.500000 0.250000 -0.0109271 +0.00000 0.500000 0.500000 -0.00807187 + +0.00000 0.500000 0.250000 -0.0109271 +0.250000 0.250000 0.500000 0.0511676 + +0.250000 0.500000 0.00000 0.0521580 +0.250000 0.250000 0.250000 0.0539741 +0.00000 0.500000 0.250000 -0.0109271 +0.250000 0.500000 0.00000 0.0521580 +0.00000 0.500000 0.00000 -0.0103118 +0.00000 0.500000 0.250000 -0.0109271 + +0.00000 0.500000 0.00000 -0.0103118 +0.250000 0.250000 0.250000 0.0539741 + +0.500000 0.250000 0.250000 0.238578 +0.250000 0.250000 0.500000 0.0511676 +0.250000 0.500000 0.250000 0.0518325 +0.500000 0.250000 0.250000 0.238578 +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.500000 0.250000 0.0518325 + +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.250000 0.500000 0.0511676 + +0.500000 0.250000 0.250000 0.238578 +0.250000 0.500000 0.00000 0.0521580 +0.250000 0.250000 0.250000 0.0539741 +0.500000 0.250000 0.250000 0.238578 +0.250000 0.500000 0.250000 0.0518325 +0.250000 0.250000 0.250000 0.0539741 + +0.250000 0.500000 0.250000 0.0518325 +0.250000 0.500000 0.00000 0.0521580 + +0.250000 0.500000 0.250000 0.0518325 +0.00000 0.500000 0.250000 -0.0109271 +0.250000 0.500000 0.00000 0.0521580 +0.250000 0.500000 0.250000 0.0518325 +0.250000 0.250000 0.250000 0.0539741 +0.250000 0.500000 0.00000 0.0521580 + +0.250000 0.250000 0.250000 0.0539741 +0.00000 0.500000 0.250000 -0.0109271 + +0.250000 0.250000 0.500000 0.0511676 +0.250000 0.250000 0.250000 0.0539741 +0.00000 0.500000 0.250000 -0.0109271 +0.250000 0.250000 0.500000 0.0511676 +0.250000 0.500000 0.250000 0.0518325 +0.00000 0.500000 0.250000 -0.0109271 + +0.250000 0.500000 0.250000 0.0518325 +0.250000 0.250000 0.250000 0.0539741 diff --git a/tests/simplex/refinement_04.cc b/tests/simplex/refinement_04.cc index 30bc0aecc8..4ea27d7c38 100644 --- a/tests/simplex/refinement_04.cc +++ b/tests/simplex/refinement_04.cc @@ -1,17 +1,16 @@ -// --------------------------------------------------------------------- +// ------------------------------------------------------------------------ // -// Copyright (C) 2020 - 2021 by the deal.II authors +// SPDX-License-Identifier: LGPL-2.1-or-later +// Copyright (C) 2021 by the deal.II authors // // This file is part of the deal.II library. // -// The deal.II library is free software; you can use it, redistribute -// it, and/or modify it under the terms of the GNU Lesser General -// Public License as published by the Free Software Foundation; either -// version 2.1 of the License, or (at your option) any later version. -// The full text of the license can be found in the file LICENSE.md at -// the top level directory of deal.II. +// Part of the source code is dual licensed under Apache-2.0 WITH +// LLVM-exception OR LGPL-2.1-or-later. Detailed license information +// governing the source code and code contributions can be found in +// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II. // -// --------------------------------------------------------------------- +// ------------------------------------------------------------------------ @@ -50,13 +49,13 @@ test(const unsigned int version) if (version == 0 || version == 1) GridGenerator::reference_cell(tria, ReferenceCells::Tetrahedron); - else if (version == 2) + else if (version == 2 || version == 3) GridGenerator::subdivided_hyper_cube_with_simplices(tria, 1); print(tria, "tria.0.vtk"); if (version == 0 || version == 2) tria.refine_global(1); - else if (version == 1) + else if (version == 1 || version == 3) tria.refine_global(2); print(tria, "tria.1.vtk"); @@ -69,4 +68,5 @@ main() test(0); // coarse grid: 1 tet, n_refinements: 1 test(1); // coarse grid: 1 tet, n_refinements: 2 test(2); // coarse grid: 5 tets (cube), n_refinements: 1 + test(3); // coarse grid: 5 tets (cube), n_refinements: 2 } diff --git a/tests/simplex/refinement_04.output b/tests/simplex/refinement_04.output index f9bfb8d951..a304f8493b 100644 --- a/tests/simplex/refinement_04.output +++ b/tests/simplex/refinement_04.output @@ -188,10 +188,10 @@ CELLS 64 320 4 34 5 33 23 4 32 33 7 22 4 29 23 22 8 -4 34 33 32 23 -4 34 29 23 32 -4 32 22 29 23 -4 33 23 22 32 +4 34 33 32 22 +4 34 29 23 22 +4 32 22 29 34 +4 33 23 22 34 4 8 25 23 22 4 25 9 24 27 4 23 24 5 33 @@ -355,3 +355,501 @@ LOOKUP_TABLE default -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 +# vtk DataFile Version 3.0 +Triangulation generated with deal.II +ASCII +DATASET UNSTRUCTURED_GRID +POINTS 8 double +0.00000 0.00000 0.00000 +1.00000 0.00000 0.00000 +0.00000 1.00000 0.00000 +1.00000 1.00000 0.00000 +0.00000 0.00000 1.00000 +1.00000 0.00000 1.00000 +0.00000 1.00000 1.00000 +1.00000 1.00000 1.00000 + +CELLS 5 25 +4 0 1 2 4 +4 2 1 3 7 +4 1 4 5 7 +4 2 4 7 6 +4 1 2 4 7 + +CELL_TYPES 5 +10 10 10 10 10 + + + +CELL_DATA 5 +SCALARS MaterialID int 1 +LOOKUP_TABLE default +0 1 2 3 4 + + + +SCALARS ManifoldID int 1 +LOOKUP_TABLE default +-1 -1 -1 -1 -1 + + +# vtk DataFile Version 3.0 +Triangulation generated with deal.II +ASCII +DATASET UNSTRUCTURED_GRID +POINTS 115 double +0.00000 0.00000 0.00000 +1.00000 0.00000 0.00000 +0.00000 1.00000 0.00000 +1.00000 1.00000 0.00000 +0.00000 0.00000 1.00000 +1.00000 0.00000 1.00000 +0.00000 1.00000 1.00000 +1.00000 1.00000 1.00000 +0.500000 0.00000 0.00000 +0.00000 0.00000 0.500000 +0.500000 0.500000 0.00000 +1.00000 0.500000 0.00000 +0.500000 0.00000 0.500000 +1.00000 0.500000 0.500000 +0.00000 0.500000 0.00000 +0.00000 0.500000 0.500000 +0.00000 1.00000 0.500000 +0.500000 1.00000 0.500000 +0.500000 1.00000 0.00000 +1.00000 1.00000 0.500000 +0.500000 0.00000 1.00000 +0.00000 0.500000 1.00000 +0.500000 0.500000 1.00000 +1.00000 0.00000 0.500000 +1.00000 0.500000 1.00000 +0.500000 1.00000 1.00000 +0.250000 0.00000 0.00000 +0.750000 0.00000 0.00000 +0.00000 0.00000 0.250000 +0.00000 0.00000 0.750000 +0.750000 0.250000 0.00000 +0.250000 0.750000 0.00000 +1.00000 0.250000 0.00000 +1.00000 0.750000 0.00000 +0.750000 0.00000 0.250000 +0.250000 0.00000 0.750000 +1.00000 0.250000 0.250000 +1.00000 0.750000 0.750000 +0.00000 0.750000 0.00000 +0.00000 0.250000 0.00000 +0.00000 0.750000 0.250000 +0.00000 0.250000 0.750000 +0.00000 1.00000 0.250000 +0.00000 1.00000 0.750000 +0.250000 1.00000 0.250000 +0.750000 1.00000 0.750000 +0.750000 1.00000 0.00000 +0.250000 1.00000 0.00000 +1.00000 1.00000 0.250000 +1.00000 1.00000 0.750000 +0.250000 0.00000 1.00000 +0.750000 0.00000 1.00000 +0.00000 0.250000 1.00000 +0.00000 0.750000 1.00000 +0.250000 0.250000 1.00000 +0.750000 0.750000 1.00000 +1.00000 0.00000 0.750000 +1.00000 0.00000 0.250000 +1.00000 0.250000 1.00000 +1.00000 0.750000 1.00000 +0.750000 1.00000 1.00000 +0.250000 1.00000 1.00000 +0.500000 0.500000 0.500000 +0.250000 0.750000 0.750000 +0.750000 0.250000 0.750000 +0.750000 0.750000 0.250000 +0.250000 0.250000 0.250000 +1.00000 0.250000 0.500000 +1.00000 0.500000 0.750000 +1.00000 0.250000 0.750000 +0.500000 0.750000 1.00000 +0.250000 0.750000 1.00000 +0.250000 0.500000 1.00000 +0.750000 0.250000 1.00000 +0.750000 0.500000 1.00000 +0.500000 0.250000 1.00000 +0.500000 1.00000 0.250000 +0.750000 1.00000 0.500000 +0.750000 1.00000 0.250000 +0.250000 1.00000 0.500000 +0.250000 1.00000 0.750000 +0.500000 1.00000 0.750000 +0.250000 0.750000 0.500000 +0.500000 0.750000 0.750000 +0.250000 0.500000 0.750000 +0.00000 0.500000 0.750000 +0.00000 0.750000 0.750000 +0.00000 0.750000 0.500000 +0.00000 0.250000 0.250000 +0.00000 0.250000 0.500000 +0.00000 0.500000 0.250000 +0.750000 0.00000 0.500000 +0.750000 0.00000 0.750000 +0.500000 0.00000 0.750000 +0.500000 0.250000 0.750000 +0.750000 0.500000 0.750000 +0.750000 0.250000 0.500000 +1.00000 0.750000 0.250000 +1.00000 0.750000 0.500000 +1.00000 0.500000 0.250000 +0.750000 0.500000 0.250000 +0.750000 0.750000 0.500000 +0.500000 0.750000 0.250000 +0.250000 0.500000 0.250000 +0.250000 0.250000 0.500000 +0.500000 0.250000 0.250000 +0.500000 0.750000 0.00000 +0.750000 0.750000 0.00000 +0.750000 0.500000 0.00000 +0.250000 0.250000 0.00000 +0.250000 0.500000 0.00000 +0.500000 0.250000 0.00000 +0.500000 0.00000 0.250000 +0.250000 0.00000 0.500000 +0.250000 0.00000 0.250000 + +CELLS 320 1600 +4 0 26 39 28 +4 26 8 109 114 +4 39 109 14 88 +4 28 114 88 9 +4 26 109 39 114 +4 26 28 114 39 +4 39 88 28 114 +4 109 114 88 39 +4 8 27 111 112 +4 27 1 30 34 +4 111 30 10 105 +4 112 34 105 12 +4 27 30 111 34 +4 27 112 34 111 +4 111 105 112 34 +4 30 34 105 111 +4 14 110 38 90 +4 110 10 31 103 +4 38 31 2 40 +4 90 103 40 15 +4 110 31 38 103 +4 110 90 103 38 +4 38 40 90 103 +4 31 103 40 38 +4 9 113 89 29 +4 113 12 104 35 +4 89 104 15 41 +4 29 35 41 4 +4 113 104 89 35 +4 113 29 35 89 +4 89 41 29 35 +4 104 35 41 89 +4 8 111 109 112 +4 111 10 110 105 +4 109 110 14 66 +4 112 105 66 12 +4 111 110 109 105 +4 111 112 105 109 +4 109 66 112 105 +4 110 105 66 109 +4 8 114 112 109 +4 114 9 113 88 +4 112 113 12 66 +4 109 88 66 14 +4 114 113 112 66 +4 114 109 88 66 +4 112 66 109 114 +4 113 88 66 114 +4 14 90 88 66 +4 90 15 89 104 +4 88 89 9 113 +4 66 104 113 12 +4 90 89 88 104 +4 90 66 104 88 +4 88 113 66 104 +4 89 104 113 88 +4 10 105 103 110 +4 105 12 104 66 +4 103 104 15 90 +4 110 66 90 14 +4 105 104 103 66 +4 105 110 66 103 +4 103 90 110 66 +4 104 66 90 103 +4 2 31 47 44 +4 31 10 106 102 +4 47 106 18 76 +4 44 102 76 17 +4 31 106 47 102 +4 31 44 102 47 +4 47 76 44 102 +4 106 102 76 47 +4 10 30 108 100 +4 30 1 32 36 +4 108 32 11 99 +4 100 36 99 13 +4 30 32 108 36 +4 30 100 36 108 +4 108 99 100 36 +4 32 36 99 108 +4 18 107 46 78 +4 107 11 33 97 +4 46 33 3 48 +4 78 97 48 19 +4 107 33 46 97 +4 107 78 97 46 +4 46 48 78 97 +4 33 97 48 46 +4 17 101 77 45 +4 101 13 98 37 +4 77 98 19 49 +4 45 37 49 7 +4 101 98 77 37 +4 101 45 37 77 +4 77 49 45 37 +4 98 37 49 77 +4 10 108 106 100 +4 108 11 107 99 +4 106 107 18 65 +4 100 99 65 13 +4 108 107 106 100 +4 108 100 99 107 +4 106 65 100 107 +4 107 99 65 100 +4 10 102 100 106 +4 102 17 101 76 +4 100 101 13 65 +4 106 76 65 18 +4 102 101 100 65 +4 102 106 76 65 +4 100 65 106 102 +4 101 76 65 102 +4 18 78 76 65 +4 78 19 77 98 +4 76 77 17 101 +4 65 98 101 13 +4 78 77 76 65 +4 78 65 98 77 +4 76 101 65 77 +4 77 98 101 65 +4 11 99 97 107 +4 99 13 98 65 +4 97 98 19 78 +4 107 65 78 18 +4 99 98 97 65 +4 99 107 65 97 +4 97 78 107 65 +4 98 65 78 97 +4 1 34 57 36 +4 34 12 91 96 +4 57 91 23 67 +4 36 96 67 13 +4 34 91 57 96 +4 34 36 96 57 +4 57 67 36 96 +4 91 96 67 57 +4 12 35 93 94 +4 35 4 50 54 +4 93 50 20 75 +4 94 54 75 22 +4 35 50 93 54 +4 35 94 54 93 +4 93 75 94 54 +4 50 54 75 93 +4 23 92 56 69 +4 92 20 51 73 +4 56 51 5 58 +4 69 73 58 24 +4 92 51 56 73 +4 92 69 73 56 +4 56 58 69 73 +4 51 73 58 56 +4 13 95 68 37 +4 95 22 74 55 +4 68 74 24 59 +4 37 55 59 7 +4 95 74 68 55 +4 95 37 55 68 +4 68 59 37 55 +4 74 55 59 68 +4 12 93 91 94 +4 93 20 92 75 +4 91 92 23 64 +4 94 75 64 22 +4 93 92 91 94 +4 93 94 75 92 +4 91 64 94 92 +4 92 75 64 94 +4 12 96 94 91 +4 96 13 95 67 +4 94 95 22 64 +4 91 67 64 23 +4 96 95 94 64 +4 96 91 67 64 +4 94 64 91 96 +4 95 67 64 96 +4 23 69 67 64 +4 69 24 68 74 +4 67 68 13 95 +4 64 74 95 22 +4 69 68 67 64 +4 69 64 74 68 +4 67 95 64 68 +4 68 74 95 64 +4 20 75 73 92 +4 75 22 74 64 +4 73 74 24 69 +4 92 64 69 23 +4 75 74 73 64 +4 75 92 64 73 +4 73 69 92 64 +4 74 64 69 73 +4 2 40 44 42 +4 40 15 82 87 +4 44 82 17 79 +4 42 87 79 16 +4 40 82 44 87 +4 40 42 87 44 +4 44 79 42 87 +4 82 87 79 44 +4 15 41 84 85 +4 41 4 54 52 +4 84 54 22 72 +4 85 52 72 21 +4 41 54 84 52 +4 41 85 52 84 +4 84 72 85 52 +4 54 52 72 84 +4 17 83 45 81 +4 83 22 55 70 +4 45 55 7 60 +4 81 70 60 25 +4 83 55 45 70 +4 83 81 70 45 +4 45 60 81 70 +4 55 70 60 45 +4 16 86 80 43 +4 86 21 71 53 +4 80 71 25 61 +4 43 53 61 6 +4 86 71 80 53 +4 86 43 53 80 +4 80 61 43 53 +4 71 53 61 80 +4 15 84 82 85 +4 84 22 83 72 +4 82 83 17 63 +4 85 72 63 21 +4 84 83 82 63 +4 84 85 72 63 +4 82 63 85 84 +4 83 72 63 84 +4 15 87 85 82 +4 87 16 86 79 +4 85 86 21 63 +4 82 79 63 17 +4 87 86 85 82 +4 87 82 79 86 +4 85 63 82 86 +4 86 79 63 82 +4 17 81 79 63 +4 81 25 80 71 +4 79 80 16 86 +4 63 71 86 21 +4 81 80 79 63 +4 81 63 71 80 +4 79 86 63 80 +4 80 71 86 63 +4 22 72 70 83 +4 72 21 71 63 +4 70 71 25 81 +4 83 63 81 17 +4 72 71 70 63 +4 72 83 63 70 +4 70 81 83 63 +4 71 63 81 70 +4 1 30 34 36 +4 30 10 105 100 +4 34 105 12 96 +4 36 100 96 13 +4 30 105 34 100 +4 30 36 100 34 +4 34 96 36 100 +4 105 100 96 34 +4 10 31 103 102 +4 31 2 40 44 +4 103 40 15 82 +4 102 44 82 17 +4 31 40 103 44 +4 31 102 44 103 +4 103 82 102 44 +4 40 44 82 103 +4 12 104 35 94 +4 104 15 41 84 +4 35 41 4 54 +4 94 84 54 22 +4 104 41 35 84 +4 104 94 84 35 +4 35 54 94 84 +4 41 84 54 35 +4 13 101 95 37 +4 101 17 83 45 +4 95 83 22 55 +4 37 45 55 7 +4 101 83 95 45 +4 101 37 45 95 +4 95 55 37 45 +4 83 45 55 95 +4 10 103 105 102 +4 103 15 104 82 +4 105 104 12 62 +4 102 82 62 17 +4 103 104 105 62 +4 103 102 82 62 +4 105 62 102 103 +4 104 82 62 103 +4 10 100 102 105 +4 100 13 101 96 +4 102 101 17 62 +4 105 96 62 12 +4 100 101 102 62 +4 100 105 96 62 +4 102 62 105 100 +4 101 96 62 100 +4 12 94 96 62 +4 94 22 95 83 +4 96 95 13 101 +4 62 83 101 17 +4 94 95 96 62 +4 94 62 83 95 +4 96 101 62 95 +4 95 83 101 62 +4 15 82 84 104 +4 82 17 83 62 +4 84 83 22 94 +4 104 62 94 12 +4 82 83 84 62 +4 82 104 62 84 +4 84 94 104 62 +4 83 62 94 84 + +CELL_TYPES 320 +10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 + + + +CELL_DATA 320 +SCALARS MaterialID int 1 +LOOKUP_TABLE default +0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 + + + +SCALARS ManifoldID int 1 +LOOKUP_TABLE default +-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 + + diff --git a/tests/simplex/refinement_05.cc b/tests/simplex/refinement_05.cc new file mode 100644 index 0000000000..99d6660be2 --- /dev/null +++ b/tests/simplex/refinement_05.cc @@ -0,0 +1,61 @@ +// ------------------------------------------------------------------------ +// +// SPDX-License-Identifier: LGPL-2.1-or-later +// Copyright (C) 2021 by the deal.II authors +// +// This file is part of the deal.II library. +// +// Part of the source code is dual licensed under Apache-2.0 WITH +// LLVM-exception OR LGPL-2.1-or-later. Detailed license information +// governing the source code and code contributions can be found in +// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II. +// +// ------------------------------------------------------------------------ + + + +// Test refinement of 3D simplex mesh. +#include + +#include +#include + +#include +#include +#include + +#include "../tests.h" + + + +void +test() +{ + constexpr int dim = 3; + Triangulation tria; + GridGenerator::subdivided_hyper_cube_with_simplices(tria, 1, 0, 1); + + for (unsigned int i = 0; i < 3; i++) + { + FE_SimplexP fe(1); + MappingFE mapping(fe); + QGaussSimplex quad(3); + + const Vector aspect_ratios = + GridTools::compute_aspect_ratio_of_cells(mapping, tria, quad); + deallog << "Number of refinements: " << i << std::endl; + deallog << "aspect ratios: "; + for (unsigned int j = 0; j < aspect_ratios.size(); j++) + deallog << aspect_ratios[j] << " "; + deallog << std::endl; + if (i < 2) + tria.refine_global(1); + } +} + +int +main() +{ + initlog(); + test(); +} diff --git a/tests/simplex/refinement_05.output b/tests/simplex/refinement_05.output new file mode 100644 index 0000000000..e2d7b900c6 --- /dev/null +++ b/tests/simplex/refinement_05.output @@ -0,0 +1,7 @@ + +DEAL::Number of refinements: 0 +DEAL::aspect ratios: 1.00000 3.73205 3.73205 3.73205 2.00000 +DEAL::Number of refinements: 1 +DEAL::aspect ratios: 1.00000 1.00000 1.00000 1.00000 2.61803 4.04892 4.04892 4.04892 3.73205 3.73205 3.73205 3.73205 2.61803 4.04892 4.04892 4.04892 3.73205 3.73205 3.73205 3.73205 2.61803 4.04892 4.04892 4.04892 3.73205 3.73205 3.73205 3.73205 4.04892 2.61803 3.73205 2.61803 2.00000 2.00000 2.00000 2.00000 2.41421 2.41421 3.18614 2.41421 +DEAL::Number of refinements: 2 +DEAL::aspect ratios: 1.00000 1.00000 1.00000 1.00000 2.61803 4.04892 4.04892 4.04892 1.00000 1.00000 1.00000 1.00000 2.61803 4.04892 4.04892 4.04892 1.00000 1.00000 1.00000 1.00000 2.61803 4.04892 4.04892 4.04892 1.00000 1.00000 1.00000 1.00000 2.61803 4.04892 4.04892 4.04892 2.61803 2.61803 2.61803 2.61803 2.61803 2.61803 4.04892 4.04892 4.04892 4.04892 4.04892 4.04892 1.00000 3.73205 4.04892 4.04892 4.04892 4.04892 4.04892 4.04892 4.04892 4.04892 4.04892 2.61803 4.04892 4.04892 4.04892 4.04892 3.73205 4.04892 1.00000 4.04892 3.73205 3.73205 3.73205 3.73205 2.61803 4.04892 4.04892 4.04892 3.73205 3.73205 3.73205 3.73205 2.61803 4.04892 4.04892 4.04892 3.73205 3.73205 3.73205 3.73205 2.61803 4.04892 4.04892 4.04892 3.73205 3.73205 3.73205 3.73205 2.61803 4.04892 4.04892 4.04892 2.61803 2.61803 2.61803 2.61803 2.61803 2.61803 4.04892 4.04892 4.04892 4.04892 4.04892 4.04892 3.73205 1.00000 4.04892 4.04892 4.04892 4.04892 4.04892 4.04892 1.00000 3.73205 5.82843 3.73205 4.04892 4.04892 4.04892 4.04892 3.73205 4.04892 3.73205 4.04892 3.73205 3.73205 3.73205 3.73205 2.61803 4.04892 4.04892 4.04892 3.73205 3.73205 3.73205 3.73205 2.61803 4.04892 4.04892 4.04892 3.73205 3.73205 3.73205 3.73205 2.61803 4.04892 4.04892 4.04892 3.73205 3.73205 3.73205 3.73205 2.61803 4.04892 4.04892 4.04892 2.61803 2.61803 2.61803 2.61803 2.61803 2.61803 4.04892 4.04892 4.04892 4.04892 4.04892 4.04892 3.73205 1.00000 4.04892 4.04892 4.04892 4.04892 4.04892 4.04892 1.00000 3.73205 5.82843 3.73205 4.04892 4.04892 4.04892 4.04892 3.73205 4.04892 3.73205 4.04892 3.73205 3.73205 3.73205 3.73205 4.04892 2.61803 3.73205 2.61803 3.73205 3.73205 3.73205 3.73205 4.04892 2.61803 3.73205 2.61803 3.73205 3.73205 3.73205 3.73205 4.04892 2.61803 3.73205 2.61803 3.73205 3.73205 3.73205 3.73205 4.04892 2.61803 3.73205 2.61803 4.04892 4.04892 4.04892 4.04892 3.73205 1.00000 4.04892 4.04892 2.61803 2.61803 2.61803 2.61803 2.61803 2.61803 4.04892 4.04892 3.73205 3.73205 3.73205 3.73205 3.73205 4.04892 4.04892 3.73205 2.61803 2.61803 2.61803 2.61803 3.73205 5.82843 3.73205 3.73205 2.00000 2.00000 2.00000 2.00000 2.41421 2.41421 3.18614 2.41421 2.00000 2.00000 2.00000 2.00000 2.41421 2.41421 3.18614 2.41421 2.00000 2.00000 2.00000 2.00000 2.41421 2.41421 3.18614 2.41421 2.00000 2.00000 2.00000 2.00000 2.41421 2.41421 3.18614 2.41421 2.41421 2.41421 2.41421 2.41421 2.00000 2.00000 3.18614 3.18614 2.41421 2.41421 2.41421 2.41421 2.00000 2.00000 3.18614 3.18614 3.18614 3.18614 3.18614 3.18614 2.00000 3.18614 3.18614 2.00000 2.41421 2.41421 2.41421 2.41421 2.00000 3.18614 2.00000 3.18614