From 4608e5fbcf6e48fa1e7dbda45107d4ad50927fea Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Mon, 25 Nov 2019 19:54:07 +0100 Subject: [PATCH] Add test case --- .../find_all_active_cells_around_point_03.cc | 111 +++++++ ...nd_all_active_cells_around_point_03.output | 281 ++++++++++++++++++ 2 files changed, 392 insertions(+) create mode 100644 tests/grid/find_all_active_cells_around_point_03.cc create mode 100644 tests/grid/find_all_active_cells_around_point_03.output diff --git a/tests/grid/find_all_active_cells_around_point_03.cc b/tests/grid/find_all_active_cells_around_point_03.cc new file mode 100644 index 0000000000..aa7b561daa --- /dev/null +++ b/tests/grid/find_all_active_cells_around_point_03.cc @@ -0,0 +1,111 @@ +// --------------------------------------------------------------------- +// +// 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 +#include + +#include + +#include +#include + +#include "../tests.h" + + + +template +void +print_result(const Mapping & mapping, + const Triangulation &tria, + const Point 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 +void +print_result(const Mapping & mapping, + const Triangulation &tria, + const Point 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 +void +test(unsigned int n_ref) +{ + Triangulation tria; + GridGenerator::channel_with_cylinder(tria, 0.03, 2, 2); + tria.refine_global(n_ref); + MappingQGeneric mapping(3); + + Point 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 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); +} diff --git a/tests/grid/find_all_active_cells_around_point_03.output b/tests/grid/find_all_active_cells_around_point_03.output new file mode 100644 index 0000000000..8bd3be602e --- /dev/null +++ b/tests/grid/find_all_active_cells_around_point_03.output @@ -0,0 +1,281 @@ + +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:: -- 2.39.5