EgtGeomKernel :

- modifica alla sincronizzazione tra curve per trimming.

CurveAux.cpp

@@ -2711,6 +2711,7 @@ GetIsoPointOnSecondCurve( const ICurve* pCrvEdge1, const ICurve* pCrvEdge2, doub
    Point3d ptCurr1 ;
    Vector3d vtCurr1 ;
    pCrvEdge1->GetPointD1D2( dUCurr1, ICurve::FROM_MINUS, ptCurr1, &vtCurr1) ;
+   vtCurr1.Normalize() ;
    // --- Piano di taglio per punto a minima distanza
    IntersCurvePlane ICP( *pCrvEdge2, ptCurr1, vtCurr1) ;
    int nIndParCloser = - 1, nIndPointCloser = -1 ;
@@ -2786,29 +2787,36 @@ GetIsoPointOnSecondCurve( const ICurve* pCrvEdge1, const ICurve* pCrvEdge2, doub
 
    // Verifico se le direzioni tangenti sono tra di loro circa parallele
    const double COS_ANG_TOL = cos( 15. * DEGTORAD) ;
-   if ( vtCurr1 * vtCurr2 < COS_ANG_TOL) {
-      // Se fuori dalla tolleranza, recupero il miglior versore tangente sul secondo bordo nell'intervallo successivo di lunghezza ( 2. * dMyDist)
+   const double COS_SMALL_ANG_TOL = cos( 4. * DEGTORAD) ;
+   double dSearchLen = dMyDist / 2 ;
+   int NUM_STEP = 10 ;
+   // se poco fuori tolleranza controllo se ho un punto abbastanza vicino con la stessa tangente
+   if ( vtCurr1 * vtCurr2 < COS_SMALL_ANG_TOL) {
+      // Se tanto fuori dalla tolleranza, recupero il miglior versore tangente sul secondo bordo nell'intervallo successivo di lunghezza ( 2. * dMyDist)
+      if ( vtCurr1 * vtCurr2 < COS_ANG_TOL) {
+         dSearchLen = dMyDist ;
+         NUM_STEP = 20 ;
+      }
       pCrvEdge2->GetLengthAtPoint( ptCurr2, dLenCurr2) ;
-      double dLimInfLen2 = Clamp( dLenCurr2 - dMyDist, dLenPrev2, dLen2) ;
-      double dLimSupLen2 = Clamp( dLenCurr2 + dMyDist, dLenPrev2, dLen2) ;
+      double dLimInfLen2 = Clamp( dLenCurr2 - dSearchLen, dLenPrev2, dLen2) ;
+      double dLimSupLen2 = Clamp( dLenCurr2 + dSearchLen, dLenPrev2, dLen2) ;
       // [Controllo migliorabile, magari mendiante metodo di bisezione (?)]
-      const int NUM_STEP = 20 ;
-      double dMinCos = - 1. - EPS_ZERO ;
-      const double DEGTOL = 5. ;
+      
+      double dMin = INFINITO ;
+      const double dCoeff = 0.1 ;
       for ( int i = 0 ; i <= NUM_STEP ; ++ i) {
          double dLen = dLimInfLen2 + i * ( dLimSupLen2 - dLimInfLen2) / NUM_STEP ;
          double dUStep2 ; pCrvEdge2->GetParamAtLength( dLen, dUStep2) ;
          Point3d ptStep2 ; Vector3d vtStep2 = V_NULL ;
          pCrvEdge2->GetPointD1D2( dUStep2, ICurve::FROM_MINUS, ptStep2, &vtStep2) ; vtStep2.Normalize() ;
          double dStepCos2 = vtCurr1 * vtStep2 ;
-         double dAngTol = ( i < NUM_STEP / 2 ? ( 2. * DEGTOL) / NUM_STEP * i :
-            ( - 2. * DEGTOL) / NUM_STEP * ( i - NUM_STEP)) ;
-         double dCosTol = 1. - cos( dAngTol * DEGTORAD) ;
-         if ( dStepCos2 + dCosTol > dMinCos) {
+         double dDist = Dist( ptCurr1, ptCurr2) ;
+         double dVal = (1 - dStepCos2) + dCoeff * dDist / dMyDist ;
+         if ( dVal < dMin) {
             ptCurr2 = ptStep2 ;
             vtCurr2 = vtStep2 ;
             dUCurr2 = dUStep2 ;
-            dMinCos = dStepCos2 + dCosTol ;
+            dMin = dVal ;
          }
       }
    }