Skip to content
Snippets Groups Projects
Commit 296d52dd authored by Boris Sedji's avatar Boris Sedji
Browse files

MSH_TRI_B elementCutMesh bug fixed ?

parent 5d079489
No related branches found
No related tags found
No related merge requests found
......@@ -891,7 +891,7 @@ static void elementCutMesh(MElement *e, std::vector<const gLevelset *> &RPN,
case MSH_POLYG_ :
case MSH_POLYG_B :
{
if(eType == MSH_TRI_3) {
if((eType == MSH_TRI_3) | (eType == MSH_TRI_B)) {
DI_Triangle T(e->getVertex(0)->x(), e->getVertex(0)->y(), e->getVertex(0)->z(),
e->getVertex(1)->x(), e->getVertex(1)->y(), e->getVertex(1)->z(),
e->getVertex(2)->x(), e->getVertex(2)->y(), e->getVertex(2)->z());
......@@ -998,7 +998,7 @@ static void elementCutMesh(MElement *e, std::vector<const gLevelset *> &RPN,
int reg = getElementaryTag(triangles[nbTr]->lsTag(), elementary, newElemTags[2]);
std::vector<int> phys;
getPhysicalTag(triangles[nbTr]->lsTag(), gePhysicals, phys, newPhysTags[2]);
if(eType == MSH_TRI_3)
if((eType == MSH_TRI_3) | (eType == MSH_TRI_B))
elements[2][reg].push_back(copy);
else if(eType == MSH_QUA_4)
elements[3][reg].push_back(copy);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment