From: Marc Fehling <mafehling.git@gmail.com>
Date: Mon, 31 May 2021 17:59:25 +0000 (-0600)
Subject: Bugfix step-75: Call `get_fe()` only on locally owned cells.
X-Git-Tag: v9.4.0-rc1~1299^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d0b22b0ad9209c3d64de091e604c0bc22775f9bf;p=dealii.git

Bugfix step-75: Call `get_fe()` only on locally owned cells.
---

diff --git a/examples/step-75/step-75.cc b/examples/step-75/step-75.cc
index 1c4f469165..b946d3455a 100644
--- a/examples/step-75/step-75.cc
+++ b/examples/step-75/step-75.cc
@@ -745,17 +745,20 @@ namespace Step75
             auto cell_other = dof_handler_fine.begin_active();
             for (auto &cell : dof_handler_coarse.active_cell_iterators())
               {
-                const unsigned int next_degree =
-                  MGTransferGlobalCoarseningTools::
-                    create_next_polynomial_coarsening_degree(
-                      cell_other->get_fe().degree, mg_data.transfer.p_sequence);
-                Assert(fe_index_for_degree.find(next_degree) !=
-                         fe_index_for_degree.end(),
-                       ExcMessage("Next polynomial degree in sequence "
-                                  "does not exist in FECollection."));
-
                 if (cell->is_locally_owned())
-                  cell->set_active_fe_index(fe_index_for_degree[next_degree]);
+                  {
+                    const unsigned int next_degree =
+                      MGTransferGlobalCoarseningTools::
+                        create_next_polynomial_coarsening_degree(
+                          cell_other->get_fe().degree,
+                          mg_data.transfer.p_sequence);
+                    Assert(fe_index_for_degree.find(next_degree) !=
+                             fe_index_for_degree.end(),
+                           ExcMessage("Next polynomial degree in sequence "
+                                      "does not exist in FECollection."));
+
+                    cell->set_active_fe_index(fe_index_for_degree[next_degree]);
+                  }
                 cell_other++;
               }
           }