EgtGeomKernel :

- piccole modifiche stilistiche.

IntersLineSurfBez.cpp

@@ -22,13 +22,14 @@
 
 using namespace std ;
 
-//----------------------
-bool
+//----------------------------------------------------------------------------
+static bool
 RefineIntersNewton( const Point3d& ptL, const Vector3d& vtL, double dLen, bool bFinite,
-   const ISurfBezier* pSurfBz, Point3d& ptSP, Point3d& ptIBz) {
-   // la funzione raffina la posisione del punto ptSP, minimizzando la distanza dalla retta e restituisce il punto di intersezione ptIBz
+                    const ISurfBezier* pSurfBz, Point3d& ptSP, Point3d& ptIBz)
+{
+  // la funzione raffina la posisione del punto ptSP, minimizzando la distanza dalla retta e restituisce il punto di intersezione ptIBz
    pSurfBz->GetPointD1D2( ptSP.x / SBZ_TREG_COEFF, ptSP.y / SBZ_TREG_COEFF, ISurfBezier::FROM_MINUS, ISurfBezier::FROM_MINUS, ptIBz) ;
-   // usando un algoritmo di newton cerco di avvicinarmi il più possibile alla retta
+  // usando un algoritmo di newton cerco di avvicinarmi il più possibile alla retta
    DistPointLine dpl( ptIBz, ptL, vtL, dLen, bFinite) ;
    double dDistNew = 0, dDistPre = 0 ;
    dpl.GetDist(dDistNew) ;
@@ -36,8 +37,8 @@ RefineIntersNewton( const Point3d& ptL, const Vector3d& vtL, double dLen, bool b
    int nCount = 0 ;
    double dh = EPS_SMALL ;
    pSurfBz->GetPointD1D2( ptSP.x, ptSP.y, ISurfBezier::FROM_MINUS, ISurfBezier::FROM_MINUS, ptIBz) ;
-   // metodo di newton in più dimensioni
-   // vario sia il parametro U che il parametro V e verifico se la distanza dalla retta diminuisce per scostamenti positivi o negativi.
+  // metodo di newton in più dimensioni
+  // vario sia il parametro U che il parametro V e verifico se la distanza dalla retta diminuisce per scostamenti positivi o negativi.
    while ( dDistNew > EPS_SMALL  &&  nCount < 100) {
       dDistPre = dDistNew ;
       Point3d ptIBzNew1 ;
@@ -50,48 +51,21 @@ RefineIntersNewton( const Point3d& ptL, const Vector3d& vtL, double dLen, bool b
       DistPointLine dplNewV( ptIBzNew2, ptL, vtL, dLen, bFinite) ;
       dplNewV.GetDist( dDistNew) ;
       double dfdV = ( dDistNew - dDistPre) / dh ;
-      //// opzione 0
-      ////scelgo h1 e h2 separatamente e in modo da annullare f(x)
-      //// opzione 1
-      //// valore fisso
-      //double dr = EPS_SMALL ;
-      //if ( dDistPre > 1)
-      //   dr = 1 ;
-      //else if ( dDistPre > 0.1)
-      //   dr = 0.1 ;
-      //else if ( dDistPre > 0.01)
-      //   dr = 0.01 ;
-      //// opzione 2
-      //// valore direttamente vincolato
-      //double dr = dDistPre ;
-      //// opzione 3
-      //// valuto la deformazione locale in base allo spostamento del punto sulla bezier // non serve
-      //double dh1 = Dist( ptIBz, ptIBzNew1) ;
-      //double dh2 = Dist( ptIBz, ptIBzNew2) ;
-      // potrei valutare il nuovo spostamento in base all'ultima variazione di dDist
-      // potrei anche vedere se sto uscendo dal triangolo ( definito nello spazio parametrico)
-      // mi avvicino cercando di annullare la distanza in un colpo solo
+     // mi avvicino cercando di annullare la distanza in un colpo solo
       double dr = - dDistPre / ( dfdU + dfdV) ;
       pSurfBz->GetPointD1D2(( ptSP.x + dr * dfdU) / SBZ_TREG_COEFF, ( ptSP.y + dr * dfdV) / SBZ_TREG_COEFF, ISurfBezier::FROM_MINUS, ISurfBezier::FROM_MINUS, ptIBz) ;
       DistPointLine dplNew( ptIBz, ptL, vtL, dLen, bFinite) ;
       dplNew.GetDist( dDistNew) ;
-      ++nCount ;
+      ++ nCount ;
    }
 
-   return nCount != 99 ;
+   return ( nCount != 99) ;
 }
 
-////----------------------------------------------------------------------------
-//bool
-//RefineIntersBisec( const Point3d& ptL, const Vector3d& vtL, double dLen, bool bFinite,
-//                   const ISurfBezier* pSurfBz, Point3d& ptSP, Point3d& ptIBz) {
-//
-//}
-
 //----------------------------------------------------------------------------
-void
+static void
 UpdateInfoIntersLineSurfBz( const Point3d& ptL, const Vector3d& vtDir, int nILT, int nT, const Point3d& ptSP, const Point3d& ptIBz, double dCos,
-   const Point3d& ptSP2, const Point3d& ptIBz2, double dCos2, ILSBIVECTOR& vInfo)
+                            const Point3d& ptSP2, const Point3d& ptIBz2, double dCos2, ILSBIVECTOR& vInfo)
 {
    if ( nILT == ILTT_IN  ||  nILT == ILTT_EDGE  ||  nILT == ILTT_VERT) {
       double dU = ( ptIBz - ptL) * vtDir ;
@@ -105,7 +79,7 @@ UpdateInfoIntersLineSurfBz( const Point3d& ptL, const Vector3d& vtDir, int nILT,
 }
 
 //----------------------------------------------------------------------------
-void
+static void
 OrderInfoIntersLineSurfBz( ILSBIVECTOR& vInfo)
 {
    // se non trovati, esco 
@@ -120,7 +94,7 @@ OrderInfoIntersLineSurfBz( ILSBIVECTOR& vInfo)
 }
 
 //----------------------------------------------------------------------------
-//  Intersezione di una linea con una superficie TriMesh
+//  Intersezione di una linea con una superficie di Bezier
 //----------------------------------------------------------------------------
 bool
 IntersLineSurfBz( const Point3d& ptL, const Vector3d& vtL, double dLen, const ISurfBezier* pSurfBz,
@@ -257,4 +231,4 @@ FilterLineSurfBzInters( const ILSBIVECTOR& vInfo, INTDBLVECTOR& vInters)
       ++ j ;
    }
    return true ;
-}
+}

IntersLineSurfBez.h

@@ -1 +0,0 @@
-#pragma once