src/SMDS/SMDS_MeshTetrahedron.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_MeshTetrahedron.cxx
  25 //  Author : Jean-Michel BOULCOURT
  26 //  Module : SMESH
  27
  28 using namespace std;
  29 #include "SMDS_MeshTetrahedron.ixx"
  30 #include "Standard_ConstructionError.hxx"
  31
  32 //=======================================================================
  33 //function : SMDS_MeshTetrahedron
  34 //purpose  : Tetrahedra
  35 //=======================================================================
  36
  37 SMDS_MeshTetrahedron::SMDS_MeshTetrahedron(const Standard_Integer ID,
  38                                            const Standard_Integer idnode1,
  39                                            const Standard_Integer idnode2,
  40                                            const Standard_Integer idnode3,
  41                                            const Standard_Integer idnode4):SMDS_MeshVolume(ID,4)
  42 {
  43   SetConnections(idnode1,idnode2,idnode3,idnode4);
  44   ComputeKey();
  45 }
  46
  47
  48 //=======================================================================
  49 //function : SetConnections
  50 //purpose  : Tetrahedra
  51 //=======================================================================
  52 void SMDS_MeshTetrahedron::SetConnections(const Standard_Integer idnode1,
  53                                           const Standard_Integer idnode2,
  54                                           const Standard_Integer idnode3,
  55                                           const Standard_Integer idnode4)
  56 {
  57   Standard_Integer idmin = (idnode1 < idnode2 ? idnode1 : idnode2);
  58   idmin = (idmin < idnode3 ? idmin : idnode3);
  59   idmin = (idmin < idnode4 ? idmin : idnode4);
  60
  61   myNodes[0] = idmin;
  62   if (idmin == idnode1) {
  63     myNodes[1] = idnode2;
  64     myNodes[2] = idnode3;
  65     myNodes[3] = idnode4;
  66   } else if (idmin == idnode2) {
  67     myNodes[1] = idnode3;
  68     myNodes[2] = idnode4;
  69     myNodes[3] = idnode1;
  70   } else if (idmin == idnode3) {
  71     myNodes[1] = idnode4;
  72     myNodes[2] = idnode1;
  73     myNodes[3] = idnode2;
  74   } else {
  75     myNodes[1] = idnode1;
  76     myNodes[2] = idnode2;
  77     myNodes[3] = idnode3;
  78   }
  79
  80 }
  81
  82
  83 //=======================================================================
  84 //function : NodesOfFace
  85 //purpose  : returns the rank node in mynodes. Useful to extract faces from volume
  86 //=======================================================================
  87 Standard_Integer SMDS_MeshTetrahedron::NodesOfFace(const Standard_Integer rankface,
  88                                                    const Standard_Integer ranknode)
  89 {
  90   static Standard_Integer facenode[4][3] = {
  91     {0,1,2},
  92     {0,3,1},
  93     {1,3,2},
  94     {0,2,3}
  95   };
  96
  97   return facenode[rankface-1][ranknode-1];
  98 }
  99
 100
 101 //=======================================================================
 102 //function : NodesOfEdge
 103 //purpose  : returns the rank node in mynodes. Useful to extract edges from volume
 104 //=======================================================================
 105 Standard_Integer SMDS_MeshTetrahedron::NodesOfEdge(const Standard_Integer rankedge,
 106                                                    const Standard_Integer ranknode) const
 107 {
 108   static Standard_Integer faceedge[6][2] = {
 109     {0,1},
 110     {1,2},
 111     {0,2},
 112     {0,3},
 113     {1,3},
 114     {2,3}
 115   };
 116
 117   return faceedge[rankedge-1][ranknode-1];
 118 }
 119
 120
 121 //=======================================================================
 122 //function : GetFaceDefinedByNodes
 123 //purpose  :
 124 //=======================================================================
 125 void SMDS_MeshTetrahedron::GetFaceDefinedByNodes(const Standard_Integer rank,
 126                                                  const Standard_Address idnode,
 127                                                  Standard_Integer& nb) const
 128 {
 129   Standard_Integer *ptr;
 130   nb = 3;
 131   ptr = (Standard_Integer *)idnode;
 132   ptr[0] = myNodes[NodesOfFace(rank,1)];
 133   ptr[1] = myNodes[NodesOfFace(rank,2)];
 134   ptr[2] = myNodes[NodesOfFace(rank,3)];
 135
 136 }
 137
 138 //=======================================================================
 139 //function : GetEdgeDefinedByNodes
 140 //purpose  :
 141 //=======================================================================
 142 void SMDS_MeshTetrahedron::GetEdgeDefinedByNodes(const Standard_Integer rank,
 143                                             Standard_Integer& idnode1,
 144                                             Standard_Integer& idnode2) const
 145 {
 146   idnode1 = myNodes[NodesOfEdge(rank,1)];
 147   idnode2 = myNodes[NodesOfEdge(rank,2)];
 148 }
 149