-// Copyright (C) 2007-2014 CEA/DEN, EDF R&D, OPEN CASCADE
+// Copyright (C) 2007-2016 CEA/DEN, EDF R&D, OPEN CASCADE
//
// Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
// CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
// Module : SMESH
//
#include "SMESH_subMesh_i.hxx"
+
+#include "SMESHDS_SubMesh.hxx"
#include "SMESH_Gen_i.hxx"
#include "SMESH_Mesh_i.hxx"
+#include "SMESH_MesherHelper.hxx"
#include "SMESH_PreMeshInfo.hxx"
#include "Utils_CorbaException.hxx"
#include "Utils_ExceptHandlers.hxx"
#include <TopoDS_Iterator.hxx>
-#include <TopExp_Explorer.hxx>
using namespace std;
SMESH_subMesh_i::~SMESH_subMesh_i()
{
- MESSAGE("SMESH_subMesh_i::~SMESH_subMesh_i");
if ( _preMeshInfo ) delete _preMeshInfo;
_preMeshInfo = NULL;
}
bool getSubMeshes(::SMESH_subMesh* theSubMesh,
TListOfSubMeshes& theSubMeshList)
{
- int size = theSubMeshList.size();
-
- SMESH_Mesh* aMesh = theSubMesh->GetFather();
- SMESHDS_Mesh* aMeshDS = aMesh->GetMeshDS();
- SMESHDS_SubMesh* aSubMeshDS = theSubMesh->GetSubMeshDS();
-
- // nodes can be bound to either vertex, edge, face or solid_or_shell
- TopoDS_Shape aShape = theSubMesh->GetSubShape();
- switch ( aShape.ShapeType() )
+ size_t size = theSubMeshList.size();
+
+ // check all child sub-meshes of one complexity,
+ // if no elements found and no algo assigned, go to children of lower complexity
+
+ TopoDS_Shape shape = theSubMesh->GetSubShape();
+ TopAbs_ShapeEnum mainType = SMESH_MesherHelper::GetGroupType( shape, /*noCompound=*/true );
+ TopAbs_ShapeEnum shapeType = shape.ShapeType();
+ bool elementsFound = false;
+ bool algoFound = false;
+ SMESH_subMeshIteratorPtr smIt = theSubMesh->getDependsOnIterator(/*includeSelf=*/true,
+ /*complexFirst=*/true);
+ while ( smIt->more() )
{
- case TopAbs_SOLID: {
- // add submesh of solid itself
- aSubMeshDS = aMeshDS->MeshElements( aShape );
- if ( aSubMeshDS )
- theSubMeshList.push_back( aSubMeshDS );
- // and of the first shell
- TopExp_Explorer exp( aShape, TopAbs_SHELL );
- if ( exp.More() ) {
- aSubMeshDS = aMeshDS->MeshElements( exp.Current() );
- if ( aSubMeshDS )
- theSubMeshList.push_back( aSubMeshDS );
- }
- break;
- }
- case TopAbs_WIRE:
- case TopAbs_COMPOUND:
- case TopAbs_COMPSOLID: {
- // call getSubMeshes() for sub-shapes
- list<TopoDS_Shape> shapeList;
- shapeList.push_back( aShape );
- list<TopoDS_Shape>::iterator sh = shapeList.begin();
- for ( ; sh != shapeList.end(); ++sh ) {
- for ( TopoDS_Iterator it( *sh ); it.More(); it.Next() ) {
- if ( ::SMESH_subMesh* aSubMesh = aMesh->GetSubMeshContaining( it.Value() ))
- getSubMeshes( aSubMesh, theSubMeshList ); // add found submesh or explore deeper
- else
- // no submesh for a compound inside compound
- shapeList.push_back( it.Value() );
- }
+ ::SMESH_subMesh* sm = smIt->next();
+ if ( sm->GetSubShape().ShapeType() != shapeType )
+ {
+ if ( elementsFound || algoFound )
+ break;
+ if ( sm->GetSubShape().ShapeType() == TopAbs_VERTEX &&
+ mainType != TopAbs_VERTEX )
+ break;
}
- // return only unique submeshes
- set<SMESHDS_SubMesh*> smSet;
- TListOfSubMeshes::iterator sm = theSubMeshList.begin();
- while ( sm != theSubMeshList.end() ) {
- if ( !smSet.insert( *sm ).second )
- sm = theSubMeshList.erase( sm );
- else
- ++sm;
+ shapeType = sm->GetSubShape().ShapeType();
+ if ( !sm->IsEmpty() )
+ {
+ elementsFound = true;
+ theSubMeshList.push_back( sm->GetSubMeshDS() );
}
- break;
- }
- default:
- if ( aSubMeshDS )
- theSubMeshList.push_back( aSubMeshDS );
+ if ( sm->GetAlgo() )
+ algoFound = true;
}
+
return size < theSubMeshList.size();
}
return 0;
::SMESH_subMesh* aSubMesh = _mesh_i->_mapSubMesh[_localId];
- SMESHDS_SubMesh* aSubMeshDS = aSubMesh->GetSubMeshDS();
- int nbElems = aSubMeshDS ? aSubMeshDS->NbElements() : 0;
+ int nbElems = 0;
- // volumes are bound to shell
TListOfSubMeshes smList;
- if ( nbElems == 0 && getSubMeshes( aSubMesh, smList ))
+ if ( getSubMeshes( aSubMesh, smList ))
{
TListOfSubMeshes::iterator sm = smList.begin();
for ( ; sm != smList.end(); ++sm )
_preMeshInfo->FullLoadFromFile();
::SMESH_subMesh* aSubMesh = _mesh_i->_mapSubMesh[_localId];
- SMESHDS_SubMesh* aSubMeshDS = aSubMesh->GetSubMeshDS();
- int nbElems = aSubMeshDS ? aSubMeshDS->NbElements() : 0;
+ int nbElems = 0;
TListOfSubMeshes smList;
- if ( nbElems )
- smList.push_back( aSubMeshDS );
-
- // volumes are bound to shell
- if ( nbElems == 0 && getSubMeshes( aSubMesh, smList ))
+ if ( getSubMeshes( aSubMesh, smList ))
{
TListOfSubMeshes::iterator sm = smList.begin();
for ( ; sm != smList.end(); ++sm )
_preMeshInfo->FullLoadFromFile();
::SMESH_subMesh* aSubMesh = _mesh_i->_mapSubMesh[_localId];
- SMESHDS_SubMesh* aSubMeshDS = aSubMesh->GetSubMeshDS();
- // PAL5440, return all nodes belonging to elements of submesh
+ // PAL5440, return all nodes belonging to elements of the sub-mesh
set<int> nodeIds;
- int nbElems = aSubMeshDS ? aSubMeshDS->NbElements() : 0;
+ int nbElems = 0;
// volumes may be bound to shell instead of solid
TListOfSubMeshes smList;
- if ( nbElems == 0 && getSubMeshes( aSubMesh, smList ))
+ if ( getSubMeshes( aSubMesh, smList ))
{
TListOfSubMeshes::iterator sm = smList.begin();
for ( ; sm != smList.end(); ++sm )
nbElems += (*sm)->NbElements();
}
}
- aSubMeshDS = 0;
- }
- else
- {
- if ( nbElems )
- smList.push_back( aSubMeshDS );
- }
-
- if ( theElemType == SMESH::NODE && aSubMeshDS )
- {
- SMDS_ElemIteratorPtr eIt = aSubMeshDS->GetElements();
- if ( eIt->more() ) {
- while ( eIt->more() ) {
- const SMDS_MeshElement* anElem = eIt->next();
- SMDS_ElemIteratorPtr nIt = anElem->nodesIterator();
- while ( nIt->more() )
- nodeIds.insert( nIt->next()->GetID() );
- }
- } else {
- SMDS_NodeIteratorPtr nIt = aSubMeshDS->GetNodes();
- while ( nIt->more() )
- nodeIds.insert( nIt->next()->GetID() );
- }
}
if ( theElemType == SMESH::NODE )
TListOfSubMeshes::iterator sm = smList.begin();
for ( i = 0; sm != smList.end(); sm++ )
{
- aSubMeshDS = *sm;
- SMDS_ElemIteratorPtr anIt = aSubMeshDS->GetElements();
+ SMDS_ElemIteratorPtr anIt = (*sm)->GetElements();
while ( i < n && anIt->more() ) {
const SMDS_MeshElement* anElem = anIt->next();
if ( theElemType == SMESH::ALL || anElem->GetType() == (SMDSAbs_ElementType)theElemType )