Improve GridGenerator::channel_with_cylinder

diff --git a/doc/doxygen/images/channel_with_cylinder_2d.png b/doc/doxygen/images/channel_with_cylinder_2d.png
index e22e83b1ea3c8620172cf8d352b5591f8d576184..a0aeaf977b017d2884d9664ac0f3c68ceb2f2e5b 100644 (file)
Binary files a/doc/doxygen/images/channel_with_cylinder_2d.png and b/doc/doxygen/images/channel_with_cylinder_2d.png differ

diff --git a/doc/doxygen/images/channel_with_cylinder_2d_manifolds.png b/doc/doxygen/images/channel_with_cylinder_2d_manifolds.png
index b862e0897b30b863bea79f9d51a07c8919703934..f4d53985f67a8564e8035049c668df591d16c7f4 100644 (file)
Binary files a/doc/doxygen/images/channel_with_cylinder_2d_manifolds.png and b/doc/doxygen/images/channel_with_cylinder_2d_manifolds.png differ

diff --git a/doc/doxygen/images/channel_with_cylinder_3d.png b/doc/doxygen/images/channel_with_cylinder_3d.png
index 80deecf43e7b0ae8e08268e4bb79aed840fddfa4..a219610f6512da2d736e12b0dd883dd8b6ee167c 100644 (file)
Binary files a/doc/doxygen/images/channel_with_cylinder_3d.png and b/doc/doxygen/images/channel_with_cylinder_3d.png differ

diff --git a/include/deal.II/grid/grid_generator.h b/include/deal.II/grid/grid_generator.h
index d2844c92734a924ece67cea53de060c88482ee04..a0f5f7bb80bca50cf4f6a06554da6671b60a0b0d 100644 (file)
--- a/include/deal.II/grid/grid_generator.h
+++ b/include/deal.II/grid/grid_generator.h
@@ -407,7 +407,7 @@ namespace GridGenerator
    * </ol>
    * Since the cylinder is slightly offset from the center of the channel,
    * this geometry results in vortex shedding at moderate Reynolds
-   * numbers. Here is the grid (after one global refinement) in 2D:
+   * numbers. Here is the grid (without additional global refinement) in 2D:
    *
    * @image html channel_with_cylinder_2d.png
    *

diff --git a/source/grid/grid_generator.cc b/source/grid/grid_generator.cc
index 02ecebcb3a838bf6b89236e799035ce9e8fa50f8..4285f5c816afec289b0b3ba04dd05635613469a6 100644 (file)
--- a/source/grid/grid_generator.cc
+++ b/source/grid/grid_generator.cc
@@ -3625,6 +3625,31 @@ namespace GridGenerator
     GridGenerator::merge_triangulations(
       tria_without_cylinder, cylinder_tria, tria, vertex_tolerance, true);
 
+    // Move the vertices in the middle of the faces of cylinder_tria slightly
+    // to give a better mesh quality. We have to balance the quality of these
+    // cells with the quality of the outer cells (initially rectangles). For
+    // constant radial distance, we would place them at the distance 0.1 *
+    // sqrt(2.) from the center. In case the shell region width is more than
+    // 0.1/6., we choose to place them at 0.1 * 4./3. from the center, which
+    // ensures that the shortest edge of the outer cells is 2./3. of the
+    // original length. If the shell region width is less, we make the edge
+    // length of the inner part and outer part (in the shorter x direction)
+    // the same.
+    {
+      const double shift =
+        std::min(0.125 + shell_region_width * 0.5, 0.1 * 4. / 3.);
+      for (const auto &cell : tria.active_cell_iterators())
+        for (unsigned int v = 0; v < GeometryInfo<2>::vertices_per_cell; ++v)
+          if (cell->vertex(v).distance(Point<2>(0.1, 0.205)) < 1e-10)
+            cell->vertex(v) = Point<2>(0.2 - shift, 0.205);
+          else if (cell->vertex(v).distance(Point<2>(0.3, 0.205)) < 1e-10)
+            cell->vertex(v) = Point<2>(0.2 + shift, 0.205);
+          else if (cell->vertex(v).distance(Point<2>(0.2, 0.1025)) < 1e-10)
+            cell->vertex(v) = Point<2>(0.2, 0.2 - shift);
+          else if (cell->vertex(v).distance(Point<2>(0.2, 0.3075)) < 1e-10)
+            cell->vertex(v) = Point<2>(0.2, 0.2 + shift);
+    }
+
     // Ensure that all manifold ids on a polar cell really are set to the
     // polar manifold id:
     for (const auto &cell : tria.active_cell_iterators())

diff --git a/tests/grid/channel_flow.output b/tests/grid/channel_flow.output
index 1618d8bfbcc4ac5c8955255614278c1f5705c8ff..0638df0cc0a5c3acb4bbb7c875313b19df0477b9 100644 (file)
--- a/tests/grid/channel_flow.output
+++ b/tests/grid/channel_flow.output
@@ -74,13 +74,13 @@ DEAL:2d::0.182322, 0.157322
 DEAL:2d::0.236003, 0.163997
 DEAL:2d::0.217678, 0.157322
 DEAL:2d::manifold id: 1
-DEAL:2d::0.290000, 0.202500
+DEAL:2d::0.306667, 0.202500
 DEAL:2d::0.278284, 0.282034
-DEAL:2d::0.200000, 0.293750
+DEAL:2d::0.200000, 0.306667
 DEAL:2d::0.121716, 0.282034
-DEAL:2d::0.110000, 0.202500
+DEAL:2d::0.0933333, 0.202500
 DEAL:2d::0.121716, 0.122966
-DEAL:2d::0.200000, 0.111250
+DEAL:2d::0.200000, 0.0933333
 DEAL:2d::0.278284, 0.122966
 DEAL:2d::boundary id: 0
 DEAL:2d::0.00000, 0.0512500
@@ -597,78 +597,78 @@ DEAL:3d::0.156540, 0.181998, 0.410000
 DEAL:3d::0.181998, 0.156540, 0.410000
 DEAL:3d::0.218002, 0.156540, 0.410000
 DEAL:3d::manifold id: 1
-DEAL:3d::0.284142, 0.162733, 0.00000
-DEAL:3d::0.290000, 0.202500, 0.0512500
-DEAL:3d::0.284142, 0.242267, 0.00000
-DEAL:3d::0.239142, 0.287892, 0.00000
+DEAL:3d::0.292475, 0.162733, 0.00000
+DEAL:3d::0.306667, 0.202500, 0.0512500
+DEAL:3d::0.292475, 0.242267, 0.00000
+DEAL:3d::0.239142, 0.294350, 0.00000
 DEAL:3d::0.278284, 0.282034, 0.0512500
-DEAL:3d::0.200000, 0.293750, 0.0512500
-DEAL:3d::0.160858, 0.287892, 0.00000
+DEAL:3d::0.200000, 0.306667, 0.0512500
+DEAL:3d::0.160858, 0.294350, 0.00000
 DEAL:3d::0.121716, 0.282034, 0.0512500
-DEAL:3d::0.115858, 0.242267, 0.00000
-DEAL:3d::0.110000, 0.202500, 0.0512500
-DEAL:3d::0.115858, 0.162733, 0.00000
+DEAL:3d::0.107525, 0.242267, 0.00000
+DEAL:3d::0.0933333, 0.202500, 0.0512500
+DEAL:3d::0.107525, 0.162733, 0.00000
 DEAL:3d::0.121716, 0.122966, 0.0512500
-DEAL:3d::0.160858, 0.117108, 0.00000
-DEAL:3d::0.200000, 0.111250, 0.0512500
-DEAL:3d::0.239142, 0.117108, 0.00000
+DEAL:3d::0.160858, 0.108150, 0.00000
+DEAL:3d::0.200000, 0.0933333, 0.0512500
+DEAL:3d::0.239142, 0.108150, 0.00000
 DEAL:3d::0.278284, 0.122966, 0.0512500
-DEAL:3d::0.284142, 0.162733, 0.102500
-DEAL:3d::0.290000, 0.202500, 0.153750
-DEAL:3d::0.284142, 0.242267, 0.102500
-DEAL:3d::0.239142, 0.287892, 0.102500
+DEAL:3d::0.292475, 0.162733, 0.102500
+DEAL:3d::0.306667, 0.202500, 0.153750
+DEAL:3d::0.292475, 0.242267, 0.102500
+DEAL:3d::0.239142, 0.294350, 0.102500
 DEAL:3d::0.278284, 0.282034, 0.153750
-DEAL:3d::0.200000, 0.293750, 0.153750
-DEAL:3d::0.160858, 0.287892, 0.102500
+DEAL:3d::0.200000, 0.306667, 0.153750
+DEAL:3d::0.160858, 0.294350, 0.102500
 DEAL:3d::0.121716, 0.282034, 0.153750
-DEAL:3d::0.115858, 0.242267, 0.102500
-DEAL:3d::0.110000, 0.202500, 0.153750
-DEAL:3d::0.115858, 0.162733, 0.102500
+DEAL:3d::0.107525, 0.242267, 0.102500
+DEAL:3d::0.0933333, 0.202500, 0.153750
+DEAL:3d::0.107525, 0.162733, 0.102500
 DEAL:3d::0.121716, 0.122966, 0.153750
-DEAL:3d::0.160858, 0.117108, 0.102500
-DEAL:3d::0.200000, 0.111250, 0.153750
-DEAL:3d::0.239142, 0.117108, 0.102500
+DEAL:3d::0.160858, 0.108150, 0.102500
+DEAL:3d::0.200000, 0.0933333, 0.153750
+DEAL:3d::0.239142, 0.108150, 0.102500
 DEAL:3d::0.278284, 0.122966, 0.153750
-DEAL:3d::0.284142, 0.162733, 0.205000
-DEAL:3d::0.290000, 0.202500, 0.256250
-DEAL:3d::0.284142, 0.242267, 0.205000
-DEAL:3d::0.239142, 0.287892, 0.205000
+DEAL:3d::0.292475, 0.162733, 0.205000
+DEAL:3d::0.306667, 0.202500, 0.256250
+DEAL:3d::0.292475, 0.242267, 0.205000
+DEAL:3d::0.239142, 0.294350, 0.205000
 DEAL:3d::0.278284, 0.282034, 0.256250
-DEAL:3d::0.200000, 0.293750, 0.256250
-DEAL:3d::0.160858, 0.287892, 0.205000
+DEAL:3d::0.200000, 0.306667, 0.256250
+DEAL:3d::0.160858, 0.294350, 0.205000
 DEAL:3d::0.121716, 0.282034, 0.256250
-DEAL:3d::0.115858, 0.242267, 0.205000
-DEAL:3d::0.110000, 0.202500, 0.256250
-DEAL:3d::0.115858, 0.162733, 0.205000
+DEAL:3d::0.107525, 0.242267, 0.205000
+DEAL:3d::0.0933333, 0.202500, 0.256250
+DEAL:3d::0.107525, 0.162733, 0.205000
 DEAL:3d::0.121716, 0.122966, 0.256250
-DEAL:3d::0.160858, 0.117108, 0.205000
-DEAL:3d::0.200000, 0.111250, 0.256250
-DEAL:3d::0.239142, 0.117108, 0.205000
+DEAL:3d::0.160858, 0.108150, 0.205000
+DEAL:3d::0.200000, 0.0933333, 0.256250
+DEAL:3d::0.239142, 0.108150, 0.205000
 DEAL:3d::0.278284, 0.122966, 0.256250
-DEAL:3d::0.284142, 0.162733, 0.307500
-DEAL:3d::0.290000, 0.202500, 0.358750
-DEAL:3d::0.284142, 0.242267, 0.307500
-DEAL:3d::0.239142, 0.287892, 0.307500
+DEAL:3d::0.292475, 0.162733, 0.307500
+DEAL:3d::0.306667, 0.202500, 0.358750
+DEAL:3d::0.292475, 0.242267, 0.307500
+DEAL:3d::0.239142, 0.294350, 0.307500
 DEAL:3d::0.278284, 0.282034, 0.358750
-DEAL:3d::0.200000, 0.293750, 0.358750
-DEAL:3d::0.160858, 0.287892, 0.307500
+DEAL:3d::0.200000, 0.306667, 0.358750
+DEAL:3d::0.160858, 0.294350, 0.307500
 DEAL:3d::0.121716, 0.282034, 0.358750
-DEAL:3d::0.115858, 0.242267, 0.307500
-DEAL:3d::0.110000, 0.202500, 0.358750
-DEAL:3d::0.115858, 0.162733, 0.307500
+DEAL:3d::0.107525, 0.242267, 0.307500
+DEAL:3d::0.0933333, 0.202500, 0.358750
+DEAL:3d::0.107525, 0.162733, 0.307500
 DEAL:3d::0.121716, 0.122966, 0.358750
-DEAL:3d::0.160858, 0.117108, 0.307500
-DEAL:3d::0.200000, 0.111250, 0.358750
-DEAL:3d::0.239142, 0.117108, 0.307500
+DEAL:3d::0.160858, 0.108150, 0.307500
+DEAL:3d::0.200000, 0.0933333, 0.358750
+DEAL:3d::0.239142, 0.108150, 0.307500
 DEAL:3d::0.278284, 0.122966, 0.358750
-DEAL:3d::0.284142, 0.162733, 0.410000
-DEAL:3d::0.284142, 0.242267, 0.410000
-DEAL:3d::0.239142, 0.287892, 0.410000
-DEAL:3d::0.160858, 0.287892, 0.410000
-DEAL:3d::0.115858, 0.242267, 0.410000
-DEAL:3d::0.115858, 0.162733, 0.410000
-DEAL:3d::0.160858, 0.117108, 0.410000
-DEAL:3d::0.239142, 0.117108, 0.410000
+DEAL:3d::0.292475, 0.162733, 0.410000
+DEAL:3d::0.292475, 0.242267, 0.410000
+DEAL:3d::0.239142, 0.294350, 0.410000
+DEAL:3d::0.160858, 0.294350, 0.410000
+DEAL:3d::0.107525, 0.242267, 0.410000
+DEAL:3d::0.107525, 0.162733, 0.410000
+DEAL:3d::0.160858, 0.108150, 0.410000
+DEAL:3d::0.239142, 0.108150, 0.410000
 DEAL:3d::boundary id: 0
 DEAL:3d::0.00000, 0.0512500, 0.0512500
 DEAL:3d::0.00000, 0.153750, 0.0512500
@@ -760,8 +760,8 @@ DEAL:3d::1.95000, 0.00000, 0.0512500
 DEAL:3d::2.05000, 0.00000, 0.0512500
 DEAL:3d::2.15000, 0.00000, 0.0512500
 DEAL:3d::0.0500000, 0.0512500, 0.00000
-DEAL:3d::0.150000, 0.0512500, 0.00000
-DEAL:3d::0.250000, 0.0512500, 0.00000
+DEAL:3d::0.150000, 0.0422917, 0.00000
+DEAL:3d::0.250000, 0.0422917, 0.00000
 DEAL:3d::0.350000, 0.0512500, 0.00000
 DEAL:3d::0.450000, 0.0512500, 0.00000
 DEAL:3d::0.550000, 0.0512500, 0.00000
@@ -781,8 +781,8 @@ DEAL:3d::1.85000, 0.0512500, 0.00000
 DEAL:3d::1.95000, 0.0512500, 0.00000
 DEAL:3d::2.05000, 0.0512500, 0.00000
 DEAL:3d::2.15000, 0.0512500, 0.00000
-DEAL:3d::0.0500000, 0.153750, 0.00000
-DEAL:3d::0.350000, 0.153750, 0.00000
+DEAL:3d::0.0416667, 0.153750, 0.00000
+DEAL:3d::0.358333, 0.153750, 0.00000
 DEAL:3d::0.450000, 0.153750, 0.00000
 DEAL:3d::0.550000, 0.153750, 0.00000
 DEAL:3d::0.650000, 0.153750, 0.00000
@@ -801,8 +801,8 @@ DEAL:3d::1.85000, 0.153750, 0.00000
 DEAL:3d::1.95000, 0.153750, 0.00000
 DEAL:3d::2.05000, 0.153750, 0.00000
 DEAL:3d::2.15000, 0.153750, 0.00000
-DEAL:3d::0.0500000, 0.256250, 0.00000
-DEAL:3d::0.350000, 0.256250, 0.00000
+DEAL:3d::0.0416667, 0.256250, 0.00000
+DEAL:3d::0.358333, 0.256250, 0.00000
 DEAL:3d::0.450000, 0.256250, 0.00000
 DEAL:3d::0.550000, 0.256250, 0.00000
 DEAL:3d::0.650000, 0.256250, 0.00000
@@ -823,9 +823,9 @@ DEAL:3d::2.05000, 0.256250, 0.00000
 DEAL:3d::2.15000, 0.256250, 0.00000
 DEAL:3d::0.0500000, 0.358750, 0.00000
 DEAL:3d::0.0500000, 0.410000, 0.0512500
-DEAL:3d::0.150000, 0.358750, 0.00000
+DEAL:3d::0.150000, 0.365208, 0.00000
 DEAL:3d::0.150000, 0.410000, 0.0512500
-DEAL:3d::0.250000, 0.358750, 0.00000
+DEAL:3d::0.250000, 0.365208, 0.00000
 DEAL:3d::0.250000, 0.410000, 0.0512500
 DEAL:3d::0.350000, 0.358750, 0.00000
 DEAL:3d::0.350000, 0.410000, 0.0512500
@@ -865,14 +865,14 @@ DEAL:3d::2.05000, 0.358750, 0.00000
 DEAL:3d::2.05000, 0.410000, 0.0512500
 DEAL:3d::2.15000, 0.358750, 0.00000
 DEAL:3d::2.15000, 0.410000, 0.0512500
-DEAL:3d::0.284142, 0.162733, 0.00000
-DEAL:3d::0.284142, 0.242267, 0.00000
-DEAL:3d::0.239142, 0.287892, 0.00000
-DEAL:3d::0.160858, 0.287892, 0.00000
-DEAL:3d::0.115858, 0.242267, 0.00000
-DEAL:3d::0.115858, 0.162733, 0.00000
-DEAL:3d::0.160858, 0.117108, 0.00000
-DEAL:3d::0.239142, 0.117108, 0.00000
+DEAL:3d::0.292475, 0.162733, 0.00000
+DEAL:3d::0.292475, 0.242267, 0.00000
+DEAL:3d::0.239142, 0.294350, 0.00000
+DEAL:3d::0.160858, 0.294350, 0.00000
+DEAL:3d::0.107525, 0.242267, 0.00000
+DEAL:3d::0.107525, 0.162733, 0.00000
+DEAL:3d::0.160858, 0.108150, 0.00000
+DEAL:3d::0.239142, 0.108150, 0.00000
 DEAL:3d::0.262147, 0.174258, 0.00000
 DEAL:3d::0.262147, 0.225742, 0.00000
 DEAL:3d::0.225742, 0.262147, 0.00000
@@ -1038,8 +1038,8 @@ DEAL:3d::1.95000, 0.410000, 0.358750
 DEAL:3d::2.05000, 0.410000, 0.358750
 DEAL:3d::2.15000, 0.410000, 0.358750
 DEAL:3d::0.0500000, 0.0512500, 0.410000
-DEAL:3d::0.150000, 0.0512500, 0.410000
-DEAL:3d::0.250000, 0.0512500, 0.410000
+DEAL:3d::0.150000, 0.0422917, 0.410000
+DEAL:3d::0.250000, 0.0422917, 0.410000
 DEAL:3d::0.350000, 0.0512500, 0.410000
 DEAL:3d::0.450000, 0.0512500, 0.410000
 DEAL:3d::0.550000, 0.0512500, 0.410000
@@ -1059,8 +1059,8 @@ DEAL:3d::1.85000, 0.0512500, 0.410000
 DEAL:3d::1.95000, 0.0512500, 0.410000
 DEAL:3d::2.05000, 0.0512500, 0.410000
 DEAL:3d::2.15000, 0.0512500, 0.410000
-DEAL:3d::0.0500000, 0.153750, 0.410000
-DEAL:3d::0.350000, 0.153750, 0.410000
+DEAL:3d::0.0416667, 0.153750, 0.410000
+DEAL:3d::0.358333, 0.153750, 0.410000
 DEAL:3d::0.450000, 0.153750, 0.410000
 DEAL:3d::0.550000, 0.153750, 0.410000
 DEAL:3d::0.650000, 0.153750, 0.410000
@@ -1079,8 +1079,8 @@ DEAL:3d::1.85000, 0.153750, 0.410000
 DEAL:3d::1.95000, 0.153750, 0.410000
 DEAL:3d::2.05000, 0.153750, 0.410000
 DEAL:3d::2.15000, 0.153750, 0.410000
-DEAL:3d::0.0500000, 0.256250, 0.410000
-DEAL:3d::0.350000, 0.256250, 0.410000
+DEAL:3d::0.0416667, 0.256250, 0.410000
+DEAL:3d::0.358333, 0.256250, 0.410000
 DEAL:3d::0.450000, 0.256250, 0.410000
 DEAL:3d::0.550000, 0.256250, 0.410000
 DEAL:3d::0.650000, 0.256250, 0.410000
@@ -1100,8 +1100,8 @@ DEAL:3d::1.95000, 0.256250, 0.410000
 DEAL:3d::2.05000, 0.256250, 0.410000
 DEAL:3d::2.15000, 0.256250, 0.410000
 DEAL:3d::0.0500000, 0.358750, 0.410000
-DEAL:3d::0.150000, 0.358750, 0.410000
-DEAL:3d::0.250000, 0.358750, 0.410000
+DEAL:3d::0.150000, 0.365208, 0.410000
+DEAL:3d::0.250000, 0.365208, 0.410000
 DEAL:3d::0.350000, 0.358750, 0.410000
 DEAL:3d::0.450000, 0.358750, 0.410000
 DEAL:3d::0.550000, 0.358750, 0.410000
@@ -1121,14 +1121,14 @@ DEAL:3d::1.85000, 0.358750, 0.410000
 DEAL:3d::1.95000, 0.358750, 0.410000
 DEAL:3d::2.05000, 0.358750, 0.410000
 DEAL:3d::2.15000, 0.358750, 0.410000
-DEAL:3d::0.284142, 0.162733, 0.410000
-DEAL:3d::0.284142, 0.242267, 0.410000
-DEAL:3d::0.239142, 0.287892, 0.410000
-DEAL:3d::0.160858, 0.287892, 0.410000
-DEAL:3d::0.115858, 0.242267, 0.410000
-DEAL:3d::0.115858, 0.162733, 0.410000
-DEAL:3d::0.160858, 0.117108, 0.410000
-DEAL:3d::0.239142, 0.117108, 0.410000
+DEAL:3d::0.292475, 0.162733, 0.410000
+DEAL:3d::0.292475, 0.242267, 0.410000
+DEAL:3d::0.239142, 0.294350, 0.410000
+DEAL:3d::0.160858, 0.294350, 0.410000
+DEAL:3d::0.107525, 0.242267, 0.410000
+DEAL:3d::0.107525, 0.162733, 0.410000
+DEAL:3d::0.160858, 0.108150, 0.410000
+DEAL:3d::0.239142, 0.108150, 0.410000
 DEAL:3d::0.262147, 0.174258, 0.410000
 DEAL:3d::0.262147, 0.225742, 0.410000
 DEAL:3d::0.225742, 0.262147, 0.410000

diff --git a/tests/grid/find_all_active_cells_around_point_03.output b/tests/grid/find_all_active_cells_around_point_03.output
index 8bd3be602e851e467f77b32d32a5fd8312488d36..5324afdc98a990e63f81a45032fe53fb1496652c 100644 (file)
--- a/tests/grid/find_all_active_cells_around_point_03.output
+++ b/tests/grid/find_all_active_cells_around_point_03.output
@@ -1,281 +1,197 @@
 
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance 0.001000000000
-DEAL::Cell: 23_1:1 unit point 1.000000000 0.09756097561
-DEAL::Cell: 107_1:2 unit point 0.09756097561 1.000000000
+DEAL::Cell: 107_1:0 unit point 0.04887876126 0.7619679173
 DEAL::
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance 0.0001000000000
-DEAL::Cell: 23_1:1 unit point 1.000000000 0.09756097561
-DEAL::Cell: 107_1:2 unit point 0.09756097561 1.000000000
+DEAL::Cell: 107_1:0 unit point 0.04887876126 0.7619679173
 DEAL::
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance 1.000000000e-05
-DEAL::Cell: 23_1:1 unit point 1.000000000 0.09756097561
-DEAL::Cell: 107_1:2 unit point 0.09756097561 1.000000000
+DEAL::Cell: 107_1:0 unit point 0.04887876126 0.7619679173
 DEAL::
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance 1.000000000e-06
-DEAL::Cell: 23_1:1 unit point 1.000000000 0.09756097561
-DEAL::Cell: 107_1:2 unit point 0.09756097561 1.000000000
+DEAL::Cell: 107_1:0 unit point 0.04887876126 0.7619679173
 DEAL::
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance 1.000000000e-07
-DEAL::Cell: 23_1:1 unit point 1.000000000 0.09756097561
-DEAL::Cell: 107_1:2 unit point 0.09756097561 1.000000000
+DEAL::Cell: 107_1:0 unit point 0.04887876126 0.7619679173
 DEAL::
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance 1.000000000e-10
-DEAL::Cell: 23_1:1 unit point 1.000000000 0.09756097561
-DEAL::Cell: 107_1:2 unit point 0.09756097561 1.000000000
+DEAL::Cell: 107_1:0 unit point 0.04887876126 0.7619679173
 DEAL::
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance default 
-DEAL::Cell: 23_1:1 unit point 1.000000000 0.09756097561
-DEAL::Cell: 107_1:2 unit point 0.09756097561 1.000000000
+DEAL::Cell: 107_1:0 unit point 0.04887876126 0.7619679173
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance 0.001000000000
-DEAL::Cell: 22_1:0 unit point 1.252385729e-16 0.09756097561
-DEAL::Cell: 104_1:1 unit point 1.000000000 0.09756097561
+DEAL::Cell: 104_1:0 unit point 0.7619679173 0.04887876126
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance 0.0001000000000
-DEAL::Cell: 22_1:0 unit point 1.252385729e-16 0.09756097561
-DEAL::Cell: 104_1:1 unit point 1.000000000 0.09756097561
+DEAL::Cell: 104_1:0 unit point 0.7619679173 0.04887876126
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance 1.000000000e-05
-DEAL::Cell: 22_1:0 unit point 1.252385729e-16 0.09756097561
-DEAL::Cell: 104_1:1 unit point 1.000000000 0.09756097561
+DEAL::Cell: 104_1:0 unit point 0.7619679173 0.04887876126
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance 1.000000000e-06
-DEAL::Cell: 22_1:0 unit point 1.252385729e-16 0.09756097561
-DEAL::Cell: 104_1:1 unit point 1.000000000 0.09756097561
+DEAL::Cell: 104_1:0 unit point 0.7619679173 0.04887876126
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance 1.000000000e-07
-DEAL::Cell: 22_1:0 unit point 1.252385729e-16 0.09756097561
-DEAL::Cell: 104_1:1 unit point 1.000000000 0.09756097561
+DEAL::Cell: 104_1:0 unit point 0.7619679173 0.04887876126
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance 1.000000000e-10
-DEAL::Cell: 22_1:0 unit point 1.252385729e-16 0.09756097561
-DEAL::Cell: 104_1:1 unit point 1.000000000 0.09756097561
+DEAL::Cell: 104_1:0 unit point 0.7619679173 0.04887876126
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance default 
-DEAL::Cell: 22_1:0 unit point 1.252385729e-16 0.09756097561
-DEAL::Cell: 104_1:1 unit point 1.000000000 0.09756097561
+DEAL::Cell: 104_1:0 unit point 0.7619679173 0.04887876126
 DEAL::
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance 0.001000000000
-DEAL::Cell: 23_2:11 unit point 1.000000000 0.1951219512
-DEAL::Cell: 107_2:22 unit point 0.1951219512 1.000000000
+DEAL::Cell: 107_2:02 unit point 0.09776156208 0.5239608435
 DEAL::
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance 0.0001000000000
-DEAL::Cell: 23_2:11 unit point 1.000000000 0.1951219512
-DEAL::Cell: 107_2:22 unit point 0.1951219512 1.000000000
+DEAL::Cell: 107_2:02 unit point 0.09776156208 0.5239608435
 DEAL::
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance 1.000000000e-05
-DEAL::Cell: 23_2:11 unit point 1.000000000 0.1951219512
-DEAL::Cell: 107_2:22 unit point 0.1951219512 1.000000000
+DEAL::Cell: 107_2:02 unit point 0.09776156208 0.5239608435
 DEAL::
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance 1.000000000e-06
-DEAL::Cell: 23_2:11 unit point 1.000000000 0.1951219512
-DEAL::Cell: 107_2:22 unit point 0.1951219512 1.000000000
+DEAL::Cell: 107_2:02 unit point 0.09776156208 0.5239608435
 DEAL::
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance 1.000000000e-07
-DEAL::Cell: 23_2:11 unit point 1.000000000 0.1951219512
-DEAL::Cell: 107_2:22 unit point 0.1951219512 1.000000000
+DEAL::Cell: 107_2:02 unit point 0.09776156208 0.5239608435
 DEAL::
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance 1.000000000e-10
-DEAL::Cell: 23_2:11 unit point 1.000000000 0.1951219512
-DEAL::Cell: 107_2:22 unit point 0.1951219512 1.000000000
+DEAL::Cell: 107_2:02 unit point 0.09776156208 0.5239608435
 DEAL::
 DEAL::Testing 2D with point 0.3000000000 0.2000000000 tolerance default 
-DEAL::Cell: 23_2:11 unit point 1.000000000 0.1951219512
-DEAL::Cell: 107_2:22 unit point 0.1951219512 1.000000000
+DEAL::Cell: 107_2:02 unit point 0.09776156208 0.5239608435
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance 0.001000000000
-DEAL::Cell: 22_2:00 unit point 2.775557562e-16 0.1951219512
-DEAL::Cell: 104_2:11 unit point 1.000000000 0.1951219512
+DEAL::Cell: 104_2:01 unit point 0.5239608435 0.09776156208
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance 0.0001000000000
-DEAL::Cell: 22_2:00 unit point 2.775557562e-16 0.1951219512
-DEAL::Cell: 104_2:11 unit point 1.000000000 0.1951219512
+DEAL::Cell: 104_2:01 unit point 0.5239608435 0.09776156208
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance 1.000000000e-05
-DEAL::Cell: 22_2:00 unit point 2.775557562e-16 0.1951219512
-DEAL::Cell: 104_2:11 unit point 1.000000000 0.1951219512
+DEAL::Cell: 104_2:01 unit point 0.5239608435 0.09776156208
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance 1.000000000e-06
-DEAL::Cell: 22_2:00 unit point 2.775557562e-16 0.1951219512
-DEAL::Cell: 104_2:11 unit point 1.000000000 0.1951219512
+DEAL::Cell: 104_2:01 unit point 0.5239608435 0.09776156208
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance 1.000000000e-07
-DEAL::Cell: 22_2:00 unit point 2.775557562e-16 0.1951219512
-DEAL::Cell: 104_2:11 unit point 1.000000000 0.1951219512
+DEAL::Cell: 104_2:01 unit point 0.5239608435 0.09776156208
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance 1.000000000e-10
-DEAL::Cell: 22_2:00 unit point 2.775557562e-16 0.1951219512
-DEAL::Cell: 104_2:11 unit point 1.000000000 0.1951219512
+DEAL::Cell: 104_2:01 unit point 0.5239608435 0.09776156208
 DEAL::
 DEAL::Testing 2D with point 0.1000000000 0.2000000000 tolerance default 
-DEAL::Cell: 22_2:00 unit point 2.775557562e-16 0.1951219512
-DEAL::Cell: 104_2:11 unit point 1.000000000 0.1951219512
+DEAL::Cell: 104_2:01 unit point 0.5239608435 0.09776156208
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance 0.001000000000
-DEAL::Cell: 93_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 94_0: unit point 1.000000000 0.04878048780 0.000000000
-DEAL::Cell: 429_0: unit point 0.04878048780 1.000000000 1.000000000
-DEAL::Cell: 430_0: unit point 0.04878048781 1.000000000 1.763805258e-14
+DEAL::Cell: 429_0: unit point 0.02442406476 0.3809274846 1.000000000
+DEAL::Cell: 430_0: unit point 0.02442406476 0.3809274846 7.929039891e-15
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance 0.0001000000000
-DEAL::Cell: 93_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 94_0: unit point 1.000000000 0.04878048780 0.000000000
-DEAL::Cell: 429_0: unit point 0.04878048780 1.000000000 1.000000000
-DEAL::Cell: 430_0: unit point 0.04878048781 1.000000000 1.763805258e-14
+DEAL::Cell: 429_0: unit point 0.02442406476 0.3809274846 1.000000000
+DEAL::Cell: 430_0: unit point 0.02442406476 0.3809274846 7.929039891e-15
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-05
-DEAL::Cell: 93_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 94_0: unit point 1.000000000 0.04878048780 0.000000000
-DEAL::Cell: 429_0: unit point 0.04878048780 1.000000000 1.000000000
-DEAL::Cell: 430_0: unit point 0.04878048781 1.000000000 1.763805258e-14
+DEAL::Cell: 429_0: unit point 0.02442406476 0.3809274846 1.000000000
+DEAL::Cell: 430_0: unit point 0.02442406476 0.3809274846 7.929039891e-15
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-06
-DEAL::Cell: 93_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 94_0: unit point 1.000000000 0.04878048780 0.000000000
-DEAL::Cell: 429_0: unit point 0.04878048780 1.000000000 1.000000000
-DEAL::Cell: 430_0: unit point 0.04878048781 1.000000000 1.763805258e-14
+DEAL::Cell: 429_0: unit point 0.02442406476 0.3809274846 1.000000000
+DEAL::Cell: 430_0: unit point 0.02442406476 0.3809274846 7.929039891e-15
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-07
-DEAL::Cell: 93_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 94_0: unit point 1.000000000 0.04878048780 0.000000000
-DEAL::Cell: 429_0: unit point 0.04878048780 1.000000000 1.000000000
-DEAL::Cell: 430_0: unit point 0.04878048781 1.000000000 1.763805258e-14
+DEAL::Cell: 429_0: unit point 0.02442406476 0.3809274846 1.000000000
+DEAL::Cell: 430_0: unit point 0.02442406476 0.3809274846 7.929039891e-15
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-10
-DEAL::Cell: 93_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 94_0: unit point 1.000000000 0.04878048780 0.000000000
-DEAL::Cell: 429_0: unit point 0.04878048780 1.000000000 1.000000000
-DEAL::Cell: 430_0: unit point 0.04878048781 1.000000000 1.763805258e-14
+DEAL::Cell: 429_0: unit point 0.02442406476 0.3809274846 1.000000000
+DEAL::Cell: 430_0: unit point 0.02442406476 0.3809274846 7.929039891e-15
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance default 
-DEAL::Cell: 93_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 94_0: unit point 1.000000000 0.04878048780 0.000000000
-DEAL::Cell: 429_0: unit point 0.04878048780 1.000000000 1.000000000
-DEAL::Cell: 430_0: unit point 0.04878048781 1.000000000 1.763805258e-14
+DEAL::Cell: 429_0: unit point 0.02442406476 0.3809274846 1.000000000
+DEAL::Cell: 430_0: unit point 0.02442406476 0.3809274846 7.929039891e-15
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance 0.001000000000
-DEAL::Cell: 89_0: unit point 0.000000000 0.04878048780 1.000000000
-DEAL::Cell: 90_0: unit point 0.000000000 0.04878048780 0.000000000
-DEAL::Cell: 417_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 418_0: unit point 1.000000000 0.04878048781 -5.916644075e-14
+DEAL::Cell: 417_0: unit point 0.3809274846 0.02442406476 1.000000000
+DEAL::Cell: 418_0: unit point 0.3809274846 0.02442406476 -3.000965623e-14
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance 0.0001000000000
-DEAL::Cell: 89_0: unit point 0.000000000 0.04878048780 1.000000000
-DEAL::Cell: 90_0: unit point 0.000000000 0.04878048780 0.000000000
-DEAL::Cell: 417_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 418_0: unit point 1.000000000 0.04878048781 -5.916644075e-14
+DEAL::Cell: 417_0: unit point 0.3809274846 0.02442406476 1.000000000
+DEAL::Cell: 418_0: unit point 0.3809274846 0.02442406476 -3.000965623e-14
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-05
-DEAL::Cell: 89_0: unit point 0.000000000 0.04878048780 1.000000000
-DEAL::Cell: 90_0: unit point 0.000000000 0.04878048780 0.000000000
-DEAL::Cell: 417_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 418_0: unit point 1.000000000 0.04878048781 -5.916644075e-14
+DEAL::Cell: 417_0: unit point 0.3809274846 0.02442406476 1.000000000
+DEAL::Cell: 418_0: unit point 0.3809274846 0.02442406476 -3.000965623e-14
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-06
-DEAL::Cell: 89_0: unit point 0.000000000 0.04878048780 1.000000000
-DEAL::Cell: 90_0: unit point 0.000000000 0.04878048780 0.000000000
-DEAL::Cell: 417_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 418_0: unit point 1.000000000 0.04878048781 -5.916644075e-14
+DEAL::Cell: 417_0: unit point 0.3809274846 0.02442406476 1.000000000
+DEAL::Cell: 418_0: unit point 0.3809274846 0.02442406476 -3.000965623e-14
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-07
-DEAL::Cell: 89_0: unit point 0.000000000 0.04878048780 1.000000000
-DEAL::Cell: 90_0: unit point 0.000000000 0.04878048780 0.000000000
-DEAL::Cell: 417_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 418_0: unit point 1.000000000 0.04878048781 -5.916644075e-14
+DEAL::Cell: 417_0: unit point 0.3809274846 0.02442406476 1.000000000
+DEAL::Cell: 418_0: unit point 0.3809274846 0.02442406476 -3.000965623e-14
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-10
-DEAL::Cell: 89_0: unit point 0.000000000 0.04878048780 1.000000000
-DEAL::Cell: 90_0: unit point 0.000000000 0.04878048780 0.000000000
-DEAL::Cell: 417_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 418_0: unit point 1.000000000 0.04878048781 -5.916644075e-14
+DEAL::Cell: 417_0: unit point 0.3809274846 0.02442406476 1.000000000
+DEAL::Cell: 418_0: unit point 0.3809274846 0.02442406476 -3.000965623e-14
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance default 
-DEAL::Cell: 89_0: unit point 0.000000000 0.04878048780 1.000000000
-DEAL::Cell: 90_0: unit point 0.000000000 0.04878048780 0.000000000
-DEAL::Cell: 417_0: unit point 1.000000000 0.04878048780 1.000000000
-DEAL::Cell: 418_0: unit point 1.000000000 0.04878048781 -5.916644075e-14
+DEAL::Cell: 417_0: unit point 0.3809274846 0.02442406476 1.000000000
+DEAL::Cell: 418_0: unit point 0.3809274846 0.02442406476 -3.000965623e-14
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance 0.001000000000
-DEAL::Cell: 93_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 94_1:1 unit point 1.000000000 0.09756097561 0.000000000
-DEAL::Cell: 429_1:6 unit point 0.09756097561 1.000000000 1.000000000
-DEAL::Cell: 430_1:2 unit point 0.09756097561 1.000000000 5.807787345e-14
+DEAL::Cell: 429_1:4 unit point 0.04887876126 0.7619679173 1.000000000
+DEAL::Cell: 430_1:0 unit point 0.04887876126 0.7619679173 2.509646622e-15
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance 0.0001000000000
-DEAL::Cell: 93_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 94_1:1 unit point 1.000000000 0.09756097561 0.000000000
-DEAL::Cell: 429_1:6 unit point 0.09756097561 1.000000000 1.000000000
-DEAL::Cell: 430_1:2 unit point 0.09756097561 1.000000000 5.807787345e-14
+DEAL::Cell: 429_1:4 unit point 0.04887876126 0.7619679173 1.000000000
+DEAL::Cell: 430_1:0 unit point 0.04887876126 0.7619679173 2.509646622e-15
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-05
-DEAL::Cell: 93_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 94_1:1 unit point 1.000000000 0.09756097561 0.000000000
-DEAL::Cell: 429_1:6 unit point 0.09756097561 1.000000000 1.000000000
-DEAL::Cell: 430_1:2 unit point 0.09756097561 1.000000000 5.807787345e-14
+DEAL::Cell: 429_1:4 unit point 0.04887876126 0.7619679173 1.000000000
+DEAL::Cell: 430_1:0 unit point 0.04887876126 0.7619679173 2.509646622e-15
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-06
-DEAL::Cell: 93_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 94_1:1 unit point 1.000000000 0.09756097561 0.000000000
-DEAL::Cell: 429_1:6 unit point 0.09756097561 1.000000000 1.000000000
-DEAL::Cell: 430_1:2 unit point 0.09756097561 1.000000000 5.807787345e-14
+DEAL::Cell: 429_1:4 unit point 0.04887876126 0.7619679173 1.000000000
+DEAL::Cell: 430_1:0 unit point 0.04887876126 0.7619679173 2.509646622e-15
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-07
-DEAL::Cell: 93_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 94_1:1 unit point 1.000000000 0.09756097561 0.000000000
-DEAL::Cell: 429_1:6 unit point 0.09756097561 1.000000000 1.000000000
-DEAL::Cell: 430_1:2 unit point 0.09756097561 1.000000000 5.807787345e-14
+DEAL::Cell: 429_1:4 unit point 0.04887876126 0.7619679173 1.000000000
+DEAL::Cell: 430_1:0 unit point 0.04887876126 0.7619679173 2.509646622e-15
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-10
-DEAL::Cell: 93_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 94_1:1 unit point 1.000000000 0.09756097561 0.000000000
-DEAL::Cell: 429_1:6 unit point 0.09756097561 1.000000000 1.000000000
-DEAL::Cell: 430_1:2 unit point 0.09756097561 1.000000000 5.807787345e-14
+DEAL::Cell: 429_1:4 unit point 0.04887876126 0.7619679173 1.000000000
+DEAL::Cell: 430_1:0 unit point 0.04887876126 0.7619679173 2.509646622e-15
 DEAL::
 DEAL::Testing 3D with point 0.3000000000 0.2000000000 0.2050000000 tolerance default 
-DEAL::Cell: 93_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 94_1:1 unit point 1.000000000 0.09756097561 0.000000000
-DEAL::Cell: 429_1:6 unit point 0.09756097561 1.000000000 1.000000000
-DEAL::Cell: 430_1:2 unit point 0.09756097561 1.000000000 5.807787345e-14
+DEAL::Cell: 429_1:4 unit point 0.04887876126 0.7619679173 1.000000000
+DEAL::Cell: 430_1:0 unit point 0.04887876126 0.7619679173 2.509646622e-15
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance 0.001000000000
-DEAL::Cell: 89_1:4 unit point 1.320082255e-16 0.09756097561 1.000000000
-DEAL::Cell: 90_1:0 unit point 0.000000000 0.09756097561 0.000000000
-DEAL::Cell: 417_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 418_1:1 unit point 1.000000000 0.09756097561 -6.250888646e-14
+DEAL::Cell: 417_1:4 unit point 0.7619679173 0.04887876126 1.000000000
+DEAL::Cell: 418_1:0 unit point 0.7619679173 0.04887876126 2.674987782e-15
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance 0.0001000000000
-DEAL::Cell: 89_1:4 unit point 1.320082255e-16 0.09756097561 1.000000000
-DEAL::Cell: 90_1:0 unit point 0.000000000 0.09756097561 0.000000000
-DEAL::Cell: 417_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 418_1:1 unit point 1.000000000 0.09756097561 -6.250888646e-14
+DEAL::Cell: 417_1:4 unit point 0.7619679173 0.04887876126 1.000000000
+DEAL::Cell: 418_1:0 unit point 0.7619679173 0.04887876126 2.674987782e-15
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-05
-DEAL::Cell: 89_1:4 unit point 1.320082255e-16 0.09756097561 1.000000000
-DEAL::Cell: 90_1:0 unit point 0.000000000 0.09756097561 0.000000000
-DEAL::Cell: 417_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 418_1:1 unit point 1.000000000 0.09756097561 -6.250888646e-14
+DEAL::Cell: 417_1:4 unit point 0.7619679173 0.04887876126 1.000000000
+DEAL::Cell: 418_1:0 unit point 0.7619679173 0.04887876126 2.674987782e-15
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-06
-DEAL::Cell: 89_1:4 unit point 1.320082255e-16 0.09756097561 1.000000000
-DEAL::Cell: 90_1:0 unit point 0.000000000 0.09756097561 0.000000000
-DEAL::Cell: 417_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 418_1:1 unit point 1.000000000 0.09756097561 -6.250888646e-14
+DEAL::Cell: 417_1:4 unit point 0.7619679173 0.04887876126 1.000000000
+DEAL::Cell: 418_1:0 unit point 0.7619679173 0.04887876126 2.674987782e-15
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-07
-DEAL::Cell: 89_1:4 unit point 1.320082255e-16 0.09756097561 1.000000000
-DEAL::Cell: 90_1:0 unit point 0.000000000 0.09756097561 0.000000000
-DEAL::Cell: 417_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 418_1:1 unit point 1.000000000 0.09756097561 -6.250888646e-14
+DEAL::Cell: 417_1:4 unit point 0.7619679173 0.04887876126 1.000000000
+DEAL::Cell: 418_1:0 unit point 0.7619679173 0.04887876126 2.674987782e-15
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance 1.000000000e-10
-DEAL::Cell: 89_1:4 unit point 1.320082255e-16 0.09756097561 1.000000000
-DEAL::Cell: 90_1:0 unit point 0.000000000 0.09756097561 0.000000000
-DEAL::Cell: 417_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 418_1:1 unit point 1.000000000 0.09756097561 -6.250888646e-14
+DEAL::Cell: 417_1:4 unit point 0.7619679173 0.04887876126 1.000000000
+DEAL::Cell: 418_1:0 unit point 0.7619679173 0.04887876126 2.674987782e-15
 DEAL::
 DEAL::Testing 3D with point 0.1000000000 0.2000000000 0.2050000000 tolerance default 
-DEAL::Cell: 89_1:4 unit point 1.320082255e-16 0.09756097561 1.000000000
-DEAL::Cell: 90_1:0 unit point 0.000000000 0.09756097561 0.000000000
-DEAL::Cell: 417_1:5 unit point 1.000000000 0.09756097561 1.000000000
-DEAL::Cell: 418_1:1 unit point 1.000000000 0.09756097561 -6.250888646e-14
+DEAL::Cell: 417_1:4 unit point 0.7619679173 0.04887876126 1.000000000
+DEAL::Cell: 418_1:0 unit point 0.7619679173 0.04887876126 2.674987782e-15
 DEAL::