--- /dev/null
+//---------------------------- number_cache_03.cc ---------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2008, 2012 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//---------------------------- number_cache_03.cc ---------------------------
+
+
+// test that we correctly compute the number of lines, quads, hexes,
+// etc and store them in the number cache of the Triangulation
+// class. compare a pristine mesh with one that has first been
+// globally refined and then globally coarsened
+//
+// this is sort of a pathological case because there now is a level in
+// the hierarchy but it is completely unused and only consists of raw
+// cells
+
+#include "../tests.h"
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_boundary.h>
+#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/grid_out.h>
+#include <deal.II/grid/grid_in.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/base/logstream.h>
+
+#include <fstream>
+#include <iomanip>
+
+std::ofstream logfile("number_cache_03/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 ()
+{
+ deallog << dim << "d" << std::endl;
+
+ Triangulation<dim> tria;
+ GridGenerator::hyper_cube (tria);
+
+ output (tria);
+
+ tria.refine_global(1);
+
+ for (typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active();
+ cell != tria.end(); ++cell)
+ cell->set_coarsen_flag ();
+ tria.execute_coarsening_and_refinement ();
+
+ output (tria);
+}
+
+
+int main ()
+{
+ deallog << std::setprecision (2);
+ deallog.attach(logfile);
+ deallog.depth_console(0);
+ deallog.threshold_double(1.e-10);
+
+ test<1> ();
+ test<2> ();
+ test<3> ();
+}
+
projects / dealii.git / commitdiff