src/SMESHDS/SMESHDS_Mesh.hxx

   1 // Copyright (C) 2007-2014  CEA/DEN, EDF R&D, OPEN CASCADE
   2 //
   3 // Copyright (C) 2003-2007  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, or (at your option) any later version.
  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.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 //
  22
  23 //  SMESH SMESHDS : management of mesh data and SMESH document
  24 //  File   : SMESHDS_Mesh.hxx
  25 //  Module : SMESH
  26 //
  27 #ifndef _SMESHDS_Mesh_HeaderFile
  28 #define _SMESHDS_Mesh_HeaderFile
  29
  30 #include "SMESH_SMESHDS.hxx"
  31
  32 #include "SMDS_Mesh.hxx"
  33 #include "SMESHDS_SubMesh.hxx"
  34
  35 #include <TopTools_IndexedMapOfShape.hxx>
  36 #include <TopoDS_Shape.hxx>
  37
  38 class TopoDS_Solid ;
  39 class TopoDS_Shell ;
  40 class TopoDS_Face  ;
  41 class TopoDS_Vertex;
  42 class TopoDS_Edge  ;
  43
  44 class SMESHDS_Script;
  45 class SMESHDS_Hypothesis;
  46 class SMDS_MeshNode     ;
  47 class SMDS_MeshEdge     ;
  48 class SMDS_MeshFace     ;
  49 class SMDS_MeshVolume   ;
  50 class SMDS_Mesh0DElement;
  51 class SMDS_BallElement;
  52
  53 #include <NCollection_DataMap.hxx>
  54 #include <map>
  55 /*
  56  * Using of native hash_map isn't portable and don't work on WIN32 platform.
  57  * So this functionality implement on new NCollection_DataMap technology
  58  */
  59 #include "SMESHDS_DataMapOfShape.hxx"
  60
  61 class SMESHDS_GroupBase;
  62 class DownIdType;
  63
  64 class SMESHDS_EXPORT SMESHDS_Mesh : public SMDS_Mesh
  65 {
  66 public:
  67   SMESHDS_Mesh(int theMeshID, bool theIsEmbeddedMode);
  68   bool IsEmbeddedMode();
  69   void SetPersistentId(int id);
  70   int GetPersistentId() const;
  71
  72   void ShapeToMesh(const TopoDS_Shape & S);
  73   TopoDS_Shape ShapeToMesh() const;
  74   bool AddHypothesis(const TopoDS_Shape & SS, const SMESHDS_Hypothesis * H);
  75   bool RemoveHypothesis(const TopoDS_Shape & S, const SMESHDS_Hypothesis * H);
  76
  77   virtual SMDS_MeshNode* AddNodeWithID(double x, double y, double z, int ID);
  78   virtual SMDS_MeshNode* AddNode(double x, double y, double z);
  79
  80   virtual SMDS_Mesh0DElement* Add0DElementWithID(int nodeID, int ID);
  81   virtual SMDS_Mesh0DElement* Add0DElementWithID(const SMDS_MeshNode * node, int ID);
  82   virtual SMDS_Mesh0DElement* Add0DElement      (const SMDS_MeshNode * node);
  83
  84   virtual SMDS_BallElement* AddBallWithID(int n,                   double diameter, int ID);
  85   virtual SMDS_BallElement* AddBallWithID(const SMDS_MeshNode * n, double diameter, int ID);
  86   virtual SMDS_BallElement* AddBall      (const SMDS_MeshNode * n, double diameter);
  87
  88   virtual SMDS_MeshEdge* AddEdgeWithID(int n1, int n2, int ID);
  89   virtual SMDS_MeshEdge* AddEdgeWithID(const SMDS_MeshNode * n1,
  90                                        const SMDS_MeshNode * n2,
  91                                        int ID);
  92   virtual SMDS_MeshEdge* AddEdge(const SMDS_MeshNode * n1,
  93                                  const SMDS_MeshNode * n2);
  94
  95   // 2d order edge with 3 nodes: n12 - node between n1 and n2
  96   virtual SMDS_MeshEdge* AddEdgeWithID(int n1, int n2, int n12, int ID);
  97   virtual SMDS_MeshEdge* AddEdgeWithID(const SMDS_MeshNode * n1,
  98                                        const SMDS_MeshNode * n2,
  99                                        const SMDS_MeshNode * n12,
 100                                        int ID);
 101   virtual SMDS_MeshEdge* AddEdge(const SMDS_MeshNode * n1,
 102                                  const SMDS_MeshNode * n2,
 103                                  const SMDS_MeshNode * n12);
 104   // tria 3
 105   virtual SMDS_MeshFace* AddFaceWithID(int n1, int n2, int n3, int ID);
 106   virtual SMDS_MeshFace* AddFaceWithID(const SMDS_MeshNode * n1,
 107                                        const SMDS_MeshNode * n2,
 108                                        const SMDS_MeshNode * n3,
 109                                        int ID);
 110   virtual SMDS_MeshFace* AddFace(const SMDS_MeshNode * n1,
 111                                  const SMDS_MeshNode * n2,
 112                                  const SMDS_MeshNode * n3);
 113   // quad 4
 114   virtual SMDS_MeshFace* AddFaceWithID(int n1, int n2, int n3, int n4, int ID);
 115   virtual SMDS_MeshFace* AddFaceWithID(const SMDS_MeshNode * n1,
 116                                        const SMDS_MeshNode * n2,
 117                                        const SMDS_MeshNode * n3,
 118                                        const SMDS_MeshNode * n4,
 119                                        int ID);
 120   virtual SMDS_MeshFace* AddFace(const SMDS_MeshNode * n1,
 121                                  const SMDS_MeshNode * n2,
 122                                  const SMDS_MeshNode * n3,
 123                                  const SMDS_MeshNode * n4);
 124
 125   // 2d order triangle of 6 nodes
 126   virtual SMDS_MeshFace* AddFaceWithID(int n1, int n2, int n3,
 127                                        int n12,int n23,int n31, int ID);
 128   virtual SMDS_MeshFace* AddFaceWithID(const SMDS_MeshNode * n1,
 129                                        const SMDS_MeshNode * n2,
 130                                        const SMDS_MeshNode * n3,
 131                                        const SMDS_MeshNode * n12,
 132                                        const SMDS_MeshNode * n23,
 133                                        const SMDS_MeshNode * n31,
 134                                        int ID);
 135   virtual SMDS_MeshFace* AddFace(const SMDS_MeshNode * n1,
 136                                  const SMDS_MeshNode * n2,
 137                                  const SMDS_MeshNode * n3,
 138                                  const SMDS_MeshNode * n12,
 139                                  const SMDS_MeshNode * n23,
 140                                  const SMDS_MeshNode * n31);
 141
 142   // biquadratic triangle of 7 nodes
 143   virtual SMDS_MeshFace* AddFaceWithID(int n1, int n2, int n3,
 144                                        int n12,int n23,int n31, int nCenter, int ID);
 145   virtual SMDS_MeshFace* AddFaceWithID(const SMDS_MeshNode * n1,
 146                                        const SMDS_MeshNode * n2,
 147                                        const SMDS_MeshNode * n3,
 148                                        const SMDS_MeshNode * n12,
 149                                        const SMDS_MeshNode * n23,
 150                                        const SMDS_MeshNode * n31,
 151                                        const SMDS_MeshNode * nCenter,
 152                                        int ID);
 153   virtual SMDS_MeshFace* AddFace(const SMDS_MeshNode * n1,
 154                                  const SMDS_MeshNode * n2,
 155                                  const SMDS_MeshNode * n3,
 156                                  const SMDS_MeshNode * n12,
 157                                  const SMDS_MeshNode * n23,
 158                                  const SMDS_MeshNode * n31,
 159                                  const SMDS_MeshNode * nCenter);
 160
 161   // 2d order quadrangle
 162   virtual SMDS_MeshFace* AddFaceWithID(int n1, int n2, int n3, int n4,
 163                                        int n12,int n23,int n34,int n41, int ID);
 164   virtual SMDS_MeshFace* AddFaceWithID(const SMDS_MeshNode * n1,
 165                                        const SMDS_MeshNode * n2,
 166                                        const SMDS_MeshNode * n3,
 167                                        const SMDS_MeshNode * n4,
 168                                        const SMDS_MeshNode * n12,
 169                                        const SMDS_MeshNode * n23,
 170                                        const SMDS_MeshNode * n34,
 171                                        const SMDS_MeshNode * n41,
 172                                        int ID);
 173   virtual SMDS_MeshFace* AddFace(const SMDS_MeshNode * n1,
 174                                  const SMDS_MeshNode * n2,
 175                                  const SMDS_MeshNode * n3,
 176                                  const SMDS_MeshNode * n4,
 177                                  const SMDS_MeshNode * n12,
 178                                  const SMDS_MeshNode * n23,
 179                                  const SMDS_MeshNode * n34,
 180                                  const SMDS_MeshNode * n41);
 181
 182   // biquadratic quadrangle of 9 nodes
 183   virtual SMDS_MeshFace* AddFaceWithID(int n1, int n2, int n3, int n4,
 184                                        int n12,int n23,int n34,int n41, int nCenter, int ID);
 185   virtual SMDS_MeshFace* AddFaceWithID(const SMDS_MeshNode * n1,
 186                                        const SMDS_MeshNode * n2,
 187                                        const SMDS_MeshNode * n3,
 188                                        const SMDS_MeshNode * n4,
 189                                        const SMDS_MeshNode * n12,
 190                                        const SMDS_MeshNode * n23,
 191                                        const SMDS_MeshNode * n34,
 192                                        const SMDS_MeshNode * n41,
 193                                        const SMDS_MeshNode * nCenter,
 194                                        int ID);
 195   virtual SMDS_MeshFace* AddFace(const SMDS_MeshNode * n1,
 196                                  const SMDS_MeshNode * n2,
 197                                  const SMDS_MeshNode * n3,
 198                                  const SMDS_MeshNode * n4,
 199                                  const SMDS_MeshNode * n12,
 200                                  const SMDS_MeshNode * n23,
 201                                  const SMDS_MeshNode * n34,
 202                                  const SMDS_MeshNode * n41,
 203                                  const SMDS_MeshNode * nCenter);
 204   // tetra 4
 205   virtual SMDS_MeshVolume* AddVolumeWithID(int n1, int n2, int n3, int n4, int ID);
 206   virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
 207                                            const SMDS_MeshNode * n2,
 208                                            const SMDS_MeshNode * n3,
 209                                            const SMDS_MeshNode * n4,
 210                                            int ID);
 211   virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
 212                                      const SMDS_MeshNode * n2,
 213                                      const SMDS_MeshNode * n3,
 214                                      const SMDS_MeshNode * n4);
 215   // pyra 5
 216   virtual SMDS_MeshVolume* AddVolumeWithID(int n1, int n2, int n3, int n4, int n5, int ID);
 217   virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
 218                                            const SMDS_MeshNode * n2,
 219                                            const SMDS_MeshNode * n3,
 220                                            const SMDS_MeshNode * n4,
 221                                            const SMDS_MeshNode * n5,
 222                                            int ID);
 223   virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
 224                                      const SMDS_MeshNode * n2,
 225                                      const SMDS_MeshNode * n3,
 226                                      const SMDS_MeshNode * n4,
 227                                      const SMDS_MeshNode * n5);
 228   // penta 6
 229   virtual SMDS_MeshVolume* AddVolumeWithID(int n1, int n2, int n3, int n4, int n5, int n6, int ID);
 230   virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
 231                                            const SMDS_MeshNode * n2,
 232                                            const SMDS_MeshNode * n3,
 233                                            const SMDS_MeshNode * n4,
 234                                            const SMDS_MeshNode * n5,
 235                                            const SMDS_MeshNode * n6,
 236                                            int ID);
 237   virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
 238                                      const SMDS_MeshNode * n2,
 239                                      const SMDS_MeshNode * n3,
 240                                      const SMDS_MeshNode * n4,
 241                                      const SMDS_MeshNode * n5,
 242                                      const SMDS_MeshNode * n6);
 243   // hexa 8
 244   virtual SMDS_MeshVolume* AddVolumeWithID(int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8, int ID);
 245   virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
 246                                            const SMDS_MeshNode * n2,
 247                                            const SMDS_MeshNode * n3,
 248                                            const SMDS_MeshNode * n4,
 249                                            const SMDS_MeshNode * n5,
 250                                            const SMDS_MeshNode * n6,
 251                                            const SMDS_MeshNode * n7,
 252                                            const SMDS_MeshNode * n8,
 253                                            int ID);
 254   virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
 255                                      const SMDS_MeshNode * n2,
 256                                      const SMDS_MeshNode * n3,
 257                                      const SMDS_MeshNode * n4,
 258                                      const SMDS_MeshNode * n5,
 259                                      const SMDS_MeshNode * n6,
 260                                      const SMDS_MeshNode * n7,
 261                                      const SMDS_MeshNode * n8);
 262   // hexagonal prism of 12 nodes
 263   virtual SMDS_MeshVolume* AddVolumeWithID(int n1, int n2, int n3, int n4, int n5, int n6,
 264                                            int n7, int n8, int n9, int n10, int n11, int n12, int ID);
 265   virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
 266                                            const SMDS_MeshNode * n2,
 267                                            const SMDS_MeshNode * n3,
 268                                            const SMDS_MeshNode * n4,
 269                                            const SMDS_MeshNode * n5,
 270                                            const SMDS_MeshNode * n6,
 271                                            const SMDS_MeshNode * n7,
 272                                            const SMDS_MeshNode * n8,
 273                                            const SMDS_MeshNode * n9,
 274                                            const SMDS_MeshNode * n10,
 275                                            const SMDS_MeshNode * n11,
 276                                            const SMDS_MeshNode * n12,
 277                                            int ID);
 278   virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
 279                                      const SMDS_MeshNode * n2,
 280                                      const SMDS_MeshNode * n3,
 281                                      const SMDS_MeshNode * n4,
 282                                      const SMDS_MeshNode * n5,
 283                                      const SMDS_MeshNode * n6,
 284                                      const SMDS_MeshNode * n7,
 285                                      const SMDS_MeshNode * n8,
 286                                      const SMDS_MeshNode * n9,
 287                                      const SMDS_MeshNode * n10,
 288                                      const SMDS_MeshNode * n11,
 289                                      const SMDS_MeshNode * n12);
 290
 291   // 2d order tetrahedron of 10 nodes
 292   virtual SMDS_MeshVolume* AddVolumeWithID(int n1, int n2, int n3, int n4,
 293                                            int n12,int n23,int n31,
 294                                            int n14,int n24,int n34, int ID);
 295   virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
 296                                            const SMDS_MeshNode * n2,
 297                                            const SMDS_MeshNode * n3,
 298                                            const SMDS_MeshNode * n4,
 299                                            const SMDS_MeshNode * n12,
 300                                            const SMDS_MeshNode * n23,
 301                                            const SMDS_MeshNode * n31,
 302                                            const SMDS_MeshNode * n14,
 303                                            const SMDS_MeshNode * n24,
 304                                            const SMDS_MeshNode * n34,
 305                                            int ID);
 306   virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
 307                                      const SMDS_MeshNode * n2,
 308                                      const SMDS_MeshNode * n3,
 309                                      const SMDS_MeshNode * n4,
 310                                      const SMDS_MeshNode * n12,
 311                                      const SMDS_MeshNode * n23,
 312                                      const SMDS_MeshNode * n31,
 313                                      const SMDS_MeshNode * n14,
 314                                      const SMDS_MeshNode * n24,
 315                                      const SMDS_MeshNode * n34);
 316
 317   // 2d order pyramid of 13 nodes
 318   virtual SMDS_MeshVolume* AddVolumeWithID(int n1, int n2, int n3, int n4, int n5,
 319                                            int n12,int n23,int n34,int n41,
 320                                            int n15,int n25,int n35,int n45,
 321                                            int ID);
 322   virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
 323                                            const SMDS_MeshNode * n2,
 324                                            const SMDS_MeshNode * n3,
 325                                            const SMDS_MeshNode * n4,
 326                                            const SMDS_MeshNode * n5,
 327                                            const SMDS_MeshNode * n12,
 328                                            const SMDS_MeshNode * n23,
 329                                            const SMDS_MeshNode * n34,
 330                                            const SMDS_MeshNode * n41,
 331                                            const SMDS_MeshNode * n15,
 332                                            const SMDS_MeshNode * n25,
 333                                            const SMDS_MeshNode * n35,
 334                                            const SMDS_MeshNode * n45,
 335                                            int ID);
 336   virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
 337                                      const SMDS_MeshNode * n2,
 338                                      const SMDS_MeshNode * n3,
 339                                      const SMDS_MeshNode * n4,
 340                                      const SMDS_MeshNode * n5,
 341                                      const SMDS_MeshNode * n12,
 342                                      const SMDS_MeshNode * n23,
 343                                      const SMDS_MeshNode * n34,
 344                                      const SMDS_MeshNode * n41,
 345                                      const SMDS_MeshNode * n15,
 346                                      const SMDS_MeshNode * n25,
 347                                      const SMDS_MeshNode * n35,
 348                                      const SMDS_MeshNode * n45);
 349
 350   // 2d order Pentahedron with 15 nodes
 351   virtual SMDS_MeshVolume* AddVolumeWithID(int n1, int n2, int n3,
 352                                            int n4, int n5, int n6,
 353                                            int n12,int n23,int n31,
 354                                            int n45,int n56,int n64,
 355                                            int n14,int n25,int n36,
 356                                            int ID);
 357   virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
 358                                            const SMDS_MeshNode * n2,
 359                                            const SMDS_MeshNode * n3,
 360                                            const SMDS_MeshNode * n4,
 361                                            const SMDS_MeshNode * n5,
 362                                            const SMDS_MeshNode * n6,
 363                                            const SMDS_MeshNode * n12,
 364                                            const SMDS_MeshNode * n23,
 365                                            const SMDS_MeshNode * n31,
 366                                            const SMDS_MeshNode * n45,
 367                                            const SMDS_MeshNode * n56,
 368                                            const SMDS_MeshNode * n64,
 369                                            const SMDS_MeshNode * n14,
 370                                            const SMDS_MeshNode * n25,
 371                                            const SMDS_MeshNode * n36,
 372                                            int ID);
 373   virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
 374                                      const SMDS_MeshNode * n2,
 375                                      const SMDS_MeshNode * n3,
 376                                      const SMDS_MeshNode * n4,
 377                                      const SMDS_MeshNode * n5,
 378                                      const SMDS_MeshNode * n6,
 379                                      const SMDS_MeshNode * n12,
 380                                      const SMDS_MeshNode * n23,
 381                                      const SMDS_MeshNode * n31,
 382                                      const SMDS_MeshNode * n45,
 383                                      const SMDS_MeshNode * n56,
 384                                      const SMDS_MeshNode * n64,
 385                                      const SMDS_MeshNode * n14,
 386                                      const SMDS_MeshNode * n25,
 387                                      const SMDS_MeshNode * n36);
 388
 389   // 2d order Hexahedrons with 20 nodes
 390   virtual SMDS_MeshVolume* AddVolumeWithID(int n1, int n2, int n3, int n4,
 391                                            int n5, int n6, int n7, int n8,
 392                                            int n12,int n23,int n34,int n41,
 393                                            int n56,int n67,int n78,int n85,
 394                                            int n15,int n26,int n37,int n48,
 395                                            int ID);
 396   virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
 397                                            const SMDS_MeshNode * n2,
 398                                            const SMDS_MeshNode * n3,
 399                                            const SMDS_MeshNode * n4,
 400                                            const SMDS_MeshNode * n5,
 401                                            const SMDS_MeshNode * n6,
 402                                            const SMDS_MeshNode * n7,
 403                                            const SMDS_MeshNode * n8,
 404                                            const SMDS_MeshNode * n12,
 405                                            const SMDS_MeshNode * n23,
 406                                            const SMDS_MeshNode * n34,
 407                                            const SMDS_MeshNode * n41,
 408                                            const SMDS_MeshNode * n56,
 409                                            const SMDS_MeshNode * n67,
 410                                            const SMDS_MeshNode * n78,
 411                                            const SMDS_MeshNode * n85,
 412                                            const SMDS_MeshNode * n15,
 413                                            const SMDS_MeshNode * n26,
 414                                            const SMDS_MeshNode * n37,
 415                                            const SMDS_MeshNode * n48,
 416                                            int ID);
 417   virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
 418                                      const SMDS_MeshNode * n2,
 419                                      const SMDS_MeshNode * n3,
 420                                      const SMDS_MeshNode * n4,
 421                                      const SMDS_MeshNode * n5,
 422                                      const SMDS_MeshNode * n6,
 423                                      const SMDS_MeshNode * n7,
 424                                      const SMDS_MeshNode * n8,
 425                                      const SMDS_MeshNode * n12,
 426                                      const SMDS_MeshNode * n23,
 427                                      const SMDS_MeshNode * n34,
 428                                      const SMDS_MeshNode * n41,
 429                                      const SMDS_MeshNode * n56,
 430                                      const SMDS_MeshNode * n67,
 431                                      const SMDS_MeshNode * n78,
 432                                      const SMDS_MeshNode * n85,
 433                                      const SMDS_MeshNode * n15,
 434                                      const SMDS_MeshNode * n26,
 435                                      const SMDS_MeshNode * n37,
 436                                      const SMDS_MeshNode * n48);
 437
 438   // 2d order Hexahedrons with 27 nodes
 439   virtual SMDS_MeshVolume* AddVolumeWithID(int n1, int n2, int n3, int n4,
 440                                            int n5, int n6, int n7, int n8,
 441                                            int n12,int n23,int n34,int n41,
 442                                            int n56,int n67,int n78,int n85,
 443                                            int n15,int n26,int n37,int n48,
 444                                            int n1234,int n1256,int n2367,int n3478,
 445                                            int n1458,int n5678,int nCenter,
 446                                            int ID);
 447   virtual SMDS_MeshVolume* AddVolumeWithID(const SMDS_MeshNode * n1,
 448                                            const SMDS_MeshNode * n2,
 449                                            const SMDS_MeshNode * n3,
 450                                            const SMDS_MeshNode * n4,
 451                                            const SMDS_MeshNode * n5,
 452                                            const SMDS_MeshNode * n6,
 453                                            const SMDS_MeshNode * n7,
 454                                            const SMDS_MeshNode * n8,
 455                                            const SMDS_MeshNode * n12,
 456                                            const SMDS_MeshNode * n23,
 457                                            const SMDS_MeshNode * n34,
 458                                            const SMDS_MeshNode * n41,
 459                                            const SMDS_MeshNode * n56,
 460                                            const SMDS_MeshNode * n67,
 461                                            const SMDS_MeshNode * n78,
 462                                            const SMDS_MeshNode * n85,
 463                                            const SMDS_MeshNode * n15,
 464                                            const SMDS_MeshNode * n26,
 465                                            const SMDS_MeshNode * n37,
 466                                            const SMDS_MeshNode * n48,
 467                                            const SMDS_MeshNode * n1234,
 468                                            const SMDS_MeshNode * n1256,
 469                                            const SMDS_MeshNode * n2367,
 470                                            const SMDS_MeshNode * n3478,
 471                                            const SMDS_MeshNode * n1458,
 472                                            const SMDS_MeshNode * n5678,
 473                                            const SMDS_MeshNode * nCenter,
 474                                            int ID);
 475   virtual SMDS_MeshVolume* AddVolume(const SMDS_MeshNode * n1,
 476                                      const SMDS_MeshNode * n2,
 477                                      const SMDS_MeshNode * n3,
 478                                      const SMDS_MeshNode * n4,
 479                                      const SMDS_MeshNode * n5,
 480                                      const SMDS_MeshNode * n6,
 481                                      const SMDS_MeshNode * n7,
 482                                      const SMDS_MeshNode * n8,
 483                                      const SMDS_MeshNode * n12,
 484                                      const SMDS_MeshNode * n23,
 485                                      const SMDS_MeshNode * n34,
 486                                      const SMDS_MeshNode * n41,
 487                                      const SMDS_MeshNode * n56,
 488                                      const SMDS_MeshNode * n67,
 489                                      const SMDS_MeshNode * n78,
 490                                      const SMDS_MeshNode * n85,
 491                                      const SMDS_MeshNode * n15,
 492                                      const SMDS_MeshNode * n26,
 493                                      const SMDS_MeshNode * n37,
 494                                      const SMDS_MeshNode * n48,
 495                                      const SMDS_MeshNode * n1234,
 496                                      const SMDS_MeshNode * n1256,
 497                                      const SMDS_MeshNode * n2367,
 498                                      const SMDS_MeshNode * n3478,
 499                                      const SMDS_MeshNode * n1458,
 500                                      const SMDS_MeshNode * n5678,
 501                                      const SMDS_MeshNode * nCenter);
 502
 503   virtual SMDS_MeshFace* AddPolygonalFaceWithID (const std::vector<int>& nodes_ids,
 504                                                  const int               ID);
 505
 506   virtual SMDS_MeshFace* AddPolygonalFaceWithID (const std::vector<const SMDS_MeshNode*>& nodes,
 507                                                  const int                                ID);
 508
 509   virtual SMDS_MeshFace* AddPolygonalFace (const std::vector<const SMDS_MeshNode*>& nodes);
 510
 511   virtual SMDS_MeshVolume* AddPolyhedralVolumeWithID
 512                            (const std::vector<int>& nodes_ids,
 513                             const std::vector<int>& quantities,
 514                             const int               ID);
 515
 516   virtual SMDS_MeshVolume* AddPolyhedralVolumeWithID
 517                            (const std::vector<const SMDS_MeshNode*>& nodes,
 518                             const std::vector<int>&                  quantities,
 519                             const int                                ID);
 520
 521   virtual SMDS_MeshVolume* AddPolyhedralVolume
 522                            (const std::vector<const SMDS_MeshNode*>& nodes,
 523                             const std::vector<int>&                  quantities);
 524
 525   virtual void MoveNode(const SMDS_MeshNode *, double x, double y, double z);
 526   virtual void RemoveNode(const SMDS_MeshNode *);
 527   void RemoveElement(const SMDS_MeshElement *);
 528
 529   /*! Remove only the given element/node and only if it is free.
 530    *  Methods do not work for meshes with descendants.
 531    *  Implemented for fast cleaning of meshes.
 532    */
 533   void RemoveFreeNode   (const SMDS_MeshNode *,    SMESHDS_SubMesh *, bool fromGroups=true);
 534   void RemoveFreeElement(const SMDS_MeshElement *, SMESHDS_SubMesh *, bool fromGroups=true);
 535
 536   void ClearMesh();
 537
 538   bool ChangeElementNodes(const SMDS_MeshElement * elem,
 539                           const SMDS_MeshNode    * nodes[],
 540                           const int                nbnodes);
 541   bool ChangePolygonNodes(const SMDS_MeshElement * elem,
 542                           std::vector<const SMDS_MeshNode*> nodes);
 543   bool ChangePolyhedronNodes(const SMDS_MeshElement * elem,
 544                              std::vector<const SMDS_MeshNode*> nodes,
 545                              std::vector<int>                  quantities);
 546   bool ModifyCellNodes(int smdsVolId, std::map<int,int> localClonedNodeIds);
 547   void Renumber (const bool isNodes, const int startID=1, const int deltaID=1);
 548
 549   void SetNodeInVolume(const SMDS_MeshNode * aNode, const TopoDS_Shell & S);
 550   void SetNodeInVolume(const SMDS_MeshNode * aNode, const TopoDS_Solid & S);
 551   void SetNodeOnFace  (const SMDS_MeshNode * aNode, const TopoDS_Face& S, double u=0.,double v=0.);
 552   void SetNodeOnEdge  (const SMDS_MeshNode * aNode, const TopoDS_Edge& S, double u=0.);
 553   void SetNodeOnVertex(const SMDS_MeshNode * aNode, const TopoDS_Vertex & S);
 554   void UnSetNodeOnShape(const SMDS_MeshNode * aNode);
 555   void SetMeshElementOnShape  (const SMDS_MeshElement * anElt, const TopoDS_Shape & S);
 556   void UnSetMeshElementOnShape(const SMDS_MeshElement * anElt, const TopoDS_Shape & S);
 557   void SetNodeInVolume(const SMDS_MeshNode * aNode, int Index);
 558   void SetNodeOnFace  (const SMDS_MeshNode * aNode, int Index, double u=0., double v=0.);
 559   void SetNodeOnEdge  (const SMDS_MeshNode * aNode, int Index, double u=0.);
 560   void SetNodeOnVertex(const SMDS_MeshNode * aNode, int Index);
 561   void SetMeshElementOnShape(const SMDS_MeshElement * anElt, int Index);
 562   bool HasMeshElements(const TopoDS_Shape & S) const;
 563   SMESHDS_SubMesh * MeshElements(const TopoDS_Shape & S) const;
 564   SMESHDS_SubMesh * MeshElements(const int Index) const;
 565   std::list<int> SubMeshIndices() const;
 566   const std::map<int,SMESHDS_SubMesh*>& SubMeshes() const
 567   { return myShapeIndexToSubMesh; }
 568   const TopoDS_Shape& GetCurrentSubShape() const { return myCurSubShape; }
 569
 570   bool HasHypothesis(const TopoDS_Shape & S);
 571   const std::list<const SMESHDS_Hypothesis*>& GetHypothesis(const TopoDS_Shape & S) const;
 572   bool IsUsedHypothesis(const SMESHDS_Hypothesis * H) const;
 573   SMESHDS_Script * GetScript();
 574   void ClearScript();
 575   int ShapeToIndex(const TopoDS_Shape & aShape) const;
 576   const TopoDS_Shape& IndexToShape(int ShapeIndex) const;
 577   int MaxShapeIndex() const { return myIndexToShape.Extent(); }
 578   int MaxSubMeshIndex() const;
 579
 580   SMESHDS_SubMesh * NewSubMesh(int Index);
 581   int AddCompoundSubmesh(const TopoDS_Shape& S, TopAbs_ShapeEnum type = TopAbs_SHAPE);
 582
 583   // Groups. SMESHDS_Mesh is not an owner of groups
 584   void AddGroup (SMESHDS_GroupBase* theGroup)      { myGroups.insert(theGroup); }
 585   void RemoveGroup (SMESHDS_GroupBase* theGroup)   { myGroups.erase(theGroup); }
 586   int GetNbGroups() const                      { return myGroups.size(); }
 587   const std::set<SMESHDS_GroupBase*>& GetGroups() const { return myGroups; }
 588
 589   bool IsGroupOfSubShapes (const TopoDS_Shape& aSubShape) const;
 590
 591   virtual void compactMesh();
 592   void CleanDownWardConnectivity();
 593   void BuildDownWardConnectivity(bool withEdges);
 594
 595   ~SMESHDS_Mesh();
 596
 597 private:
 598   void addNodeToSubmesh( const SMDS_MeshNode* aNode, int Index )
 599   {
 600     //Update or build submesh
 601     std::map<int,SMESHDS_SubMesh*>::iterator it = myShapeIndexToSubMesh.find( Index );
 602     if ( it == myShapeIndexToSubMesh.end() )
 603       it = myShapeIndexToSubMesh.insert( std::make_pair(Index, new SMESHDS_SubMesh(this, Index) )).first;
 604     it->second->AddNode( aNode ); // add aNode to submesh
 605   }
 606
 607   typedef std::list<const SMESHDS_Hypothesis*> THypList;
 608
 609   typedef NCollection_DataMap< TopoDS_Shape, THypList, SMESHDS_Hasher > ShapeToHypothesis;
 610
 611   ShapeToHypothesis          myShapeToHypothesis;
 612
 613   int                        myMeshID, myPersistentID;
 614   TopoDS_Shape               myShape;
 615
 616   typedef std::map<int,SMESHDS_SubMesh*> TShapeIndexToSubMesh;
 617   TShapeIndexToSubMesh myShapeIndexToSubMesh;
 618
 619   TopTools_IndexedMapOfShape myIndexToShape;
 620
 621   typedef std::set<SMESHDS_GroupBase*> TGroups;
 622   TGroups myGroups;
 623
 624   SMESHDS_Script*            myScript;
 625   bool                       myIsEmbeddedMode;
 626
 627   // optimize addition of nodes/elements to submeshes by, SetNodeInVolume() etc:
 628   // avoid search of submeshes in maps
 629   bool add( const SMDS_MeshElement* elem, SMESHDS_SubMesh* subMesh );
 630   SMESHDS_SubMesh* getSubmesh( const TopoDS_Shape & shape);
 631   SMESHDS_SubMesh* getSubmesh( const int            Index );
 632   int                        myCurSubID;
 633   TopoDS_Shape               myCurSubShape;
 634   SMESHDS_SubMesh*           myCurSubMesh;
 635 };
 636
 637
 638 #endif