Don't overallocate memory.

diff --git a/source/grid/grid_tools_cache.cc b/source/grid/grid_tools_cache.cc
index 4caa9a671a308efd33aa47277f7a78cf3cf39d30..d6e873cad3d3e29d608e261eebef0b4746fc152a 100644 (file)
--- a/source/grid/grid_tools_cache.cc
+++ b/source/grid/grid_tools_cache.cc
@@ -155,7 +155,12 @@ namespace GridTools
           BoundingBox<spacedim>,
           typename Triangulation<dim, spacedim>::active_cell_iterator>>
           boxes;
-        boxes.reserve(tria->n_active_cells());
+        if (const parallel::TriangulationBase<dim, spacedim> *parallel_tria =
+              dynamic_cast<const parallel::TriangulationBase<dim, spacedim> *>(
+                &*tria))
+          boxes.reserve(parallel_tria->n_locally_owned_active_cells());
+        else
+          boxes.reserve(tria->n_active_cells());
         for (const auto &cell : tria->active_cell_iterators() |
                                   IteratorFilters::LocallyOwnedCell())
           boxes.emplace_back(mapping->get_bounding_box(cell), cell);