1 // Created by: Peter KURNEV
2 // Copyright (c) 2010-2014 OPEN CASCADE SAS
3 // Copyright (c) 2007-2010 CEA/DEN, EDF R&D, OPEN CASCADE
4 // Copyright (c) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN, CEDRAT,
5 // EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
7 // This file is part of Open CASCADE Technology software library.
9 // This library is free software; you can redistribute it and/or modify it under
10 // the terms of the GNU Lesser General Public License version 2.1 as published
11 // by the Free Software Foundation, with special exception defined in the file
12 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
13 // distribution for complete text of the license and disclaimer of any warranty.
15 // Alternatively, this file may be used under the terms of Open CASCADE
16 // commercial license or contractual agreement.
18 #include <BOPAlgo_Builder.ixx>
20 #include <NCollection_IncAllocator.hxx>
22 #include <TopoDS_Shape.hxx>
23 #include <TopoDS_Face.hxx>
24 #include <TopoDS_Edge.hxx>
25 #include <TopoDS_Vertex.hxx>
26 #include <TopoDS_Compound.hxx>
28 #include <BRep_Tool.hxx>
29 #include <BRep_Builder.hxx>
31 #include <TopExp_Explorer.hxx>
33 #include <BOPCol_ListOfShape.hxx>
34 #include <BOPCol_ListOfInteger.hxx>
35 #include <BOPCol_MapOfInteger.hxx>
36 #include <BOPCol_DataMapOfIntegerListOfShape.hxx>
37 #include <BOPCol_DataMapOfShapeShape.hxx>
39 #include <BOPInt_Context.hxx>
41 #include <BOPDS_PaveBlock.hxx>
42 #include <BOPDS_ShapeInfo.hxx>
43 #include <BOPDS_DS.hxx>
44 #include <BOPDS_FaceInfo.hxx>
45 #include <BOPDS_MapOfPaveBlock.hxx>
46 #include <BOPDS_VectorOfInterfFF.hxx>
47 #include <BOPDS_Interf.hxx>
48 #include <BOPDS_VectorOfCurve.hxx>
49 #include <BOPDS_VectorOfPoint.hxx>
51 #include <BOPTools.hxx>
52 #include <BOPTools_AlgoTools.hxx>
53 #include <BOPTools_AlgoTools2D.hxx>
54 #include <BOPTools_AlgoTools3D.hxx>
55 #include <BOPAlgo_BuilderFace.hxx>
56 #include <BOPTools_CoupleOfShape.hxx>
57 #include <BOPTools_ListOfCoupleOfShape.hxx>
58 #include <BOPTools_MapOfSet.hxx>
59 #include <BOPTools_DataMapOfShapeSet.hxx>
60 #include <BOPAlgo_Builder_2Cnt.hxx>
63 Standard_Boolean HasPaveBlocksOnIn(const BOPDS_FaceInfo& aFI1,
64 const BOPDS_FaceInfo& aFI2);
66 void FillMap(const TopoDS_Shape& aS1,
67 const TopoDS_Shape& aS2,
68 BOPCol_IndexedDataMapOfShapeListOfShape& aDMSLS,
69 Handle(NCollection_IncAllocator)& aAllocator);
71 void MakeBlocksCnx(const BOPCol_IndexedDataMapOfShapeListOfShape& aMILI,
72 BOPCol_DataMapOfIntegerListOfShape& aMBlocks,
73 Handle(NCollection_IncAllocator)& aAllocator);
76 //=======================================================================
77 //function : FillImagesFaces
79 //=======================================================================
80 void BOPAlgo_Builder::FillImagesFaces()
85 FillSameDomainFaces();
88 //=======================================================================
89 //function : BuildSplitFaces
91 //=======================================================================
92 void BOPAlgo_Builder::BuildSplitFaces()
94 Standard_Boolean bHasFaceInfo, bIsClosed, bIsDegenerated, bToReverse;
95 Standard_Integer i, j, k, aNbS, aNbPBIn, aNbPBOn, aNbPBSc, aNbAV, nSp;
97 TopoDS_Face aFF, aFSD;
99 TopAbs_Orientation anOriF, anOriE;
100 TopExp_Explorer aExp;
101 BOPCol_ListIteratorOfListOfShape aIt;
102 BOPCol_ListOfInteger aLIAV;
103 BOPCol_MapOfShape aMFence;
104 Handle(NCollection_BaseAllocator) aAllocator;
105 BOPCol_ListOfShape aLFIm(myAllocator);
106 BOPCol_MapIteratorOfMapOfShape aItMS;
107 BOPAlgo_VectorOfBuilderFace aVBF;
111 //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~scope f
112 aAllocator=new NCollection_IncAllocator();
114 BOPCol_ListOfShape aLE(aAllocator);
115 BOPCol_MapOfShape aMDE(100, aAllocator);
117 aNbS=myDS->NbSourceShapes();
119 for (i=0; i<aNbS; ++i) {
120 const BOPDS_ShapeInfo& aSI=myDS->ShapeInfo(i);
121 if (aSI.ShapeType()!=TopAbs_FACE) {
125 const TopoDS_Face& aF=(*(TopoDS_Face*)(&aSI.Shape()));
127 bHasFaceInfo=myDS->HasFaceInfo(i);
132 const BOPDS_FaceInfo& aFI=myDS->FaceInfo(i);
134 const BOPDS_IndexedMapOfPaveBlock& aMPBIn=aFI.PaveBlocksIn();
135 const BOPDS_IndexedMapOfPaveBlock& aMPBOn=aFI.PaveBlocksOn();
136 const BOPDS_IndexedMapOfPaveBlock& aMPBSc=aFI.PaveBlocksSc();
138 myDS->AloneVertices(i, aLIAV);
140 aNbPBIn=aMPBIn.Extent();
141 aNbPBOn=aMPBOn.Extent();
142 aNbPBSc=aMPBSc.Extent();
143 aNbAV=aLIAV.Extent();
144 if (!aNbPBIn && !aNbPBOn && !aNbPBSc && !aNbAV) { // not compete
150 anOriF=aF.Orientation();
152 aFF.Orientation(TopAbs_FORWARD);
156 // 1. Fill the egdes set for the face aFF -> LE
160 // 1.1 Bounding edges
161 aExp.Init(aFF, TopAbs_EDGE);
162 for (; aExp.More(); aExp.Next()) {
163 const TopoDS_Edge& aE=(*(TopoDS_Edge*)(&aExp.Current()));
164 anOriE=aE.Orientation();
165 bIsDegenerated=BRep_Tool::Degenerated(aE);
166 bIsClosed=BRep_Tool::IsClosed(aE, aF);
168 if (!myImages.IsBound(aE)) {
169 if (anOriE==TopAbs_INTERNAL) {
171 aEE.Orientation(TopAbs_FORWARD);
173 aEE.Orientation(TopAbs_REVERSED);
181 const BOPCol_ListOfShape& aLIE=myImages.Find(aE);
182 aIt.Initialize(aLIE);
183 for (; aIt.More(); aIt.Next()) {
184 aSp=(*(TopoDS_Edge*)(&aIt.Value()));
185 if (bIsDegenerated) {
186 aSp.Orientation(anOriE);
191 if (anOriE==TopAbs_INTERNAL) {
192 aSp.Orientation(TopAbs_FORWARD);
194 aSp.Orientation(TopAbs_REVERSED);
200 if (aMFence.Add(aSp)) {
201 if (!BRep_Tool::IsClosed(aSp, aF)){
202 BOPTools_AlgoTools3D::DoSplitSEAMOnFace(aSp, aF);
205 aSp.Orientation(TopAbs_FORWARD);
207 aSp.Orientation(TopAbs_REVERSED);
209 }// if (aMFence.Add(aSp))
213 aSp.Orientation(anOriE);
214 bToReverse=BOPTools_AlgoTools::IsSplitToReverse(aSp, aE, myContext);
219 }// for (; aIt.More(); aIt.Next()) {
221 }// for (; aExp.More(); aExp.Next()) {
225 for (j=1; j<=aNbPBIn; ++j) {
226 const Handle(BOPDS_PaveBlock)& aPB=aMPBIn(j);
228 aSp=(*(TopoDS_Edge*)(&myDS->Shape(nSp)));
230 aSp.Orientation(TopAbs_FORWARD);
232 aSp.Orientation(TopAbs_REVERSED);
238 for (j=1; j<=aNbPBSc; ++j) {
239 const Handle(BOPDS_PaveBlock)& aPB=aMPBSc(j);
241 aSp=(*(TopoDS_Edge*)(&myDS->Shape(nSp)));
243 aSp.Orientation(TopAbs_FORWARD);
245 aSp.Orientation(TopAbs_REVERSED);
249 BOPTools_AlgoTools2D::BuildPCurveForEdgesOnPlane (aLE, aFF);
251 // 3 Build split faces
252 BOPAlgo_BuilderFace& aBF=aVBF.Append1();
255 aBF.SetRunParallel(myRunParallel);
257 }// for (i=0; i<aNbS; ++i) {
261 //===================================================
262 BOPAlgo_BuilderFaceCnt::Perform(myRunParallel, aVBF);
263 //===================================================
265 for (k=0; k<(Standard_Integer)aNbBF; ++k) {
268 BOPAlgo_BuilderFace& aBF=aVBF(k);
269 TopoDS_Face aF=aBF.Face();
270 anOriF=aBF.Orientation();
271 aF.Orientation(anOriF);
273 const BOPCol_ListOfShape& aLFR=aBF.Areas();
274 aIt.Initialize(aLFR);
275 for (; aIt.More(); aIt.Next()) {
276 TopoDS_Shape& aFR=aIt.ChangeValue();
277 if (anOriF==TopAbs_REVERSED) {
278 aFR.Orientation(TopAbs_REVERSED);
280 //aFR.Orientation(anOriF);
284 mySplits.Bind(aF, aLFIm);
285 }// for (k=0; k<aNbBF; ++k) {
287 aAllocator.Nullify();
288 //~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~scope t
290 //=======================================================================
291 //function : FillSameDomainFaces
293 //=======================================================================
294 void BOPAlgo_Builder::FillSameDomainFaces()
296 Standard_Boolean bFlag;
297 Standard_Integer i, j, k, aNbFFs, aNbCurves, aNbPoints, nF1, nF2, aNbS;
298 Handle(NCollection_IncAllocator) aAllocator;
299 BOPCol_ListIteratorOfListOfShape aItF;
300 BOPCol_MapOfShape aMFence;
301 BOPAlgo_IndexedDataMapOfSetInteger aIDMSS;
302 BOPAlgo_VectorOfVectorOfShape aVVS;
306 const BOPDS_VectorOfInterfFF& aFFs=myDS->InterfFF();
308 aNbFFs=aFFs.Extent();
313 for (i=0; i<aNbFFs; ++i) {
314 const BOPDS_InterfFF& aFF=aFFs(i);
315 aFF.Indices(nF1, nF2);
317 const BOPDS_VectorOfCurve& aCurves=aFF.Curves();
318 aNbCurves=aCurves.Extent();
321 bFlag=Standard_False;
322 for (j=0; j<aNbCurves; ++j) {
323 const BOPDS_Curve& aNC=aCurves.Value(j);
335 const BOPDS_VectorOfPoint& aPoints=aFF.Points();
336 aNbPoints=aPoints.Extent();
341 if (!myDS->HasFaceInfo(nF1) || !myDS->HasFaceInfo(nF2) ) {
345 const BOPDS_FaceInfo& aFI1=myDS->FaceInfo(nF1);
346 const BOPDS_FaceInfo& aFI2=myDS->FaceInfo(nF2);
348 const TopoDS_Shape& aF1=myDS->Shape(nF1);
349 const TopoDS_Shape& aF2=myDS->Shape(nF2);
351 bFlag=HasPaveBlocksOnIn(aFI1, aFI2);
352 bFlag=bFlag && (mySplits.IsBound(aF1) && mySplits.IsBound(aF2));
355 for (k=0; k<2; ++k) {
356 const TopoDS_Shape& aF=(!k) ? aF1 : aF2;
357 const BOPCol_ListOfShape& aLF=mySplits.Find(aF);
359 aItF.Initialize(aLF);
360 for (; aItF.More(); aItF.Next()) {
361 const TopoDS_Shape& aFx=aItF.Value();
363 if (aMFence.Add(aFx)) {
368 if (!aIDMSS.Contains(aSTx)) {
369 BOPAlgo_VectorOfShape& aVS=aVVS.Append1();
373 aIDMSS.Add (aSTx, j);
376 j=aIDMSS.ChangeFromKey(aSTx);
377 BOPAlgo_VectorOfShape& aVS=aVVS(j);
385 BOPTools_Set aST1, aST2;
390 if (aST1.IsEqual(aST2)) {
391 if (!aIDMSS.Contains(aST1)) {
392 BOPAlgo_VectorOfShape& aVS=aVVS.Append1();
393 if (aMFence.Add(aF1)) {
396 if (aMFence.Add(aF2)) {
401 aIDMSS.Add (aST1, k);
404 k=aIDMSS.ChangeFromKey(aST1);
405 BOPAlgo_VectorOfShape& aVS=aVVS(k);
406 if (aMFence.Add(aF1)) {
409 if (aMFence.Add(aF2)) {
413 }//if (aST1.IsEqual(aST2)) {
414 }// else {// if (!bFlag)
416 }// for (i=0; i<aNbFFs; ++i) {
420 Standard_Boolean bFlagSD;
421 Standard_Integer aNbVPSB, aNbVVS, aNbF, aNbF1;
422 BOPAlgo_VectorOfPairOfShapeBoolean aVPSB;
424 aNbVVS=aVVS.Extent();
425 for (i=0; i<aNbVVS; ++i) {
426 const BOPAlgo_VectorOfShape& aVS=aVVS(i);
433 for (j=0; j<aNbF1; ++j) {
434 const TopoDS_Shape& aFj=aVS(j);
435 for (k=j+1; k<aNbF; ++k) {
436 const TopoDS_Shape& aFk=aVS(k);
437 BOPAlgo_PairOfShapeBoolean& aPSB=aVPSB.Append1();
443 //====================================================
444 BOPAlgo_BuilderSDFaceCnt::Perform(myRunParallel, aVPSB);
445 //====================================================
446 aAllocator=new NCollection_IncAllocator();
447 BOPCol_IndexedDataMapOfShapeListOfShape aDMSLS(100, aAllocator);
448 BOPCol_DataMapOfIntegerListOfShape aMBlocks(100, aAllocator);
450 aNbVPSB=aVPSB.Extent();
451 for (i=0; i<aNbVPSB; ++i) {
452 BOPAlgo_PairOfShapeBoolean& aPSB=aVPSB(i);
455 const TopoDS_Shape& aFj=aPSB.Shape1();
456 const TopoDS_Shape& aFk=aPSB.Shape2();
457 FillMap(aFj, aFk, aDMSLS, aAllocator);
463 MakeBlocksCnx(aDMSLS, aMBlocks, aAllocator);
465 // 3. Fill same domain faces map -> aMSDF
466 aNbS = aMBlocks.Extent();
467 for (i=0; i<aNbS; ++i) {
468 const BOPCol_ListOfShape& aLSD=aMBlocks.Find(i);
469 if (aLSD.IsEmpty()) {
473 const TopoDS_Shape& aFSD1=aLSD.First();
474 aItF.Initialize(aLSD);
475 for (; aItF.More(); aItF.Next()) {
476 const TopoDS_Shape& aFSD=aItF.Value();
477 myShapesSD.Bind(aFSD, aFSD1);
479 // If the face has no splits but are SD face,
480 // it is considered as splitted face
481 if (!mySplits.IsBound(aFSD)) {
482 BOPCol_ListOfShape aLS;
484 mySplits.Bind(aFSD, aLS);
490 aAllocator.Nullify();
492 //=======================================================================
493 // function: FillImagesFaces1
495 //=======================================================================
496 void BOPAlgo_Builder::FillImagesFaces1()
498 Standard_Integer i, aNbS, iSense;
500 BOPCol_ListOfInteger aLIAV;
501 BOPCol_ListOfShape aLFIm;
502 BOPCol_ListIteratorOfListOfShape aItLS;
504 aNbS=myDS->NbSourceShapes();
505 for (i=0; i<aNbS; ++i) {
506 const BOPDS_ShapeInfo& aSI=myDS->ShapeInfo(i);
507 if (aSI.ShapeType()!=TopAbs_FACE) {
511 const TopoDS_Face& aF=(*(TopoDS_Face*)(&aSI.Shape()));
513 if (!mySplits.IsBound(aF)) {
518 myDS->AloneVertices(i, aLIAV);
521 const BOPCol_ListOfShape& aLSp=mySplits.Find(aF);
522 aItLS.Initialize(aLSp);
523 for (; aItLS.More(); aItLS.Next()) {
524 const TopoDS_Face& aFSp=(*(TopoDS_Face*)(&aItLS.Value()));
525 if (!myShapesSD.IsBound(aFSp)) {
529 aFSD=(*(TopoDS_Face*)(&myShapesSD.Find(aFSp)));
530 iSense=BOPTools_AlgoTools::Sense(aFSp, aFSD);
538 FillInternalVertices(aLFIm, aLIAV);
540 myImages.Bind(aF, aLFIm);
543 aItLS.Initialize(aLFIm);
544 for (; aItLS.More(); aItLS.Next()) {
545 const TopoDS_Face& aFSp=(*(TopoDS_Face*)(&aItLS.Value()));
546 myOrigins.Bind(aFSp, aF);
548 }// for (i=0; i<aNbS; ++i) {
550 //=======================================================================
551 // function: FillInternalVertices
553 //=======================================================================
554 void BOPAlgo_Builder::FillInternalVertices(BOPCol_ListOfShape& aLFIm,
555 BOPCol_ListOfInteger& aLIAV)
557 Standard_Integer nV, iFlag;
558 Standard_Real aU1, aU2;
561 BOPCol_ListIteratorOfListOfInteger aItV;
562 BOPCol_ListIteratorOfListOfShape aItF;
564 aItV.Initialize(aLIAV);
565 for (; aItV.More(); aItV.Next()) {
567 aV=(*(TopoDS_Vertex*)(&myDS->Shape(nV)));
568 aV.Orientation(TopAbs_INTERNAL);
570 aItF.Initialize(aLFIm);
571 for (; aItF.More(); aItF.Next()) {
572 TopoDS_Face& aF=(*(TopoDS_Face*)(&aItF.Value()));
573 iFlag=myContext->ComputeVF(aV, aF, aU1, aU2);
581 //=======================================================================
582 //function : MakeBlocksCnx
584 //=======================================================================
585 void MakeBlocksCnx(const BOPCol_IndexedDataMapOfShapeListOfShape& aMILI,
586 BOPCol_DataMapOfIntegerListOfShape& aMBlocks,
587 Handle(NCollection_IncAllocator)& aAllocator)
589 Standard_Integer aNbV, aNbVS, aNbVP, aNbEC, k, i, j;
590 BOPCol_ListIteratorOfListOfShape aItLI;
592 BOPCol_MapOfShape aMVS(100, aAllocator);
593 BOPCol_IndexedMapOfShape aMEC(100, aAllocator);
594 BOPCol_IndexedMapOfShape aMVP(100, aAllocator);
595 BOPCol_IndexedMapOfShape aMVAdd(100, aAllocator);
599 for (k=0,i=1; i<=aNbV; ++i) {
605 const TopoDS_Shape& nV=aMILI.FindKey(i);
606 if (aMVS.Contains(nV)){
618 for (j=1; j<=aNbVP; ++j) {
619 const TopoDS_Shape& nVP=aMVP(j);
620 const BOPCol_ListOfShape& aLV=aMILI.FindFromKey(nVP);
621 aItLI.Initialize(aLV);
622 for (; aItLI.More(); aItLI.Next()) {
623 const TopoDS_Shape& nVx=aItLI.Value();
624 if (aMEC.Contains(nVx)) {
634 aNbVP=aMVAdd.Extent();
636 break; // from while(1)
640 for (j=1; j<=aNbVP; ++j) {
646 BOPCol_ListOfShape aLIx(aAllocator);
648 aNbEC = aMEC.Extent();
649 for (j=1; j<=aNbEC; ++j) {
650 const TopoDS_Shape& nVx=aMEC(j);
654 aMBlocks.Bind(k, aLIx);
656 }//for (k=0,i=1; i<=aNbV; ++i)
663 //=======================================================================
666 //=======================================================================
667 void FillMap(const TopoDS_Shape& aS1,
668 const TopoDS_Shape& aS2,
669 BOPCol_IndexedDataMapOfShapeListOfShape& aDMSLS,
670 Handle(NCollection_IncAllocator)& aAllocator)
672 if (aDMSLS.Contains(aS1)) {
673 BOPCol_ListOfShape& aLS=aDMSLS.ChangeFromKey(aS1);
677 BOPCol_ListOfShape aLS(aAllocator);
679 aDMSLS.Add(aS1, aLS);
682 if (aDMSLS.Contains(aS2)) {
683 BOPCol_ListOfShape& aLS=aDMSLS.ChangeFromKey(aS2);
687 BOPCol_ListOfShape aLS(aAllocator);
689 aDMSLS.Add(aS2, aLS);
692 //=======================================================================
693 //function :HasPaveBlocksOnIn
695 //=======================================================================
696 Standard_Boolean HasPaveBlocksOnIn(const BOPDS_FaceInfo& aFI1,
697 const BOPDS_FaceInfo& aFI2)
699 Standard_Boolean bRet;
700 BOPDS_MapIteratorOfMapOfPaveBlock aItMPB;
703 const BOPDS_IndexedMapOfPaveBlock& aMPBOn1=aFI1.PaveBlocksOn();
704 const BOPDS_IndexedMapOfPaveBlock& aMPBIn1=aFI1.PaveBlocksIn();
706 const BOPDS_IndexedMapOfPaveBlock& aMPBOn2=aFI2.PaveBlocksOn();
707 aItMPB.Initialize(aMPBOn2);
708 for (; aItMPB.More(); aItMPB.Next()) {
709 const Handle(BOPDS_PaveBlock)& aPB=aItMPB.Value();
710 bRet=aMPBOn1.Contains(aPB) || aMPBIn1.Contains(aPB);
716 const BOPDS_IndexedMapOfPaveBlock& aMPBIn2=aFI2.PaveBlocksIn();
717 aItMPB.Initialize(aMPBIn2);
718 for (; aItMPB.More(); aItMPB.Next()) {
719 const Handle(BOPDS_PaveBlock)& aPB=aItMPB.Value();
720 bRet=aMPBOn1.Contains(aPB) || aMPBIn1.Contains(aPB);
732 BOPCol_ListIteratorOfListOfShape aItx;
734 aBBx.MakeCompound(aCx);
736 aItx.Initialize(aLE);
737 for (; aItx.More(); aItx.Next()) {
738 const TopoDS_Shape& aEx=aItx.Value();