0028240: Avoid redundant search for span index in evaluation of BSpline cache

diff --git a/src/Adaptor3d/Adaptor3d_TopolTool.cxx b/src/Adaptor3d/Adaptor3d_TopolTool.cxx
index fcf9ae978f0513dbfddbe3a93d60bf5b1731624b..a42d633196f06ce95568b7073753d801f0930b6b 100644 (file)
--- a/src/Adaptor3d/Adaptor3d_TopolTool.cxx
+++ b/src/Adaptor3d/Adaptor3d_TopolTool.cxx
@@ -1090,10 +1090,6 @@ void Adaptor3d_TopolTool::BSplSamplePnts(const Standard_Real theDefl,
   Standard_Real tol = Max(0.01*aDefl2, 1.e-9);
   Standard_Integer l;
 
-  // Calculations of B-spline values will be made using adaptor,
-  // because it caches the data for performance
-  GeomAdaptor_Surface aBSplAdaptor(aBS);
-
   anUFlg(1) = Standard_True;
   anUFlg(nbsu) = Standard_True;
   //myNbSamplesU = 2; 
@@ -1109,12 +1105,10 @@ void Adaptor3d_TopolTool::BSplSamplePnts(const Standard_Real theDefl,
       }
 
       t2 = anUPars(j);
-//      gp_Pnt p1 = aBS->Value(t2, t1);
-      gp_Pnt p1 = aBSplAdaptor.Value(t2, t1);
+      gp_Pnt p1 = myS->Value(t2, t1);
       for(k = j+2; k <= nbsu; ++k) {
        t2 = anUPars(k);
-//     gp_Pnt p2 = aBS->Value(t2, t1);
-       gp_Pnt p2 = aBSplAdaptor.Value(t2, t1);
+       gp_Pnt p2 = myS->Value(t2, t1);
        //gce_MakeLin MkLin(p1, p2);
        //const gp_Lin& lin = MkLin.Value();
 
@@ -1129,8 +1123,7 @@ void Adaptor3d_TopolTool::BSplSamplePnts(const Standard_Real theDefl,
            break;
          }
 
-//       gp_Pnt pp =  aBS->Value(anUPars(l), t1);
-         gp_Pnt pp =  aBSplAdaptor.Value(anUPars(l), t1);
+         gp_Pnt pp =  myS->Value(anUPars(l), t1);
          Standard_Real d = lin.SquareDistance(pp);
          
          if(d <= aDefl2) continue;
@@ -1197,12 +1190,10 @@ void Adaptor3d_TopolTool::BSplSamplePnts(const Standard_Real theDefl,
       }
 
       t2 = aVPars(j);
-//      gp_Pnt p1 = aBS->Value(t1, t2);
-      gp_Pnt p1 = aBSplAdaptor.Value(t1, t2);
+      gp_Pnt p1 = myS->Value(t1, t2);
       for(k = j+2; k <= nbsv; ++k) {
        t2 = aVPars(k);
-//     gp_Pnt p2 = aBS->Value(t1, t2);
-       gp_Pnt p2 = aBSplAdaptor.Value(t1, t2);
+       gp_Pnt p2 = myS->Value(t1, t2);
 
        if(p1.SquareDistance(p2) <= tol) continue;
        //gce_MakeLin MkLin(p1, p2);
@@ -1216,8 +1207,7 @@ void Adaptor3d_TopolTool::BSplSamplePnts(const Standard_Real theDefl,
            break;
          }
 
-//       gp_Pnt pp =  aBS->Value(t1, aVPars(l));
-         gp_Pnt pp =  aBSplAdaptor.Value(t1, aVPars(l));
+         gp_Pnt pp =  myS->Value(t1, aVPars(l));
          Standard_Real d = lin.SquareDistance(pp);
          
          if(d <= aDefl2) continue;

diff --git a/src/BSplCLib/BSplCLib.cxx b/src/BSplCLib/BSplCLib.cxx
index f6a7a298d82128f185c9b618c4b32f1513433b68..5db74e7a299e9bd3f48ee88e27028ef88de98aca 100644 (file)
--- a/src/BSplCLib/BSplCLib.cxx
+++ b/src/BSplCLib/BSplCLib.cxx
@@ -28,6 +28,7 @@
 //                   in TangExtendToConstraint; Continuity can be equal to 0
 
 #include <BSplCLib.hxx>
+#include <ElCLib.hxx>
 #include <gp_Pnt.hxx>
 #include <gp_Pnt2d.hxx>
 #include <gp_Vec.hxx>
@@ -257,15 +258,8 @@ void BSplCLib::LocateParameter
   }
   Standard_Integer Last1 = Last - 1;
   NewU = U;
-  if (IsPeriodic) {
-    Standard_Real Period = ULast - UFirst;
-
-    while (NewU > ULast )
-      NewU  -= Period;
-
-    while (NewU < UFirst)
-      NewU  += Period;
-  }
+  if (IsPeriodic && (NewU < UFirst || NewU > ULast))
+    NewU = ElCLib::InPeriod(NewU, UFirst, ULast);
   
   BSplCLib::Hunt (Knots, NewU, KnotIndex);
   

diff --git a/src/BSplCLib/BSplCLib.hxx b/src/BSplCLib/BSplCLib.hxx
index b2f94783846c924745563dfe17a319a859429ed9..0189aa72499e2d1fbb3b20367421ec0c2f25fc9f 100644 (file)
--- a/src/BSplCLib/BSplCLib.hxx
+++ b/src/BSplCLib/BSplCLib.hxx
@@ -1195,12 +1195,12 @@ public:
   //! Perform the evaluation of the Taylor expansion
   //! of the Bspline normalized between 0 and 1.
   //! Structure of result optimized for BSplCLib_Cache.
-  Standard_EXPORT static void BuildCache (const Standard_Real theParameter, const Standard_Real theSpanDomain, const Standard_Boolean thePeriodicFlag, const Standard_Integer theDegree, const TColStd_Array1OfReal& theFlatKnots, const TColgp_Array1OfPnt& thePoles, const TColStd_Array1OfReal* theWeights, TColStd_Array2OfReal& theCacheArray);
+  Standard_EXPORT static void BuildCache (const Standard_Real theParameter, const Standard_Real theSpanDomain, const Standard_Boolean thePeriodicFlag, const Standard_Integer theDegree, const Standard_Integer theSpanIndex, const TColStd_Array1OfReal& theFlatKnots, const TColgp_Array1OfPnt& thePoles, const TColStd_Array1OfReal* theWeights, TColStd_Array2OfReal& theCacheArray);
   
   //! Perform the evaluation of the Taylor expansion
   //! of the Bspline normalized between 0 and 1.
   //! Structure of result optimized for BSplCLib_Cache.
-  Standard_EXPORT static void BuildCache (const Standard_Real theParameter, const Standard_Real theSpanDomain, const Standard_Boolean thePeriodicFlag, const Standard_Integer theDegree, const TColStd_Array1OfReal& theFlatKnots, const TColgp_Array1OfPnt2d& thePoles, const TColStd_Array1OfReal* theWeights, TColStd_Array2OfReal& theCacheArray);
+  Standard_EXPORT static void BuildCache (const Standard_Real theParameter, const Standard_Real theSpanDomain, const Standard_Boolean thePeriodicFlag, const Standard_Integer theDegree, const Standard_Integer theSpanIndex, const TColStd_Array1OfReal& theFlatKnots, const TColgp_Array1OfPnt2d& thePoles, const TColStd_Array1OfReal* theWeights, TColStd_Array2OfReal& theCacheArray);
   
     static void PolesCoefficients (const TColgp_Array1OfPnt2d& Poles, TColgp_Array1OfPnt2d& CachePoles);
   

diff --git a/src/BSplCLib/BSplCLib_Cache.cxx b/src/BSplCLib/BSplCLib_Cache.cxx
index 275a7b91bafd96e334ad56a844f9478df0c25a5e..59efb2f8f13f086ac691729e1376ad3ccf1f6c0b 100644 (file)
--- a/src/BSplCLib/BSplCLib_Cache.cxx
+++ b/src/BSplCLib/BSplCLib_Cache.cxx
@@ -132,7 +132,7 @@ void BSplCLib_Cache::BuildCache(const Standard_Real&           theParameter,
 
   // Calculate new cache data
   BSplCLib::BuildCache(mySpanStart, mySpanLength, thePeriodic, theDegree, 
-                       theFlatKnots, thePoles2d, theWeights, 
+                       mySpanIndex, theFlatKnots, thePoles2d, theWeights, 
                        myPolesWeights->ChangeArray2());
 }
 
@@ -171,7 +171,7 @@ void BSplCLib_Cache::BuildCache(const Standard_Real&           theParameter,
 
   // Calculate new cache data
   BSplCLib::BuildCache(mySpanStart, mySpanLength, thePeriodic, theDegree, 
-                       theFlatKnots, thePoles, theWeights, 
+                       mySpanIndex, theFlatKnots, thePoles, theWeights, 
                        myPolesWeights->ChangeArray2());
 }
 

diff --git a/src/BSplCLib/BSplCLib_CurveComputation.gxx b/src/BSplCLib/BSplCLib_CurveComputation.gxx
index d9201e2c7700c54322b93dc0ac172834d7e12692..58ad873cf62eedcc570cdc2140a9cdd0816948d3 100644 (file)
--- a/src/BSplCLib/BSplCLib_CurveComputation.gxx
+++ b/src/BSplCLib/BSplCLib_CurveComputation.gxx
@@ -1098,13 +1098,14 @@ void BSplCLib::BuildCache(const Standard_Real          theParameter,
                           const Standard_Real          theSpanDomain,
                           const Standard_Boolean       thePeriodicFlag,
                           const Standard_Integer       theDegree,
+                          const Standard_Integer       theSpanIndex,
                           const TColStd_Array1OfReal&  theFlatKnots,
                           const Array1OfPoints&        thePoles,
                           const TColStd_Array1OfReal*  theWeights,
                                 TColStd_Array2OfReal&  theCacheArray)
 {
   Standard_Real    aParam = theParameter;
-  Standard_Integer anIndex = 0;
+  Standard_Integer anIndex = theSpanIndex;
   Standard_Integer aDimension;
   Standard_Boolean isRational;