0023442: Provide the access to the functionality of Geom2dHatch_Hacher via IntTools_C...
[occt.git] / src / IntTools / IntTools_Context.cxx
index d0837ef..56486b8 100755 (executable)
@@ -1,7 +1,22 @@
-// File:       IntTools_Context.cxx
-// Created:    Wed Apr  3 16:57:54 2002
-// Author:     Peter KURNEV
-//             <pkv@irinox>
+// Created on: 2002-04-03
+// Created by: Peter KURNEV
+// Copyright (c) 2002-2012 OPEN CASCADE SAS
+//
+// The content of this file is subject to the Open CASCADE Technology Public
+// License Version 6.5 (the "License"). You may not use the content of this file
+// except in compliance with the License. Please obtain a copy of the License
+// at http://www.opencascade.org and read it completely before using this file.
+//
+// The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
+// main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
+//
+// The Original Code and all software distributed under the License is
+// distributed on an "AS IS" basis, without warranty of any kind, and the
+// Initial Developer hereby disclaims all such warranties, including without
+// limitation, any warranties of merchantability, fitness for a particular
+// purpose or non-infringement. Please see the License for the specific terms
+// and conditions governing the rights and limitations under the License.
+
 
 
 #include <IntTools_Context.ixx>
 
 #include <IntTools_Tools.hxx>
 #include <IntTools_FClass2d.hxx>
-// 
 #include <Extrema_LocateExtPC.hxx>
 
 #include <Geom2d_Curve.hxx>
+//
+#include <Geom2dHatch_Intersector.hxx>
+#include <BRepTools.hxx>
+#include <TopExp_Explorer.hxx>
+#include <Geom2d_TrimmedCurve.hxx>
+#include <Precision.hxx>
+#include <Geom2d_Curve.hxx>
+#include <Geom2d_TrimmedCurve.hxx>
 
 //=======================================================================
 //function : 
 //purpose  : 
 //=======================================================================
-  IntTools_Context::IntTools_Context()
+IntTools_Context::IntTools_Context()
 {
 }
 //=======================================================================
@@ -44,6 +66,18 @@ IntTools_Context::~IntTools_Context()
   Standard_Address anAdr;
   Standard_Integer i, aNb;
   //
+  //modified by NIZNHY-PKV Mon Sep 24 08:01:03 2012f
+  Geom2dHatch_Hatcher* pHatcher;
+  aNb=myHatcherMap.Extent();
+  for (i=1; i<=aNb; ++i) {
+    anAdr=myHatcherMap(i);
+    pHatcher=(Geom2dHatch_Hatcher*)anAdr;
+    delete pHatcher;
+  }
+  myHatcherMap.Clear();
+  //
+  //modified by NIZNHY-PKV Mon Sep 24 08:01:07 2012t
+  //
   IntTools_FClass2d* pFClass2d;
   aNb=myFClass2dMap.Extent();
   for (i=1; i<=aNb; ++i) {
@@ -100,11 +134,78 @@ IntTools_Context::~IntTools_Context()
   }
   myProjSDataMap.Clear();
 }
+//modified by NIZNHY-PKV Mon Sep 24 07:56:45 2012f
+//=======================================================================
+//function : Hatcher
+//purpose  : 
+//=======================================================================
+Geom2dHatch_Hatcher& IntTools_Context::Hatcher(const TopoDS_Face& aF)
+{
+  Standard_Address anAdr;
+  Geom2dHatch_Hatcher* pHatcher;
+  //
+  anAdr=myHatcherMap.FindFromKey1(aF);
+  if (!anAdr) {
+    Standard_Real aTolArcIntr, aTolTangfIntr, aTolHatch2D, aTolHatch3D;
+    Standard_Real aUMin, aUMax, aVMin, aVMax, aU1, aU2, aEpsT;
+    TopAbs_Orientation aOrE;
+    Handle(Geom_Surface) aS;
+    Handle(Geom2d_Curve) aC2D;
+    Handle(Geom2d_TrimmedCurve) aCT2D;
+    TopoDS_Face aFF;
+    TopExp_Explorer aExp;
+    //
+    aTolHatch2D=1.e-8;
+    aTolHatch3D=1.e-8;
+    aTolArcIntr=1.e-10;
+    aTolTangfIntr=1.e-10;
+    aEpsT=Precision::PConfusion();
+    //
+    Geom2dHatch_Intersector aIntr(aTolArcIntr, aTolTangfIntr);
+    pHatcher=new Geom2dHatch_Hatcher(aIntr,
+                                    aTolHatch2D, aTolHatch3D,
+                                    Standard_True, Standard_False);
+    
+    //
+    aFF=aF;
+    aFF.Orientation(TopAbs_FORWARD);
+    aS=BRep_Tool::Surface(aFF);
+    //BRepTools::UVBounds(aFF, aUMin, aUMax, aVMin, aVMax);
+    //
+    aExp.Init (aFF, TopAbs_EDGE);
+    for (; aExp.More() ; aExp.Next()) {
+      const TopoDS_Edge& aE=*((TopoDS_Edge*)&aExp.Current());
+      aOrE=aE.Orientation();
+      //
+      aC2D=BRep_Tool::CurveOnSurface (aE, aFF, aU1, aU2);
+      if (aC2D.IsNull() ) {
+       continue;
+      }
+      if (fabs(aU1-aU2) < aEpsT) {
+       continue;
+      }
+      //
+      aCT2D=new Geom2d_TrimmedCurve(aC2D, aU1, aU2);
+      pHatcher->AddElement(aCT2D, aOrE);
+    }// for (; aExp.More() ; aExp.Next()) {
+    //
+    anAdr=(Standard_Address)pHatcher;
+    myHatcherMap.Add(aFF, anAdr);
+  }// if (!anAdr) {
+  //
+  else {
+    Standard_Address *pAdr=(Standard_Address *)anAdr;
+    pHatcher=(Geom2dHatch_Hatcher*)*pAdr;
+  }
+
+  return *pHatcher;
+}
+//modified by NIZNHY-PKV Mon Sep 24 07:56:45 2012t
 //=======================================================================
 //function : FClass2d
 //purpose  : 
 //=======================================================================
-  IntTools_FClass2d& IntTools_Context::FClass2d(const TopoDS_Face& aF)
+IntTools_FClass2d& IntTools_Context::FClass2d(const TopoDS_Face& aF)
 {
   Standard_Address anAdr;
   IntTools_FClass2d* pFClass2d;
@@ -133,7 +234,7 @@ IntTools_Context::~IntTools_Context()
 //function : ProjPS
 //purpose  : 
 //=======================================================================
-  GeomAPI_ProjectPointOnSurf& IntTools_Context::ProjPS(const TopoDS_Face& aF)
+GeomAPI_ProjectPointOnSurf& IntTools_Context::ProjPS(const TopoDS_Face& aF)
 {
   Standard_Address anAdr;
   GeomAPI_ProjectPointOnSurf* pProjPS;
@@ -169,7 +270,7 @@ IntTools_Context::~IntTools_Context()
 //function : ProjPC
 //purpose  : 
 //=======================================================================
-  GeomAPI_ProjectPointOnCurve& IntTools_Context::ProjPC(const TopoDS_Edge& aE)
+GeomAPI_ProjectPointOnCurve& IntTools_Context::ProjPC(const TopoDS_Edge& aE)
 {
   Standard_Address anAdr;
   GeomAPI_ProjectPointOnCurve* pProjPC;
@@ -196,7 +297,7 @@ IntTools_Context::~IntTools_Context()
 //function : ProjPT
 //purpose  : 
 //=======================================================================
-  GeomAPI_ProjectPointOnCurve& IntTools_Context::ProjPT(const Handle(Geom_Curve)& aC3D)
+GeomAPI_ProjectPointOnCurve& IntTools_Context::ProjPT(const Handle(Geom_Curve)& aC3D)
                                                        
 {
   Standard_Address anAdr;
@@ -225,7 +326,7 @@ IntTools_Context::~IntTools_Context()
 //function : SurfaceData
 //purpose  : 
 //=======================================================================
-  IntTools_SurfaceRangeLocalizeData& IntTools_Context::SurfaceData(const TopoDS_Face& aF) 
+IntTools_SurfaceRangeLocalizeData& IntTools_Context::SurfaceData(const TopoDS_Face& aF) 
 {
   Standard_Address anAdr;
   IntTools_SurfaceRangeLocalizeData* pSData;
@@ -253,7 +354,7 @@ IntTools_Context::~IntTools_Context()
 //function : SolidClassifier
 //purpose  : 
 //=======================================================================
-  BRepClass3d_SolidClassifier& IntTools_Context::SolidClassifier(const TopoDS_Solid& aSolid)
+BRepClass3d_SolidClassifier& IntTools_Context::SolidClassifier(const TopoDS_Solid& aSolid)
 {
   Standard_Address anAdr;
   BRepClass3d_SolidClassifier* pSC;
@@ -276,9 +377,9 @@ IntTools_Context::~IntTools_Context()
 //function : ComputeVE
 //purpose  : 
 //=======================================================================
-  Standard_Integer IntTools_Context::ComputeVE(const TopoDS_Vertex& aV1, 
-                                              const TopoDS_Edge&   aE2,
-                                              Standard_Real& aT)
+Standard_Integer IntTools_Context::ComputeVE(const TopoDS_Vertex& aV1, 
+                                            const TopoDS_Edge&   aE2,
+                                            Standard_Real& aT)
 {
   Standard_Boolean bToUpdate;
   Standard_Integer iFlag;
@@ -292,11 +393,11 @@ IntTools_Context::~IntTools_Context()
 //function : ComputeVE
 //purpose  : 
 //=======================================================================
-  Standard_Integer IntTools_Context::ComputeVE(const TopoDS_Vertex& aV1, 
-                                              const TopoDS_Edge&   aE2,
-                                              Standard_Real& aT,
-                                              Standard_Boolean& bToUpdateVertex, 
-                                              Standard_Real& aDist)
+Standard_Integer IntTools_Context::ComputeVE(const TopoDS_Vertex& aV1, 
+                                            const TopoDS_Edge&   aE2,
+                                            Standard_Real& aT,
+                                            Standard_Boolean& bToUpdateVertex, 
+                                            Standard_Real& aDist)
 {
   bToUpdateVertex=Standard_False;
   aDist=0.;
@@ -305,12 +406,6 @@ IntTools_Context::~IntTools_Context()
     return -1;
   }
   //
-  //modified by NIZNHY-PKV Wed Jul 13 08:30:08 2011f
-  //if (!BRep_Tool::IsGeometric(aE2)) { 
-  //  return -2;
-  //}
-  //modified by NIZNHY-PKV Wed Jul 13 08:30:13 2011t
-  //
   Standard_Real aTolV1, aTolE2, aTolSum, aTolVx;
   Standard_Integer aNbProj;
   gp_Pnt aP;
@@ -346,10 +441,10 @@ IntTools_Context::~IntTools_Context()
 //function : ComputeVS
 //purpose  : 
 //=======================================================================
-  Standard_Integer IntTools_Context::ComputeVS(const TopoDS_Vertex& aV1, 
-                                              const TopoDS_Face&   aF2,
-                                              Standard_Real& U,
-                                              Standard_Real& V)
+Standard_Integer IntTools_Context::ComputeVS(const TopoDS_Vertex& aV1, 
+                                            const TopoDS_Face&   aF2,
+                                            Standard_Real& U,
+                                            Standard_Real& V)
 {
   Standard_Real aTolV1, aTolF2, aTolSum, aDist;
   gp_Pnt aP;
@@ -390,8 +485,8 @@ IntTools_Context::~IntTools_Context()
 //function : StatePointFace
 //purpose  : 
 //=======================================================================
-  TopAbs_State IntTools_Context::StatePointFace(const TopoDS_Face& aF,
-                                               const gp_Pnt2d& aP2d)
+TopAbs_State IntTools_Context::StatePointFace(const TopoDS_Face& aF,
+                                             const gp_Pnt2d& aP2d)
 {
   TopAbs_State aState;
   IntTools_FClass2d& aClass2d=FClass2d(aF);
@@ -402,8 +497,8 @@ IntTools_Context::~IntTools_Context()
 //function : IsPointInFace
 //purpose  : 
 //=======================================================================
-  Standard_Boolean IntTools_Context::IsPointInFace(const TopoDS_Face& aF,
-                                                  const gp_Pnt2d& aP2d)
+Standard_Boolean IntTools_Context::IsPointInFace(const TopoDS_Face& aF,
+                                                const gp_Pnt2d& aP2d)
 {
   TopAbs_State aState=StatePointFace(aF, aP2d);
   if (aState==TopAbs_OUT || aState==TopAbs_ON) {
@@ -415,8 +510,8 @@ IntTools_Context::~IntTools_Context()
 //function : IsPointInOnFace
 //purpose  : 
 //=======================================================================
-  Standard_Boolean IntTools_Context::IsPointInOnFace(const TopoDS_Face& aF,
-                                                    const gp_Pnt2d& aP2d)
+Standard_Boolean IntTools_Context::IsPointInOnFace(const TopoDS_Face& aF,
+                                                  const gp_Pnt2d& aP2d)
 { 
   TopAbs_State aState=StatePointFace(aF, aP2d);
   if (aState==TopAbs_OUT) {
@@ -428,9 +523,9 @@ IntTools_Context::~IntTools_Context()
 //function : IsValidPointForFace
 //purpose  : 
 //=======================================================================
-  Standard_Boolean IntTools_Context::IsValidPointForFace(const gp_Pnt& aP,
-                                                        const TopoDS_Face& aF,
-                                                        const Standard_Real aTol) 
+Standard_Boolean IntTools_Context::IsValidPointForFace(const gp_Pnt& aP,
+                                                      const TopoDS_Face& aF,
+                                                      const Standard_Real aTol) 
 {
   Standard_Boolean bFlag;
   Standard_Real Umin, myEpsT, U, V;
@@ -458,10 +553,10 @@ IntTools_Context::~IntTools_Context()
 //function : IsValidPointForFaces
 //purpose  : 
 //=======================================================================
-  Standard_Boolean IntTools_Context::IsValidPointForFaces (const gp_Pnt& aP,
-                                                          const TopoDS_Face& aF1,
-                                                          const TopoDS_Face& aF2,
-                                                          const Standard_Real aTol) 
+Standard_Boolean IntTools_Context::IsValidPointForFaces (const gp_Pnt& aP,
+                                                        const TopoDS_Face& aF1,
+                                                        const TopoDS_Face& aF2,
+                                                        const Standard_Real aTol) 
 {
   Standard_Boolean bFlag1, bFlag2;
 
@@ -476,11 +571,11 @@ IntTools_Context::~IntTools_Context()
 //function : IsValidBlockForFace
 //purpose  : 
 //=======================================================================
-  Standard_Boolean IntTools_Context::IsValidBlockForFace (const Standard_Real aT1,
-                                                         const Standard_Real aT2,
-                                                         const IntTools_Curve& aC, 
-                                                         const TopoDS_Face& aF,
-                                                         const Standard_Real aTol) 
+Standard_Boolean IntTools_Context::IsValidBlockForFace (const Standard_Real aT1,
+                                                       const Standard_Real aT2,
+                                                       const IntTools_Curve& aC, 
+                                                       const TopoDS_Face& aF,
+                                                       const Standard_Real aTol) 
 {
   Standard_Boolean bFlag;
   Standard_Real aTInterm, aFirst, aLast;
@@ -501,12 +596,12 @@ IntTools_Context::~IntTools_Context()
 //function : IsValidBlockForFaces
 //purpose  : 
 //=======================================================================
-  Standard_Boolean IntTools_Context::IsValidBlockForFaces (const Standard_Real aT1,
-                                                          const Standard_Real aT2,
-                                                          const IntTools_Curve& aC, 
-                                                          const TopoDS_Face& aF1,
-                                                          const TopoDS_Face& aF2,
-                                                          const Standard_Real aTol) 
+Standard_Boolean IntTools_Context::IsValidBlockForFaces (const Standard_Real aT1,
+                                                        const Standard_Real aT2,
+                                                        const IntTools_Curve& aC, 
+                                                        const TopoDS_Face& aF1,
+                                                        const TopoDS_Face& aF2,
+                                                        const Standard_Real aTol) 
 {
   Standard_Boolean bFlag1, bFlag2;
   //
@@ -540,10 +635,10 @@ IntTools_Context::~IntTools_Context()
 //function : IsVertexOnLine
 //purpose  : 
 //=======================================================================
-  Standard_Boolean IntTools_Context::IsVertexOnLine (const TopoDS_Vertex& aV,
-                                                    const IntTools_Curve& aC, 
-                                                    const Standard_Real aTolC,
-                                                    Standard_Real& aT)
+Standard_Boolean IntTools_Context::IsVertexOnLine (const TopoDS_Vertex& aV,
+                                                  const IntTools_Curve& aC, 
+                                                  const Standard_Real aTolC,
+                                                  Standard_Real& aT)
 {
   Standard_Boolean bRet;
   Standard_Real aTolV;
@@ -557,11 +652,11 @@ IntTools_Context::~IntTools_Context()
 //function : IsVertexOnLine
 //purpose  : 
 //=======================================================================
-  Standard_Boolean IntTools_Context::IsVertexOnLine (const TopoDS_Vertex& aV,
-                                                    const Standard_Real aTolV,
-                                                    const IntTools_Curve& aC, 
-                                                    const Standard_Real aTolC,
-                                                    Standard_Real& aT)
+Standard_Boolean IntTools_Context::IsVertexOnLine (const TopoDS_Vertex& aV,
+                                                  const Standard_Real aTolV,
+                                                  const IntTools_Curve& aC, 
+                                                  const Standard_Real aTolC,
+                                                  Standard_Real& aT)
 {
   Standard_Real aFirst, aLast, aDist, aTolSum;
   Standard_Integer aNbProj;
@@ -682,9 +777,9 @@ IntTools_Context::~IntTools_Context()
 //function : ProjectPointOnEdge
 //purpose  : 
 //=======================================================================
-  Standard_Boolean IntTools_Context::ProjectPointOnEdge(const gp_Pnt& aP,
-                                                       const TopoDS_Edge& anEdge,
-                                                       Standard_Real& aT)
+Standard_Boolean IntTools_Context::ProjectPointOnEdge(const gp_Pnt& aP,
+                                                     const TopoDS_Edge& anEdge,
+                                                     Standard_Real& aT)
 {
   Standard_Integer aNbPoints;