]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Extend test 12615/head
authorPeter Munch <peterrmuench@gmail.com>
Sat, 31 Jul 2021 07:58:45 +0000 (09:58 +0200)
committerPeter Munch <peterrmuench@gmail.com>
Sat, 31 Jul 2021 10:34:38 +0000 (12:34 +0200)
doc/news/changes/minor/20210731Munch [new file with mode: 0644]
tests/remote_point_evaluation/mapping_03.cc
tests/remote_point_evaluation/mapping_03.with_p4est=true.mpirun=2.output

diff --git a/doc/news/changes/minor/20210731Munch b/doc/news/changes/minor/20210731Munch
new file mode 100644 (file)
index 0000000..06a0be9
--- /dev/null
@@ -0,0 +1,4 @@
+Fixed: Use tolerances consistently in RemotePointEvaluation, 
+VectorTools::point_values(), and GridTools::distributed_compute_point_locations().
+<br>
+(Peter Munch, Magdalena Schreter, 2021/07/31)
index 54bcfc811b810845b26da5c41fb7606692591b63..3c5baeae2df0e153659536ca8f0148a9d70ddaa9 100644 (file)
@@ -40,7 +40,6 @@ test(const bool enforce_unique_map)
 {
   parallel::distributed::Triangulation<dim> tria(MPI_COMM_WORLD);
   GridGenerator::subdivided_hyper_cube(tria, 2);
-  tria.refine_global(1);
 
   FE_Q<dim>       fe(1);
   DoFHandler<dim> dof_handler(tria);
@@ -54,8 +53,9 @@ test(const bool enforce_unique_map)
 
   std::vector<Point<dim>> evaluation_points;
 
-  evaluation_points.emplace_back(0.5, 0.5 + 1e-7); // should be assigned to rank
-                                                   // 0 for unique mapping
+  for (int i = 1; i < 15; ++i)
+    // should be assigned to rank 0 for unique mapping
+    evaluation_points.emplace_back(0.5, 0.5 + std::pow<double>(10.0, -i));
 
   Utilities::MPI::RemotePointEvaluation<dim> eval(1e-6, enforce_unique_map);
 
@@ -74,8 +74,8 @@ test(const bool enforce_unique_map)
     eval, dof_handler, vec, VectorTools::EvaluationFlags::insert);
 
   for (unsigned int i = 0; i < evaluation_points.size(); ++i)
-    deallog << i << " " << result_avg[i] << " " << result_min[i] << " "
-            << result_max[i] << " " << result_insert[i] << " "
+    deallog << "1e-" << (i + 1) << " " << result_avg[i] << " " << result_min[i]
+            << " " << result_max[i] << " " << result_insert[i] << " "
             << eval.get_point_ptrs()[i + 1] - eval.get_point_ptrs()[i]
             << std::endl;
 }
index 9529db7f95dbc438963c5a44bfa50ecafbd24f0e..6fc910c6f3624d4ed8afe39476dbbfe43b80f4c6 100644 (file)
@@ -1,7 +1,59 @@
 
-DEAL:0::0 1.50000 1.00000 2.00000 1.00000 4
-DEAL:0::0 1.00000 1.00000 1.00000 1.00000 1
+DEAL:0::1e-1 2.00000 2.00000 2.00000 2.00000 2
+DEAL:0::1e-2 2.00000 2.00000 2.00000 2.00000 2
+DEAL:0::1e-3 2.00000 2.00000 2.00000 2.00000 2
+DEAL:0::1e-4 2.00000 2.00000 2.00000 2.00000 2
+DEAL:0::1e-5 2.00000 2.00000 2.00000 2.00000 2
+DEAL:0::1e-6 2.00000 2.00000 2.00000 2.00000 2
+DEAL:0::1e-7 1.50000 1.00000 2.00000 1.00000 4
+DEAL:0::1e-8 1.50000 1.00000 2.00000 1.00000 4
+DEAL:0::1e-9 1.50000 1.00000 2.00000 1.00000 4
+DEAL:0::1e-10 1.50000 1.00000 2.00000 1.00000 4
+DEAL:0::1e-11 1.50000 1.00000 2.00000 1.00000 4
+DEAL:0::1e-12 1.50000 1.00000 2.00000 1.00000 4
+DEAL:0::1e-13 1.50000 1.00000 2.00000 1.00000 4
+DEAL:0::1e-14 1.50000 1.00000 2.00000 1.00000 4
+DEAL:0::1e-1 2.00000 2.00000 2.00000 2.00000 1
+DEAL:0::1e-2 2.00000 2.00000 2.00000 2.00000 1
+DEAL:0::1e-3 2.00000 2.00000 2.00000 2.00000 1
+DEAL:0::1e-4 2.00000 2.00000 2.00000 2.00000 1
+DEAL:0::1e-5 2.00000 2.00000 2.00000 2.00000 1
+DEAL:0::1e-6 2.00000 2.00000 2.00000 2.00000 1
+DEAL:0::1e-7 1.00000 1.00000 1.00000 1.00000 1
+DEAL:0::1e-8 1.00000 1.00000 1.00000 1.00000 1
+DEAL:0::1e-9 1.00000 1.00000 1.00000 1.00000 1
+DEAL:0::1e-10 1.00000 1.00000 1.00000 1.00000 1
+DEAL:0::1e-11 1.00000 1.00000 1.00000 1.00000 1
+DEAL:0::1e-12 1.00000 1.00000 1.00000 1.00000 1
+DEAL:0::1e-13 1.00000 1.00000 1.00000 1.00000 1
+DEAL:0::1e-14 1.00000 1.00000 1.00000 1.00000 1
 
-DEAL:1::0 1.50000 1.00000 2.00000 1.00000 4
-DEAL:1::0 1.00000 1.00000 1.00000 1.00000 1
+DEAL:1::1e-1 2.00000 2.00000 2.00000 2.00000 2
+DEAL:1::1e-2 2.00000 2.00000 2.00000 2.00000 2
+DEAL:1::1e-3 2.00000 2.00000 2.00000 2.00000 2
+DEAL:1::1e-4 2.00000 2.00000 2.00000 2.00000 2
+DEAL:1::1e-5 2.00000 2.00000 2.00000 2.00000 2
+DEAL:1::1e-6 2.00000 2.00000 2.00000 2.00000 2
+DEAL:1::1e-7 1.50000 1.00000 2.00000 1.00000 4
+DEAL:1::1e-8 1.50000 1.00000 2.00000 1.00000 4
+DEAL:1::1e-9 1.50000 1.00000 2.00000 1.00000 4
+DEAL:1::1e-10 1.50000 1.00000 2.00000 1.00000 4
+DEAL:1::1e-11 1.50000 1.00000 2.00000 1.00000 4
+DEAL:1::1e-12 1.50000 1.00000 2.00000 1.00000 4
+DEAL:1::1e-13 1.50000 1.00000 2.00000 1.00000 4
+DEAL:1::1e-14 1.50000 1.00000 2.00000 1.00000 4
+DEAL:1::1e-1 2.00000 2.00000 2.00000 2.00000 1
+DEAL:1::1e-2 2.00000 2.00000 2.00000 2.00000 1
+DEAL:1::1e-3 2.00000 2.00000 2.00000 2.00000 1
+DEAL:1::1e-4 2.00000 2.00000 2.00000 2.00000 1
+DEAL:1::1e-5 2.00000 2.00000 2.00000 2.00000 1
+DEAL:1::1e-6 2.00000 2.00000 2.00000 2.00000 1
+DEAL:1::1e-7 1.00000 1.00000 1.00000 1.00000 1
+DEAL:1::1e-8 1.00000 1.00000 1.00000 1.00000 1
+DEAL:1::1e-9 1.00000 1.00000 1.00000 1.00000 1
+DEAL:1::1e-10 1.00000 1.00000 1.00000 1.00000 1
+DEAL:1::1e-11 1.00000 1.00000 1.00000 1.00000 1
+DEAL:1::1e-12 1.00000 1.00000 1.00000 1.00000 1
+DEAL:1::1e-13 1.00000 1.00000 1.00000 1.00000 1
+DEAL:1::1e-14 1.00000 1.00000 1.00000 1.00000 1
 

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.