0030868: Modeling Algorithms - BRepBuilderAPI_GTransform hangs (2)

Attempt on healing the BSpline surface with the pole in which the DU and DV vectors are almost parallel by moving the pole outside.

src/BRepOffset/BRepOffset.cxx

@@ -309,6 +309,74 @@ Handle(Geom_Surface) BRepOffset::CollapseSingularities (const Handle(Geom_Surfac
       }
     }
 
+    if (aPoles.RowLength() > 2 && aPoles.ColLength() > 2)
+    {
+      // Check if the adjacent sides of the surfaces are collapsed, i.e. if normal is
+      // is impossible to compute in the 'corners'
+      // If such corner is found, move the pole outside a bit, to make the 'corner' more obvious
+
+      // We are dealing with four points:
+      // (UMin, VMin), (UMin, VMax), (UMax, VMin), (UMax, VMax)
+      // or (1, 1), (1, NbVPoles), (NbUPoles, 1), (NbUPoles, NbUPoles)
+
+      Handle(Geom_Surface)& aSurf = aCopy.IsNull() ? aBSpline : aCopy;
+
+      Standard_Real anU[2], aV[2];
+      aBSpline->Bounds (anU[0], anU[1], aV[0], aV[1]);
+
+      Standard_Integer aCornerPoleUIndex[2] = { aPoles.LowerRow(), aPoles.UpperRow() };
+      Standard_Integer aCornerPoleVIndex[2] = { aPoles.LowerCol(), aPoles.UpperCol() };
+
+      Standard_Integer aPoleInsideUIndex[2] = { aCornerPoleUIndex[0] + 1, aCornerPoleUIndex[1] - 1 };
+      Standard_Integer aPoleInsideVIndex[2] = { aCornerPoleVIndex[0] + 1, aCornerPoleVIndex[1] - 1 };
+
+      for (int iU = 0; iU < 2; ++iU)
+      {
+        for (int iV = 0; iV < 2; ++iV)
+        {
+          gp_Pnt aP;
+          gp_Vec aDU, aDV;
+          aSurf->D1 (anU[iU], aV[iV], aP, aDU, aDV);
+
+          gp_Vec aSNorm = aDU ^ aDV;
+          if (aSNorm.SquareMagnitude() < thePrecision * thePrecision)
+          {
+            if (aCopy.IsNull())
+            {
+              aCopy = Handle (Geom_BSplineSurface)::DownCast (theSurface->Copy());
+              aSurf = aCopy;
+            }
+
+            // Move the pole outside surface
+
+            // Get the nearest pole inside
+            gp_Pnt aPInside = aPoles.Value (aPoleInsideUIndex[iU], aPoleInsideVIndex[iV]);
+            gp_Vec aVecOutside (aPInside, aP);
+
+            aDU.Rotate (gp_Ax1 (aP, aDU ^ aVecOutside), M_PI_2);
+            aDV.Rotate (gp_Ax1 (aP, aDV ^ aVecOutside), M_PI_2);
+            
+            aVecOutside = (aDU + aDV).Normalized();
+
+            int k = 1;
+            while (aSNorm.SquareMagnitude() < thePrecision * thePrecision && k <= 1024)
+            {
+              k *= 2;
+              aP.Translate (k * Precision::Confusion() * aVecOutside);
+              aCopy->SetPole (aCornerPoleUIndex[iU], aCornerPoleVIndex[iV], aP);
+
+              aCopy->D1 (anU[iU], aV[iV], aP, aDU, aDV);
+              gp_Vec aN1 = aDU ^ aDV;
+
+              //if (aN1 * aSNorm < 0.0)
+              //  break;
+              aSNorm = aN1;
+            }
+          }
+        }
+      }
+    }
+
     if (! aCopy.IsNull())
       return aCopy;
   }

src/GeomInt/GeomInt_IntSS.hxx

@@ -90,9 +90,6 @@ public:
   
   Standard_EXPORT void TolFixTangents (Standard_Real& aTolCheck, Standard_Real& aTolAngCheck);
   
-  //! Set angular tolerance
-  Standard_EXPORT void SetAngularTolerance (Standard_Real theAngularTolerance);
-  
   //! converts RLine to Geom(2d)_Curve.
   Standard_EXPORT static void TreatRLine (const Handle(IntPatch_RLine)& theRL, const Handle(GeomAdaptor_Surface)& theHS1, const Handle(GeomAdaptor_Surface)& theHS2, Handle(Geom_Curve)& theC3d, Handle(Geom2d_Curve)& theC2d1, Handle(Geom2d_Curve)& theC2d2, Standard_Real& theTolReached);
   

src/GeomInt/GeomInt_IntSS.lxx

@@ -146,12 +146,3 @@ inline GeomInt_IntSS::GeomInt_IntSS ()
 {
   return myIntersector.Point(Index).Value();
 }
-
-//=======================================================================
-//function : SetAngularTolerance
-//purpose  : 
-//=======================================================================
-inline void GeomInt_IntSS::SetAngularTolerance (Standard_Real theAngularTolerance)
-{
-  myIntersector.SetAngularTolerance(theAngularTolerance);
-}

src/IntAna/IntAna_QuadQuadGeo.cxx

@@ -254,53 +254,24 @@ gp_Ax2 DirToAx2(const gp_Pnt& P,const gp_Dir& D)
 //purpose  : Empty constructor
 //=======================================================================
 IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(void)
-: done(Standard_False),
-  nbint(0),
-  typeres(IntAna_Empty),
-  pt1(0,0,0),
-  pt2(0,0,0),
-  pt3(0,0,0),
-  pt4(0,0,0),
-  param1(0),
-  param2(0),
-  param3(0),
-  param4(0),
-  param1bis(0),
-  param2bis(0),
-  myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+    : done(Standard_False),
+      nbint(0),
+      typeres(IntAna_Empty),
+      pt1(0,0,0),
+      pt2(0,0,0),
+      pt3(0,0,0),
+      pt4(0,0,0),
+      param1(0),
+      param2(0),
+      param3(0),
+      param4(0),
+      param1bis(0),
+      param2bis(0),
+      myCommonGen(Standard_False),
+      myPChar(0,0,0)
 {
   InitTolerances();
 }
-
-//=======================================================================
-//function : IntAna_QuadQuadGeo
-//purpose  : Constructor with angular tolerance
-//=======================================================================
-IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(Standard_Real theAngularTolerance, Standard_Boolean theUseAngularTolerance)
-: done(Standard_False),
-  nbint(0),
-  typeres(IntAna_Empty),
-  pt1(0,0,0),
-  pt2(0,0,0),
-  pt3(0,0,0),
-  pt4(0,0,0),
-  param1(0),
-  param2(0),
-  param3(0),
-  param4(0),
-  param1bis(0),
-  param2bis(0),
-  myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(theAngularTolerance),
-  myUseAngularTolerance(theUseAngularTolerance)
-{
-  InitTolerances();
-}
-
 //=======================================================================
 //function : InitTolerances
 //purpose  : 
@@ -308,20 +279,12 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(Standard_Real theAngularTolerance, Standa
 void IntAna_QuadQuadGeo::InitTolerances()
 {
   myEPSILON_DISTANCE               = 1.0e-14;
-  if (myUseAngularTolerance)
-  {
-    myEPSILON_ANGLE_CONE             = myAngularTolerance;
-  }
-  else
-  {
-    myEPSILON_ANGLE_CONE             = Precision::Angular();
-  }
+  myEPSILON_ANGLE_CONE             = Precision::Angular();
   myEPSILON_MINI_CIRCLE_RADIUS     = 0.01*Precision::Confusion();
   myEPSILON_CYLINDER_DELTA_RADIUS  = 1.0e-13;
   myEPSILON_CYLINDER_DELTA_DISTANCE= Precision::Confusion();
   myEPSILON_AXES_PARA              = Precision::Angular();
 }
-
 //=======================================================================
 //function : IntAna_QuadQuadGeo
 //purpose  : Pln  Pln 
@@ -344,9 +307,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(const gp_Pln& P1,
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(P1,P2,TolAng,Tol);
@@ -691,9 +652,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo( const gp_Pln& P
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(P,Co,Tolang,Tol);
@@ -909,9 +868,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(const gp_Pln& P,
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(P,S);
@@ -980,9 +937,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(const gp_Cylinder& Cyl1,
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(Cyl1,Cyl2,Tol);
@@ -1267,9 +1222,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(const gp_Cylinder& Cyl,
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(Cyl,Con,Tol);
@@ -1325,9 +1278,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(const gp_Cylinder& Cyl,
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(Cyl,Sph,Tol);
@@ -1393,9 +1344,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(const gp_Cylinder& Cyl,
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(Con1,Con2,Tol);
@@ -1851,9 +1800,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(const gp_Cylinder& Cyl,
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(Sph,Con,Tol);
@@ -1965,9 +1912,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(const gp_Cylinder& Cyl,
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(Sph1,Sph2,Tol);
@@ -2081,9 +2026,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(const gp_Pln& Pln,
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(Pln,Tor,Tol);
@@ -2195,9 +2138,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(const gp_Cylinder& Cyl,
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(Cyl,Tor,Tol);
@@ -2280,9 +2221,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(const gp_Cone& Con,
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(Con,Tor,Tol);
@@ -2420,9 +2359,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(const gp_Sphere& Sph,
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(Sph,Tor,Tol);
@@ -2519,9 +2456,7 @@ IntAna_QuadQuadGeo::IntAna_QuadQuadGeo(const gp_Torus& Tor1,
   param1bis(0),
   param2bis(0),
   myCommonGen(Standard_False),
-  myPChar(0,0,0),
-  myAngularTolerance(0.0),
-  myUseAngularTolerance(Standard_False)
+  myPChar(0,0,0)
 {
   InitTolerances();
   Perform(Tor1,Tor2,Tol);
@@ -2757,7 +2692,6 @@ const gp_Pnt& IntAna_QuadQuadGeo::PChar() const
 {
   return myPChar;
 }
-
 //=======================================================================
 //function : RefineDir
 //purpose  : 

src/IntAna/IntAna_QuadQuadGeo.hxx

@@ -61,9 +61,6 @@ public:
   //! Empty constructor.
   Standard_EXPORT IntAna_QuadQuadGeo();
   
-  //! Constructor with angular tolerance
-  Standard_EXPORT IntAna_QuadQuadGeo(Standard_Real theAngularTolerance, Standard_Boolean theUseAngularTolerance);
-  
   //! Creates the intersection between two planes.
   //! TolAng is the angular tolerance used to determine
   //! if the planes are parallel.
@@ -264,8 +261,7 @@ protected:
   Standard_Real myEPSILON_AXES_PARA;
   Standard_Boolean myCommonGen;
   gp_Pnt myPChar;
-  Standard_Real myAngularTolerance;
-  Standard_Boolean myUseAngularTolerance;
+
 
 private:
 

src/IntPatch/IntPatch_ImpImpIntersection.hxx

@@ -108,8 +108,8 @@ public:
   //! An exception is raised if Index<=0 or Index>NbLine.
     const Handle(IntPatch_Line)& Line (const Standard_Integer Index) const;
 
-  //! Sets the angular tolerance
-  void SetAngularTolerance (Standard_Real theAngularTolerance);
+
+
 
 protected:
 
@@ -128,8 +128,7 @@ private:
   IntPatch_SequenceOfPoint spnt;
   IntPatch_SequenceOfLine slin;
   IntPatch_TheSOnBounds solrst;
-  Standard_Real myAngularTolerance;
-  Standard_Boolean myUseAngularTolerance;
+
 
 };
 

src/IntPatch/IntPatch_ImpImpIntersection.lxx

@@ -71,9 +71,3 @@ inline const Handle(IntPatch_Line)& IntPatch_ImpImpIntersection::Line (const Sta
   if (!IsDone ()) { throw StdFail_NotDone(); }
   return slin(Index);
 }
-
-inline void IntPatch_ImpImpIntersection::SetAngularTolerance (Standard_Real theAngularTolerance)
-{
-  myAngularTolerance = theAngularTolerance;
-  myUseAngularTolerance = Standard_True;
-}

src/IntPatch/IntPatch_ImpImpIntersection_1.gxx

@@ -125,9 +125,7 @@ static Standard_Boolean IntCoCo(const IntSurf_Quadric&,
 				Standard_Boolean&,
 				Standard_Boolean&,
 				IntPatch_SequenceOfLine&,
-				IntPatch_SequenceOfPoint&,
-				Standard_Boolean theUseAngularTolerance = Standard_False,
-				Standard_Real theAngularTolerance = 0.0);
+				IntPatch_SequenceOfPoint&);
 
 //torus
 static Standard_Boolean IntPTo(const IntSurf_Quadric&,

src/IntPatch/IntPatch_ImpImpIntersection_2.gxx

@@ -29,8 +29,7 @@ IntPatch_ImpImpIntersection::IntPatch_ImpImpIntersection ():
 myDone(IntStatus_Fail),
 empt(Standard_True),
 tgte(Standard_False),
-oppo(Standard_False),
-myUseAngularTolerance(Standard_False)
+oppo(Standard_False)
 {
 }
 //=======================================================================
@@ -232,7 +231,7 @@ void IntPatch_ImpImpIntersection::Perform(const Handle(Adaptor3d_Surface)&  S1,
     }
     //
     case 33: { // Cone/Cone
-      if (!IntCoCo(quad1, quad2, TolTang, empt, SameSurf, multpoint, slin, spnt, myUseAngularTolerance, myAngularTolerance)) {
+      if (!IntCoCo(quad1, quad2, TolTang, empt, SameSurf, multpoint, slin, spnt)) {
         return;
       }
       bEmpty = empt;

src/IntPatch/IntPatch_ImpImpIntersection_5.gxx

@@ -27,9 +27,7 @@ Standard_Boolean IntCoCo(const IntSurf_Quadric& Quad1,
 			 Standard_Boolean& Same,
 			 Standard_Boolean& Multpoint,
 			 IntPatch_SequenceOfLine& slin,
-			 IntPatch_SequenceOfPoint& spnt,
-			 Standard_Boolean theUseAngularTolerance,
-			 Standard_Real theAngularTolerance)
+			 IntPatch_SequenceOfPoint& spnt)
 
 {
   Standard_Integer i, NbSol;
@@ -40,8 +38,7 @@ Standard_Boolean IntCoCo(const IntSurf_Quadric& Quad1,
   gp_Cone Co1(Quad1.Cone());
   gp_Cone Co2(Quad2.Cone());
   //
-  IntAna_QuadQuadGeo inter(theAngularTolerance, theUseAngularTolerance);
-  inter.Perform(Co1,Co2,Tol);
+  IntAna_QuadQuadGeo inter(Co1,Co2,Tol);
   if (!inter.IsDone()) {
     return Standard_False;
   }

src/IntPatch/IntPatch_Intersection.cxx

@@ -53,8 +53,7 @@ IntPatch_Intersection::IntPatch_Intersection ()
    myU1Start(0.0),
    myV1Start(0.0),
    myU2Start(0.0),
-   myV2Start(0.0),
-   myUseAngularTolerance(Standard_False)
+   myV2Start(0.0)
 {
 }
 
@@ -77,8 +76,7 @@ IntPatch_Intersection::IntPatch_Intersection(const Handle(Adaptor3d_Surface)&  S
    myU1Start(0.0),
    myV1Start(0.0),
    myU2Start(0.0),
-   myV2Start(0.0),
-   myUseAngularTolerance(Standard_False)
+   myV2Start(0.0)
 {
   if(myTolArc<1e-8) myTolArc=1e-8;
   if(myTolTang<1e-8) myTolTang=1e-8;
@@ -104,8 +102,7 @@ IntPatch_Intersection::IntPatch_Intersection(const Handle(Adaptor3d_Surface)&  S
    myU1Start(0.0),
    myV1Start(0.0),
    myU2Start(0.0),
-   myV2Start(0.0),
-   myUseAngularTolerance(Standard_False)
+   myV2Start(0.0)
 {
   Perform(S1,D1,TolArc,TolTang);
 }
@@ -132,12 +129,6 @@ void IntPatch_Intersection::SetTolerances(const Standard_Real TolArc,
   if(myUVMaxStep>0.5) myUVMaxStep=0.5;
 }
 
-void IntPatch_Intersection::SetAngularTolerance (Standard_Real theAngularTolerance)
-{
-  myAngularTolerance = theAngularTolerance;
-  myUseAngularTolerance = Standard_True;
-}
-
 //======================================================================
 // function: Perform
 //======================================================================
@@ -1363,13 +1354,8 @@ void IntPatch_Intersection::GeomGeomPerfom(const Handle(Adaptor3d_Surface)& theS
                                            const GeomAbs_SurfaceType theTyps2,
                                            const Standard_Boolean theIsReqToKeepRLine)
 {
-  IntPatch_ImpImpIntersection interii;
-
-  if (myUseAngularTolerance)
-  {
-    interii.SetAngularTolerance(myAngularTolerance);
-  }
-  interii.Perform(theS1, theD1, theS2, theD2, myTolArc, myTolTang, theIsReqToKeepRLine);
+  IntPatch_ImpImpIntersection interii(theS1,theD1,theS2,theD2,
+                                      myTolArc,myTolTang, theIsReqToKeepRLine);
 
   if (!interii.IsDone())
   {

src/IntPatch/IntPatch_Intersection.hxx

@@ -62,8 +62,6 @@ public:
   //! points in their respective parametric spaces.
   Standard_EXPORT void SetTolerances (const Standard_Real TolArc, const Standard_Real TolTang, const Standard_Real UVMaxStep, const Standard_Real Fleche);
   
-  Standard_EXPORT void SetAngularTolerance (Standard_Real theAngularTolerance);
-  
   //! Flag theIsReqToKeepRLine has been entered only for
   //! compatibility with TopOpeBRep package. It shall be deleted
   //! after deleting TopOpeBRep.
@@ -199,8 +197,7 @@ private:
   Standard_Real myV1Start;
   Standard_Real myU2Start;
   Standard_Real myV2Start;
-  Standard_Real myAngularTolerance;
-  Standard_Boolean myUseAngularTolerance;
+
 
 };