Add tests.

diff --git a/tests/grid/refine_and_coarsen_fixed_fraction_02.cc b/tests/grid/refine_and_coarsen_fixed_fraction_02.cc
new file mode 100644 (file)
index 0000000..7a71aae
--- /dev/null
+++ b/tests/grid/refine_and_coarsen_fixed_fraction_02.cc
@@ -0,0 +1,74 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2015 - 2017 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.md at
+// the top level directory of deal.II.
+//
+// ---------------------------------------------------------------------
+
+// When setting refine_fraction=0.6 and coarsen_fraction=0.4, these
+// numbers may add up to slightly larger than a total of 1.0 in
+// floating point arithmetic. Make sure we allow this nevertheless.
+
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_refinement.h>
+#include <deal.II/grid/tria.h>
+
+#include <deal.II/lac/vector.h>
+
+#include <limits>
+
+#include "../tests.h"
+
+using namespace dealii;
+
+int
+main(int argc, const char *argv[])
+{
+  initlog();
+
+  Triangulation<2> tria;
+  GridGenerator::hyper_cube(tria);
+  tria.refine_global(4);
+
+  Vector<float> indicator(tria.n_active_cells());
+  for (unsigned int i = 0; i != indicator.size(); ++i)
+    {
+      indicator[i] = i;
+    }
+
+  deallog << "n_active_cells: " << tria.n_active_cells() << std::endl;
+
+  GridRefinement::refine_and_coarsen_fixed_fraction(
+    tria,
+    indicator,
+    0.6,
+    1. - 0.6 + std::numeric_limits<double>::epsilon() * 5);
+  tria.execute_coarsening_and_refinement();
+
+  deallog << "n_active_cells: " << tria.n_active_cells() << std::endl;
+
+  indicator.reinit(tria.n_active_cells());
+  for (unsigned int i = 0; i != indicator.size(); ++i)
+    {
+      indicator[i] = i;
+    }
+
+  GridRefinement::refine_and_coarsen_fixed_fraction(
+    tria,
+    indicator,
+    0.4,
+    1. - 0.4 + std::numeric_limits<double>::epsilon() * 5);
+  tria.execute_coarsening_and_refinement();
+
+  deallog << "n_active_cells: " << tria.n_active_cells() << std::endl;
+
+  return 0;
+}

diff --git a/tests/grid/refine_and_coarsen_fixed_fraction_02.output b/tests/grid/refine_and_coarsen_fixed_fraction_02.output
new file mode 100644 (file)
index 0000000..49e6993
--- /dev/null
+++ b/tests/grid/refine_and_coarsen_fixed_fraction_02.output
@@ -0,0 +1,4 @@
+
+DEAL::n_active_cells: 256
+DEAL::n_active_cells: 442
+DEAL::n_active_cells: 523

diff --git a/tests/grid/refine_and_coarsen_fixed_number_02.cc b/tests/grid/refine_and_coarsen_fixed_number_02.cc
new file mode 100644 (file)
index 0000000..b93932b
--- /dev/null
+++ b/tests/grid/refine_and_coarsen_fixed_number_02.cc
@@ -0,0 +1,74 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2015 - 2017 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.md at
+// the top level directory of deal.II.
+//
+// ---------------------------------------------------------------------
+
+// When setting refine_fraction=0.6 and coarsen_fraction=0.4, these
+// numbers may add up to slightly larger than a total of 1.0 in
+// floating point arithmetic. Make sure we allow this nevertheless.
+
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_refinement.h>
+#include <deal.II/grid/tria.h>
+
+#include <deal.II/lac/vector.h>
+
+#include <limits>
+
+#include "../tests.h"
+
+using namespace dealii;
+
+int
+main(int argc, const char *argv[])
+{
+  initlog();
+
+  Triangulation<2> tria;
+  GridGenerator::hyper_cube(tria);
+  tria.refine_global(4);
+
+  Vector<float> indicator(tria.n_active_cells());
+  for (unsigned int i = 0; i != indicator.size(); ++i)
+    {
+      indicator[i] = i;
+    }
+
+  deallog << "n_active_cells: " << tria.n_active_cells() << std::endl;
+
+  GridRefinement::refine_and_coarsen_fixed_number(
+    tria,
+    indicator,
+    0.6,
+    1. - 0.6 + std::numeric_limits<double>::epsilon() * 5);
+  tria.execute_coarsening_and_refinement();
+
+  deallog << "n_active_cells: " << tria.n_active_cells() << std::endl;
+
+  indicator.reinit(tria.n_active_cells());
+  for (unsigned int i = 0; i != indicator.size(); ++i)
+    {
+      indicator[i] = i;
+    }
+
+  GridRefinement::refine_and_coarsen_fixed_number(
+    tria,
+    indicator,
+    0.4,
+    1. - 0.4 + std::numeric_limits<double>::epsilon() * 5);
+  tria.execute_coarsening_and_refinement();
+
+  deallog << "n_active_cells: " << tria.n_active_cells() << std::endl;
+
+  return 0;
+}

diff --git a/tests/grid/refine_and_coarsen_fixed_number_02.output b/tests/grid/refine_and_coarsen_fixed_number_02.output
new file mode 100644 (file)
index 0000000..07020bc
--- /dev/null
+++ b/tests/grid/refine_and_coarsen_fixed_number_02.output
@@ -0,0 +1,4 @@
+
+DEAL::n_active_cells: 256
+DEAL::n_active_cells: 664
+DEAL::n_active_cells: 1225