Salome HOME
Porting to OCCT development version: Standard_PI -> M_PI
[modules/geom.git] / src / GEOMAlgo / GEOMAlgo_FinderShapeOn1.cxx
index 6f9f6a0b45277d8fb301c663d30d7db19a375557..fe0a48d5d15b1822757e8fe61ee95ce5ef550989 100644 (file)
@@ -1,28 +1,30 @@
-// Copyright (C) 2005  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
+// Copyright (C) 2007-2011  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 
+// 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 
+//
+// 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 
+// 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
 //
-// File:       GEOMAlgo_FinderShapeOn1.cxx
-// Created:    Fri Mar  4 10:31:06 2005
-// Author:     Peter KURNEV
-//             <pkv@irinox>
-
 
+// File:        GEOMAlgo_FinderShapeOn1.cxx
+// Created:     Fri Mar  4 10:31:06 2005
+// Author:      Peter KURNEV
+//              <pkv@irinox>
+//
 #include <GEOMAlgo_FinderShapeOn1.ixx>
 #include <math.h>
 
@@ -279,7 +281,7 @@ void GEOMAlgo_FinderShapeOn1::ProcessVertices()
     //
     if (myShapeType==TopAbs_VERTEX){
       if (bIsConformState) {
-       myMSS.Add(aV, aSt);
+        myMSS.Add(aV, aSt);
       }
     }
     else if (bIsConformState || aSt==TopAbs_ON) {
@@ -319,11 +321,11 @@ void GEOMAlgo_FinderShapeOn1::ProcessEdges()
       //
       bIsConformState=myMSS.Contains(aV);
       if (!bIsConformState) {
-       break;// vertex has non-conformed state
+        break;// vertex has non-conformed state
       }
       else {
-       aSt=myMSS.FindFromKey(aV);
-       aSC.AppendState(aSt);
+        aSt=myMSS.FindFromKey(aV);
+        aSC.AppendState(aSt);
       }
     }
     //
@@ -347,7 +349,7 @@ void GEOMAlgo_FinderShapeOn1::ProcessEdges()
       //
       aType2=aGAC.GetType();
       if (aType2==GeomAbs_Line) {
-       continue;
+        continue;
       }
     }
     //
@@ -360,16 +362,16 @@ void GEOMAlgo_FinderShapeOn1::ProcessEdges()
     aIt.Initialize(aLP);
     for (iCnt=0; aIt.More(); aIt.Next(), ++iCnt) {
       if (myNbPntsMax) {
-       if (iCnt > myNbPntsMax) {
-         break;
-       }
+        if (iCnt > myNbPntsMax) {
+          break;
+        }
       }
       //
       const gp_Pnt& aP=aIt.Value();
       aSt = GetPointState( aP );
       bIsToBreak=aSC.AppendState(aSt);
       if (bIsToBreak) {
-       break;
+        break;
       }
     }
     //
@@ -378,7 +380,7 @@ void GEOMAlgo_FinderShapeOn1::ProcessEdges()
     bIsConformState=GEOMAlgo_SurfaceTools::IsConformState(aSt, myState);
     if (myShapeType==TopAbs_EDGE) {
       if (bIsConformState) {
-       myMSS.Add(aE, aSt);
+        myMSS.Add(aE, aSt);
       }
     }
     else if (bIsConformState || aSt==TopAbs_ON) {
@@ -420,7 +422,7 @@ void GEOMAlgo_FinderShapeOn1::ProcessFaces()
       aGAS.Load(aS);
       aType2=aGAS.GetType();
       if (aType2!=aType1) {
-       continue;
+        continue;
       }
     }
     //
@@ -429,11 +431,11 @@ void GEOMAlgo_FinderShapeOn1::ProcessFaces()
       const TopoDS_Shape& aE=aExp.Current();
       bIsConformState=myMSS.Contains(aE);
       if (!bIsConformState) {
-       break;// edge has non-conformed state
+        break;// edge has non-conformed state
       }
       else {
-       aSt=myMSS.FindFromKey(aE);
-       aSC.AppendState(aSt);
+        aSt=myMSS.FindFromKey(aE);
+        aSC.AppendState(aSt);
       }
     }
     //
@@ -450,16 +452,16 @@ void GEOMAlgo_FinderShapeOn1::ProcessFaces()
     aIt.Initialize(aLP);
     for (iCnt=0; aIt.More(); aIt.Next(), ++iCnt) {
       if (myNbPntsMax) {
-       if (iCnt > myNbPntsMax) {
-         break;
-       }
+        if (iCnt > myNbPntsMax) {
+          break;
+        }
       }
       //
       const gp_Pnt& aP=aIt.Value();
       aSt = GetPointState( aP );
       bIsToBreak=aSC.AppendState(aSt);
       if (bIsToBreak) {
-       break;
+        break;
       }
     }
     //
@@ -468,7 +470,7 @@ void GEOMAlgo_FinderShapeOn1::ProcessFaces()
     bIsConformState=GEOMAlgo_SurfaceTools::IsConformState(aSt, myState);
     if (myShapeType==TopAbs_FACE) {
       if (bIsConformState) {
-       myMSS.Add(aF, aSt);
+        myMSS.Add(aF, aSt);
       }
     }
     else if (bIsConformState || aSt==TopAbs_ON) {
@@ -502,11 +504,11 @@ void GEOMAlgo_FinderShapeOn1::ProcessSolids()
       const TopoDS_Shape& aF=aMF(j);
       bIsConformState=myMSS.Contains(aF);
       if (!bIsConformState) {
-       break;// face has non-conformed state
+        break;// face has non-conformed state
       }
       else {
-       aSt=myMSS.FindFromKey(aF);
-       aSC.AppendState(aSt);
+        aSt=myMSS.FindFromKey(aF);
+        aSC.AppendState(aSt);
       }
     }
     //
@@ -528,7 +530,7 @@ void GEOMAlgo_FinderShapeOn1::ProcessSolids()
 //purpose  : 
 //=======================================================================
 void GEOMAlgo_FinderShapeOn1::InnerPoints(const TopoDS_Face& aF,
-                                         GEOMAlgo_ListOfPnt& aLP)
+                                          GEOMAlgo_ListOfPnt& aLP)
 {
   myErrorStatus=0;
   //
@@ -545,7 +547,7 @@ void GEOMAlgo_FinderShapeOn1::InnerPoints(const TopoDS_Face& aF,
   aTRF=BRep_Tool::Triangulation(aF, aLoc);
   if (aTRF.IsNull()) {
     if (!GEOMAlgo_FinderShapeOn::BuildTriangulation(aF)) {
-      myErrorStatus=20; // no triangulation found
+      myWarningStatus=20; // no triangulation found
       return;
     }
     aTRF=BRep_Tool::Triangulation(aF, aLoc);
@@ -567,11 +569,11 @@ void GEOMAlgo_FinderShapeOn1::InnerPoints(const TopoDS_Face& aF,
       //
       aPK.SetIds(n[k], n[k+1]);
       if (aMPKI.IsBound(aPK)) {
-       Standard_Integer& iCntX=aMPKI.ChangeFind(aPK);
-       ++iCntX;
+        Standard_Integer& iCntX=aMPKI.ChangeFind(aPK);
+        ++iCntX;
       }
       else {
-       aMPKI.Bind(aPK, 1);
+        aMPKI.Bind(aPK, 1);
       }
     }
   }
@@ -588,8 +590,8 @@ void GEOMAlgo_FinderShapeOn1::InnerPoints(const TopoDS_Face& aF,
       aNbMax=aPK.NbMax();
       pIds=(Standard_Integer*)aPK.Key();
       for (k=1; k<3; ++k) {
-       aNx=*(pIds+aNbMax-k);
-       aMBN.Add(aNx);
+        aNx=*(pIds+aNbMax-k);
+        aMBN.Add(aNx);
       }
       */
       aNx=(Standard_Integer)aPK.Id(1);
@@ -628,49 +630,49 @@ void GEOMAlgo_FinderShapeOn1::InnerPoints(const TopoDS_Face& aF,
       aNbLinks=aMPKI.Extent();
       aIt.Initialize(aMPKI);
       for (; aIt.More(); aIt.Next()) {
-       iCnt=aIt.Value();
-       if (iCnt>1) {
-         // take the first having occured inner link
-         // and discretize it
-         const GEOMAlgo_PassKey& aPK=aIt.Key();
-         //qf
-         /*
-         aNbMax=aPK.NbMax();
-         pIds=(Standard_Integer*)aPK.Key();
-         aN1=*(pIds+aNbMax-1);
-         aN2=*(pIds+aNbMax-2);
-         */
-         //
-         aN1=(Standard_Integer)aPK.Id(1);
-         aN2=(Standard_Integer)aPK.Id(2);
-         //qt
-         aP1=aNodes(aN1).Transformed(aTrsf);
-         aP2=aNodes(aN2).Transformed(aTrsf);
-         //
-         if (aType==GeomAbs_Cylinder) {
-           Standard_Real aTolSM;
-           gp_Cylinder aCyl;
-           //
-           aTolSM=1.523e-6;//~1.-cos(0.1 deg)
-           aCyl=aGAS.Cylinder();
-           if (!GEOMAlgo_SurfaceTools::IsCoaxial(aP1, aP2, aCyl, aTolSM)) {
-             continue;
-           }
-         }
-         //
-         BRepLib_MakeEdge aBME(aP1, aP2);
-         bIsDone=aBME.IsDone();
-         if (!bIsDone) {
-           myErrorStatus=30; //can not obtain the line fron the link
-           return;
-         }
-         //
-         const TopoDS_Shape& aSx=aBME.Shape();
-         const TopoDS_Edge& aE=TopoDS::Edge(aSx);
-         //
-         InnerPoints(aE, myNbPntsMin, aLP);
-         break;
-       }// if (iCnt>1)
+        iCnt=aIt.Value();
+        if (iCnt>1) {
+          // take the first having occured inner link
+          // and discretize it
+          const GEOMAlgo_PassKey& aPK=aIt.Key();
+          //qf
+          /*
+          aNbMax=aPK.NbMax();
+          pIds=(Standard_Integer*)aPK.Key();
+          aN1=*(pIds+aNbMax-1);
+          aN2=*(pIds+aNbMax-2);
+          */
+          //
+          aN1=(Standard_Integer)aPK.Id(1);
+          aN2=(Standard_Integer)aPK.Id(2);
+          //qt
+          aP1=aNodes(aN1).Transformed(aTrsf);
+          aP2=aNodes(aN2).Transformed(aTrsf);
+          //
+          if (aType==GeomAbs_Cylinder) {
+            Standard_Real aTolSM;
+            gp_Cylinder aCyl;
+            //
+            aTolSM=1.523e-6;//~1.-cos(0.1 deg)
+            aCyl=aGAS.Cylinder();
+            if (!GEOMAlgo_SurfaceTools::IsCoaxial(aP1, aP2, aCyl, aTolSM)) {
+              continue;
+            }
+          }
+          //
+          BRepLib_MakeEdge aBME(aP1, aP2);
+          bIsDone=aBME.IsDone();
+          if (!bIsDone) {
+            myErrorStatus=30; //can not obtain the line fron the link
+            return;
+          }
+          //
+          const TopoDS_Shape& aSx=aBME.Shape();
+          const TopoDS_Edge& aE=TopoDS::Edge(aSx);
+          //
+          InnerPoints(aE, myNbPntsMin, aLP);
+          break;
+        }// if (iCnt>1)
       }// for (; aIt.More(); aIt.Next())
     }// if (aType==GeomAbs_Plane || aType==GeomAbs_Cylinder) 
   }// if (!aNb && myNbPntsMin) { 
@@ -680,7 +682,7 @@ void GEOMAlgo_FinderShapeOn1::InnerPoints(const TopoDS_Face& aF,
 //purpose  : 
 //=======================================================================
 void GEOMAlgo_FinderShapeOn1::InnerPoints(const TopoDS_Edge& aE,
-                                         GEOMAlgo_ListOfPnt& aLP)
+                                          GEOMAlgo_ListOfPnt& aLP)
 {
   myErrorStatus=0;
   //
@@ -736,8 +738,8 @@ void GEOMAlgo_FinderShapeOn1::InnerPoints(const TopoDS_Edge& aE,
 //purpose  : 
 //=======================================================================
 void GEOMAlgo_FinderShapeOn1::InnerPoints(const TopoDS_Edge& aE,
-                                         const Standard_Integer aNbPntsMin,
-                                         GEOMAlgo_ListOfPnt& aLP)
+                                          const Standard_Integer aNbPntsMin,
+                                          GEOMAlgo_ListOfPnt& aLP)
 {
   // try to fill it yourself
   Standard_Boolean bInf1, bInf2;
@@ -785,9 +787,9 @@ void GEOMAlgo_FinderShapeOn1::CheckData()
   }
   //
   if (!(myShapeType==TopAbs_VERTEX ||
-       myShapeType==TopAbs_EDGE ||
-       myShapeType==TopAbs_FACE ||
-       myShapeType==TopAbs_SOLID)) {
+        myShapeType==TopAbs_EDGE ||
+        myShapeType==TopAbs_FACE ||
+        myShapeType==TopAbs_SOLID)) {
     myErrorStatus=12; // unallowed subshape type
     return;
   }
@@ -803,8 +805,8 @@ void GEOMAlgo_FinderShapeOn1::CheckData()
   myGAS.Load(mySurface);
   aType=myGAS.GetType();
   if (!(aType==GeomAbs_Plane || 
-       aType==GeomAbs_Cylinder ||
-       aType==GeomAbs_Sphere)) {
+        aType==GeomAbs_Cylinder ||
+        aType==GeomAbs_Sphere)) {
     myErrorStatus=14; // unallowed surface type
   }
 }