From: Marc Fehling <mafehling.git@gmail.com>
Date: Fri, 30 Aug 2024 22:03:24 +0000 (-0600)
Subject: Fix pcout in step-75.
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F17641%2Fhead;p=dealii.git

Fix pcout in step-75.
---

diff --git a/examples/step-75/step-75.cc b/examples/step-75/step-75.cc
index 387219f135..520ea00bab 100644
--- a/examples/step-75/step-75.cc
+++ b/examples/step-75/step-75.cc
@@ -1138,6 +1138,11 @@ namespace Step75
   // which will then output all information. Output of local quantities is
   // limited to the first 8 processes to avoid cluttering the terminal.
   //
+  // On all other processes, the containers for the collected data remain empty.
+  // To ensure that we do not access invalid memory with the insertion operator
+  // (`<<`) on these processes, we need to check that the containers are not
+  // empty.
+  //
   // Furthermore, we would like to print the frequencies of the polynomial
   // degrees in the numerical discretization. Since this information is only
   // stored locally, we will count the finite elements on locally owned cells
@@ -1160,7 +1165,8 @@ namespace Step75
         Utilities::MPI::gather(mpi_communicator,
                                triangulation.n_locally_owned_active_cells());
       for (unsigned int i = 0; i < first_n_processes; ++i)
-        pcout << ' ' << n_active_cells_per_subdomain[i];
+        if (n_active_cells_per_subdomain.size() > 0)
+          pcout << ' ' << n_active_cells_per_subdomain[i];
       if (output_cropped)
         pcout << " ...";
       pcout << std::endl;
@@ -1175,7 +1181,8 @@ namespace Step75
         Utilities::MPI::gather(mpi_communicator,
                                dof_handler.n_locally_owned_dofs());
       for (unsigned int i = 0; i < first_n_processes; ++i)
-        pcout << ' ' << n_dofs_per_subdomain[i];
+        if (n_dofs_per_subdomain.size() > 0)
+          pcout << ' ' << n_dofs_per_subdomain[i];
       if (output_cropped)
         pcout << " ...";
       pcout << std::endl;
@@ -1192,7 +1199,8 @@ namespace Step75
             << std::endl
             << "     by partition:              ";
       for (unsigned int i = 0; i < first_n_processes; ++i)
-        pcout << ' ' << n_constraints_per_subdomain[i];
+        if (n_constraints_per_subdomain.size() > 0)
+          pcout << ' ' << n_constraints_per_subdomain[i];
       if (output_cropped)
         pcout << " ...";
       pcout << std::endl;