src/HEXABLOCK/HexBiCylinder.hxx

   1
   2 // Class : Gestion des tutaux en T (like castles in scotland)
   3
   4 // Copyright (C) 2009-2013  CEA/DEN, EDF R&D
   5 //
   6 // This library is free software; you can redistribute it and/or
   7 // modify it under the terms of the GNU Lesser General Public
   8 // License as published by the Free Software Foundation; either
   9 // version 2.1 of the License.
  10 //
  11 // This library is distributed in the hope that it will be useful,
  12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14 // Lesser General Public License for more details.
  15 //
  16 // You should have received a copy of the GNU Lesser General Public
  17 // License along with this library; if not, write to the Free Software
  18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  19 //
  20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  21 //
  22 #ifndef __BI_CYLINDER_H_
  23 #define __BI_CYLINDER_H_
  24
  25 #include "HexElements.hxx"
  26 #include <map>
  27
  28 BEGIN_NAMESPACE_HEXA
  29
  30 class BiCylinder : public Elements
  31 {
  32 public:
  33    virtual int countHexa   () { return nbr_hexas ; }
  34    virtual int countQuad   () { return nbr_quads ; }
  35    virtual int countEdge   () { return nbr_edges ; }
  36    virtual int countVertex () { return nbr_vertex ; }
  37
  38 public:
  39    BiCylinder (Document* doc);
  40
  41    virtual ~BiCylinder () {}
  42
  43    Hexa*   getHexaIJK   (int part, int nx, int ny, int nz);
  44    Quad*   getQuadIJ    (int part, int nx, int ny, int nz);
  45    Quad*   getQuadJK    (int part, int nx, int ny, int nz);
  46    Quad*   getQuadIK    (int part, int nx, int ny, int nz);
  47    Edge*   getEdgeI     (int part, int nx, int ny, int nz);
  48    Edge*   getEdgeJ     (int part, int nx, int ny, int nz);
  49    Edge*   getEdgeK     (int part, int nx, int ny, int nz);
  50    Vertex* getVertexIJK (int part, int nx, int ny, int nz);
  51
  52                      // Tutaux en T (like castles in scotland)
  53    int  crossCylinders (Cylinder* cyl1, Cylinder* cyl2);
  54
  55 private :
  56
  57    enum { NbrSlices=5, NbrVslices=NbrSlices+1, NbrCotes=4 };
  58
  59    void createLittleCyl ();
  60    void createBigCyl    ();
  61    void adjustLittleSlice (int nx, int nzs);
  62
  63    Vertex* addVertex (double px, double py, double pz,
  64                       int cyl, int ni, int nj, int nk);
  65    Edge* addEdge (Vertex* v1, Vertex* v2,
  66                   int cyl, int dir, int nx, int ny, int nz);
  67    Quad* addQuad (Edge* e1, Edge* e2, Edge* e3, Edge* e4,
  68                   int cyl, int dir, int nx, int ny, int nz);
  69    Hexa* addHexa (Quad* qa, Quad* qb, Quad* qc, Quad* qd, Quad* qe, Quad* qf,
  70                   int cyl, int nx, int ny, int nz);
  71
  72    Vertex* findVertex (double px, double py, double pz, int ni, int nj, int nk);
  73    Edge*   findEdge (Vertex* v1, Vertex* v2, int dir, int nx, int ny, int nz);
  74    Quad*   findQuad (Vertex* v1, Vertex* v2, int dir, int nx, int ny, int nz);
  75    Quad*   findQuad (Edge* e1, Edge* e2, int dir, int nx, int ny, int nz);
  76    Hexa*   findHexa (Quad* qa, Quad* qb, int nx, int ny, int nz);
  77
  78    Edge* findEdge (Vertex* v1, Vertex* v2);
  79    Quad* findQuad (Edge*   e1, Edge*   e2);
  80
  81    void assoCylinder  (int cyl, double* normal);
  82    void assoCylinders (double* snormal, double* gnormal);
  83    void assoSlice    (int cyl, int nx, int nzs, double* norm);
  84    void assoArc (int cyl, int nx, int ny, int nz, string& brep, double ray);
  85    int  assoIntersection (int nx, int nzlice, double* snorm, double* bnorm);
  86
  87
  88    void addSlice  (int cyl, int i, int k, double h, double r, bool find=false);
  89    void fillSlice (int cyl, int nia, int nka, int nib, int nkb, int nic,
  90                             int nkc, int nid, int nkd, bool med=false);
  91    void fillSmallCyl ();
  92
  93    double getAngle (int nj)             { return  (2*((nj)%4)+1)*M_PI/4; }
  94    int getKey (int cyl, int ni, int nj, int nk)
  95               { return 1000*(cyl+1) + 100*ni + 10*nj + nk; }
  96    int getKey (int cyl, int dd, int ni, int nj, int nk)
  97               { return 10000*(cyl+1) + 1000*dd + 100*ni + 10*nj + nk; }
  98 private :
  99     Cylinder* cross_cyl1;    // Petit
 100     Cylinder* cross_cyl2;    // gros
 101     Real3     cross_center;
 102
 103     map <int,int> map_vertex, map_edge, map_quad, map_hexa;
 104     map <int,int> :: iterator it_map;
 105
 106     int    nbr_hexas1, nbr_quads1, nbr_edges1, nbr_vertex1;
 107     double cross_rayext  [V_TWO];
 108     double cross_rayint  [V_TWO];
 109     double cross_hauteur [V_TWO];
 110
 111     bool   at_right, at_left;
 112 };
 113 END_NAMESPACE_HEXA
 114 #endif