]> SALOME platform Git repositories - modules/hydro.git/blob - src/HYDROData/HYDROData_Transform.cxx
Salome HOME
Access methods for object groups added.
[modules/hydro.git] / src / HYDROData / HYDROData_Transform.cxx
1 #include <HYDROData_Transform.h>
2
3 #include <math.h>
4
5 #include <Bnd_Box.hxx>
6 #include <gp_Vec.hxx>
7
8 #include <TopoDS_Shape.hxx>
9 #include <TopoDS_Iterator.hxx>
10
11 #include <BRepBndLib.hxx>
12
13 #include <TopTools_ListIteratorOfListOfShape.hxx>
14
15 #include <BOPCol_ListOfShape.hxx>
16 #include <BOPCol_MapOfShape.hxx>
17
18 #include <BOPCol_ListOfShape.hxx>
19
20 #include <BOPAlgo_BuilderShape.hxx>
21 #include <BOPAlgo_PaveFiller.hxx>
22
23 //#define DEB_SPLIT_TO_ZONES_CHECK_PARTITION 1
24 #ifdef DEB_SPLIT_TO_ZONES_CHECK_PARTITION
25 #include <BRepTools.hxx>
26 #include <TCollection_AsciiString.hxx>
27 static TCollection_AsciiString fileNameAfter("AfterTranslation");
28 #endif
29
30 void HYDROData_Transform::SetFileName(const TCollection_AsciiString& theNameAfter)
31 {
32 #ifdef DEB_SPLIT_TO_ZONES_CHECK_PARTITION
33   fileNameAfter = theNameAfter;
34 #endif
35 }
36 //=======================================================================
37 //function : HYDROData_Transform
38 //purpose  : 
39 //=======================================================================
40 HYDROData_Transform::HYDROData_Transform()
41 :
42   BOPAlgo_BuilderShape()
43 {
44   myToTransform=Standard_False;
45   myTreshold=10000.;
46   myTransform1=NULL;
47   myTransform2=NULL;
48   myBuilder=NULL;
49 }
50 //=======================================================================
51 //function : ~
52 //purpose  : 
53 //=======================================================================
54 HYDROData_Transform::~HYDROData_Transform()
55 {
56   Clear();
57 }
58 //=======================================================================
59 //function : Clear
60 //purpose  : 
61 //=======================================================================
62 void HYDROData_Transform::Clear()
63 {
64   myTreshold=10000.;
65   //
66   if (myTransform1) {
67     delete myTransform1;
68     myTransform1=NULL;
69   }
70   if (myTransform2) {
71     delete myTransform2;
72     myTransform2=NULL;
73   }
74   if (myBuilder) {
75     delete myBuilder;
76     myBuilder=NULL;
77   }
78   //
79   BOPAlgo_BuilderShape::PrepareHistory();
80 }
81 //=======================================================================
82 //function : SetArgument
83 //purpose  : 
84 //=======================================================================
85 void HYDROData_Transform::SetArgument(const TopoDS_Shape& theShape)
86 {
87   myArgument=theShape;
88 }
89 //=======================================================================
90 //function : AddArgument
91 //purpose  : 
92 //=======================================================================
93 const TopoDS_Shape& HYDROData_Transform::Argument()const
94 {
95   return myArgument;
96 }
97 //=======================================================================
98 //function : SetTreshold
99 //purpose  : 
100 //=======================================================================
101 void HYDROData_Transform::SetTreshold(const Standard_Real theTreshold)
102 {
103   myTreshold=theTreshold;
104 }
105 //=======================================================================
106 //function : Treshold
107 //purpose  : 
108 //=======================================================================
109 Standard_Real HYDROData_Transform::Treshold()const 
110 {
111   return myTreshold;
112 }
113 //=======================================================================
114 //function : SetToTransform
115 //purpose  : 
116 //=======================================================================
117 void HYDROData_Transform::SetToTransform(const Standard_Boolean theFlag)
118 {
119   myToTransform=theFlag;
120 }
121 //=======================================================================
122 //function : ToTransform
123 //purpose  : 
124 //=======================================================================
125 Standard_Boolean HYDROData_Transform::ToTransform()const
126 {
127   return myToTransform;
128 }
129 //=======================================================================
130 //function : CheckData
131 //purpose  : 
132 //=======================================================================
133 void HYDROData_Transform::CheckData()
134 {
135   if(myArgument.IsNull()) {
136     myErrorStatus=10;
137   }
138 }
139 //=======================================================================
140 //function : Prepare
141 //purpose  : 
142 //=======================================================================
143 void HYDROData_Transform::Prepare()
144 {
145   TopoDS_Shape aSx;
146   //
147   myErrorStatus=0;
148   myWarningStatus=0;
149   myShape=aSx;
150 }
151 //=======================================================================
152 //function : Perform
153 //purpose  : 
154 //=======================================================================
155 void HYDROData_Transform::Perform()
156 {
157   Prepare();
158   //
159   CheckData();
160   if(myErrorStatus) {
161     return;
162   }
163   //
164   Clear();
165   if(myErrorStatus) {
166     return;
167   }
168   //
169   ComputeTrsf();
170   if(myErrorStatus) {
171     return;
172   }
173   //
174   Transform1();
175   if(myErrorStatus) {
176     return;
177   }
178   //
179   BuildFuse();
180   if(myErrorStatus) {
181     return;
182   }
183   //
184   Transform2();
185 }
186 //=======================================================================
187 //function : Detect
188 //purpose  : 
189 //=======================================================================
190 void HYDROData_Transform::Detect()
191 {
192   Standard_Real aXmin, aYmin, aZmin, aXmax, aYmax, aZmax;
193   Standard_Real aXc, aYc, aZc;
194   Bnd_Box aBox; 
195   gp_Vec aVec;
196   //
197   myErrorStatus=0;
198   //
199   BRepBndLib::Add(myArgument, aBox);
200   //
201   aBox.Get(aXmin, aYmin, aZmin, aXmax, aYmax, aZmax);
202   //
203   aXc=0.5*(aXmin+aXmax);
204   aYc=0.5*(aYmin+aYmax);
205   aZc=0.5*(aZmin+aZmax);
206   //printf(" aVc : { %lf %lf %lf }\n", aXc, aYc, aZc);
207   //
208   myToTransform=Standard_False;
209   if (fabs(aXc)>myTreshold || 
210       fabs(aYc)>myTreshold || 
211       fabs(aZc)>myTreshold) {
212     myToTransform=!myToTransform;
213   }  
214   //
215   aVec.SetCoord(-aXc, -aYc, -aZc);
216   myTrsf1.SetTranslation(aVec);
217   //
218   aVec.SetCoord(aXc, aYc, aZc);
219   myTrsf2.SetTranslation(aVec);
220 }
221 //=======================================================================
222 //function : ComputeTrsf
223 //purpose  : 
224 //=======================================================================
225 void HYDROData_Transform::ComputeTrsf()
226 {
227   if (!myToTransform) {
228     gp_Trsf aTrsf;
229     //
230     myTrsf1=aTrsf;
231     myTrsf2=aTrsf;
232   }
233 }
234 //=======================================================================
235 //function : Transform1
236 //purpose  : 
237 //=======================================================================
238 void HYDROData_Transform::Transform1()
239
240   myErrorStatus=0;
241   //
242   if (myToTransform) {//pkv ft
243     Standard_Boolean bIsDone;
244     //
245     myErrorStatus=0;
246     //
247     myTransform1=new BRepBuilderAPI_Transform(myTrsf1);
248     //
249     myTransform1->Perform(myArgument);
250     bIsDone=myTransform1->IsDone();
251     if (!bIsDone) {
252       myErrorStatus=11;
253     }
254   }
255 }
256 //=======================================================================
257 //function : BuildFuse
258 //purpose  : 
259 //=======================================================================
260 void HYDROData_Transform::BuildFuse()
261 {
262   Standard_Integer iErr;
263   TopoDS_Iterator aIt;
264   BOPCol_ListOfShape aLC;
265   BOPAlgo_PaveFiller aPF;
266   //
267   myErrorStatus=0;
268   //
269    //modified by NIZNHY-PKV Wed Jan 15 13:24:48 2014f
270   const TopoDS_Shape& aS1=
271     (myToTransform)? myTransform1->Shape() : myArgument;
272   //const TopoDS_Shape& aS1=()Transform1->Shape();
273   //modified by NIZNHY-PKV Wed Jan 15 13:24:52 2014t
274
275 #ifdef DEB_SPLIT_TO_ZONES_CHECK_PARTITION
276   TCollection_AsciiString aNameBefore = fileNameAfter + ".brep";
277   BRepTools::Write(aS1, aNameBefore.ToCString());
278 #endif
279   //
280   aIt.Initialize(aS1);
281   for (; aIt.More(); aIt.Next()) {
282     const TopoDS_Shape& aS1x=aIt.Value();
283     aLC.Append(aS1x);
284   }
285   //
286   aPF.SetArguments(aLC);
287   //
288   aPF.Perform();
289   iErr=aPF.ErrorStatus();
290   if (iErr) {
291     myErrorStatus=20;
292     return;
293   }
294   //
295   myBuilder=new BOPAlgo_Builder;
296   //
297   aIt.Initialize(aS1);
298   for (; aIt.More(); aIt.Next()) {
299     const TopoDS_Shape& aS1x=aIt.Value();
300     myBuilder->AddArgument(aS1x);
301   }
302   //
303   myBuilder->PerformWithFiller(aPF);
304   iErr=myBuilder->ErrorStatus();
305   if (iErr) {
306     myErrorStatus=21;
307     return;
308   }
309 }
310 //=======================================================================
311 //function : Transform2
312 //purpose  : 
313 //=======================================================================
314 void HYDROData_Transform::Transform2()
315 {
316   myErrorStatus=0;
317   //
318   const TopoDS_Shape& aR1=myBuilder->Shape();
319   //
320   myShape=aR1; //pkv ft
321   if (myToTransform) { //pkv ft
322     Standard_Boolean bIsDone;
323     //
324     myTransform2=new BRepBuilderAPI_Transform(myTrsf2);
325     //
326     myTransform2->Perform(aR1);
327     bIsDone=myTransform2->IsDone();
328     if (!bIsDone) {
329       myErrorStatus=12;
330       return;
331     }
332     //
333     const TopoDS_Shape& aR2=myTransform2->Shape();
334     //
335     myShape=aR2;
336   }
337 }
338 //
339 //=======================================================================
340 //function : Modified
341 //purpose  : 
342 //=======================================================================
343 const TopTools_ListOfShape& HYDROData_Transform::Modified
344   (const TopoDS_Shape& aS)
345 {
346   TopTools_ListOfShape* pHistShapes;
347   TopTools_ListIteratorOfListOfShape aItLS;
348   //
349   pHistShapes=(TopTools_ListOfShape*)&myHistShapes;
350   pHistShapes->Clear();
351   //
352   if (myToTransform) { //pkv ft
353     const TopTools_ListOfShape& aLS1=myTransform1->Modified(aS);
354     const TopoDS_Shape& aS1=aLS1.First();
355     //
356     const TopTools_ListOfShape& aLS1B=myBuilder->Modified(aS1);
357     //
358     if (aLS1B.IsEmpty()) { //pkv ft
359       const TopTools_ListOfShape& aLS2B=myTransform2->Modified(aS1);
360       const TopoDS_Shape& aS2B=aLS2B.First();
361       pHistShapes->Append(aS2B);
362     }
363     else {
364       aItLS.Initialize(aLS1B);
365       for (; aItLS.More(); aItLS.Next()) {
366         const TopoDS_Shape& aS1B=aItLS.Value();
367         const TopTools_ListOfShape& aLS2B=myTransform2->Modified(aS1B);
368         const TopoDS_Shape& aS2B=aLS2B.First();
369         pHistShapes->Append(aS2B);
370       }
371     }
372   }
373   else {
374     const TopTools_ListOfShape& aLS1B=myBuilder->Modified(aS);
375     aItLS.Initialize(aLS1B);
376     for (; aItLS.More(); aItLS.Next()) {
377       const TopoDS_Shape& aS1B=aItLS.Value();
378       pHistShapes->Append(aS1B);
379     }
380   }
381   //
382   return myHistShapes;
383 }
384 //=======================================================================
385 //function : Generated
386 //purpose  : 
387 //=======================================================================
388 const TopTools_ListOfShape& HYDROData_Transform::Generated
389   (const TopoDS_Shape& aS)
390 {
391   TopTools_ListOfShape* pHistShapes;
392   TopTools_ListIteratorOfListOfShape aItLS;
393   //
394   pHistShapes=(TopTools_ListOfShape*)&myHistShapes;
395   pHistShapes->Clear();
396   //
397   if (myToTransform) { //pkv ft
398     const TopTools_ListOfShape& aLS1=myTransform1->Modified(aS);
399     const TopoDS_Shape& aS1=aLS1.First();
400     //
401     const TopTools_ListOfShape& aLS1B=myBuilder->Generated(aS1);
402     if (aLS1B.IsEmpty()) { //pkv ft
403       const TopTools_ListOfShape& aLS2B=myTransform2->Modified(aS1);
404       const TopoDS_Shape& aS2B=aLS2B.First();
405       pHistShapes->Append(aS2B);
406     }
407     else {
408       aItLS.Initialize(aLS1B);
409       for (; aItLS.More(); aItLS.Next()) {
410         const TopoDS_Shape& aS1B=aItLS.Value();
411         const TopTools_ListOfShape& aLS2B=myTransform2->Modified(aS1B);
412         const TopoDS_Shape& aS2B=aLS2B.First();
413         pHistShapes->Append(aS2B);
414       }
415     }
416   }
417   else {
418     const TopTools_ListOfShape& aLS1B=myBuilder->Generated(aS);
419     aItLS.Initialize(aLS1B);
420     for (; aItLS.More(); aItLS.Next()) {
421       const TopoDS_Shape& aS1B=aItLS.Value();
422       pHistShapes->Append(aS1B);
423     }
424   }
425   //
426   return myHistShapes;
427 }
428 //=======================================================================
429 //function : IsDeleted
430 //purpose  : 
431 //=======================================================================
432 Standard_Boolean HYDROData_Transform::IsDeleted(const TopoDS_Shape& aS)
433 {
434   Standard_Boolean bIsDeleted;
435   //
436   if (myToTransform) { //pkv ft
437     const TopTools_ListOfShape& aLS1=myTransform1->Modified(aS);
438     const TopoDS_Shape& aS1=aLS1.First();
439     //
440     bIsDeleted=myBuilder->IsDeleted(aS1);
441   }
442   else {
443     bIsDeleted=myBuilder->IsDeleted(aS);
444   }
445   //
446   return bIsDeleted;
447 }
448 //=======================================================================
449 //function : HasDeleted
450 //purpose  : 
451 //=======================================================================
452 Standard_Boolean HYDROData_Transform::HasDeleted()
453 {
454   return myBuilder->HasDeleted();
455 }
456 //=======================================================================
457 //function : HasGenerated
458 //purpose  : 
459 //=======================================================================
460 Standard_Boolean HYDROData_Transform::HasGenerated()
461 {
462   return myBuilder->HasGenerated();
463 }
464 //=======================================================================
465 //function : HasModified
466 //purpose  : 
467 //=======================================================================
468 Standard_Boolean HYDROData_Transform::HasModified()
469 {
470   return myBuilder->HasModified();
471 }