std::ofstream logfile("number_cache/output");
+
+template <int dim>
+void output (const Triangulation<dim> &tria)
+{
+
+ deallog << " " << tria.n_active_cells() << std::endl;
+ deallog << " " << tria.n_cells() << std::endl;
+ deallog << " " << tria.n_active_lines() << std::endl;
+ deallog << " " << tria.n_lines() << std::endl;
+ deallog << " " << tria.n_active_quads() << std::endl;
+ deallog << " " << tria.n_quads() << std::endl;
+ deallog << " " << tria.n_active_hexs() << std::endl;
+ deallog << " " << tria.n_hexs() << std::endl;
+
+ for (unsigned int i=0; i<tria.n_levels(); ++i)
+ {
+ deallog << " " << tria.n_active_cells(i) << std::endl;
+ deallog << " " << tria.n_cells(i) << std::endl;
+ if (dim == 1)
+ {
+ deallog << " " << tria.n_active_lines(i) << std::endl;
+ deallog << " " << tria.n_lines(i) << std::endl;
+ }
+
+ if (dim == 2)
+ {
+ deallog << " " << tria.n_active_quads(i) << std::endl;
+ deallog << " " << tria.n_quads(i) << std::endl;
+ }
+
+ if (dim == 3)
+ {
+ deallog << " " << tria.n_active_hexs(i) << std::endl;
+ deallog << " " << tria.n_hexs(i) << std::endl;
+ }
+ }
+}
+
+
template <int dim>
void test (const char *filename)
{
return;
}
- deallog << " " << tria.n_active_cells() << std::endl;
- deallog << " " << tria.n_cells() << std::endl;
- deallog << " " << tria.n_active_lines() << std::endl;
- deallog << " " << tria.n_lines() << std::endl;
- deallog << " " << tria.n_active_quads() << std::endl;
- deallog << " " << tria.n_quads() << std::endl;
- deallog << " " << tria.n_active_hexs() << std::endl;
- deallog << " " << tria.n_hexs() << std::endl;
+ output (tria);
// now refine a few cells and output again
deallog << " Refining..." << std::endl;
10U); ++i, ++cell)
cell->set_refine_flag ();
tria.execute_coarsening_and_refinement ();
-
- deallog << " " << tria.n_active_cells() << std::endl;
- deallog << " " << tria.n_cells() << std::endl;
- deallog << " " << tria.n_active_lines() << std::endl;
- deallog << " " << tria.n_lines() << std::endl;
- deallog << " " << tria.n_active_quads() << std::endl;
- deallog << " " << tria.n_quads() << std::endl;
- deallog << " " << tria.n_active_hexs() << std::endl;
- deallog << " " << tria.n_hexs() << std::endl;
+
+ output (tria);
}
DEAL:: 29037
DEAL:: 0
DEAL:: 0
+DEAL:: 29037
+DEAL:: 29037
+DEAL:: 29037
+DEAL:: 29037
DEAL:: Refining...
DEAL:: 29070
DEAL:: 29081
DEAL:: 29081
DEAL:: 0
DEAL:: 0
+DEAL:: 29026
+DEAL:: 29037
+DEAL:: 29026
+DEAL:: 29037
+DEAL:: 44
+DEAL:: 44
+DEAL:: 44
+DEAL:: 44
DEAL::Reading grid_in_3d/1.in
DEAL:: 16480
DEAL:: 16480
DEAL:: 52008
DEAL:: 16480
DEAL:: 16480
+DEAL:: 16480
+DEAL:: 16480
+DEAL:: 16480
+DEAL:: 16480
DEAL:: Refining...
DEAL:: 16557
DEAL:: 16568
DEAL:: 52364
DEAL:: 16557
DEAL:: 16568
+DEAL:: 16469
+DEAL:: 16480
+DEAL:: 16469
+DEAL:: 16480
+DEAL:: 88
+DEAL:: 88
+DEAL:: 88
+DEAL:: 88
DEAL::Reading grid_in_3d/2.in
DEAL:: 360
DEAL:: 360
DEAL:: 1218
DEAL:: 360
DEAL:: 360
+DEAL:: 360
+DEAL:: 360
+DEAL:: 360
+DEAL:: 360
DEAL:: Refining...
DEAL:: 437
DEAL:: 448
DEAL:: 1550
DEAL:: 437
DEAL:: 448
+DEAL:: 349
+DEAL:: 360
+DEAL:: 349
+DEAL:: 360
+DEAL:: 88
+DEAL:: 88
+DEAL:: 88
+DEAL:: 88
DEAL::Reading grid_in_3d/3.in
DEAL:: 2
DEAL:: 2
DEAL:: 11
DEAL:: 2
DEAL:: 2
+DEAL:: 2
+DEAL:: 2
+DEAL:: 2
+DEAL:: 2
DEAL:: Refining...
DEAL:: 9
DEAL:: 10
DEAL:: 47
DEAL:: 9
DEAL:: 10
+DEAL:: 1
+DEAL:: 2
+DEAL:: 1
+DEAL:: 2
+DEAL:: 8
+DEAL:: 8
+DEAL:: 8
+DEAL:: 8
DEAL::Reading grid_in_3d/4.in
DEAL:: 13824
DEAL:: 13824
DEAL:: 42816
DEAL:: 13824
DEAL:: 13824
+DEAL:: 13824
+DEAL:: 13824
+DEAL:: 13824
+DEAL:: 13824
DEAL:: Refining...
DEAL:: 13901
DEAL:: 13912
DEAL:: 43164
DEAL:: 13901
DEAL:: 13912
+DEAL:: 13813
+DEAL:: 13824
+DEAL:: 13813
+DEAL:: 13824
+DEAL:: 88
+DEAL:: 88
+DEAL:: 88
+DEAL:: 88
DEAL::Reading grid_in_3d/evil_0.in
DEAL:: 10
DEAL:: 10
DEAL:: 50
DEAL:: 10
DEAL:: 10
+DEAL:: 10
+DEAL:: 10
+DEAL:: 10
+DEAL:: 10
DEAL:: Refining...
DEAL:: 38
DEAL:: 42
DEAL:: 182
DEAL:: 38
DEAL:: 42
+DEAL:: 6
+DEAL:: 10
+DEAL:: 6
+DEAL:: 10
+DEAL:: 32
+DEAL:: 32
+DEAL:: 32
+DEAL:: 32
DEAL::Reading grid_in_3d/evil_1.in
DEAL:: 10
DEAL:: 10
DEAL:: 50
DEAL:: 10
DEAL:: 10
+DEAL:: 10
+DEAL:: 10
+DEAL:: 10
+DEAL:: 10
DEAL:: Refining...
DEAL:: 38
DEAL:: 42
DEAL:: 182
DEAL:: 38
DEAL:: 42
+DEAL:: 6
+DEAL:: 10
+DEAL:: 6
+DEAL:: 10
+DEAL:: 32
+DEAL:: 32
+DEAL:: 32
+DEAL:: 32
DEAL::Reading grid_in_3d/evil_2.in
DEAL:: 10
DEAL:: 10
DEAL:: 50
DEAL:: 10
DEAL:: 10
+DEAL:: 10
+DEAL:: 10
+DEAL:: 10
+DEAL:: 10
DEAL:: Refining...
DEAL:: 38
DEAL:: 42
DEAL:: 182
DEAL:: 38
DEAL:: 42
+DEAL:: 6
+DEAL:: 10
+DEAL:: 6
+DEAL:: 10
+DEAL:: 32
+DEAL:: 32
+DEAL:: 32
+DEAL:: 32
DEAL::Reading grid_in_3d/evil_3.in
DEAL:: 10
DEAL:: 10
DEAL:: 50
DEAL:: 10
DEAL:: 10
+DEAL:: 10
+DEAL:: 10
+DEAL:: 10
+DEAL:: 10
DEAL:: Refining...
DEAL:: 38
DEAL:: 42
DEAL:: 182
DEAL:: 38
DEAL:: 42
+DEAL:: 6
+DEAL:: 10
+DEAL:: 6
+DEAL:: 10
+DEAL:: 32
+DEAL:: 32
+DEAL:: 32
+DEAL:: 32
DEAL::Reading grid_in_3d/evil_4.in
DEAL:: 10
DEAL:: 10
DEAL:: 50
DEAL:: 10
DEAL:: 10
+DEAL:: 10
+DEAL:: 10
+DEAL:: 10
+DEAL:: 10
DEAL:: Refining...
DEAL:: 38
DEAL:: 42
DEAL:: 182
DEAL:: 38
DEAL:: 42
+DEAL:: 6
+DEAL:: 10
+DEAL:: 6
+DEAL:: 10
+DEAL:: 32
+DEAL:: 32
+DEAL:: 32
+DEAL:: 32