--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2019 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.
+//
+// ---------------------------------------------------------------------
+
+
+// find_all_active_cells_around_point around a cylinder
+
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/utilities.h>
+
+#include <deal.II/fe/mapping_q_generic.h>
+
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_tools.h>
+
+#include "../tests.h"
+
+
+
+template <int dim, int spacedim>
+void
+print_result(const Mapping<dim, spacedim> & mapping,
+ const Triangulation<dim, spacedim> &tria,
+ const Point<dim> p,
+ const double tolerance)
+{
+ deallog << "Testing " << dim << "D with point " << p << " tolerance "
+ << tolerance << std::endl;
+ auto c_p =
+ GridTools::find_all_active_cells_around_point(mapping, tria, p, tolerance);
+ for (auto i : c_p)
+ deallog << "Cell: " << i.first->id() << " unit point " << i.second
+ << std::endl;
+ deallog << std::endl;
+}
+
+
+
+template <int dim, int spacedim>
+void
+print_result(const Mapping<dim, spacedim> & mapping,
+ const Triangulation<dim, spacedim> &tria,
+ const Point<dim> p)
+{
+ deallog << "Testing " << dim << "D with point " << p << " tolerance default "
+ << std::endl;
+ auto c_p = GridTools::find_all_active_cells_around_point(mapping, tria, p);
+ for (auto i : c_p)
+ deallog << "Cell: " << i.first->id() << " unit point " << i.second
+ << std::endl;
+ deallog << std::endl;
+}
+
+
+
+template <int dim, int spacedim>
+void
+test(unsigned int n_ref)
+{
+ Triangulation<dim, spacedim> tria;
+ GridGenerator::channel_with_cylinder(tria, 0.03, 2, 2);
+ tria.refine_global(n_ref);
+ MappingQGeneric<dim, spacedim> mapping(3);
+
+ Point<dim> p1;
+ p1[0] = 0.3;
+ p1[1] = 0.2;
+ if (dim == 3)
+ p1[2] = 0.205;
+ print_result(mapping, tria, p1, 1e-3);
+ print_result(mapping, tria, p1, 1e-4);
+ print_result(mapping, tria, p1, 1e-5);
+ print_result(mapping, tria, p1, 1e-6);
+ print_result(mapping, tria, p1, 1e-7);
+ print_result(mapping, tria, p1, 1e-10);
+ print_result(mapping, tria, p1);
+ Point<dim> p2 = p1;
+ p2[0] = 0.1;
+ print_result(mapping, tria, p2, 1e-3);
+ print_result(mapping, tria, p2, 1e-4);
+ print_result(mapping, tria, p2, 1e-5);
+ print_result(mapping, tria, p2, 1e-6);
+ print_result(mapping, tria, p2, 1e-7);
+ print_result(mapping, tria, p2, 1e-10);
+ print_result(mapping, tria, p2);
+}
+
+
+
+int
+main()
+{
+ initlog();
+ deallog << std::setprecision(10);
+
+ test<2, 2>(1);
+ test<2, 2>(2);
+ test<3, 3>(0);
+ test<3, 3>(1);
+}
--- /dev/null
+
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::
+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::