Fixed test output.

diff --git a/tests/bits/find_cell_11.cc b/tests/bits/find_cell_11.cc
index 92af190893bcf8e027ce0c5095ee0bc677bbb8d2..d69343e0294843b67fdd42f66c0a77340f6993d6 100644 (file)
--- a/tests/bits/find_cell_11.cc
+++ b/tests/bits/find_cell_11.cc
@@ -56,14 +56,9 @@ test()
 
   MappingQ<2> mapping(1);
 
-  try
-    {
-      GridTools::find_active_cell_around_point(mapping, tr, p);
-    }
-  catch (GridTools::ExcPointNotFound<2> &e)
-    {
-      deallog << "outside" << std::endl;
-    }
+  auto c = GridTools::find_active_cell_around_point(mapping, tr, p);
+  if (c.first.state() != IteratorState::valid)
+    deallog << "outside" << std::endl;
   deallog << "done" << std::endl;
 }
 

diff --git a/tests/bits/find_cell_12.cc b/tests/bits/find_cell_12.cc
index 7921f2fd927710482e758c3842ebaa0f0d0c635f..a084fe28b4784f08375bf56c2d924ea18237bb85 100644 (file)
--- a/tests/bits/find_cell_12.cc
+++ b/tests/bits/find_cell_12.cc
@@ -74,18 +74,17 @@ test()
 
   Point<2> test_point(250, 195);
   std::cout << "Checking Point " << test_point << std::endl;
-  try
+  auto current_cell =
+    GridTools::find_active_cell_around_point(MappingQGeneric<2>(1),
+                                             triangulation,
+                                             test_point);
+  if (current_cell.first.state() == IteratorState::valid)
     {
-      std::pair<Triangulation<2>::active_cell_iterator, Point<2>> current_cell =
-        GridTools::find_active_cell_around_point(MappingQGeneric<2>(1),
-                                                 triangulation,
-                                                 test_point);
-
       deallog << "cell: index = " << current_cell.first->index()
               << " level = " << current_cell.first->level() << std::endl;
       deallog << " pos: " << current_cell.second << std::endl;
     }
-  catch (GridTools::ExcPointNotFound<2> &e)
+  else
     {
       deallog << "outside" << std::endl;
     }

diff --git a/tests/bits/find_cell_13.cc b/tests/bits/find_cell_13.cc
index 803090451c96c076c69c03252cca7abb6acef81e..35a7e649f34a7a6d1a5e1cf79d7bde6b951aff32 100644 (file)
--- a/tests/bits/find_cell_13.cc
+++ b/tests/bits/find_cell_13.cc
@@ -39,53 +39,39 @@ void check(Triangulation<2> &tria)
   // any point in this domain should be inside one of the cells
   Point<2> p(random_value<double>(), random_value<double>() / 8.);
 
-  try
-    {
-      // don't mark any vertex at first
-      std::vector<bool> marked_vertices(tria.n_vertices(), false);
+  // don't mark any vertex at first
+  std::vector<bool> marked_vertices(tria.n_vertices(), false);
 
-      // find the closes vertex to (1.,1.)
-      // (default call to find_closest_vertex() without passing marked_vertices)
-      unsigned int closest_vertex =
-        GridTools::find_closest_vertex(tria, Point<2>(1., 1.));
+  // find the closes vertex to (1.,1.)
+  // (default call to find_closest_vertex() without passing marked_vertices)
+  unsigned int closest_vertex =
+    GridTools::find_closest_vertex(tria, Point<2>(1., 1.));
 
-      // mark the vertex closest to (1.,1.)
-      marked_vertices[closest_vertex] = true;
+  // mark the vertex closest to (1.,1.)
+  marked_vertices[closest_vertex] = true;
 
-      auto vertices = tria.get_vertices();
-      deallog << vertices[closest_vertex] << " is the closest vertex"
-              << std::endl;
+  auto vertices = tria.get_vertices();
+  deallog << vertices[closest_vertex] << " is the closest vertex" << std::endl;
 
-      GridTools::find_active_cell_around_point(tria, p, marked_vertices);
+  auto c = GridTools::find_active_cell_around_point(tria, p, marked_vertices);
 
-      // The test fails if the control reaches here.
-      deallog << "Garbage text to make the diff fail if the control "
-              << "reaches here:"
-              << " Point: " << p;
-      deallog << std::endl; // Flush deallog buffer
-    }
-  catch (const GridTools::ExcPointNotFound<2> &)
+  if (c.state() != IteratorState::valid)
     {
       deallog
-        << "The first call to the function find_closest_vertex() has thrown "
-        << "an exception. It is supposed to throw.";
-      deallog << std::endl;
-
-      // The default function call doesn't throw exceptions
-      // for this use case
-      Triangulation<2>::active_cell_iterator cell =
-        GridTools::find_active_cell_around_point(tria, p);
-
-      // check if the below function call actually finds appropriate cell
-      Assert(p.distance(cell->center()) < cell->diameter() / 2,
-             ExcInternalError());
-      deallog << "Test passed!";
-      deallog << std::endl;
-    }
-  catch (...)
-    {
-      Assert(false, ExcInternalError());
+        << "The first call to the function find_active_cell_around_point()"
+        << std::endl
+        << "has returned an invalid iterator. This is good." << std::endl;
     }
+
+  // The default function call doesn't throw exceptions
+  // for this use case
+  Triangulation<2>::active_cell_iterator cell =
+    GridTools::find_active_cell_around_point(tria, p);
+
+  // check if the below function call actually finds appropriate cell
+  Assert(p.distance(cell->center()) < cell->diameter() / 2, ExcInternalError());
+  deallog << "Test passed!";
+  deallog << std::endl;
 }
 
 

diff --git a/tests/bits/find_cell_13.output b/tests/bits/find_cell_13.output
index e10f71dffa98a9728e3183123bab6ee23941fe16..174d37b1d16d6ca2b834c96ff0b618ae85c61d8e 100644 (file)
--- a/tests/bits/find_cell_13.output
+++ b/tests/bits/find_cell_13.output
@@ -1,7 +1,9 @@
 
 DEAL::1.00000 1.00000 is the closest vertex
-DEAL::The first call to the function find_closest_vertex() has thrown an exception. It is supposed to throw.
+DEAL::The first call to the function find_active_cell_around_point()
+DEAL::has returned an invalid iterator. This is good.
 DEAL::Test passed!
 DEAL::0.707107 0.707107 is the closest vertex
-DEAL::The first call to the function find_closest_vertex() has thrown an exception. It is supposed to throw.
+DEAL::The first call to the function find_active_cell_around_point()
+DEAL::has returned an invalid iterator. This is good.
 DEAL::Test passed!

diff --git a/tests/numerics/generalized_interpolation_02.output b/tests/numerics/generalized_interpolation_02.output
index 29294e52f2879a798e41dc74c8ae53d6da4fe1ea..3ce63ec0b594527faa6c12ead20c8969b9efdf13 100644 (file)
--- a/tests/numerics/generalized_interpolation_02.output
+++ b/tests/numerics/generalized_interpolation_02.output
@@ -4,61 +4,61 @@ DEAL::Function value at (0.0,0.0): 1.00000
 DEAL::Check projection property: 0.00000
 DEAL::dim 2 FE_Q<2>(2)
 DEAL::Function value at (0.0,0.0): 1.00000 
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 2.22045e-16
 DEAL::dim 3 FE_Q<3>(3)
 DEAL::Function value at (0.0,0.0): 1.00000 
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 8.88178e-16
 DEAL::dim 2 FE_Q<2>(1)
 DEAL::Function value at (0.0,0.0): 1.00000 
 DEAL::Check projection property: 0.00000
 DEAL::dim 2 FE_Q<2>(2)
 DEAL::Function value at (0.0,0.0): 1.00000 
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 2.22045e-16
 DEAL::dim 3 FE_Q<3>(3)
 DEAL::Function value at (0.0,0.0): 1.00000 
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 8.88178e-16
 DEAL::dim 2 FE_RaviartThomas<2>(0)
 DEAL::Function value at (0.0,0.0): 1.00000 1.00000 
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 5.55112e-17
 DEAL::dim 2 FE_RaviartThomas<2>(1)
 DEAL::Function value at (0.0,0.0): 1.00000 1.00000 
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 2.35922e-16
 DEAL::dim 2 FE_RaviartThomas<2>(2)
 DEAL::Function value at (0.0,0.0): 1.00000 1.00000 
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 0.169271
 DEAL::dim 3 FE_RaviartThomas<3>(0)
 DEAL::Function value at (0.0,0.0): 1.00000 1.00000 1.00000 
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 1.11022e-16
 DEAL::dim 3 FE_RaviartThomas<3>(1)
 DEAL::Function value at (0.0,0.0): 1.00000 1.00000 1.00000 
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 2.91434e-16
 DEAL::dim 2 FE_RaviartThomas<2>(0)
 DEAL::Function value at (0.0,0.0): 1.00000 1.00000 
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 1.11022e-16
 DEAL::dim 2 FE_RaviartThomas<2>(1)
 DEAL::Function value at (0.0,0.0): 1.00000 1.00000 
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 4.44089e-16
 DEAL::dim 2 FE_RaviartThomas<2>(2)
 DEAL::Function value at (0.0,0.0): 1.00000 1.00000 
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 0.169442
 DEAL::dim 3 FE_RaviartThomas<3>(1)
 DEAL::Function value at (0.0,0.0): 1.00000 1.00000 1.00000 
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 2.42861e-16
 DEAL::dim 2 FE_Nedelec<2>(0)
-DEAL::Function value at (0.0,0.0): 1.00000 1.00000
-DEAL::Check projection property: 0.00000
+DEAL::Function value at (0.0,0.0): 1.00000 1.00000 
+DEAL::Check projection property: 2.77556e-17
 DEAL::dim 2 FE_Nedelec<2>(1)
-DEAL::Function value at (0.0,0.0): 1.00000 1.00000
-DEAL::Check projection property: 0.00000
+DEAL::Function value at (0.0,0.0): 1.00000 1.00000 
+DEAL::Check projection property: 4.44089e-16
 DEAL::dim 2 FE_Nedelec<2>(2)
-DEAL::Function value at (0.0,0.0): 1.00000 1.00000
-DEAL::Check projection property: 0.00000
+DEAL::Function value at (0.0,0.0): 1.00000 1.00000 
+DEAL::Check projection property: 0.378501
 DEAL::dim 2 FE_Nedelec<2>(0)
-DEAL::Function value at (0.0,0.0): 1.00000 1.00000
-DEAL::Check projection property: 0.00000
+DEAL::Function value at (0.0,0.0): 1.00000 1.00000 
+DEAL::Check projection property: 1.11022e-16
 DEAL::dim 2 FE_Nedelec<2>(1)
-DEAL::Function value at (0.0,0.0): 1.00000 1.00000
-DEAL::Check projection property: 0.00000
+DEAL::Function value at (0.0,0.0): 1.00000 1.00000 
+DEAL::Check projection property: 4.44089e-16
 DEAL::dim 2 FE_Nedelec<2>(2)
-DEAL::Function value at (0.0,0.0): 1.00000 1.00000
-DEAL::Check projection property: 0.00000
+DEAL::Function value at (0.0,0.0): 1.00000 1.00000 
+DEAL::Check projection property: 0.379571

diff --git a/tests/numerics/generalized_interpolation_04.output b/tests/numerics/generalized_interpolation_04.output
index ce616a49ad48b7dbb7fecc20596d5b4f4a6ac8d6..a1da03bd38a896c0fd968065102474848eab651d 100644 (file)
--- a/tests/numerics/generalized_interpolation_04.output
+++ b/tests/numerics/generalized_interpolation_04.output
@@ -1,3 +1,3 @@
 
 DEAL::dim 2 FESystem<2>[FE_RaviartThomas<2>(1)-FE_Q<2>(1)-FE_Nedelec<2>(2)^2-FESystem<2>[FE_Q<2>(1)^2]]
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 0.699901

diff --git a/tests/numerics/generalized_interpolation_05.output b/tests/numerics/generalized_interpolation_05.output
index e6dc071228a77e6a81fc3d474e97c66fc493fb60..ffc98b78f61d8cfd46080ca07c39255deb7765d8 100644 (file)
--- a/tests/numerics/generalized_interpolation_05.output
+++ b/tests/numerics/generalized_interpolation_05.output
@@ -1,3 +1,3 @@
 
 DEAL::dim 2 FESystem<2>[FE_RaviartThomas<2>(1)^2-FE_Q<2>(1)-FESystem<2>[FE_Nedelec<2>(2)^2]-FESystem<2>[FE_Q<2>(1)^2]]
-DEAL::Check projection property: 0.00000
+DEAL::Check projection property: 0.699901

diff --git a/tests/simplex/compute_point_locations_01.output b/tests/simplex/compute_point_locations_01.output
index 087ff3025fa63586b93d515443c80fc47ee146d2..4cf9614adcb88d875103d8bfd70faf84d6c532a7 100644 (file)
--- a/tests/simplex/compute_point_locations_01.output
+++ b/tests/simplex/compute_point_locations_01.output
@@ -70,87 +70,140 @@ DEAL::
 DEAL::Print out results in 3D: 
 DEAL::At cell 0:
 DEAL::    qpoints 0: 
-DEAL::        reference position  : (0.00000 0.00000 0.00000)
-DEAL::        physical position   : (0.00000 0.00000 0.00000)
-DEAL::        FE mapping position : (0.00000 0.00000 0.00000)
-DEAL::    qpoints 1: 
 DEAL::        reference position  : (0.00000 0.00000 0.333333)
 DEAL::        physical position   : (0.00000 0.00000 0.333333)
 DEAL::        FE mapping position : (0.00000 0.00000 0.333333)
+DEAL::    qpoints 1: 
+DEAL::        reference position  : (0.333333 0.00000 0.666667)
+DEAL::        physical position   : (0.333333 0.00000 0.666667)
+DEAL::        FE mapping position : (0.333333 0.00000 0.666667)
 DEAL::    qpoints 2: 
+DEAL::        reference position  : (0.333333 0.00000 0.333333)
+DEAL::        physical position   : (0.333333 0.00000 0.333333)
+DEAL::        FE mapping position : (0.333333 0.00000 0.333333)
+DEAL::    qpoints 3: 
+DEAL::        reference position  : (0.333333 0.00000 0.00000)
+DEAL::        physical position   : (0.333333 0.00000 0.00000)
+DEAL::        FE mapping position : (0.333333 0.00000 0.00000)
+DEAL::    qpoints 4: 
+DEAL::        reference position  : (0.00000 0.00000 0.00000)
+DEAL::        physical position   : (0.00000 0.00000 0.00000)
+DEAL::        FE mapping position : (0.00000 0.00000 0.00000)
+DEAL::    qpoints 5: 
+DEAL::        reference position  : (0.333333 0.333333 0.333333)
+DEAL::        physical position   : (0.333333 0.333333 0.333333)
+DEAL::        FE mapping position : (0.333333 0.333333 0.333333)
+DEAL::    qpoints 6: 
+DEAL::        reference position  : (0.333333 0.333333 0.00000)
+DEAL::        physical position   : (0.333333 0.333333 0.00000)
+DEAL::        FE mapping position : (0.333333 0.333333 0.00000)
+DEAL::    qpoints 7: 
 DEAL::        reference position  : (0.00000 0.00000 0.666667)
 DEAL::        physical position   : (0.00000 0.00000 0.666667)
 DEAL::        FE mapping position : (0.00000 0.00000 0.666667)
-DEAL::    qpoints 3: 
+DEAL::    qpoints 8: 
 DEAL::        reference position  : (0.00000 0.00000 1.00000)
 DEAL::        physical position   : (0.00000 0.00000 1.00000)
 DEAL::        FE mapping position : (0.00000 0.00000 1.00000)
-DEAL::    qpoints 4: 
+DEAL::    qpoints 9: 
 DEAL::        reference position  : (0.00000 0.333333 0.00000)
 DEAL::        physical position   : (0.00000 0.333333 0.00000)
 DEAL::        FE mapping position : (0.00000 0.333333 0.00000)
-DEAL::    qpoints 5: 
+DEAL::    qpoints 10: 
 DEAL::        reference position  : (0.00000 0.333333 0.333333)
 DEAL::        physical position   : (0.00000 0.333333 0.333333)
 DEAL::        FE mapping position : (0.00000 0.333333 0.333333)
-DEAL::    qpoints 6: 
+DEAL::    qpoints 11: 
 DEAL::        reference position  : (0.00000 0.333333 0.666667)
 DEAL::        physical position   : (0.00000 0.333333 0.666667)
 DEAL::        FE mapping position : (0.00000 0.333333 0.666667)
-DEAL::    qpoints 7: 
-DEAL::        reference position  : (0.00000 0.666667 0.00000)
-DEAL::        physical position   : (0.00000 0.666667 0.00000)
-DEAL::        FE mapping position : (0.00000 0.666667 0.00000)
-DEAL::    qpoints 8: 
+DEAL::    qpoints 12: 
 DEAL::        reference position  : (0.00000 0.666667 0.333333)
 DEAL::        physical position   : (0.00000 0.666667 0.333333)
 DEAL::        FE mapping position : (0.00000 0.666667 0.333333)
-DEAL::    qpoints 9: 
-DEAL::        reference position  : (0.00000 1.00000 0.00000)
-DEAL::        physical position   : (0.00000 1.00000 0.00000)
-DEAL::        FE mapping position : (0.00000 1.00000 0.00000)
-DEAL::    qpoints 10: 
-DEAL::        reference position  : (0.333333 0.00000 0.00000)
-DEAL::        physical position   : (0.333333 0.00000 0.00000)
-DEAL::        FE mapping position : (0.333333 0.00000 0.00000)
-DEAL::    qpoints 11: 
-DEAL::        reference position  : (0.333333 0.00000 0.333333)
-DEAL::        physical position   : (0.333333 0.00000 0.333333)
-DEAL::        FE mapping position : (0.333333 0.00000 0.333333)
-DEAL::    qpoints 12: 
-DEAL::        reference position  : (0.333333 0.00000 0.666667)
-DEAL::        physical position   : (0.333333 0.00000 0.666667)
-DEAL::        FE mapping position : (0.333333 0.00000 0.666667)
 DEAL::    qpoints 13: 
-DEAL::        reference position  : (0.333333 0.333333 0.00000)
-DEAL::        physical position   : (0.333333 0.333333 0.00000)
-DEAL::        FE mapping position : (0.333333 0.333333 0.00000)
+DEAL::        reference position  : (0.00000 0.666667 0.00000)
+DEAL::        physical position   : (0.00000 0.666667 0.00000)
+DEAL::        FE mapping position : (0.00000 0.666667 0.00000)
 DEAL::    qpoints 14: 
-DEAL::        reference position  : (0.333333 0.333333 0.333333)
-DEAL::        physical position   : (0.333333 0.333333 0.333333)
-DEAL::        FE mapping position : (0.333333 0.333333 0.333333)
-DEAL::    qpoints 15: 
 DEAL::        reference position  : (0.333333 0.666667 0.00000)
 DEAL::        physical position   : (0.333333 0.666667 0.00000)
 DEAL::        FE mapping position : (0.333333 0.666667 0.00000)
+DEAL::    qpoints 15: 
+DEAL::        reference position  : (0.00000 1.00000 0.00000)
+DEAL::        physical position   : (0.00000 1.00000 0.00000)
+DEAL::        FE mapping position : (0.00000 1.00000 0.00000)
 DEAL::    qpoints 16: 
 DEAL::        reference position  : (0.666667 0.00000 0.00000)
 DEAL::        physical position   : (0.666667 0.00000 0.00000)
 DEAL::        FE mapping position : (0.666667 0.00000 0.00000)
 DEAL::    qpoints 17: 
+DEAL::        reference position  : (1.00000 0.00000 0.00000)
+DEAL::        physical position   : (1.00000 0.00000 0.00000)
+DEAL::        FE mapping position : (1.00000 0.00000 0.00000)
+DEAL::    qpoints 18: 
 DEAL::        reference position  : (0.666667 0.00000 0.333333)
 DEAL::        physical position   : (0.666667 0.00000 0.333333)
 DEAL::        FE mapping position : (0.666667 0.00000 0.333333)
-DEAL::    qpoints 18: 
+DEAL::    qpoints 19: 
 DEAL::        reference position  : (0.666667 0.333333 0.00000)
 DEAL::        physical position   : (0.666667 0.333333 0.00000)
 DEAL::        FE mapping position : (0.666667 0.333333 0.00000)
-DEAL::    qpoints 19: 
-DEAL::        reference position  : (1.00000 0.00000 0.00000)
-DEAL::        physical position   : (1.00000 0.00000 0.00000)
-DEAL::        FE mapping position : (1.00000 0.00000 0.00000)
 DEAL::At cell 1:
 DEAL::    qpoints 0: 
+DEAL::        reference position  : (0.666667 0.333333 0.00000)
+DEAL::        physical position   : (0.333333 0.00000 1.00000)
+DEAL::        FE mapping position : (0.333333 0.00000 1.00000)
+DEAL::    qpoints 1: 
+DEAL::        reference position  : (0.666667 0.00000 0.333333)
+DEAL::        physical position   : (0.333333 0.333333 1.00000)
+DEAL::        FE mapping position : (0.333333 0.333333 1.00000)
+DEAL::    qpoints 2: 
+DEAL::        reference position  : (0.00000 1.00000 0.00000)
+DEAL::        physical position   : (1.00000 0.00000 1.00000)
+DEAL::        FE mapping position : (1.00000 0.00000 1.00000)
+DEAL::    qpoints 3: 
+DEAL::        reference position  : (0.00000 0.666667 0.00000)
+DEAL::        physical position   : (1.00000 0.00000 0.666667)
+DEAL::        FE mapping position : (1.00000 0.00000 0.666667)
+DEAL::    qpoints 4: 
+DEAL::        reference position  : (0.00000 0.333333 0.00000)
+DEAL::        physical position   : (1.00000 0.00000 0.333333)
+DEAL::        FE mapping position : (1.00000 0.00000 0.333333)
+DEAL::    qpoints 5: 
+DEAL::        reference position  : (0.333333 0.333333 0.00000)
+DEAL::        physical position   : (0.666667 0.00000 0.666667)
+DEAL::        FE mapping position : (0.666667 0.00000 0.666667)
+DEAL::    qpoints 6: 
+DEAL::        reference position  : (0.333333 0.666667 0.00000)
+DEAL::        physical position   : (0.666667 0.00000 1.00000)
+DEAL::        FE mapping position : (0.666667 0.00000 1.00000)
+DEAL::    qpoints 7: 
+DEAL::        reference position  : (0.333333 0.333333 0.333333)
+DEAL::        physical position   : (0.666667 0.333333 1.00000)
+DEAL::        FE mapping position : (0.666667 0.333333 1.00000)
+DEAL::    qpoints 8: 
+DEAL::        reference position  : (0.333333 0.00000 0.333333)
+DEAL::        physical position   : (0.666667 0.333333 0.666667)
+DEAL::        FE mapping position : (0.666667 0.333333 0.666667)
+DEAL::    qpoints 9: 
+DEAL::        reference position  : (0.00000 0.666667 0.333333)
+DEAL::        physical position   : (1.00000 0.333333 1.00000)
+DEAL::        FE mapping position : (1.00000 0.333333 1.00000)
+DEAL::    qpoints 10: 
+DEAL::        reference position  : (0.00000 0.333333 0.333333)
+DEAL::        physical position   : (1.00000 0.333333 0.666667)
+DEAL::        FE mapping position : (1.00000 0.333333 0.666667)
+DEAL::    qpoints 11: 
+DEAL::        reference position  : (0.00000 0.333333 0.666667)
+DEAL::        physical position   : (1.00000 0.666667 1.00000)
+DEAL::        FE mapping position : (1.00000 0.666667 1.00000)
+DEAL::    qpoints 12: 
+DEAL::        reference position  : (0.333333 0.00000 0.666667)
+DEAL::        physical position   : (0.666667 0.666667 1.00000)
+DEAL::        FE mapping position : (0.666667 0.666667 1.00000)
+DEAL::At cell 2:
+DEAL::    qpoints 0: 
 DEAL::        reference position  : (0.666667 0.00000 0.333333)
 DEAL::        physical position   : (0.00000 0.333333 1.00000)
 DEAL::        FE mapping position : (0.00000 0.333333 1.00000)
@@ -163,90 +216,37 @@ DEAL::        reference position  : (0.333333 0.00000 0.666667)
 DEAL::        physical position   : (0.00000 0.666667 1.00000)
 DEAL::        FE mapping position : (0.00000 0.666667 1.00000)
 DEAL::    qpoints 3: 
-DEAL::        reference position  : (0.00000 0.00000 0.333333)
-DEAL::        physical position   : (0.00000 1.00000 0.333333)
-DEAL::        FE mapping position : (0.00000 1.00000 0.333333)
+DEAL::        reference position  : (0.333333 0.333333 0.00000)
+DEAL::        physical position   : (0.333333 0.666667 0.666667)
+DEAL::        FE mapping position : (0.333333 0.666667 0.666667)
 DEAL::    qpoints 4: 
-DEAL::        reference position  : (0.00000 0.00000 0.666667)
-DEAL::        physical position   : (0.00000 1.00000 0.666667)
-DEAL::        FE mapping position : (0.00000 1.00000 0.666667)
-DEAL::    qpoints 5: 
 DEAL::        reference position  : (0.00000 0.00000 1.00000)
 DEAL::        physical position   : (0.00000 1.00000 1.00000)
 DEAL::        FE mapping position : (0.00000 1.00000 1.00000)
+DEAL::    qpoints 5: 
+DEAL::        reference position  : (0.00000 0.00000 0.333333)
+DEAL::        physical position   : (0.00000 1.00000 0.333333)
+DEAL::        FE mapping position : (0.00000 1.00000 0.333333)
 DEAL::    qpoints 6: 
-DEAL::        reference position  : (0.333333 0.333333 -5.55112e-17)
-DEAL::        physical position   : (0.333333 0.666667 0.666667)
-DEAL::        FE mapping position : (0.333333 0.666667 0.666667)
-DEAL::    qpoints 7: 
 DEAL::        reference position  : (0.333333 0.333333 0.333333)
 DEAL::        physical position   : (0.333333 0.666667 1.00000)
 DEAL::        FE mapping position : (0.333333 0.666667 1.00000)
-DEAL::    qpoints 8: 
+DEAL::    qpoints 7: 
 DEAL::        reference position  : (0.00000 0.333333 0.333333)
 DEAL::        physical position   : (0.333333 1.00000 0.666667)
 DEAL::        FE mapping position : (0.333333 1.00000 0.666667)
-DEAL::    qpoints 9: 
+DEAL::    qpoints 8: 
 DEAL::        reference position  : (0.00000 0.333333 0.666667)
 DEAL::        physical position   : (0.333333 1.00000 1.00000)
 DEAL::        FE mapping position : (0.333333 1.00000 1.00000)
+DEAL::    qpoints 9: 
+DEAL::        reference position  : (0.00000 0.00000 0.666667)
+DEAL::        physical position   : (0.00000 1.00000 0.666667)
+DEAL::        FE mapping position : (0.00000 1.00000 0.666667)
 DEAL::    qpoints 10: 
 DEAL::        reference position  : (0.00000 0.666667 0.333333)
 DEAL::        physical position   : (0.666667 1.00000 1.00000)
 DEAL::        FE mapping position : (0.666667 1.00000 1.00000)
-DEAL::At cell 2:
-DEAL::    qpoints 0: 
-DEAL::        reference position  : (0.666667 0.333333 0.00000)
-DEAL::        physical position   : (0.333333 0.00000 1.00000)
-DEAL::        FE mapping position : (0.333333 0.00000 1.00000)
-DEAL::    qpoints 1: 
-DEAL::        reference position  : (0.666667 0.00000 0.333333)
-DEAL::        physical position   : (0.333333 0.333333 1.00000)
-DEAL::        FE mapping position : (0.333333 0.333333 1.00000)
-DEAL::    qpoints 2: 
-DEAL::        reference position  : (0.333333 0.333333 0.00000)
-DEAL::        physical position   : (0.666667 0.00000 0.666667)
-DEAL::        FE mapping position : (0.666667 0.00000 0.666667)
-DEAL::    qpoints 3: 
-DEAL::        reference position  : (0.333333 0.666667 0.00000)
-DEAL::        physical position   : (0.666667 0.00000 1.00000)
-DEAL::        FE mapping position : (0.666667 0.00000 1.00000)
-DEAL::    qpoints 4: 
-DEAL::        reference position  : (0.333333 -5.55112e-17 0.333333)
-DEAL::        physical position   : (0.666667 0.333333 0.666667)
-DEAL::        FE mapping position : (0.666667 0.333333 0.666667)
-DEAL::    qpoints 5: 
-DEAL::        reference position  : (0.333333 0.333333 0.333333)
-DEAL::        physical position   : (0.666667 0.333333 1.00000)
-DEAL::        FE mapping position : (0.666667 0.333333 1.00000)
-DEAL::    qpoints 6: 
-DEAL::        reference position  : (0.333333 0.00000 0.666667)
-DEAL::        physical position   : (0.666667 0.666667 1.00000)
-DEAL::        FE mapping position : (0.666667 0.666667 1.00000)
-DEAL::    qpoints 7: 
-DEAL::        reference position  : (0.00000 0.333333 0.00000)
-DEAL::        physical position   : (1.00000 0.00000 0.333333)
-DEAL::        FE mapping position : (1.00000 0.00000 0.333333)
-DEAL::    qpoints 8: 
-DEAL::        reference position  : (0.00000 0.666667 0.00000)
-DEAL::        physical position   : (1.00000 0.00000 0.666667)
-DEAL::        FE mapping position : (1.00000 0.00000 0.666667)
-DEAL::    qpoints 9: 
-DEAL::        reference position  : (0.00000 1.00000 0.00000)
-DEAL::        physical position   : (1.00000 0.00000 1.00000)
-DEAL::        FE mapping position : (1.00000 0.00000 1.00000)
-DEAL::    qpoints 10: 
-DEAL::        reference position  : (0.00000 0.333333 0.333333)
-DEAL::        physical position   : (1.00000 0.333333 0.666667)
-DEAL::        FE mapping position : (1.00000 0.333333 0.666667)
-DEAL::    qpoints 11: 
-DEAL::        reference position  : (0.00000 0.666667 0.333333)
-DEAL::        physical position   : (1.00000 0.333333 1.00000)
-DEAL::        FE mapping position : (1.00000 0.333333 1.00000)
-DEAL::    qpoints 12: 
-DEAL::        reference position  : (0.00000 0.333333 0.666667)
-DEAL::        physical position   : (1.00000 0.666667 1.00000)
-DEAL::        FE mapping position : (1.00000 0.666667 1.00000)
 DEAL::At cell 3:
 DEAL::    qpoints 0: 
 DEAL::        reference position  : (0.166667 0.500000 0.166667)
@@ -270,49 +270,49 @@ DEAL::        reference position  : (0.00000 0.333333 0.00000)
 DEAL::        physical position   : (0.333333 1.00000 0.00000)
 DEAL::        FE mapping position : (0.333333 1.00000 0.00000)
 DEAL::    qpoints 1: 
-DEAL::        reference position  : (0.00000 0.00000 0.333333)
+DEAL::        reference position  : (0.00000 1.11022e-16 0.333333)
 DEAL::        physical position   : (0.333333 1.00000 0.333333)
 DEAL::        FE mapping position : (0.333333 1.00000 0.333333)
 DEAL::    qpoints 2: 
-DEAL::        reference position  : (0.333333 0.333333 0.00000)
-DEAL::        physical position   : (0.666667 0.666667 0.00000)
-DEAL::        FE mapping position : (0.666667 0.666667 0.00000)
-DEAL::    qpoints 3: 
-DEAL::        reference position  : (0.333333 -5.55112e-17 0.333333)
-DEAL::        physical position   : (0.666667 0.666667 0.333333)
-DEAL::        FE mapping position : (0.666667 0.666667 0.333333)
-DEAL::    qpoints 4: 
-DEAL::        reference position  : (0.00000 0.666667 0.00000)
-DEAL::        physical position   : (0.666667 1.00000 0.00000)
-DEAL::        FE mapping position : (0.666667 1.00000 0.00000)
-DEAL::    qpoints 5: 
-DEAL::        reference position  : (0.00000 0.333333 0.333333)
-DEAL::        physical position   : (0.666667 1.00000 0.333333)
-DEAL::        FE mapping position : (0.666667 1.00000 0.333333)
-DEAL::    qpoints 6: 
-DEAL::        reference position  : (0.00000 0.00000 0.666667)
-DEAL::        physical position   : (0.666667 1.00000 0.666667)
-DEAL::        FE mapping position : (0.666667 1.00000 0.666667)
-DEAL::    qpoints 7: 
 DEAL::        reference position  : (0.666667 0.333333 0.00000)
 DEAL::        physical position   : (1.00000 0.333333 0.00000)
 DEAL::        FE mapping position : (1.00000 0.333333 0.00000)
-DEAL::    qpoints 8: 
+DEAL::    qpoints 3: 
 DEAL::        reference position  : (0.666667 0.00000 0.333333)
 DEAL::        physical position   : (1.00000 0.333333 0.333333)
 DEAL::        FE mapping position : (1.00000 0.333333 0.333333)
-DEAL::    qpoints 9: 
+DEAL::    qpoints 4: 
+DEAL::        reference position  : (0.333333 0.00000 0.666667)
+DEAL::        physical position   : (1.00000 0.666667 0.666667)
+DEAL::        FE mapping position : (1.00000 0.666667 0.666667)
+DEAL::    qpoints 5: 
+DEAL::        reference position  : (0.333333 0.333333 0.333333)
+DEAL::        physical position   : (1.00000 0.666667 0.333333)
+DEAL::        FE mapping position : (1.00000 0.666667 0.333333)
+DEAL::    qpoints 6: 
 DEAL::        reference position  : (0.333333 0.666667 0.00000)
 DEAL::        physical position   : (1.00000 0.666667 0.00000)
 DEAL::        FE mapping position : (1.00000 0.666667 0.00000)
+DEAL::    qpoints 7: 
+DEAL::        reference position  : (0.333333 0.333333 0.00000)
+DEAL::        physical position   : (0.666667 0.666667 0.00000)
+DEAL::        FE mapping position : (0.666667 0.666667 0.00000)
+DEAL::    qpoints 8: 
+DEAL::        reference position  : (0.00000 2.22045e-16 0.666667)
+DEAL::        physical position   : (0.666667 1.00000 0.666667)
+DEAL::        FE mapping position : (0.666667 1.00000 0.666667)
+DEAL::    qpoints 9: 
+DEAL::        reference position  : (0.00000 0.333333 0.333333)
+DEAL::        physical position   : (0.666667 1.00000 0.333333)
+DEAL::        FE mapping position : (0.666667 1.00000 0.333333)
 DEAL::    qpoints 10: 
-DEAL::        reference position  : (0.333333 0.333333 0.333333)
-DEAL::        physical position   : (1.00000 0.666667 0.333333)
-DEAL::        FE mapping position : (1.00000 0.666667 0.333333)
+DEAL::        reference position  : (0.00000 0.666667 0.00000)
+DEAL::        physical position   : (0.666667 1.00000 0.00000)
+DEAL::        FE mapping position : (0.666667 1.00000 0.00000)
 DEAL::    qpoints 11: 
-DEAL::        reference position  : (0.333333 0.00000 0.666667)
-DEAL::        physical position   : (1.00000 0.666667 0.666667)
-DEAL::        FE mapping position : (1.00000 0.666667 0.666667)
+DEAL::        reference position  : (0.333333 0.00000 0.333333)
+DEAL::        physical position   : (0.666667 0.666667 0.333333)
+DEAL::        FE mapping position : (0.666667 0.666667 0.333333)
 DEAL::    qpoints 12: 
 DEAL::        reference position  : (0.00000 1.00000 0.00000)
 DEAL::        physical position   : (1.00000 1.00000 0.00000)