From: marcfehling <marc.fehling@gmx.net>
Date: Thu, 4 Jul 2019 13:24:31 +0000 (+0200)
Subject: Communicate thresholds.
X-Git-Tag: v9.2.0-rc1~1407^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F8325%2Fhead;p=dealii.git

Communicate thresholds.
---

diff --git a/source/hp/refinement.cc b/source/hp/refinement.cc
index acb946f2c8..82eb819177 100644
--- a/source/hp/refinement.cc
+++ b/source/hp/refinement.cc
@@ -14,6 +14,10 @@
 // ---------------------------------------------------------------------
 
 
+#include <deal.II/base/mpi.h>
+
+#include <deal.II/distributed/tria_base.h>
+
 #include <deal.II/grid/grid_refinement.h>
 
 #include <deal.II/hp/dof_handler.h>
@@ -129,6 +133,24 @@ namespace hp
               }
           }
 
+      if (const parallel::Triangulation<dim, spacedim> *parallel_tria =
+            dynamic_cast<const parallel::Triangulation<dim, spacedim> *>(
+              &dof_handler.get_triangulation()))
+        {
+          max_smoothness_refine =
+            Utilities::MPI::max(max_smoothness_refine,
+                                parallel_tria->get_communicator());
+          min_smoothness_refine =
+            Utilities::MPI::min(min_smoothness_refine,
+                                parallel_tria->get_communicator());
+          max_smoothness_coarsen =
+            Utilities::MPI::max(max_smoothness_coarsen,
+                                parallel_tria->get_communicator());
+          min_smoothness_coarsen =
+            Utilities::MPI::min(min_smoothness_coarsen,
+                                parallel_tria->get_communicator());
+        }
+
       // Absent any better strategies, we will set the threshold by linear
       // interpolation for both classes of cells individually.
       const Number threshold_smoothness_refine =