src/SMDS/SMDS_MeshPrism.cxx

   1 //  SMESH SMDS : implementaion of Salome mesh data structure
   2 //
   3 //  Copyright (C) 2003  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.opencascade.org/SALOME/ or email : webmaster.salome@opencascade.org
  21 //
  22 //
  23 //
  24 //  File   : SMDS_MeshPrism.cxx
  25 //  Author : Jean-Michel BOULCOURT
  26 //  Module : SMESH
  27
  28 using namespace std;
  29 #include "SMDS_MeshPrism.ixx"
  30 #include "Standard_ConstructionError.hxx"
  31
  32
  33 //=======================================================================
  34 //function : SMDS_MeshPrism
  35 //purpose  : Prism
  36 //=======================================================================
  37
  38 SMDS_MeshPrism::SMDS_MeshPrism(const Standard_Integer ID,
  39                                const Standard_Integer idnode1,
  40                                const Standard_Integer idnode2,
  41                                const Standard_Integer idnode3,
  42                                const Standard_Integer idnode4,
  43                                const Standard_Integer idnode5,
  44                                const Standard_Integer idnode6):SMDS_MeshVolume(ID,6)
  45 {
  46   SetConnections(idnode1,idnode2,idnode3,idnode4,idnode5,idnode6);
  47   ComputeKey();
  48 }
  49
  50
  51 //=======================================================================
  52 //function : SetConnections
  53 //purpose  : Prism
  54 //=======================================================================
  55 void SMDS_MeshPrism::SetConnections(const Standard_Integer idnode1,
  56                                     const Standard_Integer idnode2,
  57                                     const Standard_Integer idnode3,
  58                                     const Standard_Integer idnode4,
  59                                     const Standard_Integer idnode5,
  60                                     const Standard_Integer idnode6)
  61 {
  62   Standard_Integer idmin = (idnode1 < idnode2 ? idnode1 : idnode2);
  63   idmin = (idmin < idnode3 ? idmin : idnode3);
  64   idmin = (idmin < idnode4 ? idmin : idnode4);
  65   idmin = (idmin < idnode5 ? idmin : idnode5);
  66   idmin = (idmin < idnode6 ? idmin : idnode6);
  67
  68   myNodes[0] = idmin;
  69   if (idmin == idnode1) { // 1 2 3 4 5 6
  70     myNodes[1] = idnode2;
  71     myNodes[2] = idnode3;
  72     myNodes[3] = idnode4;
  73     myNodes[4] = idnode5;
  74     myNodes[5] = idnode6;
  75   } else if (idmin == idnode2) { // 2 3 1 5 6 4
  76     myNodes[1] = idnode3;
  77     myNodes[2] = idnode1;
  78     myNodes[3] = idnode5;
  79     myNodes[4] = idnode6;
  80     myNodes[5] = idnode4;
  81   } else if (idmin == idnode3) { // 3 1 2 6 4 5
  82     myNodes[1] = idnode1;
  83     myNodes[2] = idnode2;
  84     myNodes[3] = idnode6;
  85     myNodes[4] = idnode4;
  86     myNodes[5] = idnode5;
  87   } else if (idmin == idnode4) { // 4 5 6 1 2 3
  88     myNodes[1] = idnode5;
  89     myNodes[2] = idnode6;
  90     myNodes[3] = idnode1;
  91     myNodes[4] = idnode2;
  92     myNodes[5] = idnode3;
  93   } else if (idmin == idnode5) { // 5 6 4 2 3 1
  94     myNodes[1] = idnode6;
  95     myNodes[2] = idnode4;
  96     myNodes[3] = idnode2;
  97     myNodes[4] = idnode3;
  98     myNodes[5] = idnode1;
  99   } else {                       // 6 4 5 3 1 2
 100     myNodes[1] = idnode4;
 101     myNodes[2] = idnode5;
 102     myNodes[3] = idnode3;
 103     myNodes[4] = idnode1;
 104     myNodes[5] = idnode2;
 105   }
 106
 107 }
 108
 109
 110 //=======================================================================
 111 //function : NodesOfFace
 112 //purpose  : returns the rank node in mynodes. Useful to extract faces from volume
 113 //=======================================================================
 114 Standard_Integer SMDS_MeshPrism::NodesOfFace(const Standard_Integer rankface,
 115                                                const Standard_Integer ranknode)
 116 {
 117   static Standard_Integer facenode[5][4] = {
 118     {0,1,2,-1},
 119     {3,5,4,-1},
 120     {0,3,4,1},
 121     {1,4,5,2},
 122     {2,5,3,0}
 123   };
 124
 125   return facenode[rankface-1][ranknode-1];
 126 }
 127
 128 //=======================================================================
 129 //function : NodesOfEdge
 130 //purpose  : returns the rank node in mynodes. Useful to extract edges from volume
 131 //=======================================================================
 132 Standard_Integer SMDS_MeshPrism::NodesOfEdge(const Standard_Integer rankedge,
 133                                                const Standard_Integer ranknode) const
 134 {
 135   static Standard_Integer faceedge[9][2] = {
 136     {0,1},
 137     {1,2},
 138     {0,2},
 139     {3,5},
 140     {4,5},
 141     {3,4},
 142     {0,3},
 143     {1,4},
 144     {2,5}
 145
 146   };
 147
 148   return faceedge[rankedge-1][ranknode-1];
 149 }
 150
 151 //=======================================================================
 152 //function : GetFaceDefinedByNodes
 153 //purpose  :
 154 //=======================================================================
 155 void SMDS_MeshPrism::GetFaceDefinedByNodes(const Standard_Integer rank,
 156                                            const Standard_Address idnode,
 157                                            Standard_Integer& nb) const
 158 {
 159   Standard_Integer *ptr;
 160   ptr = (Standard_Integer *)idnode;
 161   ptr[0] = myNodes[NodesOfFace(rank,1)];
 162   ptr[1] = myNodes[NodesOfFace(rank,2)];
 163   ptr[2] = myNodes[NodesOfFace(rank,3)];
 164   ptr[3] = myNodes[NodesOfFace(rank,4)];
 165   nb = (NodesOfFace(rank,4) == -1 ? 3 : 4);
 166   if (nb == 4)
 167     ptr[3] = myNodes[NodesOfFace(rank,4)];
 168
 169 }
 170
 171 //=======================================================================
 172 //function : GetEdgeDefinedByNodes
 173 //purpose  :
 174 //=======================================================================
 175 void SMDS_MeshPrism::GetEdgeDefinedByNodes(const Standard_Integer rank,
 176                                             Standard_Integer& idnode1,
 177                                             Standard_Integer& idnode2) const
 178 {
 179   idnode1 = myNodes[NodesOfEdge(rank,1)];
 180   idnode2 = myNodes[NodesOfEdge(rank,2)];
 181 }
 182