--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2015 - 2016 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+
+// Test GridTools::build_triangulation_from_patch ()
+
+
+
+#include "../tests.h"
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/tensor.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_tools.h>
+
+#include <fstream>
+
+
+template<int dim>
+void test()
+{
+ Triangulation<dim> triangulation (Triangulation<dim>::limit_level_difference_at_vertices);
+
+ GridGenerator::hyper_cube(triangulation);
+ triangulation.refine_global (1);
+
+
+ unsigned int index = 0;
+ for (typename Triangulation<dim>::active_cell_iterator
+ cell = triangulation.begin_active();
+ cell != triangulation.end(); ++cell, ++index)
+ {
+ std::vector<typename Triangulation<dim>::active_cell_iterator> patch_cells
+ = GridTools::get_patch_around_cell<Triangulation<dim> > (cell);
+
+ Triangulation<dim> local_triangulation(Triangulation<dim>::limit_level_difference_at_vertices);
+ std::map<typename Triangulation<dim>::active_cell_iterator ,
+ typename Triangulation<dim>::active_cell_iterator> patch_to_global_tria_map;
+
+ GridTools::build_triangulation_from_patch <Triangulation<dim> >
+ (patch_cells,local_triangulation,patch_to_global_tria_map);
+
+ deallog << "patch_cells " << cell << ": ";
+ for (unsigned int i=0; i<patch_cells.size(); ++i)
+ deallog << patch_cells[i] << ' ';
+ deallog << std::endl;
+
+ deallog << "local_triangulation " << cell << ": ";
+ for (typename Triangulation<dim>::active_cell_iterator
+ tria_cell = local_triangulation.begin_active();
+ tria_cell != local_triangulation.end(); ++tria_cell)
+ deallog << tria_cell << ' ';
+ deallog << std::endl;
+
+ }
+
+}
+
+
+int main()
+{
+ initlog();
+ deallog.threshold_double(1.e-10);
+
+ deallog.push("1d");
+ test<1>();
+ deallog.pop();
+
+ deallog.push("2d");
+ test<2>();
+ deallog.pop();
+
+ deallog.push("3d");
+ test<3>();
+ deallog.pop();
+}
--- /dev/null
+
+DEAL:1d::patch_cells 1.0: 1.0 1.1
+DEAL:1d::local_triangulation 1.0: 0.0 0.1
+DEAL:1d::patch_cells 1.1: 1.1 1.0
+DEAL:1d::local_triangulation 1.1: 0.0 0.1
+DEAL:2d::patch_cells 1.0: 1.0 1.1 1.2
+DEAL:2d::local_triangulation 1.0: 0.0 0.1 0.2
+DEAL:2d::patch_cells 1.1: 1.1 1.0 1.3
+DEAL:2d::local_triangulation 1.1: 0.0 0.1 0.2
+DEAL:2d::patch_cells 1.2: 1.2 1.3 1.0
+DEAL:2d::local_triangulation 1.2: 0.0 0.1 0.2
+DEAL:2d::patch_cells 1.3: 1.3 1.2 1.1
+DEAL:2d::local_triangulation 1.3: 0.0 0.1 0.2
+DEAL:3d::patch_cells 1.0: 1.0 1.1 1.2 1.4
+DEAL:3d::local_triangulation 1.0: 0.0 0.1 0.2 0.3
+DEAL:3d::patch_cells 1.1: 1.1 1.0 1.3 1.5
+DEAL:3d::local_triangulation 1.1: 0.0 0.1 0.2 0.3
+DEAL:3d::patch_cells 1.2: 1.2 1.3 1.0 1.6
+DEAL:3d::local_triangulation 1.2: 0.0 0.1 0.2 0.3
+DEAL:3d::patch_cells 1.3: 1.3 1.2 1.1 1.7
+DEAL:3d::local_triangulation 1.3: 0.0 0.1 0.2 0.3
+DEAL:3d::patch_cells 1.4: 1.4 1.5 1.6 1.0
+DEAL:3d::local_triangulation 1.4: 0.0 0.1 0.2 0.3
+DEAL:3d::patch_cells 1.5: 1.5 1.4 1.7 1.1
+DEAL:3d::local_triangulation 1.5: 0.0 0.1 0.2 0.3
+DEAL:3d::patch_cells 1.6: 1.6 1.7 1.4 1.2
+DEAL:3d::local_triangulation 1.6: 0.0 0.1 0.2 0.3
+DEAL:3d::patch_cells 1.7: 1.7 1.6 1.5 1.3
+DEAL:3d::local_triangulation 1.7: 0.0 0.1 0.2 0.3
--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2015 - 2016 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+
+// Test GridTools::build_triangulation_from_patch ()
+
+
+
+#include "../tests.h"
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/tensor.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_tools.h>
+
+#include <fstream>
+
+
+template<int dim>
+void test()
+{
+ Triangulation<dim> triangulation (Triangulation<dim>::limit_level_difference_at_vertices);
+
+ GridGenerator::hyper_cube(triangulation);
+ triangulation.refine_global (1);
+ triangulation.begin_active()->set_refine_flag();
+ triangulation.execute_coarsening_and_refinement();
+
+
+ unsigned int index = 0;
+ for (typename Triangulation<dim>::active_cell_iterator
+ cell = triangulation.begin_active();
+ cell != triangulation.end(); ++cell, ++index)
+ {
+ std::vector<typename Triangulation<dim>::active_cell_iterator> patch_cells
+ = GridTools::get_patch_around_cell<Triangulation<dim> > (cell);
+
+ Triangulation<dim> local_triangulation(Triangulation<dim>::limit_level_difference_at_vertices);
+ std::map<typename Triangulation<dim>::active_cell_iterator , typename Triangulation<dim>::active_cell_iterator> patch_to_global_tria_map;
+
+ GridTools::build_triangulation_from_patch <Triangulation<dim> >
+ (patch_cells,local_triangulation,patch_to_global_tria_map);
+
+ deallog << "patch_cells " << cell << ": ";
+ for (unsigned int i=0; i<patch_cells.size(); ++i)
+ deallog << patch_cells[i] << ' ';
+ deallog << std::endl;
+
+ deallog << "local_triangulation " << cell << ": ";
+ for (typename Triangulation<dim>::active_cell_iterator
+ tria_cell = local_triangulation.begin_active();
+ tria_cell != local_triangulation.end(); ++tria_cell)
+ deallog << tria_cell << ' ';
+ deallog << std::endl;
+
+ }
+
+
+}
+
+
+int main()
+{
+ initlog();
+ deallog.threshold_double(1.e-10);
+
+ deallog.push("1d");
+ test<1>();
+ deallog.pop();
+
+ deallog.push("2d");
+ test<2>();
+ deallog.pop();
+
+ deallog.push("3d");
+ test<3>();
+ deallog.pop();
+}
--- /dev/null
+
+DEAL:1d::patch_cells 1.1: 1.1 2.1
+DEAL:1d::local_triangulation 1.1: 0.1 1.0 1.1
+DEAL:1d::patch_cells 2.0: 2.0 2.1
+DEAL:1d::local_triangulation 2.0: 0.0 0.1
+DEAL:1d::patch_cells 2.1: 2.1 2.0 1.1
+DEAL:1d::local_triangulation 2.1: 0.1 1.0 1.1
+DEAL:2d::patch_cells 1.1: 1.1 2.1 2.3 1.3
+DEAL:2d::local_triangulation 1.1: 0.1 0.2 1.0 1.1 1.2 1.3
+DEAL:2d::patch_cells 1.2: 1.2 1.3 2.2 2.3
+DEAL:2d::local_triangulation 1.2: 0.1 0.2 1.0 1.1 1.2 1.3
+DEAL:2d::patch_cells 1.3: 1.3 1.2 1.1
+DEAL:2d::local_triangulation 1.3: 0.0 0.1 0.2
+DEAL:2d::patch_cells 2.0: 2.0 2.1 2.2
+DEAL:2d::local_triangulation 2.0: 0.0 0.1 0.2
+DEAL:2d::patch_cells 2.1: 2.1 2.0 1.1 2.3
+DEAL:2d::local_triangulation 2.1: 0.1 1.0 1.1 1.2 1.3
+DEAL:2d::patch_cells 2.2: 2.2 2.3 2.0 1.2
+DEAL:2d::local_triangulation 2.2: 0.1 1.0 1.1 1.2 1.3
+DEAL:2d::patch_cells 2.3: 2.3 2.2 1.1 2.1 1.2
+DEAL:2d::local_triangulation 2.3: 0.1 0.2 1.0 1.1 1.2 1.3
+DEAL:3d::patch_cells 1.1: 1.1 2.1 2.3 2.5 2.7 1.3 1.5
+DEAL:3d::local_triangulation 1.1: 0.1 0.2 0.3 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7
+DEAL:3d::patch_cells 1.2: 1.2 1.3 2.2 2.6 2.3 2.7 1.6
+DEAL:3d::local_triangulation 1.2: 0.1 0.2 0.3 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7
+DEAL:3d::patch_cells 1.3: 1.3 1.2 1.1 1.7
+DEAL:3d::local_triangulation 1.3: 0.0 0.1 0.2 0.3
+DEAL:3d::patch_cells 1.4: 1.4 1.5 1.6 2.4 2.5 2.6 2.7
+DEAL:3d::local_triangulation 1.4: 0.1 0.2 0.3 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7
+DEAL:3d::patch_cells 1.5: 1.5 1.4 1.7 1.1
+DEAL:3d::local_triangulation 1.5: 0.0 0.1 0.2 0.3
+DEAL:3d::patch_cells 1.6: 1.6 1.7 1.4 1.2
+DEAL:3d::local_triangulation 1.6: 0.0 0.1 0.2 0.3
+DEAL:3d::patch_cells 1.7: 1.7 1.6 1.5 1.3
+DEAL:3d::local_triangulation 1.7: 0.0 0.1 0.2 0.3
+DEAL:3d::patch_cells 2.0: 2.0 2.1 2.2 2.4
+DEAL:3d::local_triangulation 2.0: 0.0 0.1 0.2 0.3
+DEAL:3d::patch_cells 2.1: 2.1 2.0 1.1 2.3 2.5
+DEAL:3d::local_triangulation 2.1: 0.1 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7
+DEAL:3d::patch_cells 2.2: 2.2 2.3 2.0 1.2 2.6
+DEAL:3d::local_triangulation 2.2: 0.1 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7
+DEAL:3d::patch_cells 2.3: 2.3 2.2 1.1 2.1 1.2 2.7
+DEAL:3d::local_triangulation 2.3: 0.1 0.2 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7
+DEAL:3d::patch_cells 2.4: 2.4 2.5 2.6 2.0 1.4
+DEAL:3d::local_triangulation 2.4: 0.1 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7
+DEAL:3d::patch_cells 2.5: 2.5 2.4 1.1 2.7 2.1 1.4
+DEAL:3d::local_triangulation 2.5: 0.1 0.2 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7
+DEAL:3d::patch_cells 2.6: 2.6 2.7 2.4 1.2 2.2 1.4
+DEAL:3d::local_triangulation 2.6: 0.1 0.2 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7
+DEAL:3d::patch_cells 2.7: 2.7 2.6 1.1 2.5 1.2 2.3 1.4
+DEAL:3d::local_triangulation 2.7: 0.1 0.2 0.3 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7
--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2015 - 2016 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// Test GridTools::get_cells_at_coarsest_common_level()
+
+
+#include "../tests.h"
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/tensor.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_tools.h>
+
+#include <fstream>
+
+
+
+
+template<int dim>
+void test()
+{
+ Triangulation<dim> triangulation (Triangulation<dim>::limit_level_difference_at_vertices);
+
+ GridGenerator::hyper_cube(triangulation);
+ triangulation.refine_global (2);
+
+ const unsigned int n_refinements[] = { 0, 4, 3, 2 };
+ for (unsigned int i=0; i<n_refinements[dim]; ++i)
+ {
+ // refine one-fifth of cells randomly
+ std::vector<bool> flags (triangulation.n_active_cells(), false);
+ for (unsigned int k=0; k<flags.size()/5 + 1; ++k)
+ flags[Testing::rand() % flags.size()] = true;
+ // make sure there's at least one that
+ // will be refined
+ flags[0] = true;
+
+ // refine triangulation
+ unsigned int index=0;
+ for (typename Triangulation<dim>::active_cell_iterator
+ cell = triangulation.begin_active();
+ cell != triangulation.end(); ++cell, ++index)
+ if (flags[index])
+ cell->set_refine_flag();
+ AssertThrow (index == triangulation.n_active_cells(), ExcInternalError());
+
+ // flag all other cells for coarsening
+ // (this should ensure that at least
+ // some of them will actually be
+ // coarsened)
+ index=0;
+ for (typename Triangulation<dim>::active_cell_iterator
+ cell = triangulation.begin_active();
+ cell != triangulation.end(); ++cell, ++index)
+ if (!flags[index])
+ cell->set_coarsen_flag();
+
+ triangulation.execute_coarsening_and_refinement ();
+ }
+
+ // now extract patches and print every fifth of them
+ unsigned int index = 0;
+ for (typename Triangulation<dim>::active_cell_iterator
+ cell = triangulation.begin_active();
+ cell != triangulation.end(); ++cell, ++index)
+ {
+ std::vector<typename Triangulation<dim>::active_cell_iterator> patch_cells
+ = GridTools::get_patch_around_cell<Triangulation<dim> > (cell);
+
+ std::vector<typename Triangulation<dim>::cell_iterator> coarse_cells
+ = GridTools::get_cells_at_coarsest_common_level <Triangulation<dim> > (patch_cells);
+
+ if (index % 5 == 0)
+ {
+ deallog << "coarse_ cells " << cell << ": ";
+ for (unsigned int i=0; i<coarse_cells.size(); ++i)
+ deallog << coarse_cells[i] << ' ';
+ deallog << std::endl;
+ }
+ }
+
+
+}
+
+
+int main()
+{
+ initlog();
+ deallog.threshold_double(1.e-10);
+
+ deallog.push("1d");
+ test<1>();
+ deallog.pop();
+
+ deallog.push("2d");
+ test<2>();
+ deallog.pop();
+
+ deallog.push("3d");
+ test<3>();
+ deallog.pop();
+}
--- /dev/null
+
+DEAL:1d::coarse_ cells 2.3: 2.2 2.3
+DEAL:1d::coarse_ cells 4.2: 4.1 4.2 4.3
+DEAL:1d::coarse_ cells 5.1: 4.1 4.2
+DEAL:2d::coarse_ cells 2.1: 2.0 2.1 2.3 2.4
+DEAL:2d::coarse_ cells 3.2: 2.0 2.2
+DEAL:2d::coarse_ cells 3.7: 3.5 3.6 3.7 3.10 3.25
+DEAL:2d::coarse_ cells 3.12: 2.2 2.8
+DEAL:2d::coarse_ cells 3.20: 3.13 3.20 3.21 3.22 3.38
+DEAL:2d::coarse_ cells 3.27: 3.18 3.25 3.26 3.27 3.29
+DEAL:2d::coarse_ cells 3.35: 3.33 3.34 3.35
+DEAL:2d::coarse_ cells 3.40: 2.1 2.4
+DEAL:2d::coarse_ cells 3.45: 3.15 3.44 3.45 3.47 3.48
+DEAL:2d::coarse_ cells 3.50: 3.47 3.48 3.50 3.51
+DEAL:2d::coarse_ cells 4.3: 4.1 4.2 4.3 4.6 4.25
+DEAL:2d::coarse_ cells 4.10: 3.14 3.44
+DEAL:2d::coarse_ cells 4.15: 3.24 3.25 3.26
+DEAL:2d::coarse_ cells 4.20: 3.26 3.28 3.49
+DEAL:2d::coarse_ cells 4.25: 4.3 4.24 4.25 4.27 4.28
+DEAL:2d::coarse_ cells 4.30: 3.32 3.33 3.35
+DEAL:2d::coarse_ cells 5.3: 4.6 4.7 4.28
+DEAL:3d::coarse_ cells 2.11: 2.9 2.10 2.11 2.15 2.25
+DEAL:3d::coarse_ cells 2.45: 2.41 2.44 2.45 2.47
+DEAL:3d::coarse_ cells 2.60: 2.46 2.53 2.56 2.60 2.61 2.62
+DEAL:3d::coarse_ cells 3.2: 3.0 3.2 3.3 3.6 3.104
+DEAL:3d::coarse_ cells 3.10: 3.8 3.10 3.11 3.14 3.102 3.123 3.136
+DEAL:3d::coarse_ cells 3.16: 3.16 3.17 3.18 3.20 3.145
+DEAL:3d::coarse_ cells 3.22: 2.8 2.9 2.11 2.13
+DEAL:3d::coarse_ cells 3.27: 3.23 3.25 3.26 3.27 3.31 3.177
+DEAL:3d::coarse_ cells 3.37: 3.33 3.36 3.37 3.39 3.159 3.228 3.241
+DEAL:3d::coarse_ cells 3.42: 3.40 3.42 3.43 3.46 3.195
+DEAL:3d::coarse_ cells 3.49: 3.48 3.49 3.51 3.53 3.179 3.229
+DEAL:3d::coarse_ cells 3.55: 3.51 3.53 3.54 3.55 3.261 3.379
+DEAL:3d::coarse_ cells 3.61: 2.33 2.35 2.39 2.42
+DEAL:3d::coarse_ cells 3.68: 3.64 3.68 3.69 3.70 3.293
+DEAL:3d::coarse_ cells 3.73: 3.72 3.73 3.75 3.77 3.285 3.291 3.296
+DEAL:3d::coarse_ cells 3.79: 3.75 3.77 3.78 3.79 3.302 3.349
+DEAL:3d::coarse_ cells 3.84: 3.62 3.80 3.84 3.85 3.86 3.341 3.352
+DEAL:3d::coarse_ cells 3.89: 2.22 2.48 2.50 2.51
+DEAL:3d::coarse_ cells 3.95: 2.50 2.51 2.54
+DEAL:3d::coarse_ cells 3.100: 3.5 3.8 3.96 3.100 3.101 3.102
+DEAL:3d::coarse_ cells 3.105: 3.3 3.104 3.105 3.107 3.109 3.112
+DEAL:3d::coarse_ cells 3.110: 2.2 2.6 2.16
+DEAL:3d::coarse_ cells 3.115: 3.113 3.114 3.115 3.119 3.154 3.185
+DEAL:3d::coarse_ cells 3.120: 3.4 3.120 3.121 3.122 3.124
+DEAL:3d::coarse_ cells 3.125: 3.12 3.121 3.124 3.125 3.127 3.265
+DEAL:3d::coarse_ cells 3.130: 3.110 3.128 3.130 3.131 3.134 3.200
+DEAL:3d::coarse_ cells 3.135: 3.131 3.133 3.134 3.135 3.142 3.205 3.283
+DEAL:3d::coarse_ cells 3.140: 3.14 3.56 3.133 3.136 3.140 3.141 3.142
+DEAL:3d::coarse_ cells 3.145: 3.16 3.144 3.145 3.147 3.149
+DEAL:3d::coarse_ cells 3.150: 3.103 3.146 3.148 3.150 3.151 3.156 3.162
+DEAL:3d::coarse_ cells 3.155: 2.10 2.11 2.24
+DEAL:3d::coarse_ cells 3.160: 3.9 3.148 3.160 3.161 3.162 3.164
+DEAL:3d::coarse_ cells 3.165: 3.28 3.161 3.164 3.165 3.167 3.305
+DEAL:3d::coarse_ cells 3.170: 3.139 3.158 3.168 3.170 3.171 3.174 3.240
+DEAL:3d::coarse_ cells 3.175: 2.14 2.15 2.28 2.42
+DEAL:3d::coarse_ cells 3.180: 3.30 3.173 3.176 3.180 3.181 3.182 3.320
+DEAL:3d::coarse_ cells 3.185: 3.32 3.115 3.184 3.185 3.187 3.189
+DEAL:3d::coarse_ cells 3.190: 2.16 2.17 2.19 2.21
+DEAL:3d::coarse_ cells 3.195: 3.42 3.193 3.194 3.195 3.199
+DEAL:3d::coarse_ cells 3.200: 2.6 2.16 2.20
+DEAL:3d::coarse_ cells 3.205: 3.135 3.201 3.204 3.205 3.207 3.212 3.337
+DEAL:3d::coarse_ cells 3.210: 3.190 3.203 3.208 3.210 3.211 3.214 3.216
+DEAL:3d::coarse_ cells 3.215: 3.83 3.211 3.213 3.214 3.215 3.221 3.246
+DEAL:3d::coarse_ cells 3.220: 2.21 2.22 2.23 2.51
+DEAL:3d::coarse_ cells 3.225: 2.11 2.25
+DEAL:3d::coarse_ cells 3.230: 3.39 3.50 3.226 3.228 3.230 3.231 3.236
+DEAL:3d::coarse_ cells 3.235: 3.233 3.234 3.235 3.239
+DEAL:3d::coarse_ cells 3.240: 3.36 3.170 3.209 3.240 3.241 3.242 3.244
+DEAL:3d::coarse_ cells 3.245: 3.52 3.175 3.241 3.244 3.245 3.247 3.369
+DEAL:3d::coarse_ cells 3.250: 3.46 3.219 3.248 3.250 3.251 3.254
+DEAL:3d::coarse_ cells 3.255: 3.251 3.253 3.254 3.255 3.262 3.387
+DEAL:3d::coarse_ cells 3.260: 3.54 3.253 3.256 3.260 3.261 3.262 3.392
+DEAL:3d::coarse_ cells 3.265: 3.125 3.264 3.265 3.267 3.269 3.272
+DEAL:3d::coarse_ cells 3.270: 3.266 3.268 3.270 3.271 3.284 3.290
+DEAL:3d::coarse_ cells 3.275: 3.15 3.57 3.273 3.274 3.275 3.279 3.306
+DEAL:3d::coarse_ cells 3.280: 3.132 3.266 3.280 3.281 3.282 3.284
+DEAL:3d::coarse_ cells 3.285: 3.60 3.73 3.271 3.281 3.284 3.285 3.287
+DEAL:3d::coarse_ cells 3.290: 3.72 3.270 3.288 3.290 3.291 3.294
+DEAL:3d::coarse_ cells 3.295: 3.70 3.77 3.291 3.293 3.294 3.295
+DEAL:3d::coarse_ cells 3.300: 3.70 3.77 3.296 3.300 3.301 3.302
+DEAL:3d::coarse_ cells 3.305: 3.165 3.304 3.305 3.307 3.309 3.312
+DEAL:3d::coarse_ cells 3.310: 2.33 2.40 2.42 2.44
+DEAL:3d::coarse_ cells 3.315: 3.31 3.313 3.314 3.315 3.319 3.321
+DEAL:3d::coarse_ cells 3.320: 2.15 2.41 2.42 2.43
+DEAL:3d::coarse_ cells 3.325: 2.41 2.43 2.47
+DEAL:3d::coarse_ cells 3.330: 2.37 2.40 2.44 2.46
+DEAL:3d::coarse_ cells 3.335: 2.44 2.45 2.46
+DEAL:3d::coarse_ cells 3.340: 3.286 3.336 3.340 3.341 3.342 3.344
+DEAL:3d::coarse_ cells 3.345: 3.75 3.341 3.344 3.345 3.347 3.349 3.352
+DEAL:3d::coarse_ cells 3.350: 3.346 3.348 3.350 3.351 3.364
+DEAL:3d::coarse_ cells 3.355: 2.49 2.53 2.55 2.60
+DEAL:3d::coarse_ cells 3.360: 3.92 3.346 3.360 3.361 3.362 3.364
+DEAL:3d::coarse_ cells 3.365: 2.52 2.54 2.55
+DEAL:3d::coarse_ cells 3.370: 3.83 3.246 3.368 3.370 3.371 3.374 3.384
+DEAL:3d::coarse_ cells 3.375: 2.56 2.57 2.58 2.60
+DEAL:3d::coarse_ cells 3.380: 2.43 2.56 2.57 2.61
+DEAL:3d::coarse_ cells 3.385: 3.253 3.371 3.384 3.385 3.387 3.389 3.392
+DEAL:3d::coarse_ cells 3.390: 2.51 2.58 2.62
+DEAL:3d::coarse_ cells 3.395: 3.263 3.393 3.394 3.395 3.399
+DEAL:3d::coarse_ cells 3.400: 2.55 2.58 2.60 2.62
+DEAL:3d::coarse_ cells 3.405: 2.60 2.62 2.63
+DEAL:3d::coarse_ cells 4.2: 3.0 3.4 3.6
+DEAL:3d::coarse_ cells 4.7: 3.4 3.5 3.6 3.120
+DEAL:3d::coarse_ cells 4.12: 3.4 3.6 3.122
+DEAL:3d::coarse_ cells 4.17: 3.8 3.9 3.100
+DEAL:3d::coarse_ cells 4.22: 3.8 3.10 3.12 3.121
+DEAL:3d::coarse_ cells 4.27: 3.9 3.13 3.15 3.164
+DEAL:3d::coarse_ cells 4.32: 3.16 3.20 3.149
+DEAL:3d::coarse_ cells 4.37: 3.20 3.21 3.24
+DEAL:3d::coarse_ cells 4.42: 3.25 3.28 3.29 3.31
+DEAL:3d::coarse_ cells 4.47: 3.29 3.31 3.313
+DEAL:3d::coarse_ cells 4.52: 3.29 3.30 3.31 3.315
+DEAL:3d::coarse_ cells 4.57: 3.32 3.33 3.154
+DEAL:3d::coarse_ cells 4.62: 3.32 3.34 3.36 3.185
+DEAL:3d::coarse_ cells 4.67: 3.35 3.41 3.226
+DEAL:3d::coarse_ cells 4.72: 3.32 3.36 3.158 3.189
+DEAL:3d::coarse_ cells 4.77: 3.36 3.37 3.158 3.240
+DEAL:3d::coarse_ cells 4.82: 3.42 3.43
+DEAL:3d::coarse_ cells 4.87: 3.43 3.47 3.234
+DEAL:3d::coarse_ cells 4.92: 3.44 3.46 3.199 3.250
+DEAL:3d::coarse_ cells 4.97: 3.50 3.52 3.54 3.55
+DEAL:3d::coarse_ cells 4.102: 3.54 3.247 3.260 3.378
+DEAL:3d::coarse_ cells 4.107: 3.58 3.59 3.80 3.142
+DEAL:3d::coarse_ cells 4.112: 3.58 3.60 3.62 3.287
+DEAL:3d::coarse_ cells 4.117: 3.60 3.62 3.63 3.298
+DEAL:3d::coarse_ cells 4.122: 3.64 3.66 3.276 3.289
+DEAL:3d::coarse_ cells 4.127: 3.64 3.65 3.66 3.68
+DEAL:3d::coarse_ cells 4.132: 3.73 3.74 3.75 3.79
+DEAL:3d::coarse_ cells 4.137: 3.90 3.92 3.94 3.95
+DEAL:3d::coarse_ cells 4.142: 3.94 3.362
--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2015 - 2016 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// Test GridTools::get_cells_at_coarsest_common_level()
+
+
+#include "../tests.h"
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/tensor.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_tools.h>
+
+#include <fstream>
+
+
+
+
+template<int dim>
+void test()
+{
+ Triangulation<dim> triangulation (Triangulation<dim>::limit_level_difference_at_vertices);
+
+ GridGenerator::hyper_cube(triangulation);
+ triangulation.refine_global (2);
+
+
+ // now extract patches and print every fifth of them
+ unsigned int index = 0;
+ for (typename Triangulation<dim>::active_cell_iterator
+ cell = triangulation.begin_active();
+ cell != triangulation.end(); ++cell, ++index)
+ {
+ std::vector<typename Triangulation<dim>::active_cell_iterator> patch_cells
+ = GridTools::get_patch_around_cell<Triangulation<dim> > (cell);
+
+ std::vector<typename Triangulation<dim>::cell_iterator> coarse_cells
+ = GridTools::get_cells_at_coarsest_common_level <Triangulation<dim> > (patch_cells);
+
+ // in uniform refinemenet(without any hanging nodes), we expect the size of vector of patch_cells for each cell be
+ // equal to the vector of coarse_cells around that cell
+ Assert(patch_cells.size()==coarse_cells.size(), ExcInternalError());
+
+ // we sort both vectors to ensure that we are comparing the correct cells to each other
+ std::sort(patch_cells.begin(), patch_cells.end());
+ std::sort(coarse_cells.begin(), coarse_cells.end());
+
+ for (unsigned int i=0; i<patch_cells.size(); ++i)
+ Assert (patch_cells[i]==coarse_cells[i] ,ExcInternalError());
+
+ deallog << "coarse_ cells " << cell << ": ";
+ for (unsigned int i=0; i<coarse_cells.size(); ++i)
+ deallog << coarse_cells[i] << ' ';
+ deallog << std::endl;
+
+ }
+
+ deallog<< "OKay"<< std::endl;
+
+}
+
+
+int main()
+{
+ initlog();
+ deallog.threshold_double(1.e-10);
+
+ deallog.push("1d");
+ test<1>();
+ deallog.pop();
+
+ deallog.push("2d");
+ test<2>();
+ deallog.pop();
+
+ deallog.push("3d");
+ test<3>();
+ deallog.pop();
+}
--- /dev/null
+
+DEAL:1d::coarse_ cells 2.0: 2.0 2.1
+DEAL:1d::coarse_ cells 2.1: 2.0 2.1 2.2
+DEAL:1d::coarse_ cells 2.2: 2.1 2.2 2.3
+DEAL:1d::coarse_ cells 2.3: 2.2 2.3
+DEAL:1d::OKay
+DEAL:2d::coarse_ cells 2.0: 2.0 2.1 2.2
+DEAL:2d::coarse_ cells 2.1: 2.0 2.1 2.3 2.4
+DEAL:2d::coarse_ cells 2.2: 2.0 2.2 2.3 2.8
+DEAL:2d::coarse_ cells 2.3: 2.1 2.2 2.3 2.6 2.9
+DEAL:2d::coarse_ cells 2.4: 2.1 2.4 2.5 2.6
+DEAL:2d::coarse_ cells 2.5: 2.4 2.5 2.7
+DEAL:2d::coarse_ cells 2.6: 2.3 2.4 2.6 2.7 2.12
+DEAL:2d::coarse_ cells 2.7: 2.5 2.6 2.7 2.13
+DEAL:2d::coarse_ cells 2.8: 2.2 2.8 2.9 2.10
+DEAL:2d::coarse_ cells 2.9: 2.3 2.8 2.9 2.11 2.12
+DEAL:2d::coarse_ cells 2.10: 2.8 2.10 2.11
+DEAL:2d::coarse_ cells 2.11: 2.9 2.10 2.11 2.14
+DEAL:2d::coarse_ cells 2.12: 2.6 2.9 2.12 2.13 2.14
+DEAL:2d::coarse_ cells 2.13: 2.7 2.12 2.13 2.15
+DEAL:2d::coarse_ cells 2.14: 2.11 2.12 2.14 2.15
+DEAL:2d::coarse_ cells 2.15: 2.13 2.14 2.15
+DEAL:2d::OKay
+DEAL:3d::coarse_ cells 2.0: 2.0 2.1 2.2 2.4
+DEAL:3d::coarse_ cells 2.1: 2.0 2.1 2.3 2.5 2.8
+DEAL:3d::coarse_ cells 2.2: 2.0 2.2 2.3 2.6 2.16
+DEAL:3d::coarse_ cells 2.3: 2.1 2.2 2.3 2.7 2.10 2.17
+DEAL:3d::coarse_ cells 2.4: 2.0 2.4 2.5 2.6 2.32
+DEAL:3d::coarse_ cells 2.5: 2.1 2.4 2.5 2.7 2.12 2.33
+DEAL:3d::coarse_ cells 2.6: 2.2 2.4 2.6 2.7 2.20 2.34
+DEAL:3d::coarse_ cells 2.7: 2.3 2.5 2.6 2.7 2.14 2.21 2.35
+DEAL:3d::coarse_ cells 2.8: 2.1 2.8 2.9 2.10 2.12
+DEAL:3d::coarse_ cells 2.9: 2.8 2.9 2.11 2.13
+DEAL:3d::coarse_ cells 2.10: 2.3 2.8 2.10 2.11 2.14 2.24
+DEAL:3d::coarse_ cells 2.11: 2.9 2.10 2.11 2.15 2.25
+DEAL:3d::coarse_ cells 2.12: 2.5 2.8 2.12 2.13 2.14 2.40
+DEAL:3d::coarse_ cells 2.13: 2.9 2.12 2.13 2.15 2.41
+DEAL:3d::coarse_ cells 2.14: 2.7 2.10 2.12 2.14 2.15 2.28 2.42
+DEAL:3d::coarse_ cells 2.15: 2.11 2.13 2.14 2.15 2.29 2.43
+DEAL:3d::coarse_ cells 2.16: 2.2 2.16 2.17 2.18 2.20
+DEAL:3d::coarse_ cells 2.17: 2.3 2.16 2.17 2.19 2.21 2.24
+DEAL:3d::coarse_ cells 2.18: 2.16 2.18 2.19 2.22
+DEAL:3d::coarse_ cells 2.19: 2.17 2.18 2.19 2.23 2.26
+DEAL:3d::coarse_ cells 2.20: 2.6 2.16 2.20 2.21 2.22 2.48
+DEAL:3d::coarse_ cells 2.21: 2.7 2.17 2.20 2.21 2.23 2.28 2.49
+DEAL:3d::coarse_ cells 2.22: 2.18 2.20 2.22 2.23 2.50
+DEAL:3d::coarse_ cells 2.23: 2.19 2.21 2.22 2.23 2.30 2.51
+DEAL:3d::coarse_ cells 2.24: 2.10 2.17 2.24 2.25 2.26 2.28
+DEAL:3d::coarse_ cells 2.25: 2.11 2.24 2.25 2.27 2.29
+DEAL:3d::coarse_ cells 2.26: 2.19 2.24 2.26 2.27 2.30
+DEAL:3d::coarse_ cells 2.27: 2.25 2.26 2.27 2.31
+DEAL:3d::coarse_ cells 2.28: 2.14 2.21 2.24 2.28 2.29 2.30 2.56
+DEAL:3d::coarse_ cells 2.29: 2.15 2.25 2.28 2.29 2.31 2.57
+DEAL:3d::coarse_ cells 2.30: 2.23 2.26 2.28 2.30 2.31 2.58
+DEAL:3d::coarse_ cells 2.31: 2.27 2.29 2.30 2.31 2.59
+DEAL:3d::coarse_ cells 2.32: 2.4 2.32 2.33 2.34 2.36
+DEAL:3d::coarse_ cells 2.33: 2.5 2.32 2.33 2.35 2.37 2.40
+DEAL:3d::coarse_ cells 2.34: 2.6 2.32 2.34 2.35 2.38 2.48
+DEAL:3d::coarse_ cells 2.35: 2.7 2.33 2.34 2.35 2.39 2.42 2.49
+DEAL:3d::coarse_ cells 2.36: 2.32 2.36 2.37 2.38
+DEAL:3d::coarse_ cells 2.37: 2.33 2.36 2.37 2.39 2.44
+DEAL:3d::coarse_ cells 2.38: 2.34 2.36 2.38 2.39 2.52
+DEAL:3d::coarse_ cells 2.39: 2.35 2.37 2.38 2.39 2.46 2.53
+DEAL:3d::coarse_ cells 2.40: 2.12 2.33 2.40 2.41 2.42 2.44
+DEAL:3d::coarse_ cells 2.41: 2.13 2.40 2.41 2.43 2.45
+DEAL:3d::coarse_ cells 2.42: 2.14 2.35 2.40 2.42 2.43 2.46 2.56
+DEAL:3d::coarse_ cells 2.43: 2.15 2.41 2.42 2.43 2.47 2.57
+DEAL:3d::coarse_ cells 2.44: 2.37 2.40 2.44 2.45 2.46
+DEAL:3d::coarse_ cells 2.45: 2.41 2.44 2.45 2.47
+DEAL:3d::coarse_ cells 2.46: 2.39 2.42 2.44 2.46 2.47 2.60
+DEAL:3d::coarse_ cells 2.47: 2.43 2.45 2.46 2.47 2.61
+DEAL:3d::coarse_ cells 2.48: 2.20 2.34 2.48 2.49 2.50 2.52
+DEAL:3d::coarse_ cells 2.49: 2.21 2.35 2.48 2.49 2.51 2.53 2.56
+DEAL:3d::coarse_ cells 2.50: 2.22 2.48 2.50 2.51 2.54
+DEAL:3d::coarse_ cells 2.51: 2.23 2.49 2.50 2.51 2.55 2.58
+DEAL:3d::coarse_ cells 2.52: 2.38 2.48 2.52 2.53 2.54
+DEAL:3d::coarse_ cells 2.53: 2.39 2.49 2.52 2.53 2.55 2.60
+DEAL:3d::coarse_ cells 2.54: 2.50 2.52 2.54 2.55
+DEAL:3d::coarse_ cells 2.55: 2.51 2.53 2.54 2.55 2.62
+DEAL:3d::coarse_ cells 2.56: 2.28 2.42 2.49 2.56 2.57 2.58 2.60
+DEAL:3d::coarse_ cells 2.57: 2.29 2.43 2.56 2.57 2.59 2.61
+DEAL:3d::coarse_ cells 2.58: 2.30 2.51 2.56 2.58 2.59 2.62
+DEAL:3d::coarse_ cells 2.59: 2.31 2.57 2.58 2.59 2.63
+DEAL:3d::coarse_ cells 2.60: 2.46 2.53 2.56 2.60 2.61 2.62
+DEAL:3d::coarse_ cells 2.61: 2.47 2.57 2.60 2.61 2.63
+DEAL:3d::coarse_ cells 2.62: 2.55 2.58 2.60 2.62 2.63
+DEAL:3d::coarse_ cells 2.63: 2.59 2.61 2.62 2.63
+DEAL:3d::OKay
--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2015 - 2016 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+
+// Test GridTools::get_cells_at_coarsest_common_level()
+
+
+
+#include "../tests.h"
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/tensor.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_tools.h>
+
+#include <fstream>
+
+
+
+
+template<int dim>
+void test()
+{
+ Triangulation<dim> triangulation (Triangulation<dim>::limit_level_difference_at_vertices);
+
+ GridGenerator::hyper_cube(triangulation);
+ triangulation.refine_global (1);
+ // Adaptive refinement ... refine the first cell
+ triangulation.begin_active()->set_refine_flag();
+ triangulation.execute_coarsening_and_refinement();
+
+ // now extract patches and print every fifth of them
+ unsigned int index = 0;
+ for (typename Triangulation<dim>::active_cell_iterator
+ cell = triangulation.begin_active();
+ cell != triangulation.end(); ++cell, ++index)
+ {
+ std::vector<typename Triangulation<dim>::active_cell_iterator> patch_cells
+ = GridTools::get_patch_around_cell<Triangulation<dim> > (cell);
+
+ std::vector<typename Triangulation<dim>::cell_iterator> coarse_cells
+ = GridTools::get_cells_at_coarsest_common_level <Triangulation<dim> > (patch_cells);
+
+ deallog << "coarse_ cells " << cell << ": ";
+ for (unsigned int i=0; i<coarse_cells.size(); ++i)
+ deallog << coarse_cells[i] << ' ';
+ deallog << std::endl;
+
+ }
+
+
+}
+
+
+int main()
+{
+ initlog();
+ deallog.threshold_double(1.e-10);
+
+ deallog.push("1d");
+ test<1>();
+ deallog.pop();
+
+ deallog.push("2d");
+ test<2>();
+ deallog.pop();
+
+ deallog.push("3d");
+ test<3>();
+ deallog.pop();
+}
--- /dev/null
+
+DEAL:1d::coarse_ cells 1.1: 1.0 1.1
+DEAL:1d::coarse_ cells 2.0: 2.0 2.1
+DEAL:1d::coarse_ cells 2.1: 1.0 1.1
+DEAL:2d::coarse_ cells 1.1: 1.0 1.1 1.3
+DEAL:2d::coarse_ cells 1.2: 1.0 1.2 1.3
+DEAL:2d::coarse_ cells 1.3: 1.1 1.2 1.3
+DEAL:2d::coarse_ cells 2.0: 2.0 2.1 2.2
+DEAL:2d::coarse_ cells 2.1: 1.0 1.1
+DEAL:2d::coarse_ cells 2.2: 1.0 1.2
+DEAL:2d::coarse_ cells 2.3: 1.0 1.1 1.2
+DEAL:3d::coarse_ cells 1.1: 1.0 1.1 1.3 1.5
+DEAL:3d::coarse_ cells 1.2: 1.0 1.2 1.3 1.6
+DEAL:3d::coarse_ cells 1.3: 1.1 1.2 1.3 1.7
+DEAL:3d::coarse_ cells 1.4: 1.0 1.4 1.5 1.6
+DEAL:3d::coarse_ cells 1.5: 1.1 1.4 1.5 1.7
+DEAL:3d::coarse_ cells 1.6: 1.2 1.4 1.6 1.7
+DEAL:3d::coarse_ cells 1.7: 1.3 1.5 1.6 1.7
+DEAL:3d::coarse_ cells 2.0: 2.0 2.1 2.2 2.4
+DEAL:3d::coarse_ cells 2.1: 1.0 1.1
+DEAL:3d::coarse_ cells 2.2: 1.0 1.2
+DEAL:3d::coarse_ cells 2.3: 1.0 1.1 1.2
+DEAL:3d::coarse_ cells 2.4: 1.0 1.4
+DEAL:3d::coarse_ cells 2.5: 1.0 1.1 1.4
+DEAL:3d::coarse_ cells 2.6: 1.0 1.2 1.4
+DEAL:3d::coarse_ cells 2.7: 1.0 1.1 1.2 1.4
projects / dealii.git / commitdiff