Commit c04036fe authored by Christophe Geuzaine's avatar Christophe Geuzaine

Merge remote-tracking branch 'origin/KeepTrackOfSteinerVerticesOfVolumes'

parents 9a7bf826 0d7b2715
Pipeline #2604 passed with stage
in 64 minutes and 47 seconds
......@@ -1356,8 +1356,19 @@ void insertVerticesInRegion(GRegion *gr, int maxVert, bool _classify)
// has been found
Msg::Info("Found region %d", myGRegion->tag());
for(std::list<MTet4 *>::iterator it2 = theRegion.begin();
it2 != theRegion.end(); ++it2)
it2 != theRegion.end(); ++it2){
(*it2)->setOnWhat(myGRegion);
// Make sure that Steiner points will end up in the right region
std::vector<MVertex *> vertices;
(*it2)->tet()->getVertices(vertices);
for(std::vector<MVertex *>::iterator itv = vertices.begin(); itv != vertices.end(); ++itv){
if ((*itv)->onWhat() != NULL && (*itv)->onWhat()->dim() == 3 && (*itv)->onWhat() != myGRegion){
myGRegion->addMeshVertex((*itv));
(*itv)->setEntity(myGRegion);
}
}
}
}
else {
// the tets are in the void
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment