if ( isOk ) // the non-poly elem remains valid after sticking nodes
{
- if ( nbNodes != nbUniqueNodes )
+ if ( nbNodes != nbUniqueNodes ||
+ !aMesh->ChangeElementNodes( elem, & curNodes[0], nbNodes ))
{
elemType.Init( elem ).SetID( elem->GetID() );
if ( elem != newElem )
ReplaceElemInGroups( elem, newElem, aMesh );
}
- else
- {
- aMesh->ChangeElementNodes( elem, & curNodes[0], nbNodes );
- }
}
else {
// Remove invalid regular element or invalid polygon
else
newFaces.Append( myLastCreatedElems(i) );
}
+ // get segments adjacent to merged nodes
+ TListOfListOfNodes::iterator groupIt = nodeGroupsToMerge.begin();
+ for ( ; groupIt != nodeGroupsToMerge.end(); groupIt++ )
+ {
+ const list<const SMDS_MeshNode*>& nodes = *groupIt;
+ SMDS_ElemIteratorPtr segIt = nodes.front()->GetInverseElementIterator( SMDSAbs_Edge );
+ while ( segIt->more() )
+ segments.insert( segIt->next() );
+ }
+
// find coincident
TListOfListOfElementsID equalGroups;
if ( !segments.empty() )