Fix GroupCommunicator issue in scalarGrad3D (#288)

pecos · Jul 19, 2024 · c5a587a · c5a587a
1 parent 732a166
commit c5a587a
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/src/utils.cpp b/src/utils.cpp
@@ -905,8 +905,9 @@ void scalarGrad3D(ParGridFunction &u, ParGridFunction &gu) {
   gcomm.Bcast(zones_per_vdof);
 
   // Accumulate for all vdofs.
-  gcomm.Reduce<double>(gu.GetData(), GroupCommunicator::Sum);
-  gcomm.Bcast<double>(gu.GetData());
+  GroupCommunicator &gcomm_g = gu.ParFESpace()->GroupComm();
+  gcomm_g.Reduce<double>(gu.GetData(), GroupCommunicator::Sum);
+  gcomm_g.Bcast<double>(gu.GetData());
 
   // Compute means.
   for (int dir = 0; dir < dim; dir++) {