// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
//
-// See http://www.opencascade.org/SALOME/ or email : webmaster.salome@opencascade.org
+// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
//
//
//
// Module : SMESH
// $Header:
+#include "SMESHDS_SubMesh.hxx"
+
+#include "utilities.h"
+#include "SMDS_SetIterator.hxx"
+
using namespace std;
-#include "SMESHDS_SubMesh.ixx"
-#include "SMDS_MapIteratorOfExtendedMap.hxx"
//=======================================================================
-//function : SMESHDS_SubMesh
+//function : AddElement
//purpose :
//=======================================================================
-SMESHDS_SubMesh::SMESHDS_SubMesh(const Handle(SMDS_Mesh)& M) : myMesh(M)
+void SMESHDS_SubMesh::AddElement(const SMDS_MeshElement * ME)
{
- myListOfEltIDIsUpdate = Standard_False;
- myListOfNodeIDIsUpdate = Standard_False;
+ if ( !IsComplexSubmesh() )
+ myElements.insert(ME);
}
//=======================================================================
-//function : AddElement
+//function : RemoveElement
//purpose :
//=======================================================================
-void SMESHDS_SubMesh::AddElement (const Handle(SMDS_MeshElement)& ME)
+bool SMESHDS_SubMesh::RemoveElement(const SMDS_MeshElement * ME)
{
- myElements.Add(ME);
- myListOfEltIDIsUpdate = Standard_False;
+ if ( !IsComplexSubmesh() && NbElements() )
+ return myElements.erase(ME);
+
+ return false;
}
//=======================================================================
-//function : RemoveElement
+//function : AddNode
//purpose :
//=======================================================================
-void SMESHDS_SubMesh::RemoveElement(const Handle(SMDS_MeshElement)& ME)
+void SMESHDS_SubMesh::AddNode(const SMDS_MeshNode * N)
{
- myElements.Remove(ME);
- myListOfEltIDIsUpdate = Standard_False;
+ if ( !IsComplexSubmesh() )
+ myNodes.insert(N);
}
+
//=======================================================================
-//function : AddNode
+//function : RemoveNode
//purpose :
//=======================================================================
-void SMESHDS_SubMesh::AddNode (const Handle(SMDS_MeshNode)& N)
+
+bool SMESHDS_SubMesh::RemoveNode(const SMDS_MeshNode * N)
{
- myNodes.Add(N);
- myListOfNodeIDIsUpdate = Standard_False;
+ if ( !IsComplexSubmesh() && NbNodes() )
+ return myNodes.erase(N);
+
+ return false;
}
//=======================================================================
-//function : RemoveNode
+//function : NbElements
//purpose :
//=======================================================================
-void SMESHDS_SubMesh::RemoveNode (const Handle(SMDS_MeshNode)& N)
+int SMESHDS_SubMesh::NbElements() const
{
- myNodes.Remove(N);
- myListOfNodeIDIsUpdate = Standard_False;
+ if ( !IsComplexSubmesh() )
+ return myElements.size();
+
+ int nbElems = 0;
+#ifndef WNT
+ set<const SMESHDS_SubMesh*>::iterator it = mySubMeshes.begin();
+#else
+ set<const SMESHDS_SubMesh*>::const_iterator it = mySubMeshes.begin();
+#endif
+ for ( ; it != mySubMeshes.end(); it++ )
+ nbElems += (*it)->NbElements();
+
+ return nbElems;
}
//=======================================================================
-//function : NbElements
+//function : NbNodes
//purpose :
//=======================================================================
-Standard_Integer SMESHDS_SubMesh::NbElements()
+
+int SMESHDS_SubMesh::NbNodes() const
{
- return myElements.Extent();
+ if ( !IsComplexSubmesh() )
+ return myNodes.size();
+
+ int nbElems = 0;
+#ifndef WNT
+ set<const SMESHDS_SubMesh*>::iterator it = mySubMeshes.begin();
+#else
+ set<const SMESHDS_SubMesh*>::const_iterator it = mySubMeshes.begin();
+#endif
+ for ( ; it != mySubMeshes.end(); it++ )
+ nbElems += (*it)->NbNodes();
+
+ return nbElems;
}
+// =====================
+// class MySetIterator
+// =====================
+
+template<typename T> class MySetIterator:public SMDS_Iterator<const T*>
+{
+ typedef const set<const T*> TSet;
+ typename TSet::const_iterator myIt;
+ TSet& mySet;
+
+ public:
+ MySetIterator(const set<const T*>& s):mySet(s), myIt(s.begin())
+ {
+ }
+
+ bool more()
+ {
+ return myIt!=mySet.end();
+ }
+ const T* next()
+ {
+ const T* t=*myIt;
+ myIt++;
+ return t;
+ }
+};
+
+// =====================
+// class MyIterator
+// =====================
+
+template<typename VALUE> class MyIterator : public SMDS_Iterator<VALUE>
+{
+ public:
+ MyIterator (const set<const SMESHDS_SubMesh*>& theSubMeshes)
+ : mySubMeshes( theSubMeshes ), mySubIt( theSubMeshes.begin() ), myMore(false)
+ {}
+ bool more()
+ {
+ while (( !myElemIt.get() || !myElemIt->more() ) &&
+ mySubIt != mySubMeshes.end())
+ {
+ myElemIt = getElements(*mySubIt);
+ mySubIt++;
+ }
+ myMore = myElemIt.get() && myElemIt->more();
+ return myMore;
+ }
+ VALUE next()
+ {
+ VALUE elem = 0;
+ if ( myMore )
+ elem = myElemIt->next();
+ return elem;
+ }
+ protected:
+ virtual boost::shared_ptr< SMDS_Iterator<VALUE> >
+ getElements(const SMESHDS_SubMesh*) const = 0;
+
+ private:
+ bool myMore;
+ const set<const SMESHDS_SubMesh*>& mySubMeshes;
+ set<const SMESHDS_SubMesh*>::const_iterator mySubIt;
+ boost::shared_ptr< SMDS_Iterator<VALUE> > myElemIt;
+};
+
+// =====================
+// class MyElemIterator
+// =====================
+
+class MyElemIterator: public MyIterator<const SMDS_MeshElement*>
+{
+ public:
+ MyElemIterator (const set<const SMESHDS_SubMesh*>& theSubMeshes)
+ :MyIterator<const SMDS_MeshElement*>( theSubMeshes ) {}
+ SMDS_ElemIteratorPtr getElements(const SMESHDS_SubMesh* theSubMesh) const
+ { return theSubMesh->GetElements(); }
+};
+
+// =====================
+// class MyNodeIterator
+// =====================
+
+class MyNodeIterator: public MyIterator<const SMDS_MeshNode*>
+{
+ public:
+ MyNodeIterator (const set<const SMESHDS_SubMesh*>& theSubMeshes)
+ :MyIterator<const SMDS_MeshNode*>( theSubMeshes ) {}
+ SMDS_NodeIteratorPtr getElements(const SMESHDS_SubMesh* theSubMesh) const
+ { return theSubMesh->GetNodes(); }
+};
+
//=======================================================================
//function : GetElements
//purpose :
//=======================================================================
-const SMDS_MapOfMeshElement& SMESHDS_SubMesh::GetElements()
+
+SMDS_ElemIteratorPtr SMESHDS_SubMesh::GetElements() const
{
- return myElements;
+ if ( IsComplexSubmesh() )
+ return SMDS_ElemIteratorPtr( new MyElemIterator( mySubMeshes ));
+
+ return SMDS_ElemIteratorPtr(new MySetIterator<SMDS_MeshElement>(myElements));
}
+
//=======================================================================
-//function : NbNodes
+//function : GetNodes
//purpose :
//=======================================================================
-Standard_Integer SMESHDS_SubMesh::NbNodes()
+
+SMDS_NodeIteratorPtr SMESHDS_SubMesh::GetNodes() const
{
- return myNodes.Extent();
+ if ( IsComplexSubmesh() )
+ return SMDS_NodeIteratorPtr( new MyNodeIterator( mySubMeshes ));
+
+ return SMDS_NodeIteratorPtr(new MySetIterator<SMDS_MeshNode>(myNodes));
}
//=======================================================================
-//function : GetNodes
+//function : Contains
+//purpose : check if elem or node is in
+//=======================================================================
+
+bool SMESHDS_SubMesh::Contains(const SMDS_MeshElement * ME) const
+{
+ // DO NOT TRY TO FIND A REMOVED ELEMENT !!
+ if ( !ME )
+ return false;
+
+ if ( IsComplexSubmesh() )
+ {
+ set<const SMESHDS_SubMesh*>::const_iterator aSubIt = mySubMeshes.begin();
+ for ( ; aSubIt != mySubMeshes.end(); aSubIt++ )
+ if ( (*aSubIt)->Contains( ME ))
+ return true;
+ return false;
+ }
+
+ if ( ME->GetType() == SMDSAbs_Node )
+ {
+ const SMDS_MeshNode* n = static_cast<const SMDS_MeshNode*>( ME );
+ return ( myNodes.find( n ) != myNodes.end() );
+ }
+
+ return ( myElements.find( ME ) != myElements.end() );
+}
+
+//=======================================================================
+//function : AddSubMesh
//purpose :
//=======================================================================
-const SMDS_MapOfMeshElement& SMESHDS_SubMesh::GetNodes()
+
+void SMESHDS_SubMesh::AddSubMesh( const SMESHDS_SubMesh* theSubMesh )
{
- return myNodes;
+ ASSERT( theSubMesh );
+ mySubMeshes.insert( theSubMesh );
}
//=======================================================================
-//function : GetIDElements
+//function : RemoveSubMesh
//purpose :
//=======================================================================
-const TColStd_ListOfInteger& SMESHDS_SubMesh::GetIDElements()
+
+bool SMESHDS_SubMesh::RemoveSubMesh( const SMESHDS_SubMesh* theSubMesh )
{
- if (!myListOfEltIDIsUpdate) {
- myListOfEltID.Clear();
- for (SMDS_MapIteratorOfExtendedMap it(myElements); it.More(); it.Next()) {
- myListOfEltID.Append(it.Key()->GetID());
- }
- myListOfEltIDIsUpdate = Standard_True;
- }
- return myListOfEltID;
+ return mySubMeshes.erase( theSubMesh );
}
//=======================================================================
-//function : GetIDNodes
+//function : ContainsSubMesh
//purpose :
//=======================================================================
-const TColStd_ListOfInteger& SMESHDS_SubMesh::GetIDNodes()
+
+bool SMESHDS_SubMesh::ContainsSubMesh( const SMESHDS_SubMesh* theSubMesh ) const
{
- if (!myListOfNodeIDIsUpdate) {
- myListOfNodeID.Clear();
- for (SMDS_MapIteratorOfExtendedMap it(myNodes); it.More(); it.Next()) {
- myListOfNodeID.Append(it.Key()->GetID());
- }
- myListOfNodeIDIsUpdate = Standard_True;
+ return mySubMeshes.find( theSubMesh ) != mySubMeshes.end();
+}
+
+//=======================================================================
+//function : GetSubMeshIterator
+//purpose :
+//=======================================================================
+
+SMESHDS_SubMeshIteratorPtr SMESHDS_SubMesh::GetSubMeshIterator() const
+{
+ typedef set<const SMESHDS_SubMesh*>::const_iterator TIterator;
+ return SMESHDS_SubMeshIteratorPtr
+ ( new SMDS_SetIterator< const SMESHDS_SubMesh*, TIterator >( mySubMeshes.begin(),
+ mySubMeshes.end()));
+}
+
+//=======================================================================
+//function : Clear
+//purpose : remove the contents
+//=======================================================================
+
+void SMESHDS_SubMesh::Clear()
+{
+ myElements.clear();
+ myNodes.clear();
+ SMESHDS_SubMeshIteratorPtr sub = GetSubMeshIterator();
+ while ( sub->more() ) {
+ if ( SMESHDS_SubMesh* sm = (SMESHDS_SubMesh*) sub->next())
+ sm->Clear();
}
- return myListOfNodeID;
}