src/SMDS/SMDS_VolumeOfNodes.cxx

   1 //  Copyright (C) 2007-2010  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.
  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 //
  25 #ifdef _MSC_VER
  26 #pragma warning(disable:4786)
  27 #endif
  28
  29 #include "SMDS_VolumeOfNodes.hxx"
  30 #include "SMDS_MeshNode.hxx"
  31 #include "SMDS_SetIterator.hxx"
  32 #include "SMDS_VolumeTool.hxx"
  33 #include "utilities.h"
  34
  35 #include <vector>
  36
  37 using namespace std;
  38
  39 ///////////////////////////////////////////////////////////////////////////////
  40 /// Create an hexahedron. node 1,2,3,4 and 5,6,7,8 are quadrangle and
  41 /// 5,1 and 7,3 are an edges.
  42 ///////////////////////////////////////////////////////////////////////////////
  43 SMDS_VolumeOfNodes::SMDS_VolumeOfNodes(
  44                 const SMDS_MeshNode * node1,
  45                 const SMDS_MeshNode * node2,
  46                 const SMDS_MeshNode * node3,
  47                 const SMDS_MeshNode * node4,
  48                 const SMDS_MeshNode * node5,
  49                 const SMDS_MeshNode * node6,
  50                 const SMDS_MeshNode * node7,
  51                 const SMDS_MeshNode * node8)
  52 {
  53         myNbNodes = 8;
  54         myNodes = new const SMDS_MeshNode* [myNbNodes];
  55         myNodes[0]=node1;
  56         myNodes[1]=node2;
  57         myNodes[2]=node3;
  58         myNodes[3]=node4;
  59         myNodes[4]=node5;
  60         myNodes[5]=node6;
  61         myNodes[6]=node7;
  62         myNodes[7]=node8;
  63 }
  64
  65 SMDS_VolumeOfNodes::SMDS_VolumeOfNodes(
  66                 const SMDS_MeshNode * node1,
  67                 const SMDS_MeshNode * node2,
  68                 const SMDS_MeshNode * node3,
  69                 const SMDS_MeshNode * node4)
  70 {
  71         myNbNodes = 4;
  72         myNodes = new const SMDS_MeshNode* [myNbNodes];
  73         myNodes[0]=node1;
  74         myNodes[1]=node2;
  75         myNodes[2]=node3;
  76         myNodes[3]=node4;
  77 }
  78
  79 SMDS_VolumeOfNodes::SMDS_VolumeOfNodes(
  80                 const SMDS_MeshNode * node1,
  81                 const SMDS_MeshNode * node2,
  82                 const SMDS_MeshNode * node3,
  83                 const SMDS_MeshNode * node4,
  84                 const SMDS_MeshNode * node5)
  85 {
  86         myNbNodes = 5;
  87         myNodes = new const SMDS_MeshNode* [myNbNodes];
  88         myNodes[0]=node1;
  89         myNodes[1]=node2;
  90         myNodes[2]=node3;
  91         myNodes[3]=node4;
  92         myNodes[4]=node5;
  93 }
  94
  95 SMDS_VolumeOfNodes::SMDS_VolumeOfNodes(
  96                 const SMDS_MeshNode * node1,
  97                 const SMDS_MeshNode * node2,
  98                 const SMDS_MeshNode * node3,
  99                 const SMDS_MeshNode * node4,
 100                 const SMDS_MeshNode * node5,
 101                 const SMDS_MeshNode * node6)
 102 {
 103         myNbNodes = 6;
 104         myNodes = new const SMDS_MeshNode* [myNbNodes];
 105         myNodes[0]=node1;
 106         myNodes[1]=node2;
 107         myNodes[2]=node3;
 108         myNodes[3]=node4;
 109         myNodes[4]=node5;
 110         myNodes[5]=node6;
 111 }
 112
 113 bool SMDS_VolumeOfNodes::ChangeNodes(const SMDS_MeshNode* nodes[],
 114                                      const int            nbNodes)
 115 {
 116   if (nbNodes < 4 || nbNodes > 8 || nbNodes == 7)
 117     return false;
 118
 119   delete [] myNodes;
 120   myNbNodes = nbNodes;
 121   myNodes = new const SMDS_MeshNode* [myNbNodes];
 122   for ( int i = 0; i < nbNodes; i++ )
 123     myNodes[ i ] = nodes [ i ];
 124
 125   return true;
 126 }
 127
 128 SMDS_VolumeOfNodes::~SMDS_VolumeOfNodes()
 129 {
 130   if (myNodes != NULL) {
 131     delete [] myNodes;
 132     myNodes = NULL;
 133   }
 134 }
 135
 136 //=======================================================================
 137 //function : Print
 138 //purpose  :
 139 //=======================================================================
 140
 141 void SMDS_VolumeOfNodes::Print(ostream & OS) const
 142 {
 143         OS << "volume <" << GetID() << "> : ";
 144         int i;
 145         for (i = 0; i < NbNodes()-1; ++i) OS << myNodes[i] << ",";
 146         OS << myNodes[NbNodes()-1]<< ") " << endl;
 147 }
 148
 149 int SMDS_VolumeOfNodes::NbFaces() const
 150 {
 151         switch(NbNodes())
 152         {
 153         case 4: return 4;
 154         case 5: return 5;
 155         case 6: return 5;
 156         case 8: return 6;
 157         default: MESSAGE("invalid number of nodes");
 158         }
 159         return 0;
 160 }
 161
 162 int SMDS_VolumeOfNodes::NbNodes() const
 163 {
 164         return myNbNodes;
 165 }
 166
 167 int SMDS_VolumeOfNodes::NbEdges() const
 168 {
 169         switch(NbNodes())
 170         {
 171         case 4: return 6;
 172         case 5: return 8;
 173         case 6: return 9;
 174         case 8: return 12;
 175         default: MESSAGE("invalid number of nodes");
 176         }
 177         return 0;
 178 }
 179
 180 /// ===================================================================
 181 /*!
 182  * \brief Iterator on node of volume
 183  */
 184 /// ===================================================================
 185
 186 class SMDS_VolumeOfNodes_MyIterator:public SMDS_NodeArrayElemIterator
 187 {
 188  public:
 189   SMDS_VolumeOfNodes_MyIterator(const SMDS_MeshNode* const* s, int l):
 190     SMDS_NodeArrayElemIterator( s, & s[ l ]) {}
 191 };
 192
 193 /// ===================================================================
 194 /*!
 195  * \brief Iterator on faces or edges of volume
 196  */
 197 /// ===================================================================
 198
 199 class _MySubIterator : public SMDS_ElemIterator
 200 {
 201   vector< const SMDS_MeshElement* > myElems;
 202   int myIndex;
 203 public:
 204   _MySubIterator(const SMDS_VolumeOfNodes* vol, SMDSAbs_ElementType type):myIndex(0) {
 205     SMDS_VolumeTool vTool(vol);
 206     if (type == SMDSAbs_Face)
 207       vTool.GetAllExistingFaces( myElems );
 208     else
 209       vTool.GetAllExistingFaces( myElems );
 210   }
 211   /// Return true if and only if there are other object in this iterator
 212   virtual bool more() { return myIndex < myElems.size(); }
 213
 214   /// Return the current object and step to the next one
 215   virtual const SMDS_MeshElement* next() { return myElems[ myIndex++ ]; }
 216 };
 217
 218 SMDS_ElemIteratorPtr SMDS_VolumeOfNodes::elementsIterator(SMDSAbs_ElementType type) const
 219 {
 220   switch(type)
 221   {
 222   case SMDSAbs_Volume:
 223     return SMDS_MeshElement::elementsIterator(SMDSAbs_Volume);
 224   case SMDSAbs_Node:
 225     return SMDS_ElemIteratorPtr(new SMDS_VolumeOfNodes_MyIterator(myNodes,myNbNodes));
 226   case SMDSAbs_Face:
 227     return SMDS_ElemIteratorPtr(new _MySubIterator(this,SMDSAbs_Face));
 228   case SMDSAbs_Edge:
 229     return SMDS_ElemIteratorPtr(new _MySubIterator(this,SMDSAbs_Edge));
 230   default:
 231     MESSAGE("ERROR : Iterator not implemented");
 232     return SMDS_ElemIteratorPtr((SMDS_ElemIterator*)NULL);
 233   }
 234 }
 235
 236 SMDSAbs_ElementType SMDS_VolumeOfNodes::GetType() const
 237 {
 238         return SMDSAbs_Volume;
 239 }
 240
 241 /*!
 242  * \brief Return node by its index
 243  * \param ind - node index
 244  * \retval const SMDS_MeshNode* - the node
 245  */
 246 const SMDS_MeshNode* SMDS_VolumeOfNodes::GetNode(const int ind) const
 247 {
 248   return myNodes[ ind ];
 249 }
 250
 251 SMDSAbs_EntityType SMDS_VolumeOfNodes::GetEntityType() const
 252 {
 253   SMDSAbs_EntityType aType = SMDSEntity_Tetra;
 254   switch(myNbNodes)
 255   {
 256   case 4: aType = SMDSEntity_Tetra;   break;
 257   case 5: aType = SMDSEntity_Pyramid; break;
 258   case 6: aType = SMDSEntity_Penta;   break;
 259   case 8:
 260   default: aType = SMDSEntity_Hexa;    break;
 261   }
 262   return aType;
 263 }