EgtGeomKernel 1.6b4 :

- SurfTM aggiunta GetTriangleBoundaryEdges per classificare edge di poligoni.

SurfTriMesh.cpp

@@ -32,9 +32,10 @@ GEOOBJ_REGISTER( SRF_TRIMESH, NGE_S_TRM, SurfTriMesh) ;
 
 //----------------------------------------------------------------------------
 SurfTriMesh::SurfTriMesh( void)
-   : m_nStatus( TO_VERIFY), m_dLinTol( STM_STD_LIN_TOL), m_dSmoothAng( STM_STD_SMOOTH_ANG),
-     m_bClosed( false), m_nTempProp()
+   : m_nStatus( TO_VERIFY), m_dLinTol( STM_STD_LIN_TOL), m_dBoundaryAng( STM_STD_BOUNDARY_ANG),
+     m_dSmoothAng( STM_STD_SMOOTH_ANG), m_bClosed( false), m_nTempProp()
 {
+   m_dCosBndAng = cos( m_dBoundaryAng * DEGTORAD) ;
    m_dCosSmAng = cos( m_dSmoothAng * DEGTORAD) ;
 }
 
@@ -278,12 +279,31 @@ SurfTriMesh::GetNextTriangle( int nId, Triangle3d& Tria) const
    return nId ;  
 }
 
+//----------------------------------------------------------------------------
+bool
+SurfTriMesh::GetTriangleBoundaryEdges( int nId, TriFlags3d& TFlags) const 
+{
+  // verifico esistenza del triangolo 
+   if ( nId < 0  ||  nId >= GetTriangleSize()  ||  m_vTria[nId].nIdVert[0] == SVT_DEL)
+      return false ;
+
+  // verifico i tre lati
+   for ( int i = 0 ; i < 3 ; ++ i) {
+     // indice triangolo adiacente al lato
+      int nT = m_vTria[nId].nIdAdjac[i] ;
+     // se non c'è triangolo adiacente o se forma un angolo oltre il limite, il lato è un contorno
+      TFlags.bFlag[i] = ( nT == SVT_NULL  ||  m_vTria[nId].vtN * m_vTria[nT].vtN < m_dCosBndAng) ;
+   }
+
+   return true ;
+}
+
 //----------------------------------------------------------------------------
 bool
 SurfTriMesh::GetTriangleSmoothNormals( int nId, TriNormals3d& TNrms) const 
 {
   // verifico esistenza del triangolo 
-   if ( nId >= GetTriangleSize()  ||  m_vTria[nId].nIdVert[0] == SVT_DEL)
+   if ( nId < 0  ||  nId >= GetTriangleSize()  ||  m_vTria[nId].nIdVert[0] == SVT_DEL)
       return false ;
 
   // recupero le normali di ciascun vertice