diff --git a/srcs/BEM/boundaryElements.cpp b/srcs/BEM/boundaryElements.cpp
index 54dd80afb64e5bfd426d9f1e85d4fd627bdfc525..7ca2ff14515c146eae9242da51cb8d4a073fa9b5 100644
--- a/srcs/BEM/boundaryElements.cpp
+++ b/srcs/BEM/boundaryElements.cpp
@@ -77,7 +77,7 @@ std::size_t fillElementVector(std::vector<elementStruct> &elementVector,
         // tags and nodes tags) in gmsh.
         std::vector<int> entityTags;
         gmsh::model::getEntitiesForPhysicalGroup(physicalGroupDim, physicalGroupTag, entityTags);
-        for (std::size_t i = 0; i < entityTags.size(); i++)
+        for (std::size_t entity = 0; entity < entityTags.size(); entity++)
         {
             std::vector<int> elementTypes = {1, 8}; // Hardcoded types of 2/3-nodes 1D elements.
             for(std::size_t j = 0; j < elementTypes.size(); j++)
@@ -89,14 +89,14 @@ std::size_t fillElementVector(std::vector<elementStruct> &elementVector,
                 // Get the element tags and the node tags of each element.
                 std::vector<std::size_t> elementTags;
                 std::vector<std::size_t> elementNodeTags;
-                gmsh::model::mesh::getElementsByType(elementTypes[j], elementTags, elementNodeTags, entityTags[i]);
+                gmsh::model::mesh::getElementsByType(elementTypes[j], elementTags, elementNodeTags, entityTags[entity]);
 
                 // Get the coordinates of each node (extracted above).
                 std::vector<std::size_t> nodeTags;
                 std::vector<double> nodeCoords;
                 std::vector<double> parametricNodeCoords;
                 gmsh::model::mesh::getNodesByElementType(elementTypes[j], nodeTags, nodeCoords, parametricNodeCoords,
-                                                         entityTags[i], false);
+                                                         entityTags[entity], false);
 
                 // Constructs a map associating the index of each node tag in the "nodeCoords" vector.
                 std::map<std::size_t, std::size_t> nodeIndices;