New test.

git-svn-id: https://svn.dealii.org/trunk@27088 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/tests/deal.II/number_cache_03.cc b/tests/deal.II/number_cache_03.cc
new file mode 100644 (file)
index 0000000..f2bf47b
--- /dev/null
+++ b/tests/deal.II/number_cache_03.cc
@@ -0,0 +1,113 @@
+//----------------------------  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> ();
+}
+

diff --git a/tests/deal.II/number_cache_03/cmp/generic b/tests/deal.II/number_cache_03/cmp/generic
new file mode 100644 (file)
index 0000000..e974b9c
--- /dev/null
+++ b/tests/deal.II/number_cache_03/cmp/generic
@@ -0,0 +1,76 @@
+
+DEAL::1d
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  0
+DEAL::  0
+DEAL::  0
+DEAL::  0
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  0
+DEAL::  0
+DEAL::  0
+DEAL::  0
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::2d
+DEAL::  1
+DEAL::  1
+DEAL::  4
+DEAL::  4
+DEAL::  1
+DEAL::  1
+DEAL::  0
+DEAL::  0
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  4
+DEAL::  4
+DEAL::  1
+DEAL::  1
+DEAL::  0
+DEAL::  0
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::3d
+DEAL::  1
+DEAL::  1
+DEAL::  12
+DEAL::  12
+DEAL::  6
+DEAL::  6
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  12
+DEAL::  12
+DEAL::  6
+DEAL::  6
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1
+DEAL::  1