src/SMDS/SMDS_MeshPrism.cxx

   1 using namespace std;
   2 // File:        SMDS_MeshPrism.cxx
   3 // Created:     Wed Jan 23 17:02:34 2002
   4 // Author:      Jean-Michel BOULCOURT
   5 //              <jmb@coulox.paris1.matra-dtv.fr>
   6
   7
   8 #include "SMDS_MeshPrism.ixx"
   9 #include "Standard_ConstructionError.hxx"
  10
  11
  12 //=======================================================================
  13 //function : SMDS_MeshPrism
  14 //purpose  : Prism
  15 //=======================================================================
  16
  17 SMDS_MeshPrism::SMDS_MeshPrism(const Standard_Integer ID,
  18                                const Standard_Integer idnode1,
  19                                const Standard_Integer idnode2,
  20                                const Standard_Integer idnode3,
  21                                const Standard_Integer idnode4,
  22                                const Standard_Integer idnode5,
  23                                const Standard_Integer idnode6):SMDS_MeshVolume(ID,6)
  24 {
  25   SetConnections(idnode1,idnode2,idnode3,idnode4,idnode5,idnode6);
  26   ComputeKey();
  27 }
  28
  29
  30 //=======================================================================
  31 //function : SetConnections
  32 //purpose  : Prism
  33 //=======================================================================
  34 void SMDS_MeshPrism::SetConnections(const Standard_Integer idnode1,
  35                                     const Standard_Integer idnode2,
  36                                     const Standard_Integer idnode3,
  37                                     const Standard_Integer idnode4,
  38                                     const Standard_Integer idnode5,
  39                                     const Standard_Integer idnode6)
  40 {
  41   Standard_Integer idmin = (idnode1 < idnode2 ? idnode1 : idnode2);
  42   idmin = (idmin < idnode3 ? idmin : idnode3);
  43   idmin = (idmin < idnode4 ? idmin : idnode4);
  44   idmin = (idmin < idnode5 ? idmin : idnode5);
  45   idmin = (idmin < idnode6 ? idmin : idnode6);
  46
  47   myNodes[0] = idmin;
  48   if (idmin == idnode1) { // 1 2 3 4 5 6
  49     myNodes[1] = idnode2;
  50     myNodes[2] = idnode3;
  51     myNodes[3] = idnode4;
  52     myNodes[4] = idnode5;
  53     myNodes[5] = idnode6;
  54   } else if (idmin == idnode2) { // 2 3 1 5 6 4
  55     myNodes[1] = idnode3;
  56     myNodes[2] = idnode1;
  57     myNodes[3] = idnode5;
  58     myNodes[4] = idnode6;
  59     myNodes[5] = idnode4;
  60   } else if (idmin == idnode3) { // 3 1 2 6 4 5
  61     myNodes[1] = idnode1;
  62     myNodes[2] = idnode2;
  63     myNodes[3] = idnode6;
  64     myNodes[4] = idnode4;
  65     myNodes[5] = idnode5;
  66   } else if (idmin == idnode4) { // 4 5 6 1 2 3
  67     myNodes[1] = idnode5;
  68     myNodes[2] = idnode6;
  69     myNodes[3] = idnode1;
  70     myNodes[4] = idnode2;
  71     myNodes[5] = idnode3;
  72   } else if (idmin == idnode5) { // 5 6 4 2 3 1
  73     myNodes[1] = idnode6;
  74     myNodes[2] = idnode4;
  75     myNodes[3] = idnode2;
  76     myNodes[4] = idnode3;
  77     myNodes[5] = idnode1;
  78   } else {                       // 6 4 5 3 1 2
  79     myNodes[1] = idnode4;
  80     myNodes[2] = idnode5;
  81     myNodes[3] = idnode3;
  82     myNodes[4] = idnode1;
  83     myNodes[5] = idnode2;
  84   }
  85
  86 }
  87
  88
  89 //=======================================================================
  90 //function : NodesOfFace
  91 //purpose  : returns the rank node in mynodes. Useful to extract faces from volume
  92 //=======================================================================
  93 Standard_Integer SMDS_MeshPrism::NodesOfFace(const Standard_Integer rankface,
  94                                                const Standard_Integer ranknode)
  95 {
  96   static Standard_Integer facenode[5][4] = {
  97     {0,1,2,-1},
  98     {3,5,4,-1},
  99     {0,3,4,1},
 100     {1,4,5,2},
 101     {2,5,3,0}
 102   };
 103
 104   return facenode[rankface-1][ranknode-1];
 105 }
 106
 107 //=======================================================================
 108 //function : NodesOfEdge
 109 //purpose  : returns the rank node in mynodes. Useful to extract edges from volume
 110 //=======================================================================
 111 Standard_Integer SMDS_MeshPrism::NodesOfEdge(const Standard_Integer rankedge,
 112                                                const Standard_Integer ranknode) const
 113 {
 114   static Standard_Integer faceedge[9][2] = {
 115     {0,1},
 116     {1,2},
 117     {0,2},
 118     {3,5},
 119     {4,5},
 120     {3,4},
 121     {0,3},
 122     {1,4},
 123     {2,5}
 124
 125   };
 126
 127   return faceedge[rankedge-1][ranknode-1];
 128 }
 129
 130 //=======================================================================
 131 //function : GetFaceDefinedByNodes
 132 //purpose  :
 133 //=======================================================================
 134 void SMDS_MeshPrism::GetFaceDefinedByNodes(const Standard_Integer rank,
 135                                            const Standard_Address idnode,
 136                                            Standard_Integer& nb) const
 137 {
 138   Standard_Integer *ptr;
 139   ptr = (Standard_Integer *)idnode;
 140   ptr[0] = myNodes[NodesOfFace(rank,1)];
 141   ptr[1] = myNodes[NodesOfFace(rank,2)];
 142   ptr[2] = myNodes[NodesOfFace(rank,3)];
 143   ptr[3] = myNodes[NodesOfFace(rank,4)];
 144   nb = (NodesOfFace(rank,4) == -1 ? 3 : 4);
 145   if (nb == 4)
 146     ptr[3] = myNodes[NodesOfFace(rank,4)];
 147
 148 }
 149
 150 //=======================================================================
 151 //function : GetEdgeDefinedByNodes
 152 //purpose  :
 153 //=======================================================================
 154 void SMDS_MeshPrism::GetEdgeDefinedByNodes(const Standard_Integer rank,
 155                                             Standard_Integer& idnode1,
 156                                             Standard_Integer& idnode2) const
 157 {
 158   idnode1 = myNodes[NodesOfEdge(rank,1)];
 159   idnode2 = myNodes[NodesOfEdge(rank,2)];
 160 }
 161