4 // Copyright (C) 2009-2013 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.
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
26 #include "Hex_defines.hxx"
27 #include "HexEltBase.hxx"
40 class HEXABLOCKENGINE_EXPORT Document : public EltBase
42 // Fonctions utilisateur
44 Document (cpchar name, Hex* dad=NULL);
46 int setName (const char* name);
47 void setLevel (int niv);
48 int save (const char* ficxml); // Genere le xml dans un seul fichier
49 int saveVtk (cpchar nomfic);
53 void setTolerance (double tol) { doc_tolerance = std::max (tol, 0.0); }
54 double getTolerance () { return doc_tolerance; }
56 // Creation d'elements
57 Vertex* addVertex (double x=0.0, double y=0.0, double z=0.0);
59 Edge* addEdge (Vertex* va, Vertex* vb);
60 Edge* addEdgeVector (Vertex* va, Vector* vec);
62 Quad* addQuad (Edge* v1, Edge* v2, Edge* v3, Edge* v4);
63 Quad* addQuadVertices (Vertex* v1, Vertex* v2, Vertex* v3, Vertex* v4);
65 Hexa* addHexa (Quad* qa, Quad* qb, Quad* qc, Quad* qd, Quad* qe, Quad* qf);
66 Hexa* addHexaVertices (Vertex* v1, Vertex* v2, Vertex* v3, Vertex* v4,
67 Vertex* v5, Vertex* v6, Vertex* v7, Vertex* v8);
69 Hexa* addHexa2Quads (Quad* q1, Quad* q2);
70 Hexa* addHexa3Quads (Quad* q1, Quad* q2, Quad* q3);
71 Hexa* addHexa4Quads (Quad* q1, Quad* q2, Quad* q3, Quad* q4);
72 Hexa* addHexa5Quads (Quad* q1, Quad* q2, Quad* q3, Quad* q4, Quad* q5);
74 Vector* addVector (double dx=0.0, double dy=0.0, double dz=0.0);
75 Vector* addVectorVertices (Vertex* va, Vertex* vb);
76 Law* addLaw (const char* name, int nbnodes);
77 int addLaws (double lgmoy, bool usemax=true);
79 Group* addHexaGroup (cpchar name);
80 Group* addQuadGroup (cpchar name);
81 Group* addQuadNodeGroup (cpchar name);
82 Group* addHexaNodeGroup (cpchar name);
83 Group* addEdgeGroup (cpchar name);
84 Group* addEdgeNodeGroup (cpchar name);
85 Group* addVertexNodeGroup (cpchar name);
88 // ----------------- Transformations
89 Elements* makeTranslation (Elements* elts, Vector* trans);
90 Elements* makeScale (Elements* elts, Vertex* ver, double k);
91 Elements* makeRotation (Elements* elts, Vertex* ver, Vector* vec,
93 Elements* makeSymmetryPoint (Elements* elts, Vertex* ver);
94 Elements* makeSymmetryLine (Elements* elts, Vertex* ver, Vector* vec);
95 Elements* makeSymmetryPlane (Elements* elts, Vertex* ver, Vector* vec);
97 int performTranslation (Elements* elts, Vector* trans);
98 int performScale (Elements* elts, Vertex* ver, double k);
99 int performRotation (Elements* elts, Vertex* ver, Vector* vec, double angle);
100 int performSymmetryPoint (Elements* elts, Vertex* ver);
101 int performSymmetryLine (Elements* elts, Vertex* ver, Vector* vec);
102 int performSymmetryPlane (Elements* elts, Vertex* ver, Vector* vec);
104 // ----------------- Modifications
105 Elements* disconnectQuad (Hexa* maille, Quad* face);
106 Elements* disconnectEdge (Hexa* maille, Edge* arete);
107 Elements* disconnectVertex (Hexa* maille, Vertex* noeud);
108 Elements* disconnectEdges (Hexas thexas, Edges edges);
110 Elements* replace (Quads pattern, Vertex* p1, Vertex* c1,
111 Vertex* p2, Vertex* c2, Vertex* p3, Vertex* c3);
113 int mergeVertices (Vertex* v1, Vertex* v2);
114 int mergeEdges (Edge* e1, Edge* e2, Vertex* v1, Vertex* v2);
115 int mergeQuads (Quad* q1, Quad* q2, Vertex* v1, Vertex* v2,
116 Vertex* v3, Vertex* v4);
118 // ----------------- Associations
119 void clearAssociation ();
120 int associateOpenedLine (Edges mline, NewShapes gline, IntVector tabid,
121 double pstart, double pend);
123 int associateClosedLine (Vertex* mfirst, Edges mline,
124 NewShapes gline, IntVector tabid,
125 double pstart, bool inv);
126 // ----------------- count
131 int countVector () { return doc_vector.size(); }
133 int countGroup () { return (int) doc_group.size(); }
134 int countLaw () { return nbr_laws ; }
135 int countPropagation ();
136 int countShape () { return doc_tab_shape.size() ; }
138 int countUsedHexa ();
139 int countUsedQuad ();
140 int countUsedEdge ();
141 int countUsedVertex ();
143 // ----------------- get par indice
144 Hexa* getHexa (int nro);
145 Quad* getQuad (int nro);
146 Edge* getEdge (int nro);
147 Vertex* getVertex (int nro);
149 Hexa* getUsedHexa (int nro);
150 Quad* getUsedQuad (int nro);
151 Edge* getUsedEdge (int nro);
152 Vertex* getUsedVertex (int nro);
154 Vector* getVector (int nro);
155 NewShape* getShape (int nro);
156 Group* getGroup (int nro);
157 Law* getLaw (int nro);
158 Propagation* getPropagation (int nro);
160 // ----------------- find
161 Vertex* findVertex (double vx, double vy, double vz);
162 Edge* findEdge (Vertex* va, Vertex* vb);
163 Quad* findQuad (Vertex* va, Vertex* vb);
164 Quad* findQuad (Edge* ea, Edge* eb);
165 Hexa* findHexa (Vertex* va, Vertex* vb);
167 Group* findGroup (cpchar name);
168 Law* findLaw (cpchar name);
169 Propagation* findPropagation (Edge* arete);
171 // ----------------- Suppressions
172 int removeHexa (Hexa* maille);
173 int removeQuad (Quad* maille);
174 int removeConnectedHexa (Hexa* maille);
175 int removeElements (Elements* bloc);
176 int removeGroup (Group* grp);
177 int removeLaw (Law* lau);
180 // ----------------- Creations Hexav6
181 Elements* makeCartesianTop (int nx, int ny, int nz);
182 Elements* makeCartesianUni (Vertex* ori, Vector* vx, Vector* vy, Vector* vz,
183 double lx, double ly, double lz,
184 int nx, int ny, int nz);
185 Elements* makeCartesian (Vertex* ori, Vector* vx, Vector* vy, Vector* vz,
186 RealVector tlx, RealVector tly, RealVector tlz);
188 Elements* makeCylinderTop (int nr, int na, int nh);
189 Elements* makeCylinderUni (Vertex* orig, Vector* vx, Vector* vz,
190 double rint, double rext, double ang, double haut,
191 int nr, int na, int nh);
192 Elements* makeCylinder (Vertex* orig, Vector* vx, Vector* vz,
193 RealVector tray, RealVector tang, RealVector thaut);
195 Elements* makePipeTop (int nr, int na, int nh);
196 Elements* makePipeUni (Vertex* orig, Vector* vx, Vector* vz,
197 double rint, double rext, double angle, double haut,
198 int nr, int na, int nh);
199 Elements* makePipe (Vertex* orig, Vector* vx, Vector* vz,
200 RealVector tray, RealVector tang, RealVector thaut);
202 Elements* makeSphericalTop (int nbre, int crit=0);
203 Elements* makeSphericalUni (Vertex* centre, Vector* vx, Vector* vz,
204 double rayon, int nbre, int crit=0);
205 Elements* makeSpherical (Vertex* centre, Vector* vx, Vector* vz,
206 RealVector rayon, int crit=0);
207 Elements* makeSphereTop (int nr, int na, int nh);
208 Elements* makeSphereUni (Vertex* centre, Vector* vx, Vector* vz,
209 double rtrou, double rext, double ang,
210 Vertex* vplan, int nr, int na, int nh);
211 Elements* makeSphere (Vertex* centre, Vector* vx, Vector* vz,
212 RealVector tray, RealVector tang, RealVector thaut);
214 Elements* makeRindTop (int nr, int na, int nh);
215 Elements* makeRindUni (Vertex* centre, Vector* vx, Vector* vz,
216 double raytrou, double rint, double rext,
217 double ang, Vertex* vplan, int nr, int na, int nh);
218 Elements* makeRind (Vertex* centre, Vector* vx, Vector* vz,
219 RealVector tray, RealVector tang, RealVector thaut);
221 BiCylinder* makeCylinders (Vertex* ori1, Vector* z1, double r1, double h1,
222 Vertex* ori2, Vector* z2, double r2, double h2);
223 BiCylinder* makePipes (Vertex* ori1, Vector* z1, double rint1, double rex1,
224 double h1, Vertex* ori2, Vector* z2, double rint2,
225 double rex2, double h2);
227 Elements* extrudeQuadTop (Quad* start, int nbre);
228 Elements* extrudeQuadUni (Quad* start, Vector* dv, double len, int nbre);
229 Elements* extrudeQuad (Quad* start, Vector* dv, RealVector tlen);
231 Elements* extrudeQuadsTop (Quads start, int nbre);
232 Elements* extrudeQuadsUni (Quads start, Vector* axis, double len, int nbre);
233 Elements* extrudeQuads (Quads start, Vector* axis, RealVector tlen);
235 Elements* revolutionQuadUni (Quad* start, Vertex* center, Vector* axis,
236 double angle, int nbre);
237 Elements* revolutionQuad (Quad* start, Vertex* center, Vector* axis,
239 Elements* revolutionQuadsUni (Quads start, Vertex* center, Vector* axis,
240 double angle, int nbre);
241 Elements* revolutionQuads (Quads start, Vertex* center, Vector* axis,
244 Elements* joinQuadUni (Quad* start, Quad* dest, Vertex* v1, Vertex* v2,
245 Vertex* v3, Vertex* v4, int nb);
247 Elements* joinQuadsUni (Quads start, Quad* dest, Vertex* v1, Vertex* v2,
248 Vertex* v3, Vertex* v4, int nb);
250 Elements* joinQuad (Quad* start, Quad* dest, Vertex* va1, Vertex* vb1,
251 Vertex* va2, Vertex* vb2, RealVector& tlen);
253 Elements* joinQuads (Quads start, Quad* dest, Vertex* va1, Vertex* vb1,
254 Vertex* va2, Vertex* vb2, RealVector& tlen);
256 Elements* cutUni (Edge* eddge, int nbre);
257 Elements* cut (Edge* eddge, RealVector& tlen);
259 Group* addGroup (cpchar name, EnumGroup kind);
260 // Fonctions globales
261 bool isSaved () { return doc_saved ; }
263 int appendXml (pfile fstudy); // Genere le xml dans un fichier ouvert
264 cpchar getXml (); // Genere le xml et rend un flux
266 int getLevel () { return doc_db ; }
268 Vertex* findVertex (int id);
270 int closeQuads (Quad* q1, Quad* q2);
272 // ---------------------------------------------------
274 Law* addLaw (Law* law); // lo-add-lololo
276 // ---------------------------------------------------
278 int checkAssociations ();
279 // addShape version Python
280 NewShape* addShape (long forme, const char* name);
283 // addShape HorsPython
284 NewShape* addShape (TopoDS_Shape& forme, const char* name);
287 void dumpPropagation ();
289 Document* copyDocument ();
291 NewShape* getCloud () { return doc_cloud ; }
292 NewShape* addShape (const char* name, EnumShape type);
293 NewShape* findShape (rcstring name);
295 int associateLine (Vertex* mfirst, Edges& mline,
296 NewShapes& gline, IntVector& tabid,
297 double pstart, double pend, bool inv);
299 int clearAssociation (EnumElt type);
300 void clearAssoVertices ();
301 void clearAssoEdges ();
302 void clearAssoQuads ();
307 int loadXml (cpchar name);
308 int setXml (cpchar flux);
309 int setXml (cpchar flux, int& position);
311 EltBase* getFirstEltBase (EnumElt type) { return doc_first_elt [type]; }
312 EltBase* getLastEltBase (EnumElt type) { return doc_last_elt [type]; }
313 int getNbrElt (EnumElt type) { return doc_nbr_elt [type]; }
314 void setDeprecated (int level=1);
316 void markAll (int marque, int type=-1);
317 int saveVtk0 (cpchar nomfic);
318 int saveVtk (cpchar radical, int &nro);
320 void putError (cpchar mess, cpchar info1=NULL, cpchar info2=NULL);
321 void hputError (cpchar mess, EltBase* e1, EltBase* e2=NULL);
323 void majReferences (); // M.A.J relation "utilise par"
324 void update () { if (maj_connection) majReferences (); }
325 bool debug (int niv=0) { return doc_db > niv ; }
327 void reorderFaces ();
328 void reorderQuads ();
330 // --------------------------------------------------- Evols Hexa5
331 cpchar getNextName (EnumElt type, string& name);
332 string getNextName (EnumElt type);
334 virtual char* makeVarName (char* name);
339 // --------------------------------------------------- Evols Hexa6
341 void addHexa (Hexa* elt) { doc_hexas .push_back (elt) ; }
342 void addQuad (Quad* elt) { doc_quads .push_back (elt) ; }
343 void addEdge (Edge* elt) { doc_edges .push_back (elt) ; }
344 void addVertex (Vertex* elt) { doc_vertex.push_back (elt) ; }
351 int parseXml (XmlTree& xml);
353 void replaceVertex (Vertex* v1, Vertex* v2);
354 void replaceEdge (Edge* e1, Edge* e2);
355 int replaceQuad (Quad* q1, Quad* q2);
357 int countElement (EnumElt type);
358 EltBase* getElement (EnumElt type, int nro);
360 void majPropagation ();
364 Elements* clonerElements (Elements* table, Matrix* matrice);
366 // --------------------------------------------------- Evols Hexa3
367 Hexa* addHexaQuadsAB (AnaQuads& strquads);
368 Hexa* addHexaQuadsAC (AnaQuads& strquads);
369 Hexa* addHexaQuadsACE (AnaQuads& strquads);
370 Hexa* addHexaQuadsACD (AnaQuads& strquads);
371 Hexa* addHexaQuadsABCD (AnaQuads& strquads);
372 Hexa* addHexaQuadsACDE (AnaQuads& strquads);
374 EltBase* findElement (EnumElt type, int ident);
375 // --------------------------------------------------- Evols Hexa5
377 int sortShapes (NewShapes& tshapes, IntVector& subid,
378 bool closed, bool inv, double pstart, double pend,
379 EdgeShapes& glines, IntVector& sens, double& lgtotale);
380 int sortEdges (Vertex* vfirst, Edges& mline, Edges& edges, IntVector& sens);
382 void parseShapes (XmlTree& root);
383 void parseAssociation (XmlTree* root, Quad* elt);
384 void parseAssociation (XmlTree* root, Edge* elt);
385 void parseAssociation (XmlTree* root, Vertex* elt);
386 SubShape* parseSubShape (XmlTree* root);
389 friend class EltBase;
400 EltBase* doc_first_elt [EL_MAXI];
401 EltBase* doc_last_elt [EL_MAXI];
402 int doc_nbr_elt [EL_MAXI];
403 double doc_tolerance;
405 std::vector <Law*> doc_laws;
409 int nbr_propagations;
410 bool maj_propagation;
412 std::vector <Propagation*> doc_propagation;
413 std::vector <Group*> doc_group;
414 std::vector <Vector*> doc_vector;
417 // --------------------------------------------------- Evols Hexa3
419 std::vector <Hexa*> doc_used_hexas;
420 std::vector <Quad*> doc_used_quads;
421 std::vector <Edge*> doc_used_edges;
422 std::vector <Vertex*> doc_used_vertex;
429 // --------------------------------------------------- Evols Hexa5
431 static int doc_number;
435 NewShape* doc_cloud; // Nuage de points du doc;
436 std::vector <NewShape*> doc_tab_shape;
438 // --------------------------------------------------- Evols internes Hexa6
439 std::vector <Hexa*> doc_hexas;
440 std::vector <Quad*> doc_quads;
441 std::vector <Edge*> doc_edges;
442 std::vector <Vertex*> doc_vertex;