Salome HOME
Merge from BR_V5_DEV 16Feb09
[modules/smesh.git] / src / SMESHDS / SMESHDS_SubMesh.cxx
index f38b346c55d0f5447a8eab108e953a1d165d1aba..127c5c2d99ccd100786b66dd0e3c858ab40f4574 100644 (file)
-//  SMESH SMESHDS : management of mesh data and SMESH document
+//  Copyright (C) 2007-2008  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
+//
+//  This library is free software; you can redistribute it and/or
+//  modify it under the terms of the GNU Lesser General Public
+//  License as published by the Free Software Foundation; either
+//  version 2.1 of the License.
 //
-//  Copyright (C) 2003  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
-//  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS 
-// 
-//  This library is free software; you can redistribute it and/or 
-//  modify it under the terms of the GNU Lesser General Public 
-//  License as published by the Free Software Foundation; either 
-//  version 2.1 of the License. 
-// 
-//  This library is distributed in the hope that it will be useful, 
-//  but WITHOUT ANY WARRANTY; without even the implied warranty of 
-//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU 
-//  Lesser General Public License for more details. 
-// 
-//  You should have received a copy of the GNU Lesser General Public 
-//  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 
+//  This library is distributed in the hope that it will be useful,
+//  but WITHOUT ANY WARRANTY; without even the implied warranty of
+//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+//  Lesser General Public License for more details.
 //
+//  You should have received a copy of the GNU Lesser General Public
+//  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.salome-platform.org/ or email : webmaster.salome@opencascade.com
 //
+//  SMESH SMESHDS : management of mesh data and SMESH document
 //  File   : SMESH_SubMesh.cxx
 //  Author : Yves FRICAUD, OCC
 //  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;
 }