Skip to content
Snippets Groups Projects

Update boundaryElements.cpp

Open Maltez Cavalheiro Kévin requested to merge K.Maltez-master-patch-18822 into master
1 file
+ 3
3
Compare changes
  • Side-by-side
  • Inline
@@ -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;
Loading