diff --git a/NumLib/DOF/LocalToGlobalIndexMap.cpp b/NumLib/DOF/LocalToGlobalIndexMap.cpp
index 90a3de20b10b879cd41c2ab25a0cd4e8b1709f98..4fca442cdfeb884e8c211c6d4676ffe3b168f54e 100644
--- a/NumLib/DOF/LocalToGlobalIndexMap.cpp
+++ b/NumLib/DOF/LocalToGlobalIndexMap.cpp
@@ -44,6 +44,7 @@ LocalToGlobalIndexMap::findGlobalIndicesWithElementID(
     for (ElementIterator e = first; e != last; ++e)
     {
         LineIndex indices;
+        indices.reserve((*e)->getNumberOfNodes());
 
         for (auto* n = (*e)->getNodes();
              n < (*e)->getNodes()+(*e)->getNumberOfNodes(); ++n)
@@ -56,6 +57,7 @@ LocalToGlobalIndexMap::findGlobalIndicesWithElementID(
             indices.push_back(_mesh_component_map.getGlobalIndex(l, comp_id));
         }
 
+        indices.shrink_to_fit();
         _rows((*e)->getID(), comp_id_write) = std::move(indices);
     }
 }
@@ -77,6 +79,7 @@ void LocalToGlobalIndexMap::findGlobalIndices(
     for (ElementIterator e = first; e != last; ++e, ++elem_id)
     {
         LineIndex indices;
+        indices.reserve((*e)->getNumberOfNodes());
 
         for (auto* n = (*e)->getNodes();
              n < (*e)->getNodes() + (*e)->getNumberOfNodes(); ++n)
@@ -89,6 +92,7 @@ void LocalToGlobalIndexMap::findGlobalIndices(
             indices.push_back(_mesh_component_map.getGlobalIndex(l, comp_id));
         }
 
+        indices.shrink_to_fit();
         _rows(elem_id, comp_id_write) = std::move(indices);
     }
 }