From: Timo Heister Date: Tue, 13 Oct 2015 18:01:45 +0000 (-0400) Subject: Fix subdivided_parallelepiped X-Git-Tag: v8.4.0-rc2~313^2~3 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=82df1380c5311f555499cbb0ade558f6c64824ce;p=dealii.git Fix subdivided_parallelepiped - missing reorder_cells() would generate unconnected cells - faces had wrong boundary indicators (if colorizing) --- diff --git a/source/grid/grid_generator.cc b/source/grid/grid_generator.cc index b0ce6ad61d..2173a50d25 100644 --- a/source/grid/grid_generator.cc +++ b/source/grid/grid_generator.cc @@ -983,11 +983,26 @@ namespace GridGenerator } // Create triangulation + // reorder the cells to ensure that they satisfy the convention for + // edge and face directions + GridReordering::reorder_cells(cells, true); tria.create_triangulation (points, cells, SubCellData()); // Finally assign boundary indicators according to hyper_rectangle if (colorize) - colorize_hyper_rectangle (tria); + { + typename Triangulation::active_cell_iterator + cell = tria.begin_active(), + endc = tria.end(); + for (; cell!=endc; ++cell) + { + for (unsigned int face=0; face::faces_per_cell; ++face) + { + if (cell->face(face)->at_boundary()) + cell->face(face)->set_boundary_id(face); + } + } + } } diff --git a/tests/grid/grid_parallelepiped_04.cc b/tests/grid/grid_parallelepiped_04.cc new file mode 100644 index 0000000000..d672c2634b --- /dev/null +++ b/tests/grid/grid_parallelepiped_04.cc @@ -0,0 +1,168 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2013 - 2014 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 at +// the top level of the deal.II distribution. +// +// --------------------------------------------------------------------- + +// check boundary indicators of colored subdivided_parallelepiped + +#include "../tests.h" +#include +#include +#include +#include +#include +#include +#include +#include + +#include +#include +#include +#include +#include +#include +#include +#include + + + +#include +#include + +#include +#include +#include +#include + +template +void check_parallelepiped (bool colorize, bool log, const unsigned int (&subd)[dim]) +{ + deallog << "* checking dim=" << dim + << " subd="; + for (unsigned int i=0; i (corners) [dim]; + + // build corners for this particular dim: + switch (dim) + { + case 1: + corners[0] = Point (0.5); + break; + + case 2: + corners[0] = Point (0.0, 0.5); + corners[1] = Point (0.5, 0.0); + break; + + case 3: + corners[0] = Point (0.0, 0.3, 0.5); + corners[1] = Point (0.4, 0.0, 0.5); + corners[2] = Point (0.4, 0.3, 0.0); + break; + + default: + Assert (false, ExcInternalError ()); + } + + Triangulation triangulation; + + GridGenerator::subdivided_parallelepiped (triangulation, subd, corners, colorize); + + { + + std::map boundary_count; + typename Triangulation::active_cell_iterator + cell = triangulation.begin_active(), + endc = triangulation.end(); + for (; cell!=endc; ++cell) + { + for (unsigned int face=0; face::faces_per_cell; ++face) + { + if (cell->face(face)->at_boundary()) + { + boundary_count[cell->face(face)->boundary_id()]++; + deallog << " center: " << cell->face(face)->center() + << " id: " << (int)cell->face(face)->boundary_id() + << std::endl; + } + + } + } + + deallog << " boundary indicators: "; + for (std::map::iterator it=boundary_count.begin(); + it!=boundary_count.end(); + ++it) + { + deallog << it->first << "(" << it->second << " times) "; + } + deallog << std::endl; + } + + GridOut grid_out; + + if (log) + { + FE_Q fe(2); + DoFHandler dh (triangulation); + dh.distribute_dofs (fe); + DataOut d_o; + d_o.attach_dof_handler (dh); + Vector vec (dh.n_dofs()); + ConstraintMatrix constraints; + for (unsigned int c=0; c<6; ++c) + VectorTools::interpolate_boundary_values (dh, + c, + ConstantFunction(c), + constraints); + constraints.close(); + constraints.distribute(vec); + + d_o.add_data_vector (vec, "v"); + d_o.build_patches (2); + char fname[]="0.vtk"; + static int counter=0; + fname[0]+=counter; + ++counter; + + std::ofstream ss(fname); + + // d_o.write_vtk(ss); + d_o.write_vtk(deallog.get_file_stream()); + } + +} + +int main () +{ + initlog(true); + + //check_parallelepiped<1> (false, true); + //check_parallelepiped<2> (false, true); + for (unsigned int subd=1; subd<=3; ++subd) + { + unsigned int subdivisions[3]= {subd,subd,subd}; + check_parallelepiped<3> (true, subd==2, subdivisions); + } + for (unsigned int subd=1; subd<=3; ++subd) + { + unsigned int subdivisions[3]= {1,2,subd}; + check_parallelepiped<3> (true, false, subdivisions); + } + +} diff --git a/tests/grid/grid_parallelepiped_04.output b/tests/grid/grid_parallelepiped_04.output new file mode 100644 index 0000000000..6975d358b1 --- /dev/null +++ b/tests/grid/grid_parallelepiped_04.output @@ -0,0 +1,440 @@ + +DEAL::* checking dim=3 subd=1 1 1 +DEAL:: center: 0.400000 0.150000 0.250000 id: 0 +DEAL:: center: 0.400000 0.450000 0.750000 id: 1 +DEAL:: center: 0.200000 0.300000 0.250000 id: 2 +DEAL:: center: 0.600000 0.300000 0.750000 id: 3 +DEAL:: center: 0.200000 0.150000 0.500000 id: 4 +DEAL:: center: 0.600000 0.450000 0.500000 id: 5 +DEAL:: boundary indicators: 0(1 times) 1(1 times) 2(1 times) 3(1 times) 4(1 times) 5(1 times) +DEAL::* checking dim=3 subd=2 2 2 +DEAL:: center: 0.200000 0.0750000 0.125000 id: 0 +DEAL:: center: 0.100000 0.150000 0.125000 id: 2 +DEAL:: center: 0.100000 0.0750000 0.250000 id: 4 +DEAL:: center: 0.200000 0.375000 0.625000 id: 1 +DEAL:: center: 0.100000 0.300000 0.375000 id: 2 +DEAL:: center: 0.100000 0.225000 0.500000 id: 4 +DEAL:: center: 0.400000 0.0750000 0.375000 id: 0 +DEAL:: center: 0.500000 0.150000 0.625000 id: 3 +DEAL:: center: 0.300000 0.0750000 0.500000 id: 4 +DEAL:: center: 0.400000 0.375000 0.875000 id: 1 +DEAL:: center: 0.500000 0.300000 0.875000 id: 3 +DEAL:: center: 0.300000 0.225000 0.750000 id: 4 +DEAL:: center: 0.400000 0.225000 0.125000 id: 0 +DEAL:: center: 0.300000 0.300000 0.125000 id: 2 +DEAL:: center: 0.500000 0.375000 0.250000 id: 5 +DEAL:: center: 0.400000 0.525000 0.625000 id: 1 +DEAL:: center: 0.300000 0.450000 0.375000 id: 2 +DEAL:: center: 0.500000 0.525000 0.500000 id: 5 +DEAL:: center: 0.600000 0.225000 0.375000 id: 0 +DEAL:: center: 0.700000 0.300000 0.625000 id: 3 +DEAL:: center: 0.700000 0.375000 0.500000 id: 5 +DEAL:: center: 0.600000 0.525000 0.875000 id: 1 +DEAL:: center: 0.700000 0.450000 0.875000 id: 3 +DEAL:: center: 0.700000 0.525000 0.750000 id: 5 +DEAL:: boundary indicators: 0(4 times) 1(4 times) 2(4 times) 3(4 times) 4(4 times) 5(4 times) +# vtk DataFile Version 3.0 +#This file was generated +ASCII +DATASET UNSTRUCTURED_GRID + +POINTS 216 double +0.00000 0.00000 0.00000 +0.00000 0.0750000 0.125000 +0.00000 0.150000 0.250000 +0.100000 0.00000 0.125000 +0.100000 0.0750000 0.250000 +0.100000 0.150000 0.375000 +0.200000 0.00000 0.250000 +0.200000 0.0750000 0.375000 +0.200000 0.150000 0.500000 +0.100000 0.0750000 0.00000 +0.100000 0.150000 0.125000 +0.100000 0.225000 0.250000 +0.200000 0.0750000 0.125000 +0.200000 0.150000 0.250000 +0.200000 0.225000 0.375000 +0.300000 0.0750000 0.250000 +0.300000 0.150000 0.375000 +0.300000 0.225000 0.500000 +0.200000 0.150000 0.00000 +0.200000 0.225000 0.125000 +0.200000 0.300000 0.250000 +0.300000 0.150000 0.125000 +0.300000 0.225000 0.250000 +0.300000 0.300000 0.375000 +0.400000 0.150000 0.250000 +0.400000 0.225000 0.375000 +0.400000 0.300000 0.500000 +0.00000 0.150000 0.250000 +0.00000 0.225000 0.375000 +0.00000 0.300000 0.500000 +0.100000 0.150000 0.375000 +0.100000 0.225000 0.500000 +0.100000 0.300000 0.625000 +0.200000 0.150000 0.500000 +0.200000 0.225000 0.625000 +0.200000 0.300000 0.750000 +0.100000 0.225000 0.250000 +0.100000 0.300000 0.375000 +0.100000 0.375000 0.500000 +0.200000 0.225000 0.375000 +0.200000 0.300000 0.500000 +0.200000 0.375000 0.625000 +0.300000 0.225000 0.500000 +0.300000 0.300000 0.625000 +0.300000 0.375000 0.750000 +0.200000 0.300000 0.250000 +0.200000 0.375000 0.375000 +0.200000 0.450000 0.500000 +0.300000 0.300000 0.375000 +0.300000 0.375000 0.500000 +0.300000 0.450000 0.625000 +0.400000 0.300000 0.500000 +0.400000 0.375000 0.625000 +0.400000 0.450000 0.750000 +0.200000 0.00000 0.250000 +0.200000 0.0750000 0.375000 +0.200000 0.150000 0.500000 +0.300000 0.00000 0.375000 +0.300000 0.0750000 0.500000 +0.300000 0.150000 0.625000 +0.400000 0.00000 0.500000 +0.400000 0.0750000 0.625000 +0.400000 0.150000 0.750000 +0.300000 0.0750000 0.250000 +0.300000 0.150000 0.375000 +0.300000 0.225000 0.500000 +0.400000 0.0750000 0.375000 +0.400000 0.150000 0.500000 +0.400000 0.225000 0.625000 +0.500000 0.0750000 0.500000 +0.500000 0.150000 0.625000 +0.500000 0.225000 0.750000 +0.400000 0.150000 0.250000 +0.400000 0.225000 0.375000 +0.400000 0.300000 0.500000 +0.500000 0.150000 0.375000 +0.500000 0.225000 0.500000 +0.500000 0.300000 0.625000 +0.600000 0.150000 0.500000 +0.600000 0.225000 0.625000 +0.600000 0.300000 0.750000 +0.200000 0.150000 0.500000 +0.200000 0.225000 0.625000 +0.200000 0.300000 0.750000 +0.300000 0.150000 0.625000 +0.300000 0.225000 0.750000 +0.300000 0.300000 0.875000 +0.400000 0.150000 0.750000 +0.400000 0.225000 0.875000 +0.400000 0.300000 1.00000 +0.300000 0.225000 0.500000 +0.300000 0.300000 0.625000 +0.300000 0.375000 0.750000 +0.400000 0.225000 0.625000 +0.400000 0.300000 0.750000 +0.400000 0.375000 0.875000 +0.500000 0.225000 0.750000 +0.500000 0.300000 0.875000 +0.500000 0.375000 1.00000 +0.400000 0.300000 0.500000 +0.400000 0.375000 0.625000 +0.400000 0.450000 0.750000 +0.500000 0.300000 0.625000 +0.500000 0.375000 0.750000 +0.500000 0.450000 0.875000 +0.600000 0.300000 0.750000 +0.600000 0.375000 0.875000 +0.600000 0.450000 1.00000 +0.200000 0.150000 0.00000 +0.200000 0.225000 0.125000 +0.200000 0.300000 0.250000 +0.300000 0.150000 0.125000 +0.300000 0.225000 0.250000 +0.300000 0.300000 0.375000 +0.400000 0.150000 0.250000 +0.400000 0.225000 0.375000 +0.400000 0.300000 0.500000 +0.300000 0.225000 0.00000 +0.300000 0.300000 0.125000 +0.300000 0.375000 0.250000 +0.400000 0.225000 0.125000 +0.400000 0.300000 0.250000 +0.400000 0.375000 0.375000 +0.500000 0.225000 0.250000 +0.500000 0.300000 0.375000 +0.500000 0.375000 0.500000 +0.400000 0.300000 0.00000 +0.400000 0.375000 0.125000 +0.400000 0.450000 0.250000 +0.500000 0.300000 0.125000 +0.500000 0.375000 0.250000 +0.500000 0.450000 0.375000 +0.600000 0.300000 0.250000 +0.600000 0.375000 0.375000 +0.600000 0.450000 0.500000 +0.200000 0.300000 0.250000 +0.200000 0.375000 0.375000 +0.200000 0.450000 0.500000 +0.300000 0.300000 0.375000 +0.300000 0.375000 0.500000 +0.300000 0.450000 0.625000 +0.400000 0.300000 0.500000 +0.400000 0.375000 0.625000 +0.400000 0.450000 0.750000 +0.300000 0.375000 0.250000 +0.300000 0.450000 0.375000 +0.300000 0.525000 0.500000 +0.400000 0.375000 0.375000 +0.400000 0.450000 0.500000 +0.400000 0.525000 0.625000 +0.500000 0.375000 0.500000 +0.500000 0.450000 0.625000 +0.500000 0.525000 0.750000 +0.400000 0.450000 0.250000 +0.400000 0.525000 0.375000 +0.400000 0.600000 0.500000 +0.500000 0.450000 0.375000 +0.500000 0.525000 0.500000 +0.500000 0.600000 0.625000 +0.600000 0.450000 0.500000 +0.600000 0.525000 0.625000 +0.600000 0.600000 0.750000 +0.400000 0.150000 0.250000 +0.400000 0.225000 0.375000 +0.400000 0.300000 0.500000 +0.500000 0.150000 0.375000 +0.500000 0.225000 0.500000 +0.500000 0.300000 0.625000 +0.600000 0.150000 0.500000 +0.600000 0.225000 0.625000 +0.600000 0.300000 0.750000 +0.500000 0.225000 0.250000 +0.500000 0.300000 0.375000 +0.500000 0.375000 0.500000 +0.600000 0.225000 0.375000 +0.600000 0.300000 0.500000 +0.600000 0.375000 0.625000 +0.700000 0.225000 0.500000 +0.700000 0.300000 0.625000 +0.700000 0.375000 0.750000 +0.600000 0.300000 0.250000 +0.600000 0.375000 0.375000 +0.600000 0.450000 0.500000 +0.700000 0.300000 0.375000 +0.700000 0.375000 0.500000 +0.700000 0.450000 0.625000 +0.800000 0.300000 0.500000 +0.800000 0.375000 0.625000 +0.800000 0.450000 0.750000 +0.400000 0.300000 0.500000 +0.400000 0.375000 0.625000 +0.400000 0.450000 0.750000 +0.500000 0.300000 0.625000 +0.500000 0.375000 0.750000 +0.500000 0.450000 0.875000 +0.600000 0.300000 0.750000 +0.600000 0.375000 0.875000 +0.600000 0.450000 1.00000 +0.500000 0.375000 0.500000 +0.500000 0.450000 0.625000 +0.500000 0.525000 0.750000 +0.600000 0.375000 0.625000 +0.600000 0.450000 0.750000 +0.600000 0.525000 0.875000 +0.700000 0.375000 0.750000 +0.700000 0.450000 0.875000 +0.700000 0.525000 1.00000 +0.600000 0.450000 0.500000 +0.600000 0.525000 0.625000 +0.600000 0.600000 0.750000 +0.700000 0.450000 0.625000 +0.700000 0.525000 0.750000 +0.700000 0.600000 0.875000 +0.800000 0.450000 0.750000 +0.800000 0.525000 0.875000 +0.800000 0.600000 1.00000 + +CELLS 64 576 +8 0 1 4 3 9 10 13 12 +8 1 2 5 4 10 11 14 13 +8 3 4 7 6 12 13 16 15 +8 4 5 8 7 13 14 17 16 +8 9 10 13 12 18 19 22 21 +8 10 11 14 13 19 20 23 22 +8 12 13 16 15 21 22 25 24 +8 13 14 17 16 22 23 26 25 +8 27 28 31 30 36 37 40 39 +8 28 29 32 31 37 38 41 40 +8 30 31 34 33 39 40 43 42 +8 31 32 35 34 40 41 44 43 +8 36 37 40 39 45 46 49 48 +8 37 38 41 40 46 47 50 49 +8 39 40 43 42 48 49 52 51 +8 40 41 44 43 49 50 53 52 +8 54 55 58 57 63 64 67 66 +8 55 56 59 58 64 65 68 67 +8 57 58 61 60 66 67 70 69 +8 58 59 62 61 67 68 71 70 +8 63 64 67 66 72 73 76 75 +8 64 65 68 67 73 74 77 76 +8 66 67 70 69 75 76 79 78 +8 67 68 71 70 76 77 80 79 +8 81 82 85 84 90 91 94 93 +8 82 83 86 85 91 92 95 94 +8 84 85 88 87 93 94 97 96 +8 85 86 89 88 94 95 98 97 +8 90 91 94 93 99 100 103 102 +8 91 92 95 94 100 101 104 103 +8 93 94 97 96 102 103 106 105 +8 94 95 98 97 103 104 107 106 +8 108 109 112 111 117 118 121 120 +8 109 110 113 112 118 119 122 121 +8 111 112 115 114 120 121 124 123 +8 112 113 116 115 121 122 125 124 +8 117 118 121 120 126 127 130 129 +8 118 119 122 121 127 128 131 130 +8 120 121 124 123 129 130 133 132 +8 121 122 125 124 130 131 134 133 +8 135 136 139 138 144 145 148 147 +8 136 137 140 139 145 146 149 148 +8 138 139 142 141 147 148 151 150 +8 139 140 143 142 148 149 152 151 +8 144 145 148 147 153 154 157 156 +8 145 146 149 148 154 155 158 157 +8 147 148 151 150 156 157 160 159 +8 148 149 152 151 157 158 161 160 +8 162 163 166 165 171 172 175 174 +8 163 164 167 166 172 173 176 175 +8 165 166 169 168 174 175 178 177 +8 166 167 170 169 175 176 179 178 +8 171 172 175 174 180 181 184 183 +8 172 173 176 175 181 182 185 184 +8 174 175 178 177 183 184 187 186 +8 175 176 179 178 184 185 188 187 +8 189 190 193 192 198 199 202 201 +8 190 191 194 193 199 200 203 202 +8 192 193 196 195 201 202 205 204 +8 193 194 197 196 202 203 206 205 +8 198 199 202 201 207 208 211 210 +8 199 200 203 202 208 209 212 211 +8 201 202 205 204 210 211 214 213 +8 202 203 206 205 211 212 215 214 + +CELL_TYPES 64 + 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 +POINT_DATA 216 +SCALARS v double 1 +LOOKUP_TABLE default +0.00000 2.00000 2.00000 0.00000 4.00000 4.00000 0.00000 4.00000 4.00000 0.00000 2.00000 2.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 2.00000 2.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 2.00000 2.00000 1.00000 4.00000 4.00000 1.00000 4.00000 4.00000 1.00000 2.00000 2.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 2.00000 2.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 0.00000 4.00000 4.00000 0.00000 4.00000 4.00000 0.00000 3.00000 3.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 3.00000 3.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 3.00000 3.00000 4.00000 4.00000 1.00000 4.00000 4.00000 1.00000 3.00000 3.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 3.00000 3.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 3.00000 3.00000 1.00000 0.00000 2.00000 2.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 2.00000 2.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 2.00000 2.00000 0.00000 5.00000 5.00000 0.00000 5.00000 5.00000 2.00000 2.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 2.00000 2.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 2.00000 2.00000 1.00000 5.00000 5.00000 1.00000 5.00000 5.00000 1.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 3.00000 3.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 3.00000 3.00000 0.00000 5.00000 5.00000 0.00000 5.00000 5.00000 0.00000 3.00000 3.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 3.00000 3.00000 1.00000 0.00000 0.00000 1.00000 0.00000 0.00000 1.00000 3.00000 3.00000 1.00000 5.00000 5.00000 1.00000 5.00000 5.00000 1.00000 3.00000 3.00000 1.00000 +DEAL::* checking dim=3 subd=3 3 3 +DEAL:: center: 0.133333 0.0500000 0.0833333 id: 0 +DEAL:: center: 0.0666667 0.100000 0.0833333 id: 2 +DEAL:: center: 0.0666667 0.0500000 0.166667 id: 4 +DEAL:: center: 0.0666667 0.200000 0.250000 id: 2 +DEAL:: center: 0.0666667 0.150000 0.333333 id: 4 +DEAL:: center: 0.133333 0.350000 0.583333 id: 1 +DEAL:: center: 0.0666667 0.300000 0.416667 id: 2 +DEAL:: center: 0.0666667 0.250000 0.500000 id: 4 +DEAL:: center: 0.266667 0.0500000 0.250000 id: 0 +DEAL:: center: 0.200000 0.0500000 0.333333 id: 4 +DEAL:: center: 0.200000 0.150000 0.500000 id: 4 +DEAL:: center: 0.266667 0.350000 0.750000 id: 1 +DEAL:: center: 0.200000 0.250000 0.666667 id: 4 +DEAL:: center: 0.400000 0.0500000 0.416667 id: 0 +DEAL:: center: 0.466667 0.100000 0.583333 id: 3 +DEAL:: center: 0.333333 0.0500000 0.500000 id: 4 +DEAL:: center: 0.466667 0.200000 0.750000 id: 3 +DEAL:: center: 0.333333 0.150000 0.666667 id: 4 +DEAL:: center: 0.400000 0.350000 0.916667 id: 1 +DEAL:: center: 0.466667 0.300000 0.916667 id: 3 +DEAL:: center: 0.333333 0.250000 0.833333 id: 4 +DEAL:: center: 0.266667 0.150000 0.0833333 id: 0 +DEAL:: center: 0.200000 0.200000 0.0833333 id: 2 +DEAL:: center: 0.200000 0.300000 0.250000 id: 2 +DEAL:: center: 0.266667 0.450000 0.583333 id: 1 +DEAL:: center: 0.200000 0.400000 0.416667 id: 2 +DEAL:: center: 0.400000 0.150000 0.250000 id: 0 +DEAL:: center: 0.400000 0.450000 0.750000 id: 1 +DEAL:: center: 0.533333 0.150000 0.416667 id: 0 +DEAL:: center: 0.600000 0.200000 0.583333 id: 3 +DEAL:: center: 0.600000 0.300000 0.750000 id: 3 +DEAL:: center: 0.533333 0.450000 0.916667 id: 1 +DEAL:: center: 0.600000 0.400000 0.916667 id: 3 +DEAL:: center: 0.400000 0.250000 0.0833333 id: 0 +DEAL:: center: 0.333333 0.300000 0.0833333 id: 2 +DEAL:: center: 0.466667 0.350000 0.166667 id: 5 +DEAL:: center: 0.333333 0.400000 0.250000 id: 2 +DEAL:: center: 0.466667 0.450000 0.333333 id: 5 +DEAL:: center: 0.400000 0.550000 0.583333 id: 1 +DEAL:: center: 0.333333 0.500000 0.416667 id: 2 +DEAL:: center: 0.466667 0.550000 0.500000 id: 5 +DEAL:: center: 0.533333 0.250000 0.250000 id: 0 +DEAL:: center: 0.600000 0.350000 0.333333 id: 5 +DEAL:: center: 0.600000 0.450000 0.500000 id: 5 +DEAL:: center: 0.533333 0.550000 0.750000 id: 1 +DEAL:: center: 0.600000 0.550000 0.666667 id: 5 +DEAL:: center: 0.666667 0.250000 0.416667 id: 0 +DEAL:: center: 0.733333 0.300000 0.583333 id: 3 +DEAL:: center: 0.733333 0.350000 0.500000 id: 5 +DEAL:: center: 0.733333 0.400000 0.750000 id: 3 +DEAL:: center: 0.733333 0.450000 0.666667 id: 5 +DEAL:: center: 0.666667 0.550000 0.916667 id: 1 +DEAL:: center: 0.733333 0.500000 0.916667 id: 3 +DEAL:: center: 0.733333 0.550000 0.833333 id: 5 +DEAL:: boundary indicators: 0(9 times) 1(9 times) 2(9 times) 3(9 times) 4(9 times) 5(9 times) +DEAL::* checking dim=3 subd=1 2 1 +DEAL:: center: 0.300000 0.150000 0.125000 id: 0 +DEAL:: center: 0.300000 0.450000 0.625000 id: 1 +DEAL:: center: 0.200000 0.300000 0.250000 id: 2 +DEAL:: center: 0.100000 0.150000 0.375000 id: 4 +DEAL:: center: 0.500000 0.450000 0.375000 id: 5 +DEAL:: center: 0.500000 0.150000 0.375000 id: 0 +DEAL:: center: 0.500000 0.450000 0.875000 id: 1 +DEAL:: center: 0.600000 0.300000 0.750000 id: 3 +DEAL:: center: 0.300000 0.150000 0.625000 id: 4 +DEAL:: center: 0.700000 0.450000 0.625000 id: 5 +DEAL:: boundary indicators: 0(2 times) 1(2 times) 2(1 times) 3(1 times) 4(2 times) 5(2 times) +DEAL::* checking dim=3 subd=1 2 2 +DEAL:: center: 0.200000 0.0750000 0.125000 id: 0 +DEAL:: center: 0.200000 0.375000 0.625000 id: 1 +DEAL:: center: 0.100000 0.225000 0.250000 id: 2 +DEAL:: center: 0.100000 0.150000 0.375000 id: 4 +DEAL:: center: 0.400000 0.0750000 0.375000 id: 0 +DEAL:: center: 0.400000 0.375000 0.875000 id: 1 +DEAL:: center: 0.500000 0.225000 0.750000 id: 3 +DEAL:: center: 0.300000 0.150000 0.625000 id: 4 +DEAL:: center: 0.400000 0.225000 0.125000 id: 0 +DEAL:: center: 0.400000 0.525000 0.625000 id: 1 +DEAL:: center: 0.300000 0.375000 0.250000 id: 2 +DEAL:: center: 0.500000 0.450000 0.375000 id: 5 +DEAL:: center: 0.600000 0.225000 0.375000 id: 0 +DEAL:: center: 0.600000 0.525000 0.875000 id: 1 +DEAL:: center: 0.700000 0.375000 0.750000 id: 3 +DEAL:: center: 0.700000 0.450000 0.625000 id: 5 +DEAL:: boundary indicators: 0(4 times) 1(4 times) 2(2 times) 3(2 times) 4(2 times) 5(2 times) +DEAL::* checking dim=3 subd=1 2 3 +DEAL:: center: 0.166667 0.0500000 0.125000 id: 0 +DEAL:: center: 0.166667 0.350000 0.625000 id: 1 +DEAL:: center: 0.0666667 0.200000 0.250000 id: 2 +DEAL:: center: 0.100000 0.150000 0.375000 id: 4 +DEAL:: center: 0.366667 0.0500000 0.375000 id: 0 +DEAL:: center: 0.366667 0.350000 0.875000 id: 1 +DEAL:: center: 0.466667 0.200000 0.750000 id: 3 +DEAL:: center: 0.300000 0.150000 0.625000 id: 4 +DEAL:: center: 0.300000 0.150000 0.125000 id: 0 +DEAL:: center: 0.300000 0.450000 0.625000 id: 1 +DEAL:: center: 0.200000 0.300000 0.250000 id: 2 +DEAL:: center: 0.500000 0.150000 0.375000 id: 0 +DEAL:: center: 0.500000 0.450000 0.875000 id: 1 +DEAL:: center: 0.600000 0.300000 0.750000 id: 3 +DEAL:: center: 0.433333 0.250000 0.125000 id: 0 +DEAL:: center: 0.433333 0.550000 0.625000 id: 1 +DEAL:: center: 0.333333 0.400000 0.250000 id: 2 +DEAL:: center: 0.500000 0.450000 0.375000 id: 5 +DEAL:: center: 0.633333 0.250000 0.375000 id: 0 +DEAL:: center: 0.633333 0.550000 0.875000 id: 1 +DEAL:: center: 0.733333 0.400000 0.750000 id: 3 +DEAL:: center: 0.700000 0.450000 0.625000 id: 5 +DEAL:: boundary indicators: 0(6 times) 1(6 times) 2(3 times) 3(3 times) 4(2 times) 5(2 times)