From: bangerth Date: Thu, 9 Oct 2008 21:39:34 +0000 (+0000) Subject: Fix the remaining bug in dealing with level_1_difference_at_vertices in 1d. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=60d9561897f7121458eb5df465edb867035684f1;p=dealii-svn.git Fix the remaining bug in dealing with level_1_difference_at_vertices in 1d. git-svn-id: https://svn.dealii.org/trunk@17158 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/grid/tria.cc b/deal.II/deal.II/source/grid/tria.cc index ff72bd46bc..6991acaa0e 100644 --- a/deal.II/deal.II/source/grid/tria.cc +++ b/deal.II/deal.II/source/grid/tria.cc @@ -10438,211 +10438,249 @@ void Triangulation::fix_coarsen_flags () // here makes our lives a bit easier as // well as it takes care of these cases // earlier than it would otherwise happen. - if (smooth_grid & limit_level_difference_at_vertices) + // + // the main difference to the code + // in p_c_and_r is that here we + // absolutely have to make sure + // that we get things right, + // i.e. that in particular we set + // flags right if + // limit_level_difference_at_vertices + // is set. to do so we iterate + // until the flags don't change any + // more + std::vector previous_coarsen_flags (n_active_cells()); + save_coarsen_flags (previous_coarsen_flags); + + bool continue_iterating = true; + + do { - Assert(!anisotropic_refinement, - ExcMessage("In case of anisotropic refinement the " - "limit_level_difference_at_vertices flag for " - "mesh smoothing must not be set!")); - - // store highest level one - // of the cells adjacent to - // a vertex belongs to - std::vector vertex_level (vertices.size(), 0); - active_cell_iterator cell = begin_active(), - endc = end(); - for (; cell!=endc; ++cell) + if (smooth_grid & limit_level_difference_at_vertices) { - if (cell->refine_flag_set()) - for (unsigned int vertex=0; vertex::vertices_per_cell; - ++vertex) - vertex_level[cell->vertex_index(vertex)] - = std::max (vertex_level[cell->vertex_index(vertex)], - cell->level()+1); - else if (!cell->coarsen_flag_set()) - for (unsigned int vertex=0; vertex::vertices_per_cell; - ++vertex) - vertex_level[cell->vertex_index(vertex)] - = std::max (vertex_level[cell->vertex_index(vertex)], - cell->level()); - else + Assert(!anisotropic_refinement, + ExcMessage("In case of anisotropic refinement the " + "limit_level_difference_at_vertices flag for " + "mesh smoothing must not be set!")); + + // store highest level one + // of the cells adjacent to + // a vertex belongs to + std::vector vertex_level (vertices.size(), 0); + active_cell_iterator cell = begin_active(), + endc = end(); + for (; cell!=endc; ++cell) { - // if coarsen flag is set then - // tentatively assume that the - // cell will be coarsened. this - // isn't always true (the - // coarsen flag could be - // removed again) and so we may - // make an error here - Assert (cell->coarsen_flag_set(), ExcInternalError()); - for (unsigned int vertex=0; vertex::vertices_per_cell; - ++vertex) - vertex_level[cell->vertex_index(vertex)] - = std::max (vertex_level[cell->vertex_index(vertex)], - cell->level()-1); + if (cell->refine_flag_set()) + for (unsigned int vertex=0; vertex::vertices_per_cell; + ++vertex) + vertex_level[cell->vertex_index(vertex)] + = std::max (vertex_level[cell->vertex_index(vertex)], + cell->level()+1); + else if (!cell->coarsen_flag_set()) + for (unsigned int vertex=0; vertex::vertices_per_cell; + ++vertex) + vertex_level[cell->vertex_index(vertex)] + = std::max (vertex_level[cell->vertex_index(vertex)], + cell->level()); + else + { + // if coarsen flag is + // set then tentatively + // assume that the cell + // will be + // coarsened. this + // isn't always true + // (the coarsen flag + // could be removed + // again) and so we may + // make an error + // here. we try to + // correct this by + // iterating over the + // entire process until + // we are converged + Assert (cell->coarsen_flag_set(), ExcInternalError()); + for (unsigned int vertex=0; vertex::vertices_per_cell; + ++vertex) + vertex_level[cell->vertex_index(vertex)] + = std::max (vertex_level[cell->vertex_index(vertex)], + cell->level()-1); + } } - } - // loop over all cells in reverse - // order. do so because we can then - // update the vertex levels on the - // adjacent vertices and maybe - // already flag additional cells in - // this loop - // - // note that not only may we have - // to add additional refinement - // flags, but we will also have to - // remove coarsening flags on cells - // adjacent to vertices that will - // see refinement - for (cell=last_active(); cell != endc; --cell) - if (cell->refine_flag_set() == false) - { - for (unsigned int vertex=0; - vertex::vertices_per_cell; ++vertex) - if (vertex_level[cell->vertex_index(vertex)] >= - cell->level()+1) - { - // remove coarsen flag... - cell->clear_coarsen_flag(); - - // ...and if necessary also - // refine the current cell, - // at the same time - // updating the level - // information about - // vertices - if (vertex_level[cell->vertex_index(vertex)] > + // loop over all cells in reverse + // order. do so because we can then + // update the vertex levels on the + // adjacent vertices and maybe + // already flag additional cells in + // this loop + // + // note that not only may we have + // to add additional refinement + // flags, but we will also have to + // remove coarsening flags on cells + // adjacent to vertices that will + // see refinement + for (cell=last_active(); cell != endc; --cell) + if (cell->refine_flag_set() == false) + { + for (unsigned int vertex=0; + vertex::vertices_per_cell; ++vertex) + if (vertex_level[cell->vertex_index(vertex)] >= cell->level()+1) { - cell->set_refine_flag(); + // remove coarsen flag... + cell->clear_coarsen_flag(); + + // ...and if necessary also + // refine the current cell, + // at the same time + // updating the level + // information about + // vertices + if (vertex_level[cell->vertex_index(vertex)] > + cell->level()+1) + { + cell->set_refine_flag(); - for (unsigned int v=0; v::vertices_per_cell; - ++v) - vertex_level[cell->vertex_index(v)] - = std::max (vertex_level[cell->vertex_index(v)], - cell->level()+1); - } + for (unsigned int v=0; v::vertices_per_cell; + ++v) + vertex_level[cell->vertex_index(v)] + = std::max (vertex_level[cell->vertex_index(v)], + cell->level()+1); + } - // continue and see whether - // we may, for example, go - // into the inner 'if' - // above based on a - // different vertex - } - } - } + // continue and see whether + // we may, for example, go + // into the inner 'if' + // above based on a + // different vertex + } + } + } - // loop over all cells. Flag all - // cells of which all children are - // flagged for coarsening and - // delete the childrens' - // flags. Also delete all flags of - // cells for which not all children - // of a cell are flagged. In - // effect, only those cells are - // flagged of which originally all - // children were flagged and for - // which all children are on the - // same refinement level. For - // flagging, the user flags are - // used, to avoid confusion and - // because non-active cells can't - // be flagged for coarsening - // - // In effect, all coarsen flags are - // turned into user flags of the - // mother cell if coarsening is - // possible or deleted - // otherwise. - clear_user_flags (); - // Coarsen flags of - // cells with no mother cell, - // i.e. on the coarsest level are - // deleted explicitly. - active_cell_iterator acell = begin_active(0), - end_ac = end_active(0); - for (; acell!=end_ac; ++acell) - acell->clear_coarsen_flag(); + // loop over all cells. Flag all + // cells of which all children are + // flagged for coarsening and + // delete the childrens' + // flags. Also delete all flags of + // cells for which not all children + // of a cell are flagged. In + // effect, only those cells are + // flagged of which originally all + // children were flagged and for + // which all children are on the + // same refinement level. For + // flagging, the user flags are + // used, to avoid confusion and + // because non-active cells can't + // be flagged for coarsening + // + // In effect, all coarsen flags are + // turned into user flags of the + // mother cell if coarsening is + // possible or deleted + // otherwise. + clear_user_flags (); + // Coarsen flags of + // cells with no mother cell, + // i.e. on the coarsest level are + // deleted explicitly. + active_cell_iterator acell = begin_active(0), + end_ac = end_active(0); + for (; acell!=end_ac; ++acell) + acell->clear_coarsen_flag(); - cell_iterator cell = begin(), - endc = end(); - for (; cell!=endc; ++cell) - { - // nothing to do if we are - // already on the finest level - if (cell->active()) - continue; + cell_iterator cell = begin(), + endc = end(); + for (; cell!=endc; ++cell) + { + // nothing to do if we are + // already on the finest level + if (cell->active()) + continue; - const unsigned int n_children=cell->n_children(); - unsigned int flagged_children=0; - for (unsigned int child=0; childchild(child)->active() && - cell->child(child)->coarsen_flag_set()) - { - ++flagged_children; - // clear flag since we - // don't need it anymore - cell->child(child)->clear_coarsen_flag(); - } + const unsigned int n_children=cell->n_children(); + unsigned int flagged_children=0; + for (unsigned int child=0; childchild(child)->active() && + cell->child(child)->coarsen_flag_set()) + { + ++flagged_children; + // clear flag since we + // don't need it anymore + cell->child(child)->clear_coarsen_flag(); + } - // flag this cell for - // coarsening if all children - // were flagged - if (flagged_children == n_children) - cell->set_user_flag(); - } + // flag this cell for + // coarsening if all children + // were flagged + if (flagged_children == n_children) + cell->set_user_flag(); + } - // in principle no coarsen flags - // should be set any more at this - // point + // in principle no coarsen flags + // should be set any more at this + // point #if DEBUG - for (cell=begin(); cell!=endc; ++cell) - Assert (cell->coarsen_flag_set() == false, ExcInternalError()); + for (cell=begin(); cell!=endc; ++cell) + Assert (cell->coarsen_flag_set() == false, ExcInternalError()); #endif - // now loop over all cells which have the - // user flag set. their children were - // flagged for coarsening. set the coarsen - // flag again if we are sure that none of - // the neighbors of these children are - // refined, or will be refined, since then - // we would get a two-level jump in - // refinement. on the other hand, if one of - // the children's neighbors has their user - // flag set, then we know that its children - // will go away by coarsening, and we will - // be ok. - // - // note on the other hand that we do allow - // level-2 jumps in refinement between - // neighbors in 1d, so this whole procedure - // is only necessary if we are not in 1d - // - // since we remove some coarsening/user - // flags in the process, we have to work - // from the finest level to the coarsest - // one, since we occasionally inspect user - // flags of cells on finer levels and need - // to be sure that these flags are final - for (cell=last(); cell!=endc; --cell) - if (cell->user_flag_set()) - // if allowed: flag the - // children for coarsening - if (coarsening_allowed(cell)) - for (unsigned int c=0; cn_children(); ++c) - { - Assert (cell->child(c)->refine_flag_set()==false, - ExcInternalError()); + // now loop over all cells which have the + // user flag set. their children were + // flagged for coarsening. set the coarsen + // flag again if we are sure that none of + // the neighbors of these children are + // refined, or will be refined, since then + // we would get a two-level jump in + // refinement. on the other hand, if one of + // the children's neighbors has their user + // flag set, then we know that its children + // will go away by coarsening, and we will + // be ok. + // + // note on the other hand that we do allow + // level-2 jumps in refinement between + // neighbors in 1d, so this whole procedure + // is only necessary if we are not in 1d + // + // since we remove some coarsening/user + // flags in the process, we have to work + // from the finest level to the coarsest + // one, since we occasionally inspect user + // flags of cells on finer levels and need + // to be sure that these flags are final + for (cell=last(); cell!=endc; --cell) + if (cell->user_flag_set()) + // if allowed: flag the + // children for coarsening + if (coarsening_allowed(cell)) + for (unsigned int c=0; cn_children(); ++c) + { + Assert (cell->child(c)->refine_flag_set()==false, + ExcInternalError()); - cell->child(c)->set_coarsen_flag(); - } + cell->child(c)->set_coarsen_flag(); + } - // clear all user flags again, now that we - // don't need them any more - clear_user_flags (); + // clear all user flags again, now that we + // don't need them any more + clear_user_flags (); + + + // now see if anything has + // changed in the last + // iteration of this function + std::vector current_coarsen_flags (n_active_cells()); + save_coarsen_flags (current_coarsen_flags); + + continue_iterating = (current_coarsen_flags != previous_coarsen_flags); + previous_coarsen_flags = current_coarsen_flags; + } + while (continue_iterating == true); } diff --git a/tests/deal.II/coarsening_02_1d.cc b/tests/deal.II/coarsening_02_1d.cc index 3f3467cbc8..6e581d3c0b 100644 --- a/tests/deal.II/coarsening_02_1d.cc +++ b/tests/deal.II/coarsening_02_1d.cc @@ -28,6 +28,7 @@ #include #include #include +#include #include #include @@ -75,8 +76,7 @@ void test(std::ostream& /*out*/) GridGenerator::hyper_cube(triangulation); triangulation.refine_global (2); - const unsigned int n_refinements[] = { 0, 13, 11, 6 }; - for (unsigned int i=0; i flags (triangulation.n_active_cells(), false); @@ -105,7 +105,22 @@ void test(std::ostream& /*out*/) cell != triangulation.end(); ++cell, ++index) if (!flags[index]) cell->set_coarsen_flag(); - + + GridOut go; + go.write_gnuplot (triangulation, deallog.get_file_stream()); + + for (typename Triangulation::active_cell_iterator + cell = triangulation.begin_active(); + cell != triangulation.end(); ++cell) + if (cell->refine_flag_set()) + { + deallog << cell->center() << ' ' << cell->level() << std::endl; + deallog << cell->center() << ' ' << 0 << std::endl; + deallog << std::endl; + } + + + triangulation.execute_coarsening_and_refinement (); // verify that none of the cells diff --git a/tests/deal.II/coarsening_02_1d/cmp/generic b/tests/deal.II/coarsening_02_1d/cmp/generic index 25966e9bd0..d992d0f27a 100644 --- a/tests/deal.II/coarsening_02_1d/cmp/generic +++ b/tests/deal.II/coarsening_02_1d/cmp/generic @@ -1,3 +1,2821 @@ -DEAL:2d::OK -DEAL:3d::OK +0.00000 2 0 +0.250000 2 0 + +0.250000 2 0 +0.500000 2 0 + +0.500000 2 0 +0.750000 2 0 + +0.750000 2 0 +1.00000 2 0 + +DEAL:1d::0.125000 2 +DEAL:1d::0.125000 0 +DEAL:1d:: +DEAL:1d::0.875000 2 +DEAL:1d::0.875000 0 +DEAL:1d:: +DEAL:1d::Iteration 0: cells=6 +0.250000 2 0 +0.500000 2 0 + +0.500000 2 0 +0.750000 2 0 + +0.00000 3 0 +0.125000 3 0 + +0.125000 3 0 +0.250000 3 0 + +0.750000 3 0 +0.875000 3 0 + +0.875000 3 0 +1.00000 3 0 + +DEAL:1d::0.375000 2 +DEAL:1d::0.375000 0 +DEAL:1d:: +DEAL:1d::0.187500 3 +DEAL:1d::0.187500 0 +DEAL:1d:: +DEAL:1d::0.812500 3 +DEAL:1d::0.812500 0 +DEAL:1d:: +DEAL:1d::Iteration 1: cells=10 +0.00000 3 0 +0.125000 3 0 + +0.875000 3 0 +1.00000 3 0 + +0.250000 3 0 +0.375000 3 0 + +0.375000 3 0 +0.500000 3 0 + +0.500000 3 0 +0.625000 3 0 + +0.625000 3 0 +0.750000 3 0 + +0.125000 4 0 +0.187500 4 0 + +0.187500 4 0 +0.250000 4 0 + +0.750000 4 0 +0.812500 4 0 + +0.812500 4 0 +0.875000 4 0 + +DEAL:1d::0.0625000 3 +DEAL:1d::0.0625000 0 +DEAL:1d:: +DEAL:1d::0.437500 3 +DEAL:1d::0.437500 0 +DEAL:1d:: +DEAL:1d::0.687500 3 +DEAL:1d::0.687500 0 +DEAL:1d:: +DEAL:1d::Iteration 2: cells=11 +0.125000 3 0 +0.250000 3 0 + +0.750000 3 0 +0.875000 3 0 + +0.875000 3 0 +1.00000 3 0 + +0.250000 3 0 +0.375000 3 0 + +0.500000 3 0 +0.625000 3 0 + +0.00000 4 0 +0.0625000 4 0 + +0.0625000 4 0 +0.125000 4 0 + +0.375000 4 0 +0.437500 4 0 + +0.437500 4 0 +0.500000 4 0 + +0.625000 4 0 +0.687500 4 0 + +0.687500 4 0 +0.750000 4 0 + +DEAL:1d::0.187500 3 +DEAL:1d::0.187500 0 +DEAL:1d:: +DEAL:1d::0.312500 3 +DEAL:1d::0.312500 0 +DEAL:1d:: +DEAL:1d::0.0937500 4 +DEAL:1d::0.0937500 0 +DEAL:1d:: +DEAL:1d::Iteration 3: cells=11 +0.750000 2 0 +1.00000 2 0 + +0.375000 3 0 +0.500000 3 0 + +0.500000 3 0 +0.625000 3 0 + +0.625000 3 0 +0.750000 3 0 + +0.00000 4 0 +0.0625000 4 0 + +0.125000 4 0 +0.187500 4 0 + +0.187500 4 0 +0.250000 4 0 + +0.250000 4 0 +0.312500 4 0 + +0.312500 4 0 +0.375000 4 0 + +0.0625000 5 0 +0.0937500 5 0 + +0.0937500 5 0 +0.125000 5 0 + +DEAL:1d::0.875000 2 +DEAL:1d::0.875000 0 +DEAL:1d:: +DEAL:1d::0.437500 3 +DEAL:1d::0.437500 0 +DEAL:1d:: +DEAL:1d::0.281250 4 +DEAL:1d::0.281250 0 +DEAL:1d:: +DEAL:1d::0.343750 4 +DEAL:1d::0.343750 0 +DEAL:1d:: +DEAL:1d::Iteration 4: cells=14 +0.750000 3 0 +0.875000 3 0 + +0.875000 3 0 +1.00000 3 0 + +0.500000 3 0 +0.625000 3 0 + +0.625000 3 0 +0.750000 3 0 + +0.00000 4 0 +0.0625000 4 0 + +0.0625000 4 0 +0.125000 4 0 + +0.125000 4 0 +0.187500 4 0 + +0.187500 4 0 +0.250000 4 0 + +0.375000 4 0 +0.437500 4 0 + +0.437500 4 0 +0.500000 4 0 + +0.250000 5 0 +0.281250 5 0 + +0.281250 5 0 +0.312500 5 0 + +0.312500 5 0 +0.343750 5 0 + +0.343750 5 0 +0.375000 5 0 + +DEAL:1d::0.812500 3 +DEAL:1d::0.812500 0 +DEAL:1d:: +DEAL:1d::0.156250 4 +DEAL:1d::0.156250 0 +DEAL:1d:: +DEAL:1d::0.296875 5 +DEAL:1d::0.296875 0 +DEAL:1d:: +DEAL:1d::0.359375 5 +DEAL:1d::0.359375 0 +DEAL:1d:: +DEAL:1d::Iteration 5: cells=19 +0.875000 3 0 +1.00000 3 0 + +0.500000 3 0 +0.625000 3 0 + +0.625000 3 0 +0.750000 3 0 + +0.00000 4 0 +0.0625000 4 0 + +0.0625000 4 0 +0.125000 4 0 + +0.187500 4 0 +0.250000 4 0 + +0.437500 4 0 +0.500000 4 0 + +0.750000 4 0 +0.812500 4 0 + +0.812500 4 0 +0.875000 4 0 + +0.250000 5 0 +0.281250 5 0 + +0.312500 5 0 +0.343750 5 0 + +0.125000 5 0 +0.156250 5 0 + +0.156250 5 0 +0.187500 5 0 + +0.375000 5 0 +0.406250 5 0 + +0.406250 5 0 +0.437500 5 0 + +0.281250 6 0 +0.296875 6 0 + +0.296875 6 0 +0.312500 6 0 + +0.343750 6 0 +0.359375 6 0 + +0.359375 6 0 +0.375000 6 0 + +DEAL:1d::0.937500 3 +DEAL:1d::0.937500 0 +DEAL:1d:: +DEAL:1d::0.218750 4 +DEAL:1d::0.218750 0 +DEAL:1d:: +DEAL:1d::0.781250 4 +DEAL:1d::0.781250 0 +DEAL:1d:: +DEAL:1d::0.328125 5 +DEAL:1d::0.328125 0 +DEAL:1d:: +DEAL:1d::0.351562 6 +DEAL:1d::0.351562 0 +DEAL:1d:: +DEAL:1d::Iteration 6: cells=22 +0.00000 3 0 +0.125000 3 0 + +0.500000 3 0 +0.625000 3 0 + +0.875000 4 0 +0.937500 4 0 + +0.937500 4 0 +1.00000 4 0 + +0.125000 4 0 +0.187500 4 0 + +0.437500 4 0 +0.500000 4 0 + +0.812500 4 0 +0.875000 4 0 + +0.625000 4 0 +0.687500 4 0 + +0.687500 4 0 +0.750000 4 0 + +0.250000 5 0 +0.281250 5 0 + +0.281250 5 0 +0.312500 5 0 + +0.187500 5 0 +0.218750 5 0 + +0.218750 5 0 +0.250000 5 0 + +0.375000 5 0 +0.406250 5 0 + +0.406250 5 0 +0.437500 5 0 + +0.750000 5 0 +0.781250 5 0 + +0.781250 5 0 +0.812500 5 0 + +0.312500 6 0 +0.328125 6 0 + +0.328125 6 0 +0.343750 6 0 + +0.359375 6 0 +0.375000 6 0 + +0.343750 7 0 +0.351562 7 0 + +0.351562 7 0 +0.359375 7 0 + +DEAL:1d::0.0625000 3 +DEAL:1d::0.0625000 0 +DEAL:1d:: +DEAL:1d::0.468750 4 +DEAL:1d::0.468750 0 +DEAL:1d:: +DEAL:1d::0.718750 4 +DEAL:1d::0.718750 0 +DEAL:1d:: +DEAL:1d::0.203125 5 +DEAL:1d::0.203125 0 +DEAL:1d:: +DEAL:1d::0.355469 7 +DEAL:1d::0.355469 0 +DEAL:1d:: +DEAL:1d::Iteration 7: cells=29 +0.875000 3 0 +1.00000 3 0 + +0.00000 4 0 +0.0625000 4 0 + +0.0625000 4 0 +0.125000 4 0 + +0.750000 4 0 +0.812500 4 0 + +0.812500 4 0 +0.875000 4 0 + +0.625000 4 0 +0.687500 4 0 + +0.500000 4 0 +0.562500 4 0 + +0.562500 4 0 +0.625000 4 0 + +0.250000 5 0 +0.281250 5 0 + +0.281250 5 0 +0.312500 5 0 + +0.218750 5 0 +0.250000 5 0 + +0.406250 5 0 +0.437500 5 0 + +0.125000 5 0 +0.156250 5 0 + +0.156250 5 0 +0.187500 5 0 + +0.437500 5 0 +0.468750 5 0 + +0.468750 5 0 +0.500000 5 0 + +0.687500 5 0 +0.718750 5 0 + +0.718750 5 0 +0.750000 5 0 + +0.312500 6 0 +0.328125 6 0 + +0.328125 6 0 +0.343750 6 0 + +0.187500 6 0 +0.203125 6 0 + +0.203125 6 0 +0.218750 6 0 + +0.375000 6 0 +0.390625 6 0 + +0.390625 6 0 +0.406250 6 0 + +0.343750 7 0 +0.351562 7 0 + +0.359375 7 0 +0.367188 7 0 + +0.367188 7 0 +0.375000 7 0 + +0.351562 8 0 +0.355469 8 0 + +0.355469 8 0 +0.359375 8 0 + +DEAL:1d::0.937500 3 +DEAL:1d::0.937500 0 +DEAL:1d:: +DEAL:1d::0.0312500 4 +DEAL:1d::0.0312500 0 +DEAL:1d:: +DEAL:1d::0.0937500 4 +DEAL:1d::0.0937500 0 +DEAL:1d:: +DEAL:1d::0.656250 4 +DEAL:1d::0.656250 0 +DEAL:1d:: +DEAL:1d::0.234375 5 +DEAL:1d::0.234375 0 +DEAL:1d:: +DEAL:1d::0.703125 5 +DEAL:1d::0.703125 0 +DEAL:1d:: +DEAL:1d::0.734375 5 +DEAL:1d::0.734375 0 +DEAL:1d:: +DEAL:1d::Iteration 8: cells=31 +0.125000 4 0 +0.187500 4 0 + +0.437500 4 0 +0.500000 4 0 + +0.812500 4 0 +0.875000 4 0 + +0.500000 4 0 +0.562500 4 0 + +0.562500 4 0 +0.625000 4 0 + +0.875000 4 0 +0.937500 4 0 + +0.937500 4 0 +1.00000 4 0 + +0.250000 5 0 +0.281250 5 0 + +0.281250 5 0 +0.312500 5 0 + +0.187500 5 0 +0.218750 5 0 + +0.375000 5 0 +0.406250 5 0 + +0.406250 5 0 +0.437500 5 0 + +0.00000 5 0 +0.0312500 5 0 + +0.0312500 5 0 +0.0625000 5 0 + +0.0625000 5 0 +0.0937500 5 0 + +0.0937500 5 0 +0.125000 5 0 + +0.750000 5 0 +0.781250 5 0 + +0.781250 5 0 +0.812500 5 0 + +0.625000 5 0 +0.656250 5 0 + +0.656250 5 0 +0.687500 5 0 + +0.312500 6 0 +0.328125 6 0 + +0.328125 6 0 +0.343750 6 0 + +0.359375 6 0 +0.375000 6 0 + +0.218750 6 0 +0.234375 6 0 + +0.234375 6 0 +0.250000 6 0 + +0.687500 6 0 +0.703125 6 0 + +0.703125 6 0 +0.718750 6 0 + +0.718750 6 0 +0.734375 6 0 + +0.734375 6 0 +0.750000 6 0 + +0.343750 7 0 +0.351562 7 0 + +0.351562 7 0 +0.359375 7 0 + +DEAL:1d::0.156250 4 +DEAL:1d::0.156250 0 +DEAL:1d:: +DEAL:1d::0.906250 4 +DEAL:1d::0.906250 0 +DEAL:1d:: +DEAL:1d::0.0156250 5 +DEAL:1d::0.0156250 0 +DEAL:1d:: +DEAL:1d::0.0468750 5 +DEAL:1d::0.0468750 0 +DEAL:1d:: +DEAL:1d::0.367188 6 +DEAL:1d::0.367188 0 +DEAL:1d:: +DEAL:1d::0.726562 6 +DEAL:1d::0.726562 0 +DEAL:1d:: +DEAL:1d::Iteration 9: cells=34 +0.250000 4 0 +0.312500 4 0 + +0.437500 4 0 +0.500000 4 0 + +0.812500 4 0 +0.875000 4 0 + +0.500000 4 0 +0.562500 4 0 + +0.562500 4 0 +0.625000 4 0 + +0.937500 4 0 +1.00000 4 0 + +0.125000 5 0 +0.156250 5 0 + +0.156250 5 0 +0.187500 5 0 + +0.312500 5 0 +0.343750 5 0 + +0.187500 5 0 +0.218750 5 0 + +0.218750 5 0 +0.250000 5 0 + +0.406250 5 0 +0.437500 5 0 + +0.0625000 5 0 +0.0937500 5 0 + +0.0937500 5 0 +0.125000 5 0 + +0.750000 5 0 +0.781250 5 0 + +0.781250 5 0 +0.812500 5 0 + +0.625000 5 0 +0.656250 5 0 + +0.656250 5 0 +0.687500 5 0 + +0.875000 5 0 +0.906250 5 0 + +0.906250 5 0 +0.937500 5 0 + +0.375000 6 0 +0.390625 6 0 + +0.390625 6 0 +0.406250 6 0 + +0.343750 6 0 +0.359375 6 0 + +0.00000 6 0 +0.0156250 6 0 + +0.0156250 6 0 +0.0312500 6 0 + +0.687500 6 0 +0.703125 6 0 + +0.703125 6 0 +0.718750 6 0 + +0.734375 6 0 +0.750000 6 0 + +0.0312500 6 0 +0.0468750 6 0 + +0.0468750 6 0 +0.0625000 6 0 + +0.359375 7 0 +0.367188 7 0 + +0.367188 7 0 +0.375000 7 0 + +0.718750 7 0 +0.726562 7 0 + +0.726562 7 0 +0.734375 7 0 + +DEAL:1d::0.281250 4 +DEAL:1d::0.281250 0 +DEAL:1d:: +DEAL:1d::0.531250 4 +DEAL:1d::0.531250 0 +DEAL:1d:: +DEAL:1d::0.203125 5 +DEAL:1d::0.203125 0 +DEAL:1d:: +DEAL:1d::0.234375 5 +DEAL:1d::0.234375 0 +DEAL:1d:: +DEAL:1d::0.765625 5 +DEAL:1d::0.765625 0 +DEAL:1d:: +DEAL:1d::0.921875 5 +DEAL:1d::0.921875 0 +DEAL:1d:: +DEAL:1d::0.398438 6 +DEAL:1d::0.398438 0 +DEAL:1d:: +DEAL:1d::0.00781250 6 +DEAL:1d::0.00781250 0 +DEAL:1d:: +DEAL:1d::Iteration 10: cells=39 +0.0625000 4 0 +0.125000 4 0 + +0.812500 4 0 +0.875000 4 0 + +0.625000 4 0 +0.687500 4 0 + +0.562500 4 0 +0.625000 4 0 + +0.125000 5 0 +0.156250 5 0 + +0.156250 5 0 +0.187500 5 0 + +0.312500 5 0 +0.343750 5 0 + +0.0312500 5 0 +0.0625000 5 0 + +0.250000 5 0 +0.281250 5 0 + +0.281250 5 0 +0.312500 5 0 + +0.687500 5 0 +0.718750 5 0 + +0.781250 5 0 +0.812500 5 0 + +0.437500 5 0 +0.468750 5 0 + +0.468750 5 0 +0.500000 5 0 + +0.875000 5 0 +0.906250 5 0 + +0.500000 5 0 +0.531250 5 0 + +0.531250 5 0 +0.562500 5 0 + +0.937500 5 0 +0.968750 5 0 + +0.968750 5 0 +1.00000 5 0 + +0.375000 6 0 +0.390625 6 0 + +0.343750 6 0 +0.359375 6 0 + +0.359375 6 0 +0.375000 6 0 + +0.0156250 6 0 +0.0312500 6 0 + +0.187500 6 0 +0.203125 6 0 + +0.203125 6 0 +0.218750 6 0 + +0.718750 6 0 +0.734375 6 0 + +0.734375 6 0 +0.750000 6 0 + +0.218750 6 0 +0.234375 6 0 + +0.234375 6 0 +0.250000 6 0 + +0.406250 6 0 +0.421875 6 0 + +0.421875 6 0 +0.437500 6 0 + +0.750000 6 0 +0.765625 6 0 + +0.765625 6 0 +0.781250 6 0 + +0.906250 6 0 +0.921875 6 0 + +0.921875 6 0 +0.937500 6 0 + +0.390625 7 0 +0.398438 7 0 + +0.398438 7 0 +0.406250 7 0 + +0.00000 7 0 +0.00781250 7 0 + +0.00781250 7 0 +0.0156250 7 0 + +DEAL:1d::0.0937500 4 +DEAL:1d::0.0937500 0 +DEAL:1d:: +DEAL:1d::0.171875 5 +DEAL:1d::0.171875 0 +DEAL:1d:: +DEAL:1d::0.890625 5 +DEAL:1d::0.890625 0 +DEAL:1d:: +DEAL:1d::0.515625 5 +DEAL:1d::0.515625 0 +DEAL:1d:: +DEAL:1d::0.953125 5 +DEAL:1d::0.953125 0 +DEAL:1d:: +DEAL:1d::0.742188 6 +DEAL:1d::0.742188 0 +DEAL:1d:: +DEAL:1d::0.757812 6 +DEAL:1d::0.757812 0 +DEAL:1d:: +DEAL:1d::0.00390625 7 +DEAL:1d::0.00390625 0 +DEAL:1d:: +DEAL:1d::Iteration 11: cells=41 +0.250000 4 0 +0.312500 4 0 + +0.625000 4 0 +0.687500 4 0 + +0.562500 4 0 +0.625000 4 0 + +0.125000 5 0 +0.156250 5 0 + +0.312500 5 0 +0.343750 5 0 + +0.343750 5 0 +0.375000 5 0 + +0.187500 5 0 +0.218750 5 0 + +0.218750 5 0 +0.250000 5 0 + +0.406250 5 0 +0.437500 5 0 + +0.0312500 5 0 +0.0625000 5 0 + +0.0625000 5 0 +0.0937500 5 0 + +0.0937500 5 0 +0.125000 5 0 + +0.687500 5 0 +0.718750 5 0 + +0.781250 5 0 +0.812500 5 0 + +0.437500 5 0 +0.468750 5 0 + +0.468750 5 0 +0.500000 5 0 + +0.906250 5 0 +0.937500 5 0 + +0.531250 5 0 +0.562500 5 0 + +0.968750 5 0 +1.00000 5 0 + +0.812500 5 0 +0.843750 5 0 + +0.843750 5 0 +0.875000 5 0 + +0.375000 6 0 +0.390625 6 0 + +0.390625 6 0 +0.406250 6 0 + +0.156250 6 0 +0.171875 6 0 + +0.171875 6 0 +0.187500 6 0 + +0.0156250 6 0 +0.0312500 6 0 + +0.875000 6 0 +0.890625 6 0 + +0.890625 6 0 +0.906250 6 0 + +0.718750 6 0 +0.734375 6 0 + +0.500000 6 0 +0.515625 6 0 + +0.515625 6 0 +0.531250 6 0 + +0.937500 6 0 +0.953125 6 0 + +0.953125 6 0 +0.968750 6 0 + +0.765625 6 0 +0.781250 6 0 + +0.734375 7 0 +0.742188 7 0 + +0.742188 7 0 +0.750000 7 0 + +0.00781250 7 0 +0.0156250 7 0 + +0.750000 7 0 +0.757812 7 0 + +0.757812 7 0 +0.765625 7 0 + +0.00000 8 0 +0.00390625 8 0 + +0.00390625 8 0 +0.00781250 8 0 + +DEAL:1d::0.281250 4 +DEAL:1d::0.281250 0 +DEAL:1d:: +DEAL:1d::0.656250 4 +DEAL:1d::0.656250 0 +DEAL:1d:: +DEAL:1d::0.0781250 5 +DEAL:1d::0.0781250 0 +DEAL:1d:: +DEAL:1d::0.484375 5 +DEAL:1d::0.484375 0 +DEAL:1d:: +DEAL:1d::0.398438 6 +DEAL:1d::0.398438 0 +DEAL:1d:: +DEAL:1d::0.0234375 6 +DEAL:1d::0.0234375 0 +DEAL:1d:: +DEAL:1d::0.523438 6 +DEAL:1d::0.523438 0 +DEAL:1d:: +DEAL:1d::0.773438 6 +DEAL:1d::0.773438 0 +DEAL:1d:: +DEAL:1d::0.00585938 8 +DEAL:1d::0.00585938 0 +DEAL:1d:: +DEAL:1d::Iteration 12: cells=50 +0.187500 4 0 +0.250000 4 0 + +0.125000 5 0 +0.156250 5 0 + +0.156250 5 0 +0.187500 5 0 + +0.312500 5 0 +0.343750 5 0 + +0.343750 5 0 +0.375000 5 0 + +0.250000 5 0 +0.281250 5 0 + +0.281250 5 0 +0.312500 5 0 + +0.0937500 5 0 +0.125000 5 0 + +0.687500 5 0 +0.718750 5 0 + +0.437500 5 0 +0.468750 5 0 + +0.875000 5 0 +0.906250 5 0 + +0.906250 5 0 +0.937500 5 0 + +0.937500 5 0 +0.968750 5 0 + +0.968750 5 0 +1.00000 5 0 + +0.812500 5 0 +0.843750 5 0 + +0.843750 5 0 +0.875000 5 0 + +0.625000 5 0 +0.656250 5 0 + +0.656250 5 0 +0.687500 5 0 + +0.562500 5 0 +0.593750 5 0 + +0.593750 5 0 +0.625000 5 0 + +0.375000 6 0 +0.390625 6 0 + +0.406250 6 0 +0.421875 6 0 + +0.421875 6 0 +0.437500 6 0 + +0.0312500 6 0 +0.0468750 6 0 + +0.0468750 6 0 +0.0625000 6 0 + +0.718750 6 0 +0.734375 6 0 + +0.734375 6 0 +0.750000 6 0 + +0.500000 6 0 +0.515625 6 0 + +0.0625000 6 0 +0.0781250 6 0 + +0.0781250 6 0 +0.0937500 6 0 + +0.750000 6 0 +0.765625 6 0 + +0.781250 6 0 +0.796875 6 0 + +0.796875 6 0 +0.812500 6 0 + +0.468750 6 0 +0.484375 6 0 + +0.484375 6 0 +0.500000 6 0 + +0.531250 6 0 +0.546875 6 0 + +0.546875 6 0 +0.562500 6 0 + +0.390625 7 0 +0.398438 7 0 + +0.398438 7 0 +0.406250 7 0 + +0.0156250 7 0 +0.0234375 7 0 + +0.0234375 7 0 +0.0312500 7 0 + +0.515625 7 0 +0.523438 7 0 + +0.523438 7 0 +0.531250 7 0 + +0.765625 7 0 +0.773438 7 0 + +0.773438 7 0 +0.781250 7 0 + +0.00000 8 0 +0.00390625 8 0 + +0.00781250 8 0 +0.0117188 8 0 + +0.0117188 8 0 +0.0156250 8 0 + +0.00390625 9 0 +0.00585938 9 0 + +0.00585938 9 0 +0.00781250 9 0 + +DEAL:1d::0.218750 4 +DEAL:1d::0.218750 0 +DEAL:1d:: +DEAL:1d::0.265625 5 +DEAL:1d::0.265625 0 +DEAL:1d:: +DEAL:1d::0.109375 5 +DEAL:1d::0.109375 0 +DEAL:1d:: +DEAL:1d::0.984375 5 +DEAL:1d::0.984375 0 +DEAL:1d:: +DEAL:1d::0.0546875 6 +DEAL:1d::0.0546875 0 +DEAL:1d:: +DEAL:1d::0.726562 6 +DEAL:1d::0.726562 0 +DEAL:1d:: +DEAL:1d::0.507812 6 +DEAL:1d::0.507812 0 +DEAL:1d:: +DEAL:1d::0.0859375 6 +DEAL:1d::0.0859375 0 +DEAL:1d:: +DEAL:1d::0.492188 6 +DEAL:1d::0.492188 0 +DEAL:1d:: +DEAL:1d::0.554688 6 +DEAL:1d::0.554688 0 +DEAL:1d:: +DEAL:1d::0.00195312 8 +DEAL:1d::0.00195312 0 +DEAL:1d:: +DEAL:1d::0.00976562 8 +DEAL:1d::0.00976562 0 +DEAL:1d:: +DEAL:1d::Iteration 13: cells=56 +0.812500 4 0 +0.875000 4 0 + +0.875000 4 0 +0.937500 4 0 + +0.125000 5 0 +0.156250 5 0 + +0.156250 5 0 +0.187500 5 0 + +0.312500 5 0 +0.343750 5 0 + +0.343750 5 0 +0.375000 5 0 + +0.281250 5 0 +0.312500 5 0 + +0.406250 5 0 +0.437500 5 0 + +0.781250 5 0 +0.812500 5 0 + +0.437500 5 0 +0.468750 5 0 + +0.187500 5 0 +0.218750 5 0 + +0.218750 5 0 +0.250000 5 0 + +0.937500 5 0 +0.968750 5 0 + +0.625000 5 0 +0.656250 5 0 + +0.656250 5 0 +0.687500 5 0 + +0.593750 5 0 +0.625000 5 0 + +0.375000 6 0 +0.390625 6 0 + +0.390625 6 0 +0.406250 6 0 + +0.250000 6 0 +0.265625 6 0 + +0.265625 6 0 +0.281250 6 0 + +0.0312500 6 0 +0.0468750 6 0 + +0.734375 6 0 +0.750000 6 0 + +0.515625 6 0 +0.531250 6 0 + +0.0625000 6 0 +0.0781250 6 0 + +0.750000 6 0 +0.765625 6 0 + +0.765625 6 0 +0.781250 6 0 + +0.0937500 6 0 +0.109375 6 0 + +0.109375 6 0 +0.125000 6 0 + +0.468750 6 0 +0.484375 6 0 + +0.531250 6 0 +0.546875 6 0 + +0.687500 6 0 +0.703125 6 0 + +0.703125 6 0 +0.718750 6 0 + +0.968750 6 0 +0.984375 6 0 + +0.984375 6 0 +1.00000 6 0 + +0.562500 6 0 +0.578125 6 0 + +0.578125 6 0 +0.593750 6 0 + +0.0468750 7 0 +0.0546875 7 0 + +0.0546875 7 0 +0.0625000 7 0 + +0.0156250 7 0 +0.0234375 7 0 + +0.0234375 7 0 +0.0312500 7 0 + +0.718750 7 0 +0.726562 7 0 + +0.726562 7 0 +0.734375 7 0 + +0.500000 7 0 +0.507812 7 0 + +0.507812 7 0 +0.515625 7 0 + +0.0781250 7 0 +0.0859375 7 0 + +0.0859375 7 0 +0.0937500 7 0 + +0.484375 7 0 +0.492188 7 0 + +0.492188 7 0 +0.500000 7 0 + +0.546875 7 0 +0.554688 7 0 + +0.554688 7 0 +0.562500 7 0 + +0.00390625 8 0 +0.00781250 8 0 + +0.0117188 8 0 +0.0156250 8 0 + +0.00000 9 0 +0.00195312 9 0 + +0.00195312 9 0 +0.00390625 9 0 + +0.00781250 9 0 +0.00976562 9 0 + +0.00976562 9 0 +0.0117188 9 0 + +DEAL:1d::0.843750 4 +DEAL:1d::0.843750 0 +DEAL:1d:: +DEAL:1d::0.421875 5 +DEAL:1d::0.421875 0 +DEAL:1d:: +DEAL:1d::0.953125 5 +DEAL:1d::0.953125 0 +DEAL:1d:: +DEAL:1d::0.671875 5 +DEAL:1d::0.671875 0 +DEAL:1d:: +DEAL:1d::0.398438 6 +DEAL:1d::0.398438 0 +DEAL:1d:: +DEAL:1d::0.257812 6 +DEAL:1d::0.257812 0 +DEAL:1d:: +DEAL:1d::0.0703125 6 +DEAL:1d::0.0703125 0 +DEAL:1d:: +DEAL:1d::0.757812 6 +DEAL:1d::0.757812 0 +DEAL:1d:: +DEAL:1d::0.476562 6 +DEAL:1d::0.476562 0 +DEAL:1d:: +DEAL:1d::0.695312 6 +DEAL:1d::0.695312 0 +DEAL:1d:: +DEAL:1d::0.0195312 7 +DEAL:1d::0.0195312 0 +DEAL:1d:: +DEAL:1d::0.503906 7 +DEAL:1d::0.503906 0 +DEAL:1d:: +DEAL:1d::0.0136719 8 +DEAL:1d::0.0136719 0 +DEAL:1d:: +DEAL:1d::Iteration 14: cells=62 +0.125000 4 0 +0.187500 4 0 + +0.812500 5 0 +0.843750 5 0 + +0.843750 5 0 +0.875000 5 0 + +0.312500 5 0 +0.343750 5 0 + +0.343750 5 0 +0.375000 5 0 + +0.281250 5 0 +0.312500 5 0 + +0.0937500 5 0 +0.125000 5 0 + +0.781250 5 0 +0.812500 5 0 + +0.187500 5 0 +0.218750 5 0 + +0.968750 5 0 +1.00000 5 0 + +0.875000 5 0 +0.906250 5 0 + +0.906250 5 0 +0.937500 5 0 + +0.625000 5 0 +0.656250 5 0 + +0.562500 5 0 +0.593750 5 0 + +0.593750 5 0 +0.625000 5 0 + +0.375000 6 0 +0.390625 6 0 + +0.265625 6 0 +0.281250 6 0 + +0.0312500 6 0 +0.0468750 6 0 + +0.0468750 6 0 +0.0625000 6 0 + +0.718750 6 0 +0.734375 6 0 + +0.734375 6 0 +0.750000 6 0 + +0.515625 6 0 +0.531250 6 0 + +0.0781250 6 0 +0.0937500 6 0 + +0.765625 6 0 +0.781250 6 0 + +0.406250 6 0 +0.421875 6 0 + +0.421875 6 0 +0.437500 6 0 + +0.531250 6 0 +0.546875 6 0 + +0.546875 6 0 +0.562500 6 0 + +0.703125 6 0 +0.718750 6 0 + +0.437500 6 0 +0.453125 6 0 + +0.453125 6 0 +0.468750 6 0 + +0.218750 6 0 +0.234375 6 0 + +0.234375 6 0 +0.250000 6 0 + +0.937500 6 0 +0.953125 6 0 + +0.953125 6 0 +0.968750 6 0 + +0.656250 6 0 +0.671875 6 0 + +0.671875 6 0 +0.687500 6 0 + +0.390625 7 0 +0.398438 7 0 + +0.398438 7 0 +0.406250 7 0 + +0.0234375 7 0 +0.0312500 7 0 + +0.250000 7 0 +0.257812 7 0 + +0.257812 7 0 +0.265625 7 0 + +0.507812 7 0 +0.515625 7 0 + +0.0625000 7 0 +0.0703125 7 0 + +0.0703125 7 0 +0.0781250 7 0 + +0.484375 7 0 +0.492188 7 0 + +0.492188 7 0 +0.500000 7 0 + +0.750000 7 0 +0.757812 7 0 + +0.757812 7 0 +0.765625 7 0 + +0.468750 7 0 +0.476562 7 0 + +0.476562 7 0 +0.484375 7 0 + +0.687500 7 0 +0.695312 7 0 + +0.695312 7 0 +0.703125 7 0 + +0.00000 8 0 +0.00390625 8 0 + +0.00390625 8 0 +0.00781250 8 0 + +0.00781250 8 0 +0.0117188 8 0 + +0.0156250 8 0 +0.0195312 8 0 + +0.0195312 8 0 +0.0234375 8 0 + +0.500000 8 0 +0.503906 8 0 + +0.503906 8 0 +0.507812 8 0 + +0.0117188 9 0 +0.0136719 9 0 + +0.0136719 9 0 +0.0156250 9 0 + +DEAL:1d::0.156250 4 +DEAL:1d::0.156250 0 +DEAL:1d:: +DEAL:1d::0.984375 5 +DEAL:1d::0.984375 0 +DEAL:1d:: +DEAL:1d::0.890625 5 +DEAL:1d::0.890625 0 +DEAL:1d:: +DEAL:1d::0.273438 6 +DEAL:1d::0.273438 0 +DEAL:1d:: +DEAL:1d::0.0546875 6 +DEAL:1d::0.0546875 0 +DEAL:1d:: +DEAL:1d::0.726562 6 +DEAL:1d::0.726562 0 +DEAL:1d:: +DEAL:1d::0.242188 6 +DEAL:1d::0.242188 0 +DEAL:1d:: +DEAL:1d::0.945312 6 +DEAL:1d::0.945312 0 +DEAL:1d:: +DEAL:1d::0.480469 7 +DEAL:1d::0.480469 0 +DEAL:1d:: +DEAL:1d::0.505859 8 +DEAL:1d::0.505859 0 +DEAL:1d:: +DEAL:1d::Iteration 15: cells=66 +0.812500 5 0 +0.843750 5 0 + +0.843750 5 0 +0.875000 5 0 + +0.312500 5 0 +0.343750 5 0 + +0.343750 5 0 +0.375000 5 0 + +0.406250 5 0 +0.437500 5 0 + +0.0937500 5 0 +0.125000 5 0 + +0.781250 5 0 +0.812500 5 0 + +0.187500 5 0 +0.218750 5 0 + +0.625000 5 0 +0.656250 5 0 + +0.656250 5 0 +0.687500 5 0 + +0.562500 5 0 +0.593750 5 0 + +0.593750 5 0 +0.625000 5 0 + +0.125000 5 0 +0.156250 5 0 + +0.156250 5 0 +0.187500 5 0 + +0.375000 6 0 +0.390625 6 0 + +0.390625 6 0 +0.406250 6 0 + +0.250000 6 0 +0.265625 6 0 + +0.0312500 6 0 +0.0468750 6 0 + +0.734375 6 0 +0.750000 6 0 + +0.0625000 6 0 +0.0781250 6 0 + +0.0781250 6 0 +0.0937500 6 0 + +0.750000 6 0 +0.765625 6 0 + +0.765625 6 0 +0.781250 6 0 + +0.281250 6 0 +0.296875 6 0 + +0.296875 6 0 +0.312500 6 0 + +0.531250 6 0 +0.546875 6 0 + +0.546875 6 0 +0.562500 6 0 + +0.687500 6 0 +0.703125 6 0 + +0.703125 6 0 +0.718750 6 0 + +0.437500 6 0 +0.453125 6 0 + +0.453125 6 0 +0.468750 6 0 + +0.218750 6 0 +0.234375 6 0 + +0.953125 6 0 +0.968750 6 0 + +0.968750 6 0 +0.984375 6 0 + +0.984375 6 0 +1.00000 6 0 + +0.875000 6 0 +0.890625 6 0 + +0.890625 6 0 +0.906250 6 0 + +0.906250 6 0 +0.921875 6 0 + +0.921875 6 0 +0.937500 6 0 + +0.265625 7 0 +0.273438 7 0 + +0.273438 7 0 +0.281250 7 0 + +0.00000 7 0 +0.00781250 7 0 + +0.0156250 7 0 +0.0234375 7 0 + +0.0234375 7 0 +0.0312500 7 0 + +0.0468750 7 0 +0.0546875 7 0 + +0.0546875 7 0 +0.0625000 7 0 + +0.718750 7 0 +0.726562 7 0 + +0.726562 7 0 +0.734375 7 0 + +0.484375 7 0 +0.492188 7 0 + +0.492188 7 0 +0.500000 7 0 + +0.515625 7 0 +0.523438 7 0 + +0.523438 7 0 +0.531250 7 0 + +0.468750 7 0 +0.476562 7 0 + +0.234375 7 0 +0.242188 7 0 + +0.242188 7 0 +0.250000 7 0 + +0.937500 7 0 +0.945312 7 0 + +0.945312 7 0 +0.953125 7 0 + +0.507812 8 0 +0.511719 8 0 + +0.511719 8 0 +0.515625 8 0 + +0.00781250 8 0 +0.0117188 8 0 + +0.0117188 8 0 +0.0156250 8 0 + +0.476562 8 0 +0.480469 8 0 + +0.480469 8 0 +0.484375 8 0 + +0.500000 8 0 +0.503906 8 0 + +0.503906 9 0 +0.505859 9 0 + +0.505859 9 0 +0.507812 9 0 + +DEAL:1d::0.828125 5 +DEAL:1d::0.828125 0 +DEAL:1d:: +DEAL:1d::0.109375 5 +DEAL:1d::0.109375 0 +DEAL:1d:: +DEAL:1d::0.640625 5 +DEAL:1d::0.640625 0 +DEAL:1d:: +DEAL:1d::0.257812 6 +DEAL:1d::0.257812 0 +DEAL:1d:: +DEAL:1d::0.460938 6 +DEAL:1d::0.460938 0 +DEAL:1d:: +DEAL:1d::0.277344 7 +DEAL:1d::0.277344 0 +DEAL:1d:: +DEAL:1d::0.0195312 7 +DEAL:1d::0.0195312 0 +DEAL:1d:: +DEAL:1d::0.527344 7 +DEAL:1d::0.527344 0 +DEAL:1d:: +DEAL:1d::0.238281 7 +DEAL:1d::0.238281 0 +DEAL:1d:: +DEAL:1d::0.513672 8 +DEAL:1d::0.513672 0 +DEAL:1d:: +DEAL:1d::0.00976562 8 +DEAL:1d::0.00976562 0 +DEAL:1d:: +DEAL:1d::0.0136719 8 +DEAL:1d::0.0136719 0 +DEAL:1d:: +DEAL:1d::0.482422 8 +DEAL:1d::0.482422 0 +DEAL:1d:: +DEAL:1d::0.501953 8 +DEAL:1d::0.501953 0 +DEAL:1d:: +DEAL:1d::Iteration 16: cells=77 +0.843750 5 0 +0.875000 5 0 + +0.312500 5 0 +0.343750 5 0 + +0.343750 5 0 +0.375000 5 0 + +0.375000 5 0 +0.406250 5 0 + +0.406250 5 0 +0.437500 5 0 + +0.0625000 5 0 +0.0937500 5 0 + +0.687500 5 0 +0.718750 5 0 + +0.750000 5 0 +0.781250 5 0 + +0.781250 5 0 +0.812500 5 0 + +0.968750 5 0 +1.00000 5 0 + +0.875000 5 0 +0.906250 5 0 + +0.906250 5 0 +0.937500 5 0 + +0.656250 5 0 +0.687500 5 0 + +0.562500 5 0 +0.593750 5 0 + +0.593750 5 0 +0.625000 5 0 + +0.125000 5 0 +0.156250 5 0 + +0.156250 5 0 +0.187500 5 0 + +0.812500 6 0 +0.828125 6 0 + +0.828125 6 0 +0.843750 6 0 + +0.0312500 6 0 +0.0468750 6 0 + +0.0468750 6 0 +0.0625000 6 0 + +0.718750 6 0 +0.734375 6 0 + +0.734375 6 0 +0.750000 6 0 + +0.0937500 6 0 +0.109375 6 0 + +0.109375 6 0 +0.125000 6 0 + +0.187500 6 0 +0.203125 6 0 + +0.203125 6 0 +0.218750 6 0 + +0.296875 6 0 +0.312500 6 0 + +0.546875 6 0 +0.562500 6 0 + +0.625000 6 0 +0.640625 6 0 + +0.640625 6 0 +0.656250 6 0 + +0.437500 6 0 +0.453125 6 0 + +0.937500 6 0 +0.953125 6 0 + +0.953125 6 0 +0.968750 6 0 + +0.265625 7 0 +0.273438 7 0 + +0.0234375 7 0 +0.0312500 7 0 + +0.250000 7 0 +0.257812 7 0 + +0.257812 7 0 +0.265625 7 0 + +0.281250 7 0 +0.289062 7 0 + +0.289062 7 0 +0.296875 7 0 + +0.468750 7 0 +0.476562 7 0 + +0.242188 7 0 +0.250000 7 0 + +0.531250 7 0 +0.539062 7 0 + +0.539062 7 0 +0.546875 7 0 + +0.453125 7 0 +0.460938 7 0 + +0.460938 7 0 +0.468750 7 0 + +0.218750 7 0 +0.226562 7 0 + +0.226562 7 0 +0.234375 7 0 + +0.507812 8 0 +0.511719 8 0 + +0.476562 8 0 +0.480469 8 0 + +0.503906 8 0 +0.507812 8 0 + +0.273438 8 0 +0.277344 8 0 + +0.277344 8 0 +0.281250 8 0 + +0.00000 8 0 +0.00390625 8 0 + +0.00390625 8 0 +0.00781250 8 0 + +0.0156250 8 0 +0.0195312 8 0 + +0.0195312 8 0 +0.0234375 8 0 + +0.484375 8 0 +0.488281 8 0 + +0.488281 8 0 +0.492188 8 0 + +0.492188 8 0 +0.496094 8 0 + +0.496094 8 0 +0.500000 8 0 + +0.515625 8 0 +0.519531 8 0 + +0.519531 8 0 +0.523438 8 0 + +0.523438 8 0 +0.527344 8 0 + +0.527344 8 0 +0.531250 8 0 + +0.234375 8 0 +0.238281 8 0 + +0.238281 8 0 +0.242188 8 0 + +0.511719 9 0 +0.513672 9 0 + +0.513672 9 0 +0.515625 9 0 + +0.00781250 9 0 +0.00976562 9 0 + +0.00976562 9 0 +0.0117188 9 0 + +0.0117188 9 0 +0.0136719 9 0 + +0.0136719 9 0 +0.0156250 9 0 + +0.480469 9 0 +0.482422 9 0 + +0.482422 9 0 +0.484375 9 0 + +0.500000 9 0 +0.501953 9 0 + +0.501953 9 0 +0.503906 9 0 + +DEAL:1d::0.859375 5 +DEAL:1d::0.859375 0 +DEAL:1d:: +DEAL:1d::0.328125 5 +DEAL:1d::0.328125 0 +DEAL:1d:: +DEAL:1d::0.890625 5 +DEAL:1d::0.890625 0 +DEAL:1d:: +DEAL:1d::0.117188 6 +DEAL:1d::0.117188 0 +DEAL:1d:: +DEAL:1d::0.195312 6 +DEAL:1d::0.195312 0 +DEAL:1d:: +DEAL:1d::0.210938 6 +DEAL:1d::0.210938 0 +DEAL:1d:: +DEAL:1d::0.554688 6 +DEAL:1d::0.554688 0 +DEAL:1d:: +DEAL:1d::0.945312 6 +DEAL:1d::0.945312 0 +DEAL:1d:: +DEAL:1d::0.285156 7 +DEAL:1d::0.285156 0 +DEAL:1d:: +DEAL:1d::0.246094 7 +DEAL:1d::0.246094 0 +DEAL:1d:: +DEAL:1d::0.275391 8 +DEAL:1d::0.275391 0 +DEAL:1d:: +DEAL:1d::0.0175781 8 +DEAL:1d::0.0175781 0 +DEAL:1d:: +DEAL:1d::0.494141 8 +DEAL:1d::0.494141 0 +DEAL:1d:: +DEAL:1d::0.529297 8 +DEAL:1d::0.529297 0 +DEAL:1d:: +DEAL:1d::0.502930 9 +DEAL:1d::0.502930 0 +DEAL:1d:: +DEAL:1d::Iteration 17: cells=86 +0.750000 4 0 +0.812500 4 0 + +0.812500 5 0 +0.843750 5 0 + +0.343750 5 0 +0.375000 5 0 + +0.375000 5 0 +0.406250 5 0 + +0.406250 5 0 +0.437500 5 0 + +0.0625000 5 0 +0.0937500 5 0 + +0.687500 5 0 +0.718750 5 0 + +0.718750 5 0 +0.750000 5 0 + +0.968750 5 0 +1.00000 5 0 + +0.625000 5 0 +0.656250 5 0 + +0.656250 5 0 +0.687500 5 0 + +0.593750 5 0 +0.625000 5 0 + +0.843750 6 0 +0.859375 6 0 + +0.859375 6 0 +0.875000 6 0 + +0.0312500 6 0 +0.0468750 6 0 + +0.0468750 6 0 +0.0625000 6 0 + +0.312500 6 0 +0.328125 6 0 + +0.328125 6 0 +0.343750 6 0 + +0.0937500 6 0 +0.109375 6 0 + +0.296875 6 0 +0.312500 6 0 + +0.875000 6 0 +0.890625 6 0 + +0.890625 6 0 +0.906250 6 0 + +0.437500 6 0 +0.453125 6 0 + +0.453125 6 0 +0.468750 6 0 + +0.218750 6 0 +0.234375 6 0 + +0.953125 6 0 +0.968750 6 0 + +0.906250 6 0 +0.921875 6 0 + +0.921875 6 0 +0.937500 6 0 + +0.562500 6 0 +0.578125 6 0 + +0.578125 6 0 +0.593750 6 0 + +0.125000 6 0 +0.140625 6 0 + +0.140625 6 0 +0.156250 6 0 + +0.156250 6 0 +0.171875 6 0 + +0.171875 6 0 +0.187500 6 0 + +0.00000 7 0 +0.00781250 7 0 + +0.0234375 7 0 +0.0312500 7 0 + +0.250000 7 0 +0.257812 7 0 + +0.257812 7 0 +0.265625 7 0 + +0.289062 7 0 +0.296875 7 0 + +0.515625 7 0 +0.523438 7 0 + +0.468750 7 0 +0.476562 7 0 + +0.234375 7 0 +0.242188 7 0 + +0.539062 7 0 +0.546875 7 0 + +0.109375 7 0 +0.117188 7 0 + +0.117188 7 0 +0.125000 7 0 + +0.187500 7 0 +0.195312 7 0 + +0.195312 7 0 +0.203125 7 0 + +0.203125 7 0 +0.210938 7 0 + +0.210938 7 0 +0.218750 7 0 + +0.546875 7 0 +0.554688 7 0 + +0.554688 7 0 +0.562500 7 0 + +0.937500 7 0 +0.945312 7 0 + +0.945312 7 0 +0.953125 7 0 + +0.507812 8 0 +0.511719 8 0 + +0.511719 8 0 +0.515625 8 0 + +0.00781250 8 0 +0.0117188 8 0 + +0.0117188 8 0 +0.0156250 8 0 + +0.476562 8 0 +0.480469 8 0 + +0.480469 8 0 +0.484375 8 0 + +0.277344 8 0 +0.281250 8 0 + +0.265625 8 0 +0.269531 8 0 + +0.269531 8 0 +0.273438 8 0 + +0.0195312 8 0 +0.0234375 8 0 + +0.484375 8 0 +0.488281 8 0 + +0.488281 8 0 +0.492188 8 0 + +0.496094 8 0 +0.500000 8 0 + +0.281250 8 0 +0.285156 8 0 + +0.285156 8 0 +0.289062 8 0 + +0.523438 8 0 +0.527344 8 0 + +0.242188 8 0 +0.246094 8 0 + +0.246094 8 0 +0.250000 8 0 + +0.531250 8 0 +0.535156 8 0 + +0.535156 8 0 +0.539062 8 0 + +0.503906 9 0 +0.505859 9 0 + +0.505859 9 0 +0.507812 9 0 + +0.273438 9 0 +0.275391 9 0 + +0.275391 9 0 +0.277344 9 0 + +0.0156250 9 0 +0.0175781 9 0 + +0.0175781 9 0 +0.0195312 9 0 + +0.492188 9 0 +0.494141 9 0 + +0.494141 9 0 +0.496094 9 0 + +0.500000 9 0 +0.501953 9 0 + +0.527344 9 0 +0.529297 9 0 + +0.529297 9 0 +0.531250 9 0 + +0.501953 10 0 +0.502930 10 0 + +0.502930 10 0 +0.503906 10 0 + +DEAL:1d::0.781250 4 +DEAL:1d::0.781250 0 +DEAL:1d:: +DEAL:1d::0.304688 6 +DEAL:1d::0.304688 0 +DEAL:1d:: +DEAL:1d::0.882812 6 +DEAL:1d::0.882812 0 +DEAL:1d:: +DEAL:1d::0.226562 6 +DEAL:1d::0.226562 0 +DEAL:1d:: +DEAL:1d::0.960938 6 +DEAL:1d::0.960938 0 +DEAL:1d:: +DEAL:1d::0.261719 7 +DEAL:1d::0.261719 0 +DEAL:1d:: +DEAL:1d::0.238281 7 +DEAL:1d::0.238281 0 +DEAL:1d:: +DEAL:1d::0.113281 7 +DEAL:1d::0.113281 0 +DEAL:1d:: +DEAL:1d::0.121094 7 +DEAL:1d::0.121094 0 +DEAL:1d:: +DEAL:1d::0.214844 7 +DEAL:1d::0.214844 0 +DEAL:1d:: +DEAL:1d::0.941406 7 +DEAL:1d::0.941406 0 +DEAL:1d:: +DEAL:1d::0.486328 8 +DEAL:1d::0.486328 0 +DEAL:1d:: +DEAL:1d::0.533203 8 +DEAL:1d::0.533203 0 +DEAL:1d:: +DEAL:1d::0.276367 9 +DEAL:1d::0.276367 0 +DEAL:1d:: +DEAL:1d::0.0166016 9 +DEAL:1d::0.0166016 0 +DEAL:1d:: +DEAL:1d::0.500977 9 +DEAL:1d::0.500977 0 +DEAL:1d:: +DEAL:1d::Iteration 18: cells=98 +0.625000 4 0 +0.687500 4 0 + +0.687500 4 0 +0.750000 4 0 + +0.812500 5 0 +0.843750 5 0 + +0.343750 5 0 +0.375000 5 0 + +0.375000 5 0 +0.406250 5 0 + +0.406250 5 0 +0.437500 5 0 + +0.750000 5 0 +0.781250 5 0 + +0.781250 5 0 +0.812500 5 0 + +0.562500 5 0 +0.593750 5 0 + +0.593750 5 0 +0.625000 5 0 + +0.156250 5 0 +0.187500 5 0 + +0.843750 6 0 +0.859375 6 0 + +0.859375 6 0 +0.875000 6 0 + +0.0312500 6 0 +0.0468750 6 0 + +0.0468750 6 0 +0.0625000 6 0 + +0.312500 6 0 +0.328125 6 0 + +0.328125 6 0 +0.343750 6 0 + +0.187500 6 0 +0.203125 6 0 + +0.546875 6 0 +0.562500 6 0 + +0.890625 6 0 +0.906250 6 0 + +0.437500 6 0 +0.453125 6 0 + +0.453125 6 0 +0.468750 6 0 + +0.906250 6 0 +0.921875 6 0 + +0.0625000 6 0 +0.0781250 6 0 + +0.0781250 6 0 +0.0937500 6 0 + +0.140625 6 0 +0.156250 6 0 + +0.968750 6 0 +0.984375 6 0 + +0.984375 6 0 +1.00000 6 0 + +0.00000 7 0 +0.00781250 7 0 + +0.0234375 7 0 +0.0312500 7 0 + +0.250000 7 0 +0.257812 7 0 + +0.507812 7 0 +0.515625 7 0 + +0.289062 7 0 +0.296875 7 0 + +0.515625 7 0 +0.523438 7 0 + +0.468750 7 0 +0.476562 7 0 + +0.242188 7 0 +0.250000 7 0 + +0.539062 7 0 +0.546875 7 0 + +0.0937500 7 0 +0.101562 7 0 + +0.101562 7 0 +0.109375 7 0 + +0.203125 7 0 +0.210938 7 0 + +0.296875 7 0 +0.304688 7 0 + +0.304688 7 0 +0.312500 7 0 + +0.945312 7 0 +0.953125 7 0 + +0.875000 7 0 +0.882812 7 0 + +0.882812 7 0 +0.890625 7 0 + +0.218750 7 0 +0.226562 7 0 + +0.226562 7 0 +0.234375 7 0 + +0.953125 7 0 +0.960938 7 0 + +0.960938 7 0 +0.968750 7 0 + +0.921875 7 0 +0.929688 7 0 + +0.929688 7 0 +0.937500 7 0 + +0.125000 7 0 +0.132812 7 0 + +0.132812 7 0 +0.140625 7 0 + +0.257812 8 0 +0.261719 8 0 + +0.261719 8 0 +0.265625 8 0 + +0.00781250 8 0 +0.0117188 8 0 + +0.476562 8 0 +0.480469 8 0 + +0.480469 8 0 +0.484375 8 0 + +0.503906 8 0 +0.507812 8 0 + +0.265625 8 0 +0.269531 8 0 + +0.269531 8 0 +0.273438 8 0 + +0.0195312 8 0 +0.0234375 8 0 + +0.488281 8 0 +0.492188 8 0 + +0.492188 8 0 +0.496094 8 0 + +0.281250 8 0 +0.285156 8 0 + +0.285156 8 0 +0.289062 8 0 + +0.523438 8 0 +0.527344 8 0 + +0.527344 8 0 +0.531250 8 0 + +0.234375 8 0 +0.238281 8 0 + +0.238281 8 0 +0.242188 8 0 + +0.535156 8 0 +0.539062 8 0 + +0.109375 8 0 +0.113281 8 0 + +0.113281 8 0 +0.117188 8 0 + +0.117188 8 0 +0.121094 8 0 + +0.121094 8 0 +0.125000 8 0 + +0.210938 8 0 +0.214844 8 0 + +0.214844 8 0 +0.218750 8 0 + +0.937500 8 0 +0.941406 8 0 + +0.941406 8 0 +0.945312 8 0 + +0.0117188 9 0 +0.0136719 9 0 + +0.0136719 9 0 +0.0156250 9 0 + +0.273438 9 0 +0.275391 9 0 + +0.0175781 9 0 +0.0195312 9 0 + +0.277344 9 0 +0.279297 9 0 + +0.279297 9 0 +0.281250 9 0 + +0.501953 9 0 +0.503906 9 0 + +0.484375 9 0 +0.486328 9 0 + +0.486328 9 0 +0.488281 9 0 + +0.496094 9 0 +0.498047 9 0 + +0.498047 9 0 +0.500000 9 0 + +0.531250 9 0 +0.533203 9 0 + +0.533203 9 0 +0.535156 9 0 + +0.275391 10 0 +0.276367 10 0 + +0.276367 10 0 +0.277344 10 0 + +0.0156250 10 0 +0.0166016 10 0 + +0.0166016 10 0 +0.0175781 10 0 + +0.500000 10 0 +0.500977 10 0 + +0.500977 10 0 +0.501953 10 0 + +DEAL:1d::0.656250 4 +DEAL:1d::0.656250 0 +DEAL:1d:: +DEAL:1d::0.359375 5 +DEAL:1d::0.359375 0 +DEAL:1d:: +DEAL:1d::0.421875 5 +DEAL:1d::0.421875 0 +DEAL:1d:: +DEAL:1d::0.0546875 6 +DEAL:1d::0.0546875 0 +DEAL:1d:: +DEAL:1d::0.914062 6 +DEAL:1d::0.914062 0 +DEAL:1d:: +DEAL:1d::0.0859375 6 +DEAL:1d::0.0859375 0 +DEAL:1d:: +DEAL:1d::0.148438 6 +DEAL:1d::0.148438 0 +DEAL:1d:: +DEAL:1d::0.0976562 7 +DEAL:1d::0.0976562 0 +DEAL:1d:: +DEAL:1d::0.207031 7 +DEAL:1d::0.207031 0 +DEAL:1d:: +DEAL:1d::0.886719 7 +DEAL:1d::0.886719 0 +DEAL:1d:: +DEAL:1d::0.230469 7 +DEAL:1d::0.230469 0 +DEAL:1d:: +DEAL:1d::0.259766 8 +DEAL:1d::0.259766 0 +DEAL:1d:: +DEAL:1d::0.00976562 8 +DEAL:1d::0.00976562 0 +DEAL:1d:: +DEAL:1d::0.0214844 8 +DEAL:1d::0.0214844 0 +DEAL:1d:: +DEAL:1d::0.525391 8 +DEAL:1d::0.525391 0 +DEAL:1d:: +DEAL:1d::0.216797 8 +DEAL:1d::0.216797 0 +DEAL:1d:: +DEAL:1d::0.274414 9 +DEAL:1d::0.274414 0 +DEAL:1d:: +DEAL:1d::0.497070 9 +DEAL:1d::0.497070 0 +DEAL:1d:: +DEAL:1d::0.534180 9 +DEAL:1d::0.534180 0 +DEAL:1d:: +DEAL:1d::0.0161133 10 +DEAL:1d::0.0161133 0 +DEAL:1d:: +DEAL:1d::Iteration 19: cells=117 +DEAL:1d::OK