src/SMDS/SMDS_MeshElementIDFactory.cxx

   1 // Copyright (C) 2007-2016  CEA/DEN, EDF R&D, OPEN CASCADE
   2 //
   3 // Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   5 //
   6 // This library is free software; you can redistribute it and/or
   7 // modify it under the terms of the GNU Lesser General Public
   8 // License as published by the Free Software Foundation; either
   9 // version 2.1 of the License, or (at your option) any later version.
  10 //
  11 // This library is distributed in the hope that it will be useful,
  12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14 // Lesser General Public License for more details.
  15 //
  16 // You should have received a copy of the GNU Lesser General Public
  17 // License along with this library; if not, write to the Free Software
  18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  19 //
  20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 //
  22
  23 //  SMESH SMDS : implementaion of Salome mesh data structure
  24 //  File   : SMDS_MeshElementIDFactory.cxx
  25 //  Author : Jean-Michel BOULCOURT
  26 //  Module : SMESH
  27 //
  28 #ifdef _MSC_VER
  29 #pragma warning(disable:4786)
  30 #endif
  31
  32 #include "SMDS_MeshElementIDFactory.hxx"
  33 #include "SMDS_MeshElement.hxx"
  34 #include "SMDS_Mesh.hxx"
  35
  36 #include "utilities.h"
  37
  38 #include "SMDS_UnstructuredGrid.hxx"
  39 #include <vtkCellType.h>
  40
  41 #include <climits>
  42
  43 using namespace std;
  44
  45 //=======================================================================
  46 //function : SMDS_MeshElementIDFactory
  47 //purpose  :
  48 //=======================================================================
  49 SMDS_MeshElementIDFactory::SMDS_MeshElementIDFactory():
  50   SMDS_MeshNodeIDFactory()
  51 {
  52 }
  53
  54 int SMDS_MeshElementIDFactory::SetInVtkGrid(SMDS_MeshElement * elem)
  55 {
  56   // --- retrieve nodes ID
  57
  58   SMDS_MeshCell *cell = dynamic_cast<SMDS_MeshCell*>(elem);
  59   assert(cell);
  60   vector<vtkIdType> nodeIds( elem->NbNodes() );
  61   SMDS_ElemIteratorPtr it = elem->nodesIterator();
  62   for( int i = 0; it->more(); ++i )
  63   {
  64     int nodeId = (static_cast<const SMDS_MeshNode*>(it->next()))->getVtkId();
  65     nodeIds[i] = nodeId;
  66   }
  67
  68   // --- insert cell in vtkUnstructuredGrid
  69
  70   int typ = VTK_VERTEX;
  71   int cellId = myMesh->getGrid()->InsertNextLinkedCell(typ, nodeIds.size(), &nodeIds[0]);
  72   cell->setVtkId(cellId);
  73   return cellId;
  74 }
  75
  76 //=======================================================================
  77 //function : BindID
  78 //purpose  :
  79 //=======================================================================
  80
  81 bool SMDS_MeshElementIDFactory::BindID(int ID, SMDS_MeshElement * elem)
  82 {
  83   SetInVtkGrid(elem);
  84   return myMesh->registerElement(ID, elem);
  85 }
  86
  87 //=======================================================================
  88 //function : MeshElement
  89 //purpose  :
  90 //=======================================================================
  91 SMDS_MeshElement* SMDS_MeshElementIDFactory::MeshElement(int ID)
  92 {
  93   if ( ID<1 || ID >= (int) myMesh->myCells.size() )
  94     return NULL;
  95   const SMDS_MeshElement* elem = GetMesh()->FindElement(ID);
  96   return (SMDS_MeshElement*)(elem);
  97 }
  98
  99 //=======================================================================
 100 //function : GetFreeID
 101 //purpose  :
 102 //=======================================================================
 103
 104 int SMDS_MeshElementIDFactory::GetFreeID()
 105 {
 106   int ID;
 107   do {
 108     ID = SMDS_MeshIDFactory::GetFreeID();
 109   } while ( MeshElement( ID ));
 110   return ID;
 111 }
 112
 113 //=======================================================================
 114 //function : ReleaseID
 115 //purpose  :
 116 //=======================================================================
 117 void SMDS_MeshElementIDFactory::ReleaseID(int ID, int vtkId)
 118 {
 119   if (ID < 1) // TODO check case ID == O
 120   {
 121     MESSAGE("~~~~~~~~~~~~~~ SMDS_MeshElementIDFactory::ReleaseID ID = " << ID);
 122     return;
 123   }
 124   if (vtkId >= 0)
 125   {
 126     assert(vtkId < (int)myMesh->myCellIdVtkToSmds.size());
 127     myMesh->myCellIdVtkToSmds[vtkId] = -1;
 128     myMesh->setMyModified();
 129   }
 130   SMDS_MeshIDFactory::ReleaseID(ID);
 131   if (ID == myMax)
 132     myMax = 0;
 133   if (ID == myMin)
 134     myMax = 0;
 135 }
 136
 137 //=======================================================================
 138 //function : updateMinMax
 139 //purpose  :
 140 //=======================================================================
 141
 142 void SMDS_MeshElementIDFactory::updateMinMax() const
 143 {
 144   myMin = INT_MAX;
 145   myMax = 0;
 146   for (size_t i = 0; i < myMesh->myCells.size(); i++)
 147   {
 148     if (myMesh->myCells[i])
 149     {
 150       int id = myMesh->myCells[i]->GetID();
 151       if (id > myMax)
 152         myMax = id;
 153       if (id < myMin)
 154         myMin = id;
 155     }
 156   }
 157   if (myMin == INT_MAX)
 158     myMin = 0;
 159 }
 160
 161 //=======================================================================
 162 //function : elementsIterator
 163 //purpose  : Return an iterator on elements of the factory
 164 //=======================================================================
 165
 166 SMDS_ElemIteratorPtr SMDS_MeshElementIDFactory::elementsIterator() const
 167 {
 168   return myMesh->elementsIterator(SMDSAbs_All);
 169 }
 170
 171 void SMDS_MeshElementIDFactory::Clear()
 172 {
 173   myMesh->myCellIdVtkToSmds.clear();
 174   myMin = myMax = 0;
 175   SMDS_MeshIDFactory::Clear();
 176 }