]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Adapt a few existing tests to test the new active_cell_index() function.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Fri, 10 Apr 2015 22:26:22 +0000 (17:26 -0500)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Mon, 13 Apr 2015 00:10:18 +0000 (19:10 -0500)
tests/bits/refine_and_coarsen_for_active_cell_index.cc [new file with mode: 0644]
tests/bits/refine_and_coarsen_for_active_cell_index.output [new file with mode: 0644]
tests/bits/refine_and_coarsen_for_active_cell_index_02.cc [new file with mode: 0644]
tests/bits/refine_and_coarsen_for_active_cell_index_02.output [new file with mode: 0644]
tests/bits/refine_and_coarsen_for_active_cell_index_03.cc [new file with mode: 0644]
tests/bits/refine_and_coarsen_for_active_cell_index_03.output [new file with mode: 0644]

diff --git a/tests/bits/refine_and_coarsen_for_active_cell_index.cc b/tests/bits/refine_and_coarsen_for_active_cell_index.cc
new file mode 100644 (file)
index 0000000..0bb066d
--- /dev/null
@@ -0,0 +1,128 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2010 - 2015 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// check that cell->active_cell_index() works as advertised
+
+#include "../tests.h"
+
+#include <deal.II/base/geometry_info.h>
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/quadrature_lib.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 <fstream>
+
+
+template <int dim>
+void check (const Triangulation<dim> &tria)
+{
+  unsigned int index = 0;
+  for (typename Triangulation<dim>::active_cell_iterator cell=tria.begin_active();
+      cell!=tria.end(); ++cell, ++index)
+    Assert (cell->active_cell_index() == index, ExcInternalError());
+}
+
+void do_refine (Triangulation<1> &tria)
+{
+  const int dim = 1;
+
+  tria.refine_global (2);
+  tria.begin_active()->set_refine_flag();
+  tria.execute_coarsening_and_refinement ();
+}
+
+
+void do_refine (Triangulation<2> &tria)
+{
+  const int dim = 2;
+
+  tria.refine_global (2);
+  tria.begin_active()->set_refine_flag();
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active ()->set_refine_flag (RefinementPossibilities<dim>::cut_x);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active ()->set_refine_flag (RefinementPossibilities<dim>::cut_y);
+  tria.execute_coarsening_and_refinement ();
+}
+
+
+void do_refine (Triangulation<3> &tria)
+{
+  const int dim = 3;
+
+  tria.refine_global (2);
+  tria.begin_active()->set_refine_flag();
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_x);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_y);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_z);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_xy);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_xz);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_yz);
+  tria.execute_coarsening_and_refinement ();
+}
+
+
+template <int dim>
+void check ()
+{
+  Triangulation<dim> tria;
+  GridGenerator::hyper_cube (tria);
+  do_refine (tria);
+  // refine the mesh globally and
+  // verify that the parent relation
+  // holds
+  tria.refine_global (1);
+
+  check (tria);
+
+  // coarsen the mesh globally and
+  // verify again
+  for (typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active ();
+       cell != tria.end (); ++cell)
+    cell->set_coarsen_flag ();
+
+  tria.execute_coarsening_and_refinement ();
+
+  check (tria);
+
+  deallog << "OK for " << dim << "d" << std::endl;
+}
+
+
+int main ()
+{
+  std::ofstream logfile("output");
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+  deallog.threshold_double(1.e-10);
+
+  check<1> ();
+  check<2> ();
+  check<3> ();
+}
+
+
+
diff --git a/tests/bits/refine_and_coarsen_for_active_cell_index.output b/tests/bits/refine_and_coarsen_for_active_cell_index.output
new file mode 100644 (file)
index 0000000..292df5f
--- /dev/null
@@ -0,0 +1,4 @@
+
+DEAL::OK for 1d
+DEAL::OK for 2d
+DEAL::OK for 3d
diff --git a/tests/bits/refine_and_coarsen_for_active_cell_index_02.cc b/tests/bits/refine_and_coarsen_for_active_cell_index_02.cc
new file mode 100644 (file)
index 0000000..91be163
--- /dev/null
@@ -0,0 +1,132 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2010 - 2015 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// check that cell->active_cell_index() works as advertised
+
+#include "../tests.h"
+
+#include <deal.II/base/geometry_info.h>
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/quadrature_lib.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/dofs/dof_accessor.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/grid_generator.h>
+
+#include <fstream>
+
+
+
+template <int dim>
+void check (const Triangulation<dim> &tria)
+{
+  unsigned int index = 0;
+  for (typename Triangulation<dim>::active_cell_iterator cell=tria.begin_active();
+      cell!=tria.end(); ++cell, ++index)
+    Assert (cell->active_cell_index() == index, ExcInternalError());
+}
+
+
+
+void do_refine (Triangulation<1> &tria)
+{
+  tria.refine_global (2);
+  tria.begin_active()->set_refine_flag();
+  tria.execute_coarsening_and_refinement ();
+}
+
+
+void do_refine (Triangulation<2> &tria)
+{
+  const int dim = 2;
+
+  tria.refine_global (2);
+  tria.begin_active()->set_refine_flag();
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active ()->set_refine_flag (RefinementPossibilities<dim>::cut_x);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active ()->set_refine_flag (RefinementPossibilities<dim>::cut_y);
+  tria.execute_coarsening_and_refinement ();
+}
+
+
+void do_refine (Triangulation<3> &tria)
+{
+  const int dim = 3;
+
+  tria.refine_global (2);
+  tria.begin_active()->set_refine_flag();
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_x);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_y);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_z);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_xy);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_xz);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_yz);
+  tria.execute_coarsening_and_refinement ();
+}
+
+
+template <int dim>
+void check ()
+{
+  Triangulation<dim> tria;
+  GridGenerator::hyper_cube (tria);
+  do_refine (tria);
+  // refine the mesh globally and
+  // verify that the parent relation
+  // holds
+  tria.refine_global (1);
+
+  DoFHandler<dim> dof_handler (tria);
+
+  check (tria);
+
+  // coarsen the mesh globally and
+  // verify again
+  for (typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active ();
+       cell != tria.end (); ++cell)
+    cell->set_coarsen_flag ();
+
+  tria.execute_coarsening_and_refinement ();
+
+  check (tria);
+
+  deallog << "OK for " << dim << "d" << std::endl;
+}
+
+
+int main ()
+{
+  std::ofstream logfile("output");
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+  deallog.threshold_double(1.e-10);
+
+  check<1> ();
+  check<2> ();
+  check<3> ();
+}
+
+
+
diff --git a/tests/bits/refine_and_coarsen_for_active_cell_index_02.output b/tests/bits/refine_and_coarsen_for_active_cell_index_02.output
new file mode 100644 (file)
index 0000000..292df5f
--- /dev/null
@@ -0,0 +1,4 @@
+
+DEAL::OK for 1d
+DEAL::OK for 2d
+DEAL::OK for 3d
diff --git a/tests/bits/refine_and_coarsen_for_active_cell_index_03.cc b/tests/bits/refine_and_coarsen_for_active_cell_index_03.cc
new file mode 100644 (file)
index 0000000..28779bc
--- /dev/null
@@ -0,0 +1,131 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2010 - 2015 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.
+//
+// ---------------------------------------------------------------------
+
+
+
+// check that cell->active_cell_index() works as advertised
+
+#include "../tests.h"
+
+#include <deal.II/base/geometry_info.h>
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/quadrature_lib.h>
+#include <deal.II/dofs/dof_accessor.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/grid_generator.h>
+
+#include <fstream>
+
+
+
+template <int dim>
+void check (const Triangulation<dim> &tria)
+{
+  unsigned int index = 0;
+  for (typename Triangulation<dim>::active_cell_iterator cell=tria.begin_active();
+      cell!=tria.end(); ++cell, ++index)
+    Assert (cell->active_cell_index() == index, ExcInternalError());
+}
+
+
+
+void do_refine (Triangulation<1> &tria)
+{
+  tria.refine_global (2);
+  tria.begin_active()->set_refine_flag();
+  tria.execute_coarsening_and_refinement ();
+}
+
+
+void do_refine (Triangulation<2> &tria)
+{
+  const int dim = 2;
+
+  tria.refine_global (2);
+  tria.begin_active()->set_refine_flag();
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active ()->set_refine_flag (RefinementPossibilities<dim>::cut_x);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active ()->set_refine_flag (RefinementPossibilities<dim>::cut_y);
+  tria.execute_coarsening_and_refinement ();
+}
+
+
+void do_refine (Triangulation<3> &tria)
+{
+  const int dim = 3;
+
+  tria.refine_global (2);
+  tria.begin_active()->set_refine_flag();
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_x);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_y);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_z);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_xy);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_xz);
+  tria.execute_coarsening_and_refinement ();
+  tria.begin_active()->set_refine_flag(RefinementPossibilities<dim>::cut_yz);
+  tria.execute_coarsening_and_refinement ();
+}
+
+
+template <int dim>
+void check ()
+{
+  Triangulation<dim> tria;
+  GridGenerator::hyper_cube (tria);
+  do_refine (tria);
+  // refine the mesh globally and
+  // verify that the parent relation
+  // holds
+  tria.refine_global (1);
+
+  DoFHandler<dim> dof_handler (tria);
+
+  check (tria);
+
+  // coarsen the mesh globally and
+  // verify again
+  for (typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active ();
+       cell != tria.end (); ++cell)
+    cell->set_coarsen_flag ();
+
+  tria.execute_coarsening_and_refinement ();
+
+  check (tria);
+
+  deallog << "OK for " << dim << "d" << std::endl;
+}
+
+
+int main ()
+{
+  std::ofstream logfile("output");
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+  deallog.threshold_double(1.e-10);
+
+  check<1> ();
+  check<2> ();
+  check<3> ();
+}
+
+
+
diff --git a/tests/bits/refine_and_coarsen_for_active_cell_index_03.output b/tests/bits/refine_and_coarsen_for_active_cell_index_03.output
new file mode 100644 (file)
index 0000000..292df5f
--- /dev/null
@@ -0,0 +1,4 @@
+
+DEAL::OK for 1d
+DEAL::OK for 2d
+DEAL::OK for 3d

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.