src/INTERP_KERNEL/CellModel.cxx

   1 // Copyright (C) 2007-2012  CEA/DEN, EDF R&D
   2 //
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Lesser General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2.1 of the License.
   7 //
   8 // This library is distributed in the hope that it will be useful,
   9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 // Lesser General Public License for more details.
  12 //
  13 // You should have received a copy of the GNU Lesser General Public
  14 // License along with this library; if not, write to the Free Software
  15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 //
  17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 //
  19 // Author : Anthony Geay (CEA/DEN)
  20
  21 #include "CellModel.hxx"
  22
  23 #include "InterpKernelException.hxx"
  24
  25 #include <algorithm>
  26 #include <sstream>
  27 #include <vector>
  28 #include <limits>
  29
  30 namespace INTERP_KERNEL
  31 {
  32   const char *CellModel::CELL_TYPES_REPR[]={"NORM_POINT1", "NORM_SEG2", "NORM_SEG3", "NORM_TRI3", "NORM_QUAD4",// 0->4
  33                                             "NORM_POLYGON", "NORM_TRI6", "NORM_TRI7" , "NORM_QUAD8", "NORM_QUAD9",//5->9
  34                                             "NORM_SEG4", "", "", "", "NORM_TETRA4",//10->14
  35                                             "NORM_PYRA5", "NORM_PENTA6", "", "NORM_HEXA8", "",//15->19
  36                                             "NORM_TETRA10", "", "NORM_HEXGP12", "NORM_PYRA13", "",//20->24
  37                                             "NORM_PENTA15", "", "NORM_HEXA27", "", "",//25->29
  38                                             "NORM_HEXA20", "NORM_POLYHED", "NORM_QPOLYG", "NORM_POLYL", "",//30->34
  39                                             "", "", "", "", "",//35->39
  40                                             "NORM_ERROR"};
  41
  42   std::map<NormalizedCellType,CellModel> CellModel::_map_of_unique_instance;
  43
  44   const CellModel& CellModel::GetCellModel(NormalizedCellType type)
  45   {
  46     if(_map_of_unique_instance.empty())
  47       buildUniqueInstance();
  48     const std::map<NormalizedCellType,CellModel>::iterator iter=_map_of_unique_instance.find(type);
  49     if(iter==_map_of_unique_instance.end())
  50       {
  51         std::ostringstream stream; stream << "no cellmodel for normalized type " << type;
  52         throw Exception(stream.str().c_str());
  53       }
  54     return (*iter).second;
  55   }
  56
  57   const char *CellModel::getRepr() const
  58   {
  59     return CELL_TYPES_REPR[(int)_type];
  60   }
  61
  62   /*!
  63    * This method is compatible with all types including dynamic one.
  64    */
  65   bool CellModel::isCompatibleWith(NormalizedCellType type) const
  66   {
  67     if(_type==type)
  68       return true;
  69     const CellModel& other=GetCellModel(type);
  70     if(_dim!=other.getDimension())
  71       return false;
  72     bool b1=isQuadratic();
  73     bool b2=other.isQuadratic();
  74     if((b1 && !b2) || (!b1 && b2))
  75       return false;
  76     b1=isDynamic();
  77     b2=other.isDynamic();
  78     return b1 || b2;
  79   }
  80
  81   void CellModel::buildUniqueInstance()
  82   {
  83     _map_of_unique_instance.insert(std::make_pair(NORM_POINT1,CellModel(NORM_POINT1)));
  84     _map_of_unique_instance.insert(std::make_pair(NORM_SEG2,CellModel(NORM_SEG2)));
  85     _map_of_unique_instance.insert(std::make_pair(NORM_SEG3,CellModel(NORM_SEG3)));
  86     _map_of_unique_instance.insert(std::make_pair(NORM_SEG4,CellModel(NORM_SEG4)));
  87     _map_of_unique_instance.insert(std::make_pair(NORM_TRI3,CellModel(NORM_TRI3)));
  88     _map_of_unique_instance.insert(std::make_pair(NORM_QUAD4,CellModel(NORM_QUAD4)));
  89     _map_of_unique_instance.insert(std::make_pair(NORM_TRI6,CellModel(NORM_TRI6)));
  90     _map_of_unique_instance.insert(std::make_pair(NORM_TRI7,CellModel(NORM_TRI7)));
  91     _map_of_unique_instance.insert(std::make_pair(NORM_QUAD8,CellModel(NORM_QUAD8)));
  92     _map_of_unique_instance.insert(std::make_pair(NORM_QUAD9,CellModel(NORM_QUAD9)));
  93     _map_of_unique_instance.insert(std::make_pair(NORM_TETRA4,CellModel(NORM_TETRA4)));
  94     _map_of_unique_instance.insert(std::make_pair(NORM_HEXA8,CellModel(NORM_HEXA8)));
  95     _map_of_unique_instance.insert(std::make_pair(NORM_PYRA5,CellModel(NORM_PYRA5)));
  96     _map_of_unique_instance.insert(std::make_pair(NORM_PENTA6,CellModel(NORM_PENTA6)));
  97     _map_of_unique_instance.insert(std::make_pair(NORM_TETRA10,CellModel(NORM_TETRA10)));
  98     _map_of_unique_instance.insert(std::make_pair(NORM_HEXGP12,CellModel(NORM_HEXGP12)));
  99     _map_of_unique_instance.insert(std::make_pair(NORM_PYRA13,CellModel(NORM_PYRA13)));
 100     _map_of_unique_instance.insert(std::make_pair(NORM_PENTA15,CellModel(NORM_PENTA15)));
 101     _map_of_unique_instance.insert(std::make_pair(NORM_HEXA20,CellModel(NORM_HEXA20)));
 102     _map_of_unique_instance.insert(std::make_pair(NORM_HEXA27,CellModel(NORM_HEXA27)));
 103     _map_of_unique_instance.insert(std::make_pair(NORM_POLYGON,CellModel(NORM_POLYGON)));
 104     _map_of_unique_instance.insert(std::make_pair(NORM_POLYHED,CellModel(NORM_POLYHED)));
 105     _map_of_unique_instance.insert(std::make_pair(NORM_QPOLYG,CellModel(NORM_QPOLYG)));
 106     _map_of_unique_instance.insert(std::make_pair(NORM_POLYL,CellModel(NORM_POLYL)));
 107     _map_of_unique_instance.insert(std::make_pair(NORM_ERROR,CellModel(NORM_ERROR)));
 108   }
 109
 110   CellModel::CellModel(NormalizedCellType type):_type(type)
 111   {
 112     _is_extruded=false;
 113     _quadratic=false;
 114     _dyn=false;
 115     _extruded_type=NORM_ERROR;
 116     _reverse_extruded_type=NORM_ERROR;
 117     _linear_type=NORM_ERROR;
 118     _quadratic_type=NORM_ERROR;
 119     _nb_of_little_sons=std::numeric_limits<unsigned>::max();
 120     switch(type)
 121       {
 122       case NORM_POINT1:
 123         {
 124           _nb_of_pts=1; _nb_of_sons=0; _dim=0; _extruded_type=NORM_SEG2; _is_simplex=true;
 125         }
 126         break;
 127       case NORM_SEG2:
 128         {
 129           _nb_of_pts=2; _nb_of_sons=2; _dim=1; _extruded_type=NORM_QUAD4; _quadratic_type=NORM_SEG3; _is_simplex=true; _is_extruded=true; _reverse_extruded_type=NORM_POINT1;
 130           _sons_type[0]=NORM_POINT1; _sons_type[1]=NORM_POINT1;
 131           _sons_con[0][0]=0; _nb_of_sons_con[0]=1;
 132           _sons_con[1][0]=1; _nb_of_sons_con[1]=1;
 133         }
 134         break;
 135       case NORM_SEG3:
 136         {
 137           _nb_of_pts=3; _nb_of_sons=3; _dim=1; _extruded_type=NORM_QUAD8; _linear_type=NORM_SEG2; _quadratic=true; _is_simplex=false;
 138           _sons_type[0]=NORM_POINT1; _sons_type[1]=NORM_POINT1; _sons_type[2]=NORM_POINT1;
 139           _sons_con[0][0]=0; _nb_of_sons_con[0]=1;
 140           _sons_con[1][0]=1; _nb_of_sons_con[1]=1;
 141           _sons_con[2][0]=2; _nb_of_sons_con[2]=1;
 142         }
 143         break;
 144       case NORM_SEG4:
 145         {
 146           _nb_of_pts=4; _nb_of_sons=4; _dim=1; _linear_type=NORM_SEG2; _quadratic=true; _is_simplex=false; // no _extruded_type because no cubic 2D cell
 147           _sons_type[0]=NORM_POINT1; _sons_type[1]=NORM_POINT1; _sons_type[2]=NORM_POINT1; _sons_type[3]=NORM_POINT1;
 148           _sons_con[0][0]=0; _nb_of_sons_con[0]=1;
 149           _sons_con[1][0]=1; _nb_of_sons_con[1]=1;
 150           _sons_con[2][0]=2; _nb_of_sons_con[2]=1;
 151           _sons_con[3][0]=3; _nb_of_sons_con[3]=1;
 152         }
 153         break;
 154       case NORM_TETRA4:
 155         {
 156           _nb_of_pts=4; _nb_of_sons=4; _dim=3; _quadratic_type=NORM_TETRA10; _is_simplex=true;
 157           _sons_type[0]=NORM_TRI3; _sons_type[1]=NORM_TRI3; _sons_type[2]=NORM_TRI3; _sons_type[3]=NORM_TRI3;
 158           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=2; _nb_of_sons_con[0]=3;
 159           _sons_con[1][0]=0; _sons_con[1][1]=3; _sons_con[1][2]=1; _nb_of_sons_con[1]=3;
 160           _sons_con[2][0]=1; _sons_con[2][1]=3; _sons_con[2][2]=2; _nb_of_sons_con[2]=3;
 161           _sons_con[3][0]=2; _sons_con[3][1]=3; _sons_con[3][2]=0; _nb_of_sons_con[3]=3;
 162           _little_sons_con[0][0]=0; _little_sons_con[0][1]=1;  _nb_of_little_sons=6;
 163           _little_sons_con[1][0]=1; _little_sons_con[1][1]=2;
 164           _little_sons_con[2][0]=2; _little_sons_con[2][1]=0;
 165           _little_sons_con[3][0]=0; _little_sons_con[3][1]=3;
 166           _little_sons_con[4][0]=1; _little_sons_con[4][1]=3;
 167           _little_sons_con[5][0]=2; _little_sons_con[5][1]=3;
 168         }
 169         break;
 170       case NORM_HEXA8:
 171         {
 172           _nb_of_pts=8; _nb_of_sons=6; _dim=3; _quadratic_type=NORM_HEXA20; _is_simplex=false; _is_extruded=true; _reverse_extruded_type=NORM_QUAD4;
 173           _sons_type[0]=NORM_QUAD4; _sons_type[1]=NORM_QUAD4; _sons_type[2]=NORM_QUAD4; _sons_type[3]=NORM_QUAD4; _sons_type[4]=NORM_QUAD4; _sons_type[5]=NORM_QUAD4;
 174           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=2; _sons_con[0][3]=3; _nb_of_sons_con[0]=4;
 175           _sons_con[1][0]=4; _sons_con[1][1]=7; _sons_con[1][2]=6; _sons_con[1][3]=5; _nb_of_sons_con[1]=4;
 176           _sons_con[2][0]=0; _sons_con[2][1]=4; _sons_con[2][2]=5; _sons_con[2][3]=1; _nb_of_sons_con[2]=4;
 177           _sons_con[3][0]=1; _sons_con[3][1]=5; _sons_con[3][2]=6; _sons_con[3][3]=2; _nb_of_sons_con[3]=4;
 178           _sons_con[4][0]=2; _sons_con[4][1]=6; _sons_con[4][2]=7; _sons_con[4][3]=3; _nb_of_sons_con[4]=4;
 179           _sons_con[5][0]=3; _sons_con[5][1]=7; _sons_con[5][2]=4; _sons_con[5][3]=0; _nb_of_sons_con[5]=4;
 180           _little_sons_con[0][0]=0; _little_sons_con[0][1]=1;  _nb_of_little_sons=12;
 181           _little_sons_con[1][0]=1; _little_sons_con[1][1]=2;
 182           _little_sons_con[2][0]=2; _little_sons_con[2][1]=3;
 183           _little_sons_con[3][0]=3; _little_sons_con[3][1]=0;
 184           _little_sons_con[4][0]=4; _little_sons_con[4][1]=5;
 185           _little_sons_con[5][0]=5; _little_sons_con[5][1]=6;
 186           _little_sons_con[6][0]=6; _little_sons_con[6][1]=7;
 187           _little_sons_con[7][0]=7; _little_sons_con[7][1]=4;
 188           _little_sons_con[8][0]=0; _little_sons_con[8][1]=4;
 189           _little_sons_con[9][0]=1; _little_sons_con[9][1]=5;
 190           _little_sons_con[10][0]=2; _little_sons_con[10][1]=6;
 191           _little_sons_con[11][0]=3; _little_sons_con[11][1]=7;
 192         }
 193         break;
 194       case NORM_QUAD4:
 195         {
 196           _nb_of_pts=4; _nb_of_sons=4; _dim=2; _quadratic_type=NORM_QUAD8; _is_simplex=false; _is_extruded=true;
 197           _sons_type[0]=NORM_SEG2; _sons_type[1]=NORM_SEG2; _sons_type[2]=NORM_SEG2; _sons_type[3]=NORM_SEG2;
 198           _sons_con[0][0]=0; _sons_con[0][1]=1; _nb_of_sons_con[0]=2;
 199           _sons_con[1][0]=1; _sons_con[1][1]=2; _nb_of_sons_con[1]=2;
 200           _sons_con[2][0]=2; _sons_con[2][1]=3; _nb_of_sons_con[2]=2;
 201           _sons_con[3][0]=3; _sons_con[3][1]=0; _nb_of_sons_con[3]=2; _extruded_type=NORM_HEXA8;
 202         }
 203         break;
 204       case NORM_TRI3:
 205         {
 206           _nb_of_pts=3; _nb_of_sons=3; _dim=2; _quadratic_type=NORM_TRI6; _is_simplex=true;
 207           _sons_type[0]=NORM_SEG2; _sons_type[1]=NORM_SEG2; _sons_type[2]=NORM_SEG2;
 208           _sons_con[0][0]=0; _sons_con[0][1]=1; _nb_of_sons_con[0]=2;
 209           _sons_con[1][0]=1; _sons_con[1][1]=2; _nb_of_sons_con[1]=2;
 210           _sons_con[2][0]=2; _sons_con[2][1]=0; _nb_of_sons_con[2]=2; _extruded_type=NORM_PENTA6;
 211         }
 212         break;
 213       case NORM_TRI6:
 214         {
 215           _nb_of_pts=6; _nb_of_sons=3; _dim=2; _linear_type=NORM_TRI3; _is_simplex=false;
 216           _sons_type[0]=NORM_SEG3; _sons_type[1]=NORM_SEG3; _sons_type[2]=NORM_SEG3;
 217           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=3; _nb_of_sons_con[0]=3;
 218           _sons_con[1][0]=1; _sons_con[1][1]=2; _sons_con[1][2]=4; _nb_of_sons_con[1]=3;
 219           _sons_con[2][0]=2; _sons_con[2][1]=0; _sons_con[2][2]=5; _nb_of_sons_con[2]=3; _quadratic=true; _extruded_type=NORM_PENTA15;
 220         }
 221         break;
 222       case NORM_TRI7:
 223         {
 224           _nb_of_pts=7; _nb_of_sons=3; _dim=2; _linear_type=NORM_TRI3; _is_simplex=false;
 225           _sons_type[0]=NORM_SEG3; _sons_type[1]=NORM_SEG3; _sons_type[2]=NORM_SEG3;
 226           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=3; _nb_of_sons_con[0]=3;
 227           _sons_con[1][0]=1; _sons_con[1][1]=2; _sons_con[1][2]=4; _nb_of_sons_con[1]=3;
 228           _sons_con[2][0]=2; _sons_con[2][1]=0; _sons_con[2][2]=5; _nb_of_sons_con[2]=3; _quadratic=true; //no extruded type because no penta20
 229         }
 230         break;
 231       case NORM_QUAD8:
 232         {
 233           _nb_of_pts=8; _nb_of_sons=4; _dim=2; _linear_type=NORM_QUAD4; _is_simplex=false;
 234           _sons_type[0]=NORM_SEG3; _sons_type[1]=NORM_SEG3; _sons_type[2]=NORM_SEG3; _sons_type[3]=NORM_SEG3;
 235           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=4; _nb_of_sons_con[0]=3;
 236           _sons_con[1][0]=1; _sons_con[1][1]=2; _sons_con[1][2]=5; _nb_of_sons_con[1]=3;
 237           _sons_con[2][0]=2; _sons_con[2][1]=3; _sons_con[2][2]=6; _nb_of_sons_con[2]=3;
 238           _sons_con[3][0]=3; _sons_con[3][1]=0; _sons_con[3][2]=7; _nb_of_sons_con[3]=3; _quadratic=true; _extruded_type=NORM_HEXA20;
 239         }
 240         break;
 241       case NORM_QUAD9:
 242         {
 243           _nb_of_pts=9; _nb_of_sons=4; _dim=2; _linear_type=NORM_QUAD4; _is_simplex=false;
 244           _sons_type[0]=NORM_SEG3; _sons_type[1]=NORM_SEG3; _sons_type[2]=NORM_SEG3; _sons_type[3]=NORM_SEG3;
 245           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=4; _nb_of_sons_con[0]=3;
 246           _sons_con[1][0]=1; _sons_con[1][1]=2; _sons_con[1][2]=5; _nb_of_sons_con[1]=3;
 247           _sons_con[2][0]=2; _sons_con[2][1]=3; _sons_con[2][2]=6; _nb_of_sons_con[2]=3;
 248           _sons_con[3][0]=3; _sons_con[3][1]=0; _sons_con[3][2]=7; _nb_of_sons_con[3]=3; _quadratic=true; _extruded_type=NORM_HEXA27;
 249         }
 250         break;
 251       case NORM_PYRA5:
 252         {
 253           _nb_of_pts=5; _nb_of_sons=5; _dim=3; _quadratic_type=NORM_PYRA13; _is_simplex=false;
 254           _sons_type[0]=NORM_QUAD4; _sons_type[1]=NORM_TRI3; _sons_type[2]=NORM_TRI3; _sons_type[3]=NORM_TRI3; _sons_type[4]=NORM_TRI3;
 255           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=2; _sons_con[0][3]=3; _nb_of_sons_con[0]=4;
 256           _sons_con[1][0]=0; _sons_con[1][1]=4; _sons_con[1][2]=1; _nb_of_sons_con[1]=3;
 257           _sons_con[2][0]=1; _sons_con[2][1]=4; _sons_con[2][2]=2; _nb_of_sons_con[2]=3;
 258           _sons_con[3][0]=2; _sons_con[3][1]=4; _sons_con[3][2]=3; _nb_of_sons_con[3]=3;
 259           _sons_con[4][0]=3; _sons_con[4][1]=4; _sons_con[4][2]=0; _nb_of_sons_con[4]=3;
 260           _little_sons_con[0][0]=0; _little_sons_con[0][1]=1;  _nb_of_little_sons=8;
 261           _little_sons_con[1][0]=1; _little_sons_con[1][1]=2;
 262           _little_sons_con[2][0]=2; _little_sons_con[2][1]=3;
 263           _little_sons_con[3][0]=3; _little_sons_con[3][1]=0;
 264           _little_sons_con[4][0]=0; _little_sons_con[4][1]=4;
 265           _little_sons_con[5][0]=1; _little_sons_con[5][1]=4;
 266           _little_sons_con[6][0]=2; _little_sons_con[6][1]=4;
 267           _little_sons_con[7][0]=3; _little_sons_con[7][1]=4;
 268         }
 269         break;
 270       case NORM_PENTA6:
 271         {
 272           _nb_of_pts=6; _nb_of_sons=5; _dim=3; _quadratic_type=NORM_PENTA15; _is_simplex=false; _is_extruded=true; _reverse_extruded_type=NORM_TRI3;
 273           _sons_type[0]=NORM_TRI3; _sons_type[1]=NORM_TRI3; _sons_type[2]=NORM_QUAD4; _sons_type[3]=NORM_QUAD4; _sons_type[4]=NORM_QUAD4;
 274           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=2; _nb_of_sons_con[0]=3;
 275           _sons_con[1][0]=3; _sons_con[1][1]=5; _sons_con[1][2]=4; _nb_of_sons_con[1]=3;
 276           _sons_con[2][0]=0; _sons_con[2][1]=3; _sons_con[2][2]=4; _sons_con[2][3]=1; _nb_of_sons_con[2]=4;
 277           _sons_con[3][0]=1; _sons_con[3][1]=4; _sons_con[3][2]=5; _sons_con[3][3]=2; _nb_of_sons_con[3]=4;
 278           _sons_con[4][0]=2; _sons_con[4][1]=5; _sons_con[4][2]=3; _sons_con[4][3]=0; _nb_of_sons_con[4]=4;
 279           _little_sons_con[0][0]=0; _little_sons_con[0][1]=1;  _nb_of_little_sons=9;
 280           _little_sons_con[1][0]=1; _little_sons_con[1][1]=2;
 281           _little_sons_con[2][0]=2; _little_sons_con[2][1]=0;
 282           _little_sons_con[3][0]=3; _little_sons_con[3][1]=4;
 283           _little_sons_con[4][0]=4; _little_sons_con[4][1]=5;
 284           _little_sons_con[5][0]=5; _little_sons_con[5][1]=3;
 285           _little_sons_con[6][0]=0; _little_sons_con[6][1]=3;
 286           _little_sons_con[7][0]=1; _little_sons_con[7][1]=4;
 287           _little_sons_con[8][0]=2; _little_sons_con[8][1]=5;
 288         }
 289         break;
 290       case NORM_TETRA10:
 291         {
 292           _nb_of_pts=10; _nb_of_sons=4; _dim=3; _linear_type=NORM_TETRA4; _is_simplex=false;
 293           _sons_type[0]=NORM_TRI6; _sons_type[1]=NORM_TRI6; _sons_type[2]=NORM_TRI6; _sons_type[3]=NORM_TRI6;
 294           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=2; _sons_con[0][3]=4; _sons_con[0][4]=5; _sons_con[0][5]=6; _nb_of_sons_con[0]=6;
 295           _sons_con[1][0]=0; _sons_con[1][1]=3; _sons_con[1][2]=1; _sons_con[1][3]=7; _sons_con[1][4]=8; _sons_con[1][5]=4; _nb_of_sons_con[1]=6;
 296           _sons_con[2][0]=1; _sons_con[2][1]=3; _sons_con[2][2]=2; _sons_con[2][3]=8; _sons_con[2][4]=9; _sons_con[2][5]=5; _nb_of_sons_con[2]=6;
 297           _sons_con[3][0]=2; _sons_con[3][1]=3; _sons_con[3][2]=0; _sons_con[3][3]=9; _sons_con[3][4]=7; _sons_con[3][5]=6; _nb_of_sons_con[3]=6;  _quadratic=true;
 298           _little_sons_con[0][0]=0; _little_sons_con[0][1]=1;  _little_sons_con[0][2]=4;  _nb_of_little_sons=6;
 299           _little_sons_con[1][0]=1; _little_sons_con[1][1]=2;  _little_sons_con[1][2]=5;
 300           _little_sons_con[2][0]=2; _little_sons_con[2][1]=0;  _little_sons_con[2][2]=6;
 301           _little_sons_con[3][0]=0; _little_sons_con[3][1]=3;  _little_sons_con[3][2]=7;
 302           _little_sons_con[4][0]=1; _little_sons_con[4][1]=3;  _little_sons_con[4][2]=8;
 303           _little_sons_con[5][0]=2; _little_sons_con[5][1]=3;  _little_sons_con[5][2]=9;
 304         }
 305         break;
 306       case NORM_HEXGP12:
 307         {
 308           _nb_of_pts=12; _nb_of_sons=8; _dim=3; _is_simplex=false; _is_extruded=true;
 309           _sons_type[0]=NORM_POLYGON; _sons_type[1]=NORM_POLYGON; _sons_type[2]=NORM_QUAD4; _sons_type[3]=NORM_QUAD4; _sons_type[4]=NORM_QUAD4; _sons_type[5]=NORM_QUAD4;
 310           _sons_type[6]=NORM_QUAD4; _sons_type[7]=NORM_QUAD4;
 311           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=2; _sons_con[0][3]=3; _sons_con[0][4]=4; _sons_con[0][5]=5; _nb_of_sons_con[0]=6;
 312           _sons_con[1][0]=6; _sons_con[1][1]=11; _sons_con[1][2]=10; _sons_con[1][3]=9; _sons_con[1][4]=8; _sons_con[1][5]=7; _nb_of_sons_con[1]=6;
 313           _sons_con[2][0]=0; _sons_con[2][1]=6; _sons_con[2][2]=7; _sons_con[2][3]=1; _nb_of_sons_con[2]=4;
 314           _sons_con[3][0]=1; _sons_con[3][1]=7; _sons_con[3][2]=8; _sons_con[3][3]=2; _nb_of_sons_con[3]=4;
 315           _sons_con[4][0]=2; _sons_con[4][1]=8; _sons_con[4][2]=9; _sons_con[4][3]=3; _nb_of_sons_con[4]=4;
 316           _sons_con[5][0]=3; _sons_con[5][1]=9; _sons_con[5][2]=10; _sons_con[5][3]=4; _nb_of_sons_con[5]=4;
 317           _sons_con[6][0]=4; _sons_con[6][1]=10; _sons_con[6][2]=11; _sons_con[6][3]=5; _nb_of_sons_con[6]=4;
 318           _sons_con[7][0]=5; _sons_con[7][1]=11; _sons_con[7][2]=6; _sons_con[7][3]=0; _nb_of_sons_con[7]=4;
 319         }
 320         break;
 321       case NORM_PYRA13:
 322         {
 323           _nb_of_pts=13; _nb_of_sons=5; _dim=3; _linear_type=NORM_PYRA5; _is_simplex=false;
 324           _sons_type[0]=NORM_QUAD8; _sons_type[1]=NORM_TRI6; _sons_type[2]=NORM_TRI6; _sons_type[3]=NORM_TRI6; _sons_type[4]=NORM_TRI6;
 325           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=2; _sons_con[0][3]=3; _sons_con[0][4]=5; _sons_con[0][5]=6; _sons_con[0][6]=7; _sons_con[0][7]=8; _nb_of_sons_con[0]=8;
 326           _sons_con[1][0]=0; _sons_con[1][1]=4; _sons_con[1][2]=1; _sons_con[1][3]=9; _sons_con[1][4]=10; _sons_con[1][5]=5; _nb_of_sons_con[1]=6;
 327           _sons_con[2][0]=1; _sons_con[2][1]=4; _sons_con[2][2]=2; _sons_con[2][3]=10; _sons_con[2][4]=11; _sons_con[2][5]=6; _nb_of_sons_con[2]=6;
 328           _sons_con[3][0]=2; _sons_con[3][1]=4; _sons_con[3][2]=3; _sons_con[3][3]=11; _sons_con[3][4]=12; _sons_con[3][5]=7;  _nb_of_sons_con[3]=6;
 329           _sons_con[4][0]=3; _sons_con[4][1]=4; _sons_con[4][2]=0; _sons_con[4][3]=12; _sons_con[4][4]=9; _sons_con[4][5]=8; _nb_of_sons_con[4]=6; _quadratic=true;
 330           _little_sons_con[0][0]=0; _little_sons_con[0][1]=1; _little_sons_con[0][2]=5;  _nb_of_little_sons=8;
 331           _little_sons_con[1][0]=1; _little_sons_con[1][1]=2; _little_sons_con[1][2]=6;
 332           _little_sons_con[2][0]=2; _little_sons_con[2][1]=3; _little_sons_con[2][2]=7;
 333           _little_sons_con[3][0]=3; _little_sons_con[3][1]=0; _little_sons_con[3][2]=8;
 334           _little_sons_con[4][0]=0; _little_sons_con[4][1]=4; _little_sons_con[4][2]=9;
 335           _little_sons_con[5][0]=1; _little_sons_con[5][1]=4; _little_sons_con[5][2]=10;
 336           _little_sons_con[6][0]=2; _little_sons_con[6][1]=4; _little_sons_con[6][2]=11;
 337           _little_sons_con[7][0]=3; _little_sons_con[7][1]=4; _little_sons_con[7][2]=12;
 338         }
 339         break;
 340       case NORM_PENTA15:
 341         {
 342           _nb_of_pts=15; _nb_of_sons=5; _dim=3; _linear_type=NORM_PENTA6; _is_simplex=false;
 343           _sons_type[0]=NORM_TRI6; _sons_type[1]=NORM_TRI6; _sons_type[2]=NORM_QUAD8; _sons_type[3]=NORM_QUAD8; _sons_type[4]=NORM_QUAD8;
 344           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=2; _sons_con[0][3]=6; _sons_con[0][4]=7; _sons_con[0][5]=8; _nb_of_sons_con[0]=6;
 345           _sons_con[1][0]=3; _sons_con[1][1]=5; _sons_con[1][2]=4; _sons_con[1][3]=11; _sons_con[1][4]=10; _sons_con[1][5]=9; _nb_of_sons_con[1]=6;
 346           _sons_con[2][0]=0; _sons_con[2][1]=3; _sons_con[2][2]=4; _sons_con[2][3]=1; _sons_con[2][4]=12; _sons_con[2][5]=9; _sons_con[2][6]=13; _sons_con[2][7]=6; _nb_of_sons_con[2]=8;
 347           _sons_con[3][0]=1; _sons_con[3][1]=4; _sons_con[3][2]=5; _sons_con[3][3]=2; _sons_con[3][4]=13; _sons_con[3][5]=10; _sons_con[3][6]=14; _sons_con[3][7]=7; _nb_of_sons_con[3]=8;
 348           _sons_con[4][0]=2; _sons_con[4][1]=4; _sons_con[4][2]=5; _sons_con[4][3]=0; _sons_con[4][4]=14; _sons_con[4][5]=11; _sons_con[4][6]=12; _sons_con[4][7]=8; _nb_of_sons_con[4]=8; _quadratic=true;
 349           _little_sons_con[0][0]=0; _little_sons_con[0][1]=1; _little_sons_con[0][2]=6;  _nb_of_little_sons=9;
 350           _little_sons_con[1][0]=1; _little_sons_con[1][1]=2; _little_sons_con[1][2]=7;
 351           _little_sons_con[2][0]=2; _little_sons_con[2][1]=0; _little_sons_con[2][2]=8;
 352           _little_sons_con[3][0]=3; _little_sons_con[3][1]=4; _little_sons_con[3][2]=9;
 353           _little_sons_con[4][0]=4; _little_sons_con[4][1]=5; _little_sons_con[4][2]=10;
 354           _little_sons_con[5][0]=5; _little_sons_con[5][1]=3; _little_sons_con[5][2]=11;
 355           _little_sons_con[6][0]=0; _little_sons_con[6][1]=3; _little_sons_con[6][2]=12;
 356           _little_sons_con[7][0]=1; _little_sons_con[7][1]=4; _little_sons_con[7][2]=13;
 357           _little_sons_con[8][0]=2; _little_sons_con[8][1]=5; _little_sons_con[8][2]=14;
 358         }
 359         break;
 360       case NORM_HEXA20:
 361         {
 362           _nb_of_pts=20; _nb_of_sons=6; _dim=3; _linear_type=NORM_HEXA8; _is_simplex=false;
 363           _sons_type[0]=NORM_QUAD8; _sons_type[1]=NORM_QUAD8; _sons_type[2]=NORM_QUAD8; _sons_type[3]=NORM_QUAD8; _sons_type[4]=NORM_QUAD8; _sons_type[5]=NORM_QUAD8;
 364           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=2; _sons_con[0][3]=3; _sons_con[0][4]=8; _sons_con[0][5]=9; _sons_con[0][6]=10; _sons_con[0][7]=11; _nb_of_sons_con[0]=8;
 365           _sons_con[1][0]=4; _sons_con[1][1]=7; _sons_con[1][2]=6; _sons_con[1][3]=5; _sons_con[1][4]=15; _sons_con[1][5]=14; _sons_con[1][6]=13; _sons_con[1][7]=12; _nb_of_sons_con[1]=8;
 366           _sons_con[2][0]=0; _sons_con[2][1]=4; _sons_con[2][2]=5; _sons_con[2][3]=1; _sons_con[2][4]=16; _sons_con[2][5]=12; _sons_con[2][6]=17; _sons_con[2][7]=8; _nb_of_sons_con[2]=8;
 367           _sons_con[3][0]=1; _sons_con[3][1]=5; _sons_con[3][3]=6; _sons_con[3][3]=2; _sons_con[3][4]=17; _sons_con[3][5]=13; _sons_con[3][6]=18; _sons_con[3][7]=9;_nb_of_sons_con[3]=8;
 368           _sons_con[4][0]=2; _sons_con[4][1]=6; _sons_con[4][3]=7; _sons_con[4][3]=3; _sons_con[4][4]=18; _sons_con[4][5]=14; _sons_con[4][6]=19; _sons_con[4][7]=10; _nb_of_sons_con[4]=8;
 369           _sons_con[5][0]=3; _sons_con[5][1]=7; _sons_con[5][3]=4; _sons_con[5][3]=0; _sons_con[5][4]=19; _sons_con[5][5]=15; _sons_con[5][6]=16; _sons_con[5][7]=11; _nb_of_sons_con[5]=8; _quadratic=true;
 370           _little_sons_con[0][0]=0; _little_sons_con[0][1]=1;  _little_sons_con[0][2]=8; _nb_of_little_sons=12;
 371           _little_sons_con[1][0]=1; _little_sons_con[1][1]=2;  _little_sons_con[1][2]=9;
 372           _little_sons_con[2][0]=2; _little_sons_con[2][1]=3;  _little_sons_con[2][2]=10;
 373           _little_sons_con[3][0]=3; _little_sons_con[3][1]=0;  _little_sons_con[3][2]=11;
 374           _little_sons_con[4][0]=4; _little_sons_con[4][1]=5;  _little_sons_con[4][2]=12;
 375           _little_sons_con[5][0]=5; _little_sons_con[5][1]=6;  _little_sons_con[5][2]=13;
 376           _little_sons_con[6][0]=6; _little_sons_con[6][1]=7;  _little_sons_con[6][2]=14;
 377           _little_sons_con[7][0]=7; _little_sons_con[7][1]=4;  _little_sons_con[7][2]=15;
 378           _little_sons_con[8][0]=0; _little_sons_con[8][1]=4;  _little_sons_con[8][2]=16;
 379           _little_sons_con[9][0]=1; _little_sons_con[9][1]=5;  _little_sons_con[9][2]=17;
 380           _little_sons_con[10][0]=2; _little_sons_con[10][1]=6;  _little_sons_con[10][2]=18;
 381           _little_sons_con[11][0]=3; _little_sons_con[11][1]=7;  _little_sons_con[11][2]=19;
 382         }
 383         break;
 384       case NORM_HEXA27:
 385         {
 386           _nb_of_pts=27; _nb_of_sons=6; _dim=3; _linear_type=NORM_HEXA8; _is_simplex=false;
 387           _sons_type[0]=NORM_QUAD9; _sons_type[1]=NORM_QUAD9; _sons_type[2]=NORM_QUAD9; _sons_type[3]=NORM_QUAD9; _sons_type[4]=NORM_QUAD9; _sons_type[5]=NORM_QUAD9;
 388           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=2; _sons_con[0][3]=3; _sons_con[0][4]=8; _sons_con[0][5]=9; _sons_con[0][6]=10; _sons_con[0][7]=11; _sons_con[0][8]=20; _nb_of_sons_con[0]=9;
 389           _sons_con[1][0]=4; _sons_con[1][1]=7; _sons_con[1][2]=6; _sons_con[1][3]=5; _sons_con[1][4]=15; _sons_con[1][5]=14; _sons_con[1][6]=13; _sons_con[1][7]=12; _sons_con[1][8]=25; _nb_of_sons_con[1]=9;
 390           _sons_con[2][0]=0; _sons_con[2][1]=4; _sons_con[2][2]=5; _sons_con[2][3]=1; _sons_con[2][4]=16; _sons_con[2][5]=12; _sons_con[2][6]=17; _sons_con[2][7]=8; _sons_con[2][8]=21; _nb_of_sons_con[2]=9;
 391           _sons_con[3][0]=1; _sons_con[3][1]=5; _sons_con[3][3]=6; _sons_con[3][3]=2; _sons_con[3][4]=17; _sons_con[3][5]=13; _sons_con[3][6]=18; _sons_con[3][7]=9; _sons_con[3][8]=22; _nb_of_sons_con[3]=9;
 392           _sons_con[4][0]=2; _sons_con[4][1]=6; _sons_con[4][3]=7; _sons_con[4][3]=3; _sons_con[4][4]=18; _sons_con[4][5]=14; _sons_con[4][6]=19; _sons_con[4][7]=10; _sons_con[4][8]=23; _nb_of_sons_con[4]=9;
 393           _sons_con[5][0]=3; _sons_con[5][1]=7; _sons_con[5][3]=4; _sons_con[5][3]=0; _sons_con[5][4]=19; _sons_con[5][5]=15; _sons_con[5][6]=16; _sons_con[5][7]=11; _sons_con[5][8]=24; _nb_of_sons_con[5]=9;
 394           _quadratic=true;
 395         }
 396         break;
 397       case NORM_POLYGON:
 398         {
 399           _nb_of_pts=0; _nb_of_sons=0; _dim=2; _dyn=true; _extruded_type=NORM_POLYHED; _is_simplex=false;
 400         }
 401         break;
 402       case NORM_POLYHED:
 403         {
 404           _nb_of_pts=0; _nb_of_sons=0; _dim=3; _dyn=true; _is_simplex=false;
 405         }
 406         break;
 407       case NORM_QPOLYG:
 408         {
 409           _nb_of_pts=0; _nb_of_sons=0; _dim=2; _dyn=true; _is_simplex=false; _quadratic=true;
 410         }
 411         break;
 412       case NORM_POLYL:
 413         {
 414           _nb_of_pts=0; _nb_of_sons=0; _dim=1; _dyn=true; _extruded_type=NORM_POLYGON; _is_simplex=false;
 415         }
 416       case NORM_ERROR:
 417         {
 418           _nb_of_pts=std::numeric_limits<unsigned>::max(); _nb_of_sons=std::numeric_limits<unsigned>::max(); _dim=std::numeric_limits<unsigned>::max();
 419         }
 420         break;
 421       }
 422   }
 423
 424   /*!
 425    * Equivalent to getNumberOfSons except that this method deals with dynamic type.
 426    */
 427   unsigned CellModel::getNumberOfSons2(const int *conn, int lgth) const
 428   {
 429     if(!isDynamic())
 430       return getNumberOfSons();
 431     if(_dim==2)
 432       {
 433         if(_type==NORM_POLYGON)
 434           return lgth;
 435         else
 436           return lgth/2;
 437       }
 438     else if(_dim==1)
 439       return lgth;//NORM_POLYL
 440     else
 441       return std::count(conn,conn+lgth,-1)+1;
 442   }
 443
 444   unsigned CellModel::getNumberOfEdgesIn3D(const int *conn, int lgth) const
 445   {
 446     if(!isDynamic())
 447       return _nb_of_little_sons;
 448     else//polyhedron
 449       return (lgth-std::count(conn,conn+lgth,-1))/2;
 450   }
 451
 452   /*!
 453    * Equivalent to getSonType except that this method deals with dynamic type.
 454    */
 455   NormalizedCellType CellModel::getSonType2(unsigned sonId) const
 456   {
 457     if(!isDynamic())
 458       return getSonType(sonId);
 459     if(_dim==2)
 460       {
 461         if(_type==NORM_POLYGON)
 462           return NORM_SEG2;
 463         else
 464           return NORM_SEG3;
 465       }
 466     else if(_dim==1)
 467       return NORM_ERROR;//NORM_POLYL
 468     //polyedron
 469     return NORM_POLYGON;
 470   }
 471
 472   /*!
 473    * \b WARNING this method do not manage correctly types that return true at the call of isDynamic. Use fillSonCellNodalConnectivity2 instead.
 474    */
 475   unsigned CellModel::fillSonCellNodalConnectivity(int sonId, const int *nodalConn, int *sonNodalConn) const
 476   {
 477     unsigned nbOfTurnLoop=_nb_of_sons_con[sonId];
 478     const unsigned *sonConn=_sons_con[sonId];
 479     for(unsigned i=0;i<nbOfTurnLoop;i++)
 480       sonNodalConn[i]=nodalConn[sonConn[i]];
 481     return nbOfTurnLoop;
 482   }
 483
 484   unsigned CellModel::fillSonCellNodalConnectivity2(int sonId, const int *nodalConn, int lgth, int *sonNodalConn, NormalizedCellType& typeOfSon) const
 485   {
 486     typeOfSon=getSonType2(sonId);
 487     if(!isDynamic())
 488       return fillSonCellNodalConnectivity(sonId,nodalConn,sonNodalConn);
 489     else
 490       {
 491         if(_dim==2)//polygon
 492           {
 493             if(_type==NORM_POLYGON)
 494               {
 495                 sonNodalConn[0]=nodalConn[sonId];
 496                 sonNodalConn[1]=nodalConn[(sonId+1)%lgth];
 497                 return 2;
 498               }
 499             else
 500               {
 501                 sonNodalConn[0]=nodalConn[sonId];
 502                 sonNodalConn[1]=nodalConn[(sonId+1)%(lgth/2)];
 503                 sonNodalConn[2]=nodalConn[sonId+(lgth/2)];
 504                 return 3;
 505               }
 506           }
 507         else if(_dim==3)
 508           {//polyedron
 509             const int *where=nodalConn;
 510             for(int i=0;i<sonId;i++)
 511               {
 512                 where=std::find(where,nodalConn+lgth,-1);
 513                 where++;
 514               }
 515             const int *where2=std::find(where,nodalConn+lgth,-1);
 516             std::copy(where,where2,sonNodalConn);
 517             return where2-where;
 518           }
 519         else
 520           throw INTERP_KERNEL::Exception("CellModel::fillSonCellNodalConnectivity2 : no sons on NORM_POLYL !");
 521       }
 522   }
 523
 524   unsigned CellModel::fillSonEdgesNodalConnectivity3D(int sonId, const int *nodalConn, int lgth, int *sonNodalConn, NormalizedCellType& typeOfSon) const
 525   {
 526     if(!isDynamic())
 527       {
 528         if(!isQuadratic())
 529           {
 530             typeOfSon=NORM_SEG2;
 531             sonNodalConn[0]=nodalConn[_little_sons_con[sonId][0]];
 532             sonNodalConn[1]=nodalConn[_little_sons_con[sonId][1]];
 533             return 2;
 534           }
 535         else
 536           {
 537             typeOfSon=NORM_SEG3;
 538             sonNodalConn[0]=nodalConn[_little_sons_con[sonId][0]];
 539             sonNodalConn[1]=nodalConn[_little_sons_con[sonId][1]];
 540             sonNodalConn[2]=nodalConn[_little_sons_con[sonId][2]];
 541             return 3;
 542           }
 543       }
 544     else
 545       throw INTERP_KERNEL::Exception("CellModel::fillSonEdgesNodalConnectivity3D : not implemented yet for NORM_POLYHED !");
 546   }
 547
 548   //================================================================================
 549   /*!
 550    * \brief Return number of nodes in sonId-th son of a Dynamic() cell
 551    */
 552   //================================================================================
 553
 554   unsigned CellModel::getNumberOfNodesConstituentTheSon2(unsigned sonId, const int *nodalConn, int lgth) const
 555   {
 556     if(!isDynamic())
 557       return getNumberOfNodesConstituentTheSon(sonId);
 558
 559     if(_dim==2)//polygon
 560       {
 561         if(_type==NORM_POLYGON)
 562           return 2;
 563         else
 564           return 3;
 565       }
 566     else if(_dim==3)
 567       {//polyedron
 568         const int *where=nodalConn;
 569         for(unsigned int i=0;i<sonId;i++)
 570           {
 571             where=std::find(where,nodalConn+lgth,-1);
 572             where++;
 573           }
 574         const int *where2=std::find(where,nodalConn+lgth,-1);
 575         return where2-where;
 576       }
 577     else
 578       throw INTERP_KERNEL::Exception("CellModel::getNumberOfNodesConstituentTheSon2 : no sons on NORM_POLYL !");
 579   }
 580
 581   /*!
 582    * This method retrieves if cell1 represented by 'conn1' and cell2 represented by 'conn2'
 583    * are equivalent by a permutation or not. This method expects to work on 1D or 2D (only mesh dimension where it is possible to have a spaceDim) strictly higher than meshDim.
 584    * If not an exception will be thrown.
 585    * @return True if two cells have same orientation, false if not.
 586    */
 587   bool CellModel::getOrientationStatus(unsigned lgth, const int *conn1, const int *conn2) const
 588   {
 589     if(_dim!=1 && _dim!=2)
 590       throw INTERP_KERNEL::Exception("CellModel::getOrientationStatus : invalid dimension ! Must be 1 or 2 !");
 591     if(!_quadratic)
 592       {
 593         std::vector<int> tmp(2*lgth);
 594         std::vector<int>::iterator it=std::copy(conn1,conn1+lgth,tmp.begin());
 595         std::copy(conn1,conn1+lgth,it);
 596         it=std::search(tmp.begin(),tmp.end(),conn2,conn2+lgth);
 597         if(it==tmp.begin())
 598           return true;
 599         if(it!=tmp.end())
 600           return _dim!=1;
 601         std::vector<int>::reverse_iterator it2=std::search(tmp.rbegin(),tmp.rend(),conn2,conn2+lgth);
 602         if(it2!=tmp.rend())
 603           return false;
 604         throw INTERP_KERNEL::Exception("CellModel::getOrientationStatus : Request of orientation status of non equal connectively cells !");
 605       }
 606     else
 607       {
 608         if(_dim!=1)
 609           {
 610             std::vector<int> tmp(lgth);
 611             std::vector<int>::iterator it=std::copy(conn1,conn1+lgth/2,tmp.begin());
 612             std::copy(conn1,conn1+lgth/2,it);
 613             it=std::search(tmp.begin(),tmp.end(),conn2,conn2+lgth/2);
 614             int d=std::distance(tmp.begin(),it);
 615             if(it==tmp.end())
 616               return false;
 617             it=std::copy(conn1+lgth/2,conn1+lgth,tmp.begin());
 618             std::copy(conn1+lgth/2,conn1+lgth,it);
 619             it=std::search(tmp.begin(),tmp.end(),conn2,conn2+lgth);
 620             if(it==tmp.end())
 621               return false;
 622             int d2=std::distance(tmp.begin(),it);
 623             return d==d2;
 624           }
 625         else
 626           {
 627             int p=(lgth+1)/2;
 628             std::vector<int> tmp(2*p);
 629             std::vector<int>::iterator it=std::copy(conn1,conn1+p,tmp.begin());
 630             std::copy(conn1,conn1+p,it);
 631             it=std::search(tmp.begin(),tmp.end(),conn2,conn2+p);
 632             int d=std::distance(tmp.begin(),it);
 633             if(it==tmp.end())
 634               return false;
 635             tmp.resize(2*p-2);
 636             it=std::copy(conn1+p,conn1+lgth,tmp.begin());
 637             std::copy(conn1+p,conn1+lgth,it);
 638             it=std::search(tmp.begin(),tmp.end(),conn2+p,conn2+lgth);
 639             if(it==tmp.end())
 640               return false;
 641             int d2=std::distance(tmp.begin(),it);
 642             return d==d2;
 643           }
 644       }
 645   }
 646
 647 }