]> SALOME platform Git repositories - modules/geom.git/blobdiff - src/GEOMAlgo/GEOMAlgo_ShellSolid.cxx
Salome HOME
0021672: [CEA 565] Dump Study from script
[modules/geom.git] / src / GEOMAlgo / GEOMAlgo_ShellSolid.cxx
index e7da2ac055fea44cca07d7dc80c5aa57f96380e4..c032c1111cd588a8894e2c711db88eb06dd1b37c 100644 (file)
@@ -1,10 +1,31 @@
-// File:       GEOMAlgo_ShellSolid.cxx
-// Created:    Wed Jan 12 12:49:45 2005
-// Author:     Peter KURNEV
-//             <pkv@irinox>
+// Copyright (C) 2007-2012  CEA/DEN, EDF R&D, OPEN CASCADE
+//
+// Copyright (C) 2003-2007  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
+// CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
+//
+// This library is free software; you can redistribute it and/or
+// modify it under the terms of the GNU Lesser General Public
+// License as published by the Free Software Foundation; either
+// version 2.1 of the License.
+//
+// This library is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+// Lesser General Public License for more details.
+//
+// You should have received a copy of the GNU Lesser General Public
+// License along with this library; if not, write to the Free Software
+// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
+//
+// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
+//
 
-
-#include <GEOMAlgo_ShellSolid.ixx>
+// File:        GEOMAlgo_ShellSolid.cxx
+// Created:     Wed Jan 12 12:49:45 2005
+// Author:      Peter KURNEV
+//              <pkv@irinox>
+//
+#include <GEOMAlgo_ShellSolid.hxx>
 
 #include <Standard_Failure.hxx>
 
@@ -54,7 +75,7 @@
 
 //=======================================================================
 //function : GEOMAlgo_ShellSolid
-//purpose  : 
+//purpose  :
 //=======================================================================
 GEOMAlgo_ShellSolid::GEOMAlgo_ShellSolid()
 :
@@ -63,16 +84,16 @@ GEOMAlgo_ShellSolid::GEOMAlgo_ShellSolid()
 }
 //=======================================================================
 //function : ~
-//purpose  : 
+//purpose  :
 //=======================================================================
 GEOMAlgo_ShellSolid::~GEOMAlgo_ShellSolid()
 {
 }
 //=======================================================================
-// function: 
-// purpose: 
+// function:
+// purpose:
 //=======================================================================
-void GEOMAlgo_ShellSolid::Perform() 
+void GEOMAlgo_ShellSolid::Perform()
 {
   myErrorStatus=0;
   //
@@ -103,12 +124,12 @@ void GEOMAlgo_ShellSolid::Perform()
 }
 //=======================================================================
 // function: Prepare
-// purpose: 
+// purpose:
 //=======================================================================
-void GEOMAlgo_ShellSolid::Prepare() 
+void GEOMAlgo_ShellSolid::Prepare()
 {
   const BOPTools_PaveFiller& aPaveFiller=myDSFiller->PaveFiller();
-  // 
+  //
   // 1 States
   BOPTools_SolidStateFiller aStateFiller(aPaveFiller);
   aStateFiller.Do();
@@ -126,9 +147,9 @@ void GEOMAlgo_ShellSolid::Prepare()
 }
 //=================================================================================
 // function: BuildResult
-// purpose: 
+// purpose:
 //=================================================================================
-void GEOMAlgo_ShellSolid::BuildResult() 
+void GEOMAlgo_ShellSolid::BuildResult()
 {
   Standard_Boolean bIsTouchCase;
   Standard_Integer i, j, nF1, nF2, aNbFFs, aNbS, aNbCurves, nSp, iRank1;
@@ -197,7 +218,7 @@ void GEOMAlgo_ShellSolid::BuildResult()
       const BOPTools_ListOfPaveBlock& aSectEdges=aBC.NewPaveBlocks();
       aNbS=aSectEdges.Extent();
       if (aNbS) {
-       break;
+        break;
       }
     }
     //
@@ -209,7 +230,7 @@ void GEOMAlgo_ShellSolid::BuildResult()
     for (; anExp.More(); anExp.Next()) {
       const TopoDS_Edge& aE=TopoDS::Edge(anExp.Current());
       if (BRep_Tool::Degenerated(aE)) {
-       continue;
+        continue;
       }
       //
       nE=aDS.ShapeIndex(aE, myRank);
@@ -217,53 +238,53 @@ void GEOMAlgo_ShellSolid::BuildResult()
       aNbPB=aLPB.Extent();
       //
       if (aNbPB<2) {
-       nSp=nE;
-       if (aNbPB) {
-         const BOPTools_PaveBlock& aPB=aLPB.First();
-         nSp=aPB.Edge();
-       }
-       const TopoDS_Shape& aSp=aDS.Shape(nSp);
-       //
-       aState=aDS.GetState(nSp);
-       if (aState==BooleanOperations_IN) {
-         myLSIN.Append(aF1);
-       }
-       else if (aState==BooleanOperations_OUT) {
-         myLSOUT.Append(aF1);
-       }
-       else if (aState==BooleanOperations_ON) {
-         Standard_Real aTol;
-         TopAbs_State aSt;
-         //
-         //const TopoDS_Face& aF2=TopoDS::Face(aDS.Shape((iRank1==myRank)? nF2 : nF1));
-         //aTol=BRep_Tool::Tolerance(aF2);
-         aTol=1.e-7;
-         //
-         BOPTools_Tools3D::PointNearEdge(aE, aF1, aP2D, aP3D);
-         const TopoDS_Solid& aRefSolid=(myRank==1) ? 
-           TopoDS::Solid(aDS.Tool()) : TopoDS::Solid(aDS.Object());
-         //
-         BOPTools_PaveFiller* pPF=(BOPTools_PaveFiller*)& aPaveFiller;
-         IntTools_Context& aCtx=pPF->ChangeContext();
-         //
-         BRepClass3d_SolidClassifier& aSC=aCtx.SolidClassifier(aRefSolid);
-         aSC.Perform(aP3D, aTol);
-         aSt=aSC.State();
-         if (aSt==TopAbs_IN) {
-           myLSIN.Append(aF1);
-         }
-         else if (aSt==TopAbs_OUT) {
-           myLSOUT.Append(aF1);
-         }
-       } 
-       break; 
-      } // if (aNbPB<2) { 
+        nSp=nE;
+        if (aNbPB) {
+          const BOPTools_PaveBlock& aPB=aLPB.First();
+          nSp=aPB.Edge();
+        }
+        /*const TopoDS_Shape& aSp=*/aDS.Shape(nSp);
+        //
+        aState=aDS.GetState(nSp);
+        if (aState==BooleanOperations_IN) {
+          myLSIN.Append(aF1);
+        }
+        else if (aState==BooleanOperations_OUT) {
+          myLSOUT.Append(aF1);
+        }
+        else if (aState==BooleanOperations_ON) {
+          Standard_Real aTol;
+          TopAbs_State aSt;
+          //
+          //const TopoDS_Face& aF2=TopoDS::Face(aDS.Shape((iRank1==myRank)? nF2 : nF1));
+          //aTol=BRep_Tool::Tolerance(aF2);
+          aTol=1.e-7;
+          //
+          BOPTools_Tools3D::PointNearEdge(aE, aF1, aP2D, aP3D);
+          const TopoDS_Solid& aRefSolid=(myRank==1) ?
+            TopoDS::Solid(aDS.Tool()) : TopoDS::Solid(aDS.Object());
+          //
+          BOPTools_PaveFiller* pPF=(BOPTools_PaveFiller*)& aPaveFiller;
+          const Handle(IntTools_Context)& aCtx=pPF->Context();
+          //
+          BRepClass3d_SolidClassifier& aSC=aCtx->SolidClassifier(aRefSolid);
+          aSC.Perform(aP3D, aTol);
+          aSt=aSC.State();
+          if (aSt==TopAbs_IN) {
+            myLSIN.Append(aF1);
+          }
+          else if (aSt==TopAbs_OUT) {
+            myLSOUT.Append(aF1);
+          }
+        }
+        break;
+      } // if (aNbPB<2) {
     } //for (; anExp.More(); anExp.Next())
-  } 
+  }
 }
 //=======================================================================
 // function: DetectSDFaces
-// purpose: 
+// purpose:
 //=======================================================================
 void GEOMAlgo_ShellSolid::DetectSDFaces()
 {
@@ -278,9 +299,9 @@ void GEOMAlgo_ShellSolid::DetectSDFaces()
   aNb=aFFs.Extent();
   for (i=1; i<=aNb; i++) {
     bFlag=Standard_False;
-    
+
     BOPTools_SSInterference& aFF=aFFs(i);
-    
+
     nF1=aFF.Index1();
     nF2=aFF.Index2();
     const TopoDS_Face& aF1=TopoDS::Face(aDS.Shape(nF1));
@@ -293,11 +314,11 @@ void GEOMAlgo_ShellSolid::DetectSDFaces()
     if (!aNbSps) {
       continue;
     }
-    
+
     const BOPTools_PaveBlock& aPB=aLPB.First();
     const TopoDS_Edge& aSpE=TopoDS::Edge(aDS.Shape(aPB.Edge()));
-    
-    BOPTools_Tools3D::GetNormalToFaceOnEdge (aSpE, aF1, aDNF1); 
+
+    BOPTools_Tools3D::GetNormalToFaceOnEdge (aSpE, aF1, aDNF1);
     BOPTools_Tools3D::GetNormalToFaceOnEdge (aSpE, aF2, aDNF2);
     iSenseFlag=BOPTools_Tools3D::SenseFlag (aDNF1, aDNF2);
     //
@@ -306,13 +327,13 @@ void GEOMAlgo_ShellSolid::DetectSDFaces()
     //
       TopoDS_Face aF1FWD=aF1;
       aF1FWD.Orientation (TopAbs_FORWARD);
-      
+
       BOP_WireEdgeSet aWES (aF1FWD);
       BOP_SDFWESFiller aWESFiller(nF1, nF2, *myDSFiller);
       aWESFiller.SetSenseFlag(iSenseFlag);
       aWESFiller.SetOperation(BOP_COMMON);
       aWESFiller.Do(aWES);
-      
+
       BOP_FaceBuilder aFB;
       aFB.Do(aWES);
       const TopTools_ListOfShape& aLF=aFB.NewFaces();
@@ -320,31 +341,31 @@ void GEOMAlgo_ShellSolid::DetectSDFaces()
       iZone=0;
       TopTools_ListIteratorOfListOfShape anIt(aLF);
       for (; anIt.More(); anIt.Next()) {
-       const TopoDS_Shape& aFR=anIt.Value();
+        const TopoDS_Shape& aFR=anIt.Value();
 
-       if (aFR.ShapeType()==TopAbs_FACE) {
-         const TopoDS_Face& aFaceResult=TopoDS::Face(aFR);
-         //
-         Standard_Boolean bIsValidIn2D, bNegativeFlag;
-         bIsValidIn2D=BOPTools_Tools3D::IsValidArea (aFaceResult, bNegativeFlag);
-         if (bIsValidIn2D) { 
-           //if(CheckSameDomainFaceInside(aFaceResult, aF2)) {
-           iZone=1;
-           break;
-           //}
-         }
-         //
-       }
+        if (aFR.ShapeType()==TopAbs_FACE) {
+          const TopoDS_Face& aFaceResult=TopoDS::Face(aFR);
+          //
+          Standard_Boolean bIsValidIn2D, bNegativeFlag;
+          bIsValidIn2D=BOPTools_Tools3D::IsValidArea (aFaceResult, bNegativeFlag);
+          if (bIsValidIn2D) {
+            //if(CheckSameDomainFaceInside(aFaceResult, aF2)) {
+            iZone=1;
+            break;
+            //}
+          }
+          //
+        }
       }
-      
-      if (iZone) { 
-       bFlag=Standard_True;
-       aFF.SetStatesMap(aWESFiller.StatesMap());
+
+      if (iZone) {
+        bFlag=Standard_True;
+        aFF.SetStatesMap(aWESFiller.StatesMap());
       }
-      
+
     }// if (iSenseFlag)
-  
+
   aFF.SetTangentFacesFlag(bFlag);
   aFF.SetSenseFlag (iSenseFlag);
-  }// end of for (i=1; i<=aNb; i++) 
+  }// end of for (i=1; i<=aNb; i++)
 }