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
  20 #include "CellModel.hxx"
  21
  22 #include "InterpKernelException.hxx"
  23
  24 #include <algorithm>
  25 #include <sstream>
  26 #include <vector>
  27 #include <limits>
  28
  29 namespace INTERP_KERNEL
  30 {
  31   const char *CellModel::CELL_TYPES_REPR[]={"NORM_POINT1", "NORM_SEG2", "NORM_SEG3", "NORM_TRI3", "NORM_QUAD4",// 0->4
  32                                             "NORM_POLYGON", "NORM_TRI6", "NORM_TRI7" , "NORM_QUAD8", "NORM_QUAD9",//5->9
  33                                             "NORM_SEG4", "", "", "", "NORM_TETRA4",//10->14
  34                                             "NORM_PYRA5", "NORM_PENTA6", "", "NORM_HEXA8", "",//15->19
  35                                             "NORM_TETRA10", "", "NORM_HEXGP12", "NORM_PYRA13", "",//20->24
  36                                             "NORM_PENTA15", "", "NORM_HEXA27", "", "",//25->29
  37                                             "NORM_HEXA20", "NORM_POLYHED", "NORM_QPOLYG", "NORM_POLYL", "",//30->34
  38                                             "", "", "", "", "",//35->39
  39                                             "NORM_ERROR"};
  40
  41   std::map<NormalizedCellType,CellModel> CellModel::_map_of_unique_instance;
  42
  43   const CellModel& CellModel::GetCellModel(NormalizedCellType type)
  44   {
  45     if(_map_of_unique_instance.empty())
  46       buildUniqueInstance();
  47     const std::map<NormalizedCellType,CellModel>::iterator iter=_map_of_unique_instance.find(type);
  48     if(iter==_map_of_unique_instance.end())
  49       {
  50         std::ostringstream stream; stream << "no cellmodel for normalized type " << type;
  51         throw Exception(stream.str().c_str());
  52       }
  53     return (*iter).second;
  54   }
  55
  56   const char *CellModel::getRepr() const
  57   {
  58     return CELL_TYPES_REPR[(int)_type];
  59   }
  60
  61   /*!
  62    * This method is compatible with all types including dynamic one.
  63    */
  64   bool CellModel::isCompatibleWith(NormalizedCellType type) const
  65   {
  66     if(_type==type)
  67       return true;
  68     const CellModel& other=GetCellModel(type);
  69     if(_dim!=other.getDimension())
  70       return false;
  71     bool b1=isQuadratic();
  72     bool b2=other.isQuadratic();
  73     if((b1 && !b2) || (!b1 && b2))
  74       return false;
  75     b1=isDynamic();
  76     b2=other.isDynamic();
  77     return b1 || b2;
  78   }
  79
  80   void CellModel::buildUniqueInstance()
  81   {
  82     _map_of_unique_instance.insert(std::make_pair(NORM_POINT1,CellModel(NORM_POINT1)));
  83     _map_of_unique_instance.insert(std::make_pair(NORM_SEG2,CellModel(NORM_SEG2)));
  84     _map_of_unique_instance.insert(std::make_pair(NORM_SEG3,CellModel(NORM_SEG3)));
  85     _map_of_unique_instance.insert(std::make_pair(NORM_SEG4,CellModel(NORM_SEG4)));
  86     _map_of_unique_instance.insert(std::make_pair(NORM_TRI3,CellModel(NORM_TRI3)));
  87     _map_of_unique_instance.insert(std::make_pair(NORM_QUAD4,CellModel(NORM_QUAD4)));
  88     _map_of_unique_instance.insert(std::make_pair(NORM_TRI6,CellModel(NORM_TRI6)));
  89     _map_of_unique_instance.insert(std::make_pair(NORM_TRI7,CellModel(NORM_TRI7)));
  90     _map_of_unique_instance.insert(std::make_pair(NORM_QUAD8,CellModel(NORM_QUAD8)));
  91     _map_of_unique_instance.insert(std::make_pair(NORM_QUAD9,CellModel(NORM_QUAD9)));
  92     _map_of_unique_instance.insert(std::make_pair(NORM_TETRA4,CellModel(NORM_TETRA4)));
  93     _map_of_unique_instance.insert(std::make_pair(NORM_HEXA8,CellModel(NORM_HEXA8)));
  94     _map_of_unique_instance.insert(std::make_pair(NORM_PYRA5,CellModel(NORM_PYRA5)));
  95     _map_of_unique_instance.insert(std::make_pair(NORM_PENTA6,CellModel(NORM_PENTA6)));
  96     _map_of_unique_instance.insert(std::make_pair(NORM_TETRA10,CellModel(NORM_TETRA10)));
  97     _map_of_unique_instance.insert(std::make_pair(NORM_HEXGP12,CellModel(NORM_HEXGP12)));
  98     _map_of_unique_instance.insert(std::make_pair(NORM_PYRA13,CellModel(NORM_PYRA13)));
  99     _map_of_unique_instance.insert(std::make_pair(NORM_PENTA15,CellModel(NORM_PENTA15)));
 100     _map_of_unique_instance.insert(std::make_pair(NORM_HEXA20,CellModel(NORM_HEXA20)));
 101     _map_of_unique_instance.insert(std::make_pair(NORM_HEXA27,CellModel(NORM_HEXA27)));
 102     _map_of_unique_instance.insert(std::make_pair(NORM_POLYGON,CellModel(NORM_POLYGON)));
 103     _map_of_unique_instance.insert(std::make_pair(NORM_POLYHED,CellModel(NORM_POLYHED)));
 104     _map_of_unique_instance.insert(std::make_pair(NORM_QPOLYG,CellModel(NORM_QPOLYG)));
 105     _map_of_unique_instance.insert(std::make_pair(NORM_POLYL,CellModel(NORM_POLYL)));
 106     _map_of_unique_instance.insert(std::make_pair(NORM_ERROR,CellModel(NORM_ERROR)));
 107   }
 108
 109   CellModel::CellModel(NormalizedCellType type):_type(type)
 110   {
 111     _is_extruded=false;
 112     _quadratic=false;
 113     _dyn=false;
 114     _extruded_type=NORM_ERROR;
 115     _linear_type=NORM_ERROR;
 116     _quadratic_type=NORM_ERROR;
 117     switch(type)
 118       {
 119       case NORM_POINT1:
 120         {
 121           _nb_of_pts=1; _nb_of_sons=0; _dim=0; _extruded_type=NORM_SEG2; _is_simplex=true;
 122         }
 123         break;
 124       case NORM_SEG2:
 125         {
 126           _nb_of_pts=2; _nb_of_sons=2; _dim=1; _extruded_type=NORM_QUAD4; _quadratic_type=NORM_SEG3; _is_simplex=true; _is_extruded=true;
 127           _sons_type[0]=NORM_POINT1; _sons_type[1]=NORM_POINT1;
 128           _sons_con[0][0]=0; _nb_of_sons_con[0]=1;
 129           _sons_con[1][0]=1; _nb_of_sons_con[1]=1;
 130         }
 131         break;
 132       case NORM_SEG3:
 133         {
 134           _nb_of_pts=3; _nb_of_sons=3; _dim=1; _extruded_type=NORM_QUAD8; _linear_type=NORM_SEG2; _quadratic=true; _is_simplex=false;
 135           _sons_type[0]=NORM_POINT1; _sons_type[1]=NORM_POINT1; _sons_type[2]=NORM_POINT1;
 136           _sons_con[0][0]=0; _nb_of_sons_con[0]=1;
 137           _sons_con[1][0]=1; _nb_of_sons_con[1]=1;
 138           _sons_con[2][0]=2; _nb_of_sons_con[2]=1;
 139         }
 140         break;
 141       case NORM_SEG4:
 142         {
 143           _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
 144           _sons_type[0]=NORM_POINT1; _sons_type[1]=NORM_POINT1; _sons_type[2]=NORM_POINT1; _sons_type[3]=NORM_POINT1;
 145           _sons_con[0][0]=0; _nb_of_sons_con[0]=1;
 146           _sons_con[1][0]=1; _nb_of_sons_con[1]=1;
 147           _sons_con[2][0]=2; _nb_of_sons_con[2]=1;
 148           _sons_con[3][0]=3; _nb_of_sons_con[3]=1;
 149         }
 150         break;
 151       case NORM_TETRA4:
 152         {
 153           _nb_of_pts=4; _nb_of_sons=4; _dim=3; _quadratic_type=NORM_TETRA10; _is_simplex=true;
 154           _sons_type[0]=NORM_TRI3; _sons_type[1]=NORM_TRI3; _sons_type[2]=NORM_TRI3; _sons_type[3]=NORM_TRI3;
 155           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=2; _nb_of_sons_con[0]=3;
 156           _sons_con[1][0]=0; _sons_con[1][1]=3; _sons_con[1][2]=1; _nb_of_sons_con[1]=3;
 157           _sons_con[2][0]=1; _sons_con[2][1]=3; _sons_con[2][2]=2; _nb_of_sons_con[2]=3;
 158           _sons_con[3][0]=2; _sons_con[3][1]=3; _sons_con[3][2]=0; _nb_of_sons_con[3]=3;
 159         }
 160         break;
 161       case NORM_HEXA8:
 162         {
 163           _nb_of_pts=8; _nb_of_sons=6; _dim=3; _quadratic_type=NORM_HEXA20; _is_simplex=false; _is_extruded=true;
 164           _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;
 165           _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;
 166           _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;
 167           _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;
 168           _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;
 169           _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;
 170           _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;
 171         }
 172         break;
 173       case NORM_QUAD4:
 174         {
 175           _nb_of_pts=4; _nb_of_sons=4; _dim=2; _quadratic_type=NORM_QUAD8; _is_simplex=false; _is_extruded=true;
 176           _sons_type[0]=NORM_SEG2; _sons_type[1]=NORM_SEG2; _sons_type[2]=NORM_SEG2; _sons_type[3]=NORM_SEG2;
 177           _sons_con[0][0]=0; _sons_con[0][1]=1; _nb_of_sons_con[0]=2;
 178           _sons_con[1][0]=1; _sons_con[1][1]=2; _nb_of_sons_con[1]=2;
 179           _sons_con[2][0]=2; _sons_con[2][1]=3; _nb_of_sons_con[2]=2;
 180           _sons_con[3][0]=3; _sons_con[3][1]=0; _nb_of_sons_con[3]=2; _extruded_type=NORM_HEXA8;
 181         }
 182         break;
 183       case NORM_TRI3:
 184         {
 185           _nb_of_pts=3; _nb_of_sons=3; _dim=2; _quadratic_type=NORM_TRI6; _is_simplex=true;
 186           _sons_type[0]=NORM_SEG2; _sons_type[1]=NORM_SEG2; _sons_type[2]=NORM_SEG2;
 187           _sons_con[0][0]=0; _sons_con[0][1]=1; _nb_of_sons_con[0]=2;
 188           _sons_con[1][0]=1; _sons_con[1][1]=2; _nb_of_sons_con[1]=2;
 189           _sons_con[2][0]=2; _sons_con[2][1]=0; _nb_of_sons_con[2]=2; _extruded_type=NORM_PENTA6;
 190         }
 191         break;
 192       case NORM_TRI6:
 193         {
 194           _nb_of_pts=6; _nb_of_sons=3; _dim=2; _linear_type=NORM_TRI3; _is_simplex=false;
 195           _sons_type[0]=NORM_SEG3; _sons_type[1]=NORM_SEG3; _sons_type[2]=NORM_SEG3;
 196           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=3; _nb_of_sons_con[0]=3;
 197           _sons_con[1][0]=1; _sons_con[1][1]=2; _sons_con[1][2]=4; _nb_of_sons_con[1]=3;
 198           _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;
 199         }
 200         break;
 201       case NORM_TRI7:
 202         {
 203           _nb_of_pts=7; _nb_of_sons=3; _dim=2; _linear_type=NORM_TRI3; _is_simplex=false;
 204           _sons_type[0]=NORM_SEG3; _sons_type[1]=NORM_SEG3; _sons_type[2]=NORM_SEG3;
 205           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=3; _nb_of_sons_con[0]=3;
 206           _sons_con[1][0]=1; _sons_con[1][1]=2; _sons_con[1][2]=4; _nb_of_sons_con[1]=3;
 207           _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
 208         }
 209         break;
 210       case NORM_QUAD8:
 211         {
 212           _nb_of_pts=8; _nb_of_sons=4; _dim=2; _linear_type=NORM_QUAD4; _is_simplex=false;
 213           _sons_type[0]=NORM_SEG3; _sons_type[1]=NORM_SEG3; _sons_type[2]=NORM_SEG3; _sons_type[3]=NORM_SEG3;
 214           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=4; _nb_of_sons_con[0]=3;
 215           _sons_con[1][0]=1; _sons_con[1][1]=2; _sons_con[1][2]=5; _nb_of_sons_con[1]=3;
 216           _sons_con[2][0]=2; _sons_con[2][1]=3; _sons_con[2][2]=6; _nb_of_sons_con[2]=3;
 217           _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;
 218         }
 219         break;
 220       case NORM_QUAD9:
 221         {
 222           _nb_of_pts=9; _nb_of_sons=4; _dim=2; _linear_type=NORM_QUAD4; _is_simplex=false;
 223           _sons_type[0]=NORM_SEG3; _sons_type[1]=NORM_SEG3; _sons_type[2]=NORM_SEG3; _sons_type[3]=NORM_SEG3;
 224           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=4; _nb_of_sons_con[0]=3;
 225           _sons_con[1][0]=1; _sons_con[1][1]=2; _sons_con[1][2]=5; _nb_of_sons_con[1]=3;
 226           _sons_con[2][0]=2; _sons_con[2][1]=3; _sons_con[2][2]=6; _nb_of_sons_con[2]=3;
 227           _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;
 228         }
 229         break;
 230       case NORM_PYRA5:
 231         {
 232           _nb_of_pts=5; _nb_of_sons=5; _dim=3; _quadratic_type=NORM_PYRA13; _is_simplex=false;
 233           _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;
 234           _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;
 235           _sons_con[1][0]=0; _sons_con[1][1]=4; _sons_con[1][2]=1; _nb_of_sons_con[1]=3;
 236           _sons_con[2][0]=1; _sons_con[2][1]=4; _sons_con[2][2]=2; _nb_of_sons_con[2]=3;
 237           _sons_con[3][0]=2; _sons_con[3][1]=4; _sons_con[3][2]=3; _nb_of_sons_con[3]=3;
 238           _sons_con[4][0]=3; _sons_con[4][1]=4; _sons_con[4][2]=0; _nb_of_sons_con[4]=3;
 239         }
 240         break;
 241       case NORM_PENTA6:
 242         {
 243           _nb_of_pts=6; _nb_of_sons=5; _dim=3; _quadratic_type=NORM_PENTA15; _is_simplex=false; _is_extruded=true;
 244           _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;
 245           _sons_con[0][0]=0; _sons_con[0][1]=1; _sons_con[0][2]=2; _nb_of_sons_con[0]=3;
 246           _sons_con[1][0]=3; _sons_con[1][1]=5; _sons_con[1][2]=4; _nb_of_sons_con[1]=3;
 247           _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;
 248           _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;
 249           _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;
 250         }
 251         break;
 252       case NORM_TETRA10:
 253         {
 254           _nb_of_pts=10; _nb_of_sons=4; _dim=3; _linear_type=NORM_TETRA4; _is_simplex=false;
 255           _sons_type[0]=NORM_TRI6; _sons_type[1]=NORM_TRI6; _sons_type[2]=NORM_TRI6; _sons_type[3]=NORM_TRI6;
 256           _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;
 257           _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;
 258           _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;
 259           _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;
 260         }
 261         break;
 262       case NORM_HEXGP12:
 263         {
 264           _nb_of_pts=12; _nb_of_sons=8; _dim=3; _is_simplex=false; _is_extruded=true;
 265           _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;
 266           _sons_type[6]=NORM_QUAD4; _sons_type[7]=NORM_QUAD4;
 267           _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;
 268           _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;
 269           _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;
 270           _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;
 271           _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;
 272           _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;
 273           _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;
 274           _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;
 275         }
 276         break;
 277       case NORM_PYRA13:
 278         {
 279           _nb_of_pts=13; _nb_of_sons=5; _dim=3; _linear_type=NORM_PYRA5; _is_simplex=false;
 280           _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;
 281           _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;
 282           _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;
 283           _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;
 284           _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;
 285           _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;
 286         }
 287         break;
 288       case NORM_PENTA15:
 289         {
 290           _nb_of_pts=15; _nb_of_sons=5; _dim=3; _linear_type=NORM_PENTA6; _is_simplex=false;
 291           _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;
 292           _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;
 293           _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;
 294           _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;
 295           _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;
 296           _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;
 297         }
 298         break;
 299       case NORM_HEXA20:
 300         {
 301           _nb_of_pts=20; _nb_of_sons=6; _dim=3; _linear_type=NORM_HEXA8; _is_simplex=false;
 302           _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;
 303           _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;
 304           _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;
 305           _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;
 306           _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;
 307           _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;
 308           _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;
 309         }
 310         break;
 311       case NORM_HEXA27:
 312         {
 313           _nb_of_pts=27; _nb_of_sons=6; _dim=3; _linear_type=NORM_HEXA8; _is_simplex=false;
 314           _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;
 315           _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;
 316           _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;
 317           _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;
 318           _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;
 319           _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;
 320           _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;
 321           _quadratic=true;
 322         }
 323         break;
 324       case NORM_POLYGON:
 325         {
 326           _nb_of_pts=0; _nb_of_sons=0; _dim=2; _dyn=true; _extruded_type=NORM_POLYHED; _is_simplex=false;
 327         }
 328         break;
 329       case NORM_POLYHED:
 330         {
 331           _nb_of_pts=0; _nb_of_sons=0; _dim=3; _dyn=true; _is_simplex=false;
 332         }
 333         break;
 334       case NORM_QPOLYG:
 335         {
 336           _nb_of_pts=0; _nb_of_sons=0; _dim=2; _dyn=true; _is_simplex=false; _quadratic=true;
 337         }
 338         break;
 339       case NORM_POLYL:
 340         {
 341           _nb_of_pts=0; _nb_of_sons=0; _dim=1; _dyn=true; _extruded_type=NORM_POLYGON; _is_simplex=false;
 342         }
 343       case NORM_ERROR:
 344         {
 345           _nb_of_pts=std::numeric_limits<unsigned>::max(); _nb_of_sons=std::numeric_limits<unsigned>::max(); _dim=std::numeric_limits<unsigned>::max();
 346         }
 347         break;
 348       }
 349   }
 350
 351   /*!
 352    * Equivalent to getNumberOfSons except that this method deals with dynamic type.
 353    */
 354   unsigned CellModel::getNumberOfSons2(const int *conn, int lgth) const
 355   {
 356     if(!isDynamic())
 357       return getNumberOfSons();
 358     if(_dim==2)
 359       {
 360         if(_type==NORM_POLYGON)
 361           return lgth;
 362         else
 363           return lgth/2;
 364       }
 365     else if(_dim==1)
 366       return lgth;//NORM_POLYL
 367     else
 368       return std::count(conn,conn+lgth,-1)+1;
 369   }
 370
 371   /*!
 372    * Equivalent to getSonType except that this method deals with dynamic type.
 373    */
 374   NormalizedCellType CellModel::getSonType2(unsigned sonId) const
 375   {
 376     if(!isDynamic())
 377       return getSonType(sonId);
 378     if(_dim==2)
 379       {
 380         if(_type==NORM_POLYGON)
 381           return NORM_SEG2;
 382         else
 383           return NORM_SEG3;
 384       }
 385     else if(_dim==1)
 386       return NORM_ERROR;//NORM_POLYL
 387     //polyedron
 388     return NORM_POLYGON;
 389   }
 390
 391   /*!
 392    * \b WARNING this method do not manage correctly types that return true at the call of isDynamic. Use fillSonCellNodalConnectivity2 instead.
 393    */
 394   unsigned CellModel::fillSonCellNodalConnectivity(int sonId, const int *nodalConn, int *sonNodalConn) const
 395   {
 396     unsigned nbOfTurnLoop=_nb_of_sons_con[sonId];
 397     const unsigned *sonConn=_sons_con[sonId];
 398     for(unsigned i=0;i<nbOfTurnLoop;i++)
 399       sonNodalConn[i]=nodalConn[sonConn[i]];
 400     return nbOfTurnLoop;
 401   }
 402
 403   unsigned CellModel::fillSonCellNodalConnectivity2(int sonId, const int *nodalConn, int lgth, int *sonNodalConn, NormalizedCellType& typeOfSon) const
 404   {
 405     typeOfSon=getSonType2(sonId);
 406     if(!isDynamic())
 407       return fillSonCellNodalConnectivity(sonId,nodalConn,sonNodalConn);
 408     else
 409       {
 410         if(_dim==2)//polygon
 411           {
 412             if(_type==NORM_POLYGON)
 413               {
 414                 sonNodalConn[0]=nodalConn[sonId];
 415                 sonNodalConn[1]=nodalConn[(sonId+1)%lgth];
 416                 return 2;
 417               }
 418             else
 419               {
 420                 sonNodalConn[0]=nodalConn[sonId];
 421                 sonNodalConn[1]=nodalConn[(sonId+1)%lgth];
 422                 sonNodalConn[2]=nodalConn[sonId+lgth];
 423                 return 3;
 424               }
 425           }
 426         else if(_dim==3)
 427           {//polyedron
 428             const int *where=nodalConn;
 429             for(int i=0;i<sonId;i++)
 430               {
 431                 where=std::find(where,nodalConn+lgth,-1);
 432                 where++;
 433               }
 434             const int *where2=std::find(where,nodalConn+lgth,-1);
 435             std::copy(where,where2,sonNodalConn);
 436             return where2-where;
 437           }
 438         else
 439           throw INTERP_KERNEL::Exception("CellModel::fillSonCellNodalConnectivity2 : no sons on NORM_POLYL !");
 440       }
 441   }
 442
 443   //================================================================================
 444   /*!
 445    * \brief Return number of nodes in sonId-th son of a Dynamic() cell
 446    */
 447   //================================================================================
 448
 449   unsigned CellModel::getNumberOfNodesConstituentTheSon2(unsigned sonId, const int *nodalConn, int lgth) const
 450   {
 451     if(!isDynamic())
 452       return getNumberOfNodesConstituentTheSon(sonId);
 453
 454     if(_dim==2)//polygon
 455       {
 456         if(_type==NORM_POLYGON)
 457           return 2;
 458         else
 459           return 3;
 460       }
 461     else if(_dim==3)
 462       {//polyedron
 463         const int *where=nodalConn;
 464         for(unsigned int i=0;i<sonId;i++)
 465           {
 466             where=std::find(where,nodalConn+lgth,-1);
 467             where++;
 468           }
 469         const int *where2=std::find(where,nodalConn+lgth,-1);
 470         return where2-where;
 471       }
 472     else
 473       throw INTERP_KERNEL::Exception("CellModel::getNumberOfNodesConstituentTheSon2 : no sons on NORM_POLYL !");
 474   }
 475
 476   /*!
 477    * This method retrieves if cell1 represented by 'conn1' and cell2 represented by 'conn2'
 478    * 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.
 479    * If not an exception will be thrown.
 480    * @return True if two cells have same orientation, false if not.
 481    */
 482   bool CellModel::getOrientationStatus(unsigned lgth, const int *conn1, const int *conn2) const
 483   {
 484     if(_dim!=1 && _dim!=2)
 485       throw INTERP_KERNEL::Exception("CellModel::getOrientationStatus : invalid dimension ! Must be 1 or 2 !");
 486     if(!_quadratic)
 487       {
 488         std::vector<int> tmp(2*lgth);
 489         std::vector<int>::iterator it=std::copy(conn1,conn1+lgth,tmp.begin());
 490         std::copy(conn1,conn1+lgth,it);
 491         it=std::search(tmp.begin(),tmp.end(),conn2,conn2+lgth);
 492         if(it==tmp.begin())
 493           return true;
 494         if(it!=tmp.end())
 495           return _dim!=1;
 496         std::vector<int>::reverse_iterator it2=std::search(tmp.rbegin(),tmp.rend(),conn2,conn2+lgth);
 497         if(it2!=tmp.rend())
 498           return false;
 499         throw INTERP_KERNEL::Exception("CellModel::getOrientationStatus : Request of orientation status of non equal connectively cells !");
 500       }
 501     else
 502       {
 503         if(_dim!=1)
 504           {
 505             std::vector<int> tmp(lgth);
 506             std::vector<int>::iterator it=std::copy(conn1,conn1+lgth/2,tmp.begin());
 507             std::copy(conn1,conn1+lgth/2,it);
 508             it=std::search(tmp.begin(),tmp.end(),conn2,conn2+lgth/2);
 509             int d=std::distance(tmp.begin(),it);
 510             if(it==tmp.end())
 511               return false;
 512             it=std::copy(conn1+lgth/2,conn1+lgth,tmp.begin());
 513             std::copy(conn1+lgth/2,conn1+lgth,it);
 514             it=std::search(tmp.begin(),tmp.end(),conn2,conn2+lgth);
 515             if(it==tmp.end())
 516               return false;
 517             int d2=std::distance(tmp.begin(),it);
 518             return d==d2;
 519           }
 520         else
 521           {
 522             int p=(lgth+1)/2;
 523             std::vector<int> tmp(2*p);
 524             std::vector<int>::iterator it=std::copy(conn1,conn1+p,tmp.begin());
 525             std::copy(conn1,conn1+p,it);
 526             it=std::search(tmp.begin(),tmp.end(),conn2,conn2+p);
 527             int d=std::distance(tmp.begin(),it);
 528             if(it==tmp.end())
 529               return false;
 530             tmp.resize(2*p-2);
 531             it=std::copy(conn1+p,conn1+lgth,tmp.begin());
 532             std::copy(conn1+p,conn1+lgth,it);
 533             it=std::search(tmp.begin(),tmp.end(),conn2+p,conn2+lgth);
 534             if(it==tmp.end())
 535               return false;
 536             int d2=std::distance(tmp.begin(),it);
 537             return d==d2;
 538           }
 539       }
 540   }
 541
 542 }