From 3e837f72b08eea03857b7344986d8c4ac97d6c54 Mon Sep 17 00:00:00 2001
From: Christophe Geuzaine <cgeuzaine@ulg.ac.be>
Date: Mon, 12 Oct 2015 18:10:22 +0000
Subject: [PATCH] crack all elements from all entities except the crack
entities (instead of just cracking elements from entities of dim+1)
---
Plugin/Crack.cpp | 7 ++++++-
1 file changed, 6 insertions(+), 1 deletion(-)
diff --git a/Plugin/Crack.cpp b/Plugin/Crack.cpp
index a8d76a4b62..403087080e 100644
--- a/Plugin/Crack.cpp
+++ b/Plugin/Crack.cpp
@@ -111,6 +111,10 @@ PView *GMSH_CrackPlugin::execute(PView *view)
}
}
+ std::set<GEntity*> crackEntities;
+ crackEntities.insert(entities.begin(), entities.end());
+ crackEntities.insert(openEntities.begin(), openEntities.end());
+
// get crack elements
std::vector<MElement*> crackElements;
for(unsigned int i = 0; i < entities.size(); i++)
@@ -196,7 +200,8 @@ PView *GMSH_CrackPlugin::execute(PView *view)
std::vector<GEntity*> allentities;
m->getEntities(allentities);
for(unsigned int ent = 0; ent < allentities.size(); ent++){
- if(allentities[ent]->dim() != dim + 1) continue;
+ //if(allentities[ent]->dim() != dim + 1) continue;
+ if(crackEntities.find(allentities[ent]) != crackEntities.end()) continue;
for(unsigned int i = 0; i < allentities[ent]->getNumMeshElements(); i++){
MElement *e = allentities[ent]->getMeshElement(i);
for(int j = 0; j < e->getNumVertices(); j++){
--
GitLab