]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Add test case 9094/head
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Mon, 25 Nov 2019 18:54:07 +0000 (19:54 +0100)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Mon, 25 Nov 2019 18:58:16 +0000 (19:58 +0100)
tests/grid/find_all_active_cells_around_point_03.cc [new file with mode: 0644]
tests/grid/find_all_active_cells_around_point_03.output [new file with mode: 0644]

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 (file)
index 0000000..aa7b561
--- /dev/null
@@ -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 <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);
+}
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 (file)
index 0000000..8bd3be6
--- /dev/null
@@ -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::

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.