2 // Class : Gestion des tutaux en T (like castles in scotland)
4 // Copyright (C) 2009-2019 CEA/DEN, EDF R&D
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.
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.
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
20 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
22 #ifndef __BI_CYLINDER_H_
23 #define __BI_CYLINDER_H_
25 #include "HexElements.hxx"
30 class HexaExport BiCylinder : public Elements
33 BiCylinder (Document* doc);
35 Hexa* getHexaIJK (int part, int nx, int ny, int nz);
36 Quad* getQuadIJ (int part, int nx, int ny, int nz);
37 Quad* getQuadJK (int part, int nx, int ny, int nz);
38 Quad* getQuadIK (int part, int nx, int ny, int nz);
39 Edge* getEdgeI (int part, int nx, int ny, int nz);
40 Edge* getEdgeJ (int part, int nx, int ny, int nz);
41 Edge* getEdgeK (int part, int nx, int ny, int nz);
42 Vertex* getVertexIJK (int part, int nx, int ny, int nz);
46 virtual ~BiCylinder () {}
47 // Tutaux en T (like castles in scotland)
48 int makeCylinders (Vertex* ori1, double* vz1, double ray1, double h1,
49 Vertex* ori2, double* vz2, double ray2, double h2);
51 int makePipes (Vertex* ori1, double* z1, double ri1, double rx1, double h1,
52 Vertex* ori2, double* z2, double ri2, double rx2, double h2,
55 virtual int countHexa () { return nbr_hexas ; }
56 virtual int countQuad () { return nbr_quads ; }
57 virtual int countEdge () { return nbr_edges ; }
58 virtual int countVertex () { return nbr_vertex ; }
60 // int crossCylinders (Cylinder* cyl1, Cylinder* cyl2); // perime
64 enum { NbrSlices=5, NbrVslices=NbrSlices+1, NbrCotes=4 };
66 void createLittlePipe ();
67 void createBigPipe ();
68 void adjustLittleSlice (int nx, int nzs);
69 void adjustLittleSquare (int nk);
71 Vertex* addVertex (double px, double py, double pz,
72 int cyl, int ni, int nj, int nk);
73 Edge* addEdge (Vertex* v1, Vertex* v2,
74 int cyl, int dir, int nx, int ny, int nz);
75 Quad* addQuad (Edge* e1, Edge* e2, Edge* e3, Edge* e4,
76 int cyl, int dir, int nx, int ny, int nz);
77 Hexa* addHexa (Quad* qa, Quad* qb, Quad* qc, Quad* qd, Quad* qe, Quad* qf,
78 int cyl, int nx, int ny, int nz);
80 Vertex* findVertex (double px, double py, double pz, int ni, int nj, int nk);
81 Edge* findEdge (Vertex* v1, Vertex* v2, int dir, int nx, int ny, int nz);
82 Quad* findQuad (Vertex* v1, Vertex* v2, int dir, int nx, int ny, int nz);
83 Quad* findQuad (Edge* e1, Edge* e2, int dir, int nx, int ny, int nz);
84 Hexa* findHexa (Quad* qa, Quad* qb, int nx, int ny, int nz);
86 Edge* findEdge (Vertex* v1, Vertex* v2); // A replacer dans HexElements
88 void assoCylinders ();
89 void assoSlice (int cyl, int nx, int nzs, double* norm);
90 void assoArc (int cyl, int nx, int ny, int nz, int subid);
91 void assoArc (int cyl, int ny, Vertex* v1, Vertex* v2, int subid);
92 int assoIntersection (int nx, int nzlice);
95 void addSlice (int cyl, int i, int k, double h, double r, bool find=false);
96 void addCarre (int cyl, int k, double h, double r, bool find=false);
97 void addCube (int cyl, int k, int i=0, int k1=-1);
98 void fillSlice (int cyl, int nia, int nka, int nib, int nkb, int nic,
99 int nkc, int nid, int nkd, bool med=false);
101 double getAngle (int nj) { return (2*((nj)%4)+1)*M_PI/4; }
102 int getKey (int cyl, int ni, int nj, int nk)
103 { return 1000*(cyl+1) + 100*ni + 10*nj + nk; }
104 int getKey (int cyl, int dd, int ni, int nj, int nk)
105 { return 10000*(cyl+1) + 1000*dd + 100*ni + 10*nj + nk; }
108 Real3 cross_oribig, cross_dirbig;
109 Real3 cross_orismall, cross_dirsmall;
111 std::map <int,int> map_vertex, map_edge, map_quad, map_hexa;
112 std::map <int,int> :: iterator it_map;
114 int nbr_hexas1, nbr_quads1, nbr_edges1, nbr_vertex1;
115 double cross_rayext [V_TWO];
116 double cross_rayint [V_TWO];
117 double cross_hauteur [V_TWO];
119 bool at_right, at_left;