ForkMaster/netgenplugin
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fix SIGSEGV of computing a quadratic sub-mesh

Browse Source
This commit is contained in:
eap 2018-07-19 15:58:11 +03:00
parent d195b60384
commit 671f63a4da
1 changed files with 11 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
src/NETGENPlugin/NETGENPlugin_Mesher.cxx
View File

@ -325,7 +325,7 @@ void NETGENPlugin_Mesher::SetParameters(const NETGENPlugin_Hypothesis* hyp)
aSObj->UnRegister(); aSObj->UnRegister();
} }
TopoDS_Shape S = smeshGen_i->GeomObjectToShape( aGeomObj.in() ); TopoDS_Shape S = smeshGen_i->GeomObjectToShape( aGeomObj.in() );
::SetLocalSize(S, val); ::SetLocalSize(S, val);
} }
} }
} }
@ -3310,9 +3310,19 @@ bool NETGENPlugin_Mesher::Compute()
FillSMesh( occgeo, *_ngMesh, initState, *_mesh, nodeVec, comment, &quadHelper ); FillSMesh( occgeo, *_ngMesh, initState, *_mesh, nodeVec, comment, &quadHelper );
if ( quadHelper.GetIsQuadratic() ) // remove free nodes if ( quadHelper.GetIsQuadratic() ) // remove free nodes
{
for ( size_t i = 0; i < nodeVec.size(); ++i ) for ( size_t i = 0; i < nodeVec.size(); ++i )
if ( nodeVec[i] && nodeVec[i]->NbInverseElements() == 0 ) if ( nodeVec[i] && nodeVec[i]->NbInverseElements() == 0 )
{
_mesh->GetMeshDS()->RemoveFreeNode( nodeVec[i], 0, /*fromGroups=*/false ); _mesh->GetMeshDS()->RemoveFreeNode( nodeVec[i], 0, /*fromGroups=*/false );
nodeVec[i]=0;
}
for ( size_t i = nodeVec.size()-1; i > 0; --i ) // remove trailing removed nodes
if ( !nodeVec[i] )
nodeVec.resize( i );
else
break;
}
} }
SMESH_ComputeErrorPtr readErr = ReadErrors(nodeVec); SMESH_ComputeErrorPtr readErr = ReadErrors(nodeVec);
if ( readErr && readErr->HasBadElems() ) if ( readErr && readErr->HasBadElems() )