- in Lapjoint -> MakePocket aggiunto messaggio in caso si rovini il pezzo successivo

Finish Release: 2.5i2

Compile.bat

@@ -1,5 +1,5 @@
 
-REM Compilazione degli script Beam Egaltech                         2023.02.27
+REM Compilazione degli script Beam Egaltech                         2023.09.26
 REM Per togliere info di debug aggiungere flag -s prima del nome del file di input
 
 REM Compilazione 32 e 64 bit
@@ -9,6 +9,7 @@ REM Compilazione 32 e 64 bit
 \EgtProg\Dll32\luac54 -o bin\LuaLibs\DiceCut.lua                       LuaLibs\DiceCut.lua
 \EgtProg\Dll32\luac54 -o bin\LuaLibs\FaceByPocket.lua                  LuaLibs\FaceByPocket.lua
 \EgtProg\Dll32\luac54 -o bin\LuaLibs\FacesBySaw.lua                    LuaLibs\FacesBySaw.lua
+\EgtProg\Dll32\luac54 -o bin\LuaLibs\FeatureTopology.lua               LuaLibs\FeatureTopology.lua
 \EgtProg\Dll32\luac54 -o bin\LuaLibs\MachiningLib.lua                  LuaLibs\MachiningLib.lua
 \EgtProg\Dll32\luac54 -o bin\LuaLibs\ProcessBlockHausFront.lua         LuaLibs\ProcessBlockHausFront.lua
 \EgtProg\Dll32\luac54 -o bin\LuaLibs\ProcessChamfer.lua                LuaLibs\ProcessChamfer.lua

LuaLibs/BeamExec.lua

@@ -46,6 +46,9 @@
 -- 2023/02/20 Ora le mortase a coda di rondine laterali sono sempre fatte prima dei tagli longitudinali. 
 -- 2023/03/31 Corretto ordinamento per fori di coda da lasciare in coda.
 -- 2023/07/31 Corretto errore nelle mortase in doppio.
+-- 2023/09/13 Aggiunta ClassifyTopology per la classificazione topologica delle feature. In CollectFeatures aggiunta la raccolta preliminare di alcune informazioni.
+-- 2022/09/26 Funzione IsFeatureCuttingEntireSection spostata in BeamLib
+-- 2022/09/26 In ClassifyTopology aggiunto passaggio del parametro nRawId
 
 -- Tabella per definizione modulo
 local BeamExec = {}
@@ -102,6 +105,7 @@ _G.package.loaded.ProcessFreeContour = nil
 _G.package.loaded.ProcessDecor = nil
 local ML = require( 'MachiningLib')
 local BL = require( 'BeamLib')
+local Topology = require( 'FeatureTopology')
 local DC = require( 'DiceCut')
 local Fbs = require( 'FacesBySaw')
 local Hcut= require( 'ProcessHeadCut')
@@ -247,12 +251,6 @@ local function IsTailFeature( Proc, b3Raw, dCurrOvmH, dCurrOvmT)
    return false
 end
 
--------------------------------------------------------------------------------------------------------------
--- restituisce vero se la feature con box b3Proc taglia l'intera sezione della barra, rappresentata dalle sue dimensioni W e H
-local function IsFeatureCuttingEntireSection( b3Proc, dRawW, dRawH)
-   return ((abs(b3Proc:getDimY() - dRawW) < 10 * GEO.EPS_SMALL or b3Proc:getDimY() > dRawW) and (abs(b3Proc:getDimZ() - dRawH) < 10 * GEO.EPS_SMALL or b3Proc:getDimZ() > dRawH))
-end
-
 -------------------------------------------------------------------------------------------------------------
 local function CollectFeatures( PartId, b3Raw, dCurrOvmH, dCurrOvmT)
    local dRawW = b3Raw:getDimY()
@@ -280,6 +278,7 @@ local function CollectFeatures( PartId, b3Raw, dCurrOvmH, dCurrOvmT)
             local nAddMainId = EgtGetInfo( ProcId, 'MAINID', 'i')
             if nGrp and nPrc and nDo == 1 then
                local Proc = {}
+               Proc.PartId = PartId
                Proc.Id = ProcId
                Proc.Grp = nGrp
                Proc.Prc = nPrc
@@ -297,10 +296,19 @@ local function CollectFeatures( PartId, b3Raw, dCurrOvmH, dCurrOvmT)
                   Proc.MainId = Proc.Id + nAddMainId
                end
                Proc.Box = EgtGetBBoxGlob( ProcId, GDB_BB.STANDARD)
+               if b3Raw then
+                  -- recupero l'elenco delle facce della parte interessate dalla feature
+                  Proc.AffectedFaces = BL.GetProcessAffectedFaces( Proc)
+                  -- recupero informazioni sulle facce della feature
+                  Proc.Face = {}
+                  for i = 1, Proc.Fct do
+                     Proc.Face[i] = { Id = i - 1, VtN = EgtSurfTmFacetNormVersor( Proc.Id, i - 1, GDB_ID.ROOT ), Elevation = BL.GetFaceElevation( Proc.Id, i - 1, PartId)}
+                  end
+               end
                if Proc.Box and not Proc.Box:isEmpty() then 
                   Proc.Head = IsHeadFeature( Proc, b3Raw, dCurrOvmH)
                   Proc.Tail, Proc.AdvTail = IsTailFeature( Proc, b3Raw, dCurrOvmH, dCurrOvmT)
-                  if Proc.Fct == 1 and IsFeatureCuttingEntireSection( Proc.Box , dRawW, dRawH) and ( Proc.Head or Proc.Tail) and Proc.Prc ~= 340 and Proc.Prc ~= 350 then
+                  if Proc.Fct == 1 and BL.IsFeatureCuttingEntireSection( Proc.Box, dRawW, dRawH) and ( Proc.Head or Proc.Tail) and Proc.Prc ~= 340 and Proc.Prc ~= 350 then
                      if Proc.Head and Proc.Box:getCenter():getX() < dMachineBeforeIntersectingDrillingsXHead then
                         dMachineBeforeIntersectingDrillingsXHead = Proc.Box:getCenter():getX()
                         nMachineBeforeIntersectingDrillingsIdHead = Proc.Id
@@ -324,6 +332,7 @@ local function CollectFeatures( PartId, b3Raw, dCurrOvmH, dCurrOvmT)
                         Proc.Tail = Drill.IsTailFeature( Proc, b3Raw, dCurrOvmH, dCurrOvmT)
                         -- definisco dati seconda parte
                         local Proc2 = {}
+                        Proc2.PartId = PartId
                         Proc2.Id = ProcId
                         Proc2.Grp = nGrp
                         Proc2.Prc = nPrc
@@ -343,7 +352,7 @@ local function CollectFeatures( PartId, b3Raw, dCurrOvmH, dCurrOvmT)
                      end
                   -- se BlockHaus HalfLap
                   elseif Proc.Prc == 37 then
-                     local nFacInd = BL.GetFaceWithMostAdj( Proc.Id, PartId)
+                     local nFacInd = BL.GetFaceWithMostAdj( Proc, PartId)
                      if nFacInd then
                         local vtN = EgtSurfTmFacetNormVersor( Proc.Id, nFacInd, GDB_ID.ROOT)
                         if vtN then
@@ -387,7 +396,7 @@ local function AnalyzeHeadFeatures( b3Solid, vProc, dRawW, dRawH)
       end
       if Proc.Head and Proc.Id ~= nReplacedFeatureId and Proc.Prc ~= 340 then
          -- controllo se la feature taglia l'intera sezione; in caso positivo la finitura non è necessaria
-         bHeadFinishingNeeded = not IsFeatureCuttingEntireSection( Proc.Box, dRawW, dRawH)
+         bHeadFinishingNeeded = not BL.IsFeatureCuttingEntireSection( Proc.Box, dRawW, dRawH)
       end
    end
    return bHeadFinishingNeeded, nReplacedFeatureId
@@ -410,7 +419,7 @@ local function AnalyzeTailFeatures( b3Solid, vProc, dRawW, dRawH)
       end
       if Proc.Tail and Proc.Id ~= nReplacedFeatureId and Proc.Prc ~= 350 then
          -- controllo se la feature taglia l'intera sezione; in caso positivo la finitura non è necessaria
-         bTailFinishingNeeded = not IsFeatureCuttingEntireSection( Proc.Box, dRawW, dRawH)
+         bTailFinishingNeeded = not BL.IsFeatureCuttingEntireSection( Proc.Box, dRawW, dRawH)
       end
    end
    return bTailFinishingNeeded, nReplacedFeatureId
@@ -651,11 +660,11 @@ local function PrintFeatures( vProc, b3Raw)
    EgtOutLog( ' RawBox=' .. tostring( b3Raw))
    for i = 1, #vProc do
       local Proc = vProc[i]
-      local sOut = string.format( ' Id=%3d Grp=%1d Prc=%3d TC=%2d/%d Flg=%2d Down=%s Side=%s Head=%s Tail=%s Fcse=%1d,%1d Diam=%.2f Fct=%2d Box=%s',
+      local sOut = string.format( ' Id=%3d Grp=%1d Prc=%3d TC=%2d/%d Flg=%2d Down=%s Side=%s Head=%s Tail=%s Fcse=%1d,%1d Diam=%.2f Fct=%2d Box=%s TopoName=%s',
                                   Proc.Id, Proc.Grp, Proc.Prc, Proc.TaskId, Proc.CutId,
                                   Proc.Flg, EgtIf( Proc.Down, 'T', 'F'), EgtIf( Proc.Side, 'T', 'F'),
                                   EgtIf( Proc.Head, 'T', 'F'), EgtIf( Proc.Tail, 'T', EgtIf( Proc.AdvTail, 'A', 'F')),
-                                  Proc.Fcs, Proc.Fce, Proc.Diam, Proc.Fct, tostring( Proc.Box))
+                                  Proc.Fcs, Proc.Fce, Proc.Diam, Proc.Fct, tostring( Proc.Box), Proc.TopologyLongName or '')
       -- info speciali per Block Haus Half Lap
       if Proc.Prc == 37 then
          local sSpec = string.format( ' N=%s Hd=%s', tostring( Proc.vtN or V_NULL()), EgtIf( Proc.HeadDir, 'T', 'F'))
@@ -1137,6 +1146,20 @@ local function ClassifyFeatures( vProc, b3Raw, Stats)
    return bAllOk, bSomeDown, bSomeSide, bSplitRot
 end
 
+-------------------------------------------------------------------------------------------------------------
+local function ClassifyTopology( vProc, nRawId)
+   local b3Raw = EgtGetRawPartBBox( nRawId)
+   local nRecognized = 0
+   for i = 1, #vProc do
+      local Proc = vProc[i]
+      if Topology.Classify( Proc, b3Raw) then
+         nRecognized = nRecognized + 1
+      end
+   end
+   
+   return nRecognized
+end
+
 -------------------------------------------------------------------------------------------------------------
 local function AddFeatureMachining( Proc, nPhase, nRawId, nPartId, dCurrOvmH, bNeedHCut, b3Raw, nOrd, sDownOrSideOrStd, bPreMove, vtMove, dCurrOvmT)
    local bOk = true
@@ -1816,6 +1839,8 @@ function BeamExec.ProcessFeatures()
       if BD.TWO_EQUAL_HEADS or BD.DOWN_HEAD then
          SetMirroredFeatures( vProc, b3Raw)
       end
+      -- classifico topologicamente le feature
+      ClassifyTopology( vProc, nRawId)
       -- le ordino lungo X
       OrderFeatures( vProc, b3Raw)
       -- le classifico

LuaLibs/BeamLib.lua

@@ -22,6 +22,10 @@
 -- 2023/01/31 Aggiunta funzione ConvertToClosedCurve, precedentemente parte di ProcessMortise.Make
 -- 2023/02/22 Modifiche a SetOpenSide, aggiunte ChangeOrOpenStart e CurveWithOnlyStraightLines.
 -- 2023/06/12 In ChangeOrOpenStart corretta ricerca segmento più lungo.
+-- 2023/09/13 Aggiunte funzioni Is3EdgesApprox e GetProcessAffectedFaces.
+-- 2023/09/25 In GetFaceWithMostAdj aggiunta verifica sottosquadro anche per facce non adiacenti.
+-- 2023/09/26 In GetFaceWithMostAdj gestito primo parametro anche come Proc; gestito caso strip con facce tutte in sottosquadro.
+-- 2023/09/26 Spostata qui funzione IsFeatureCuttingEntireSection da BeamExec.
 
 -- Tabella per definizione modulo
 local BeamLib = {}
@@ -509,7 +513,8 @@ function BeamLib.GetFaceElevationFromPointDir( nSurfId, nPartId, ptC, vtN, nIdGe
 end
 
 ---------------------------------------------------------------------
-function BeamLib.GetFaceWithMostAdj( nSurfId, nPartId, bCompare3Fc, dCosSideAng)
+function BeamLib.GetFaceWithMostAdj( Proc, nPartId, bCompare3Fc, dCosSideAng)
+   local nSurfId = Proc.Id
    -- verifica che la superficie non sia chiusa e quindi non lavorabile
    if EgtSurfIsClosed( nSurfId) then
       return
@@ -521,14 +526,17 @@ function BeamLib.GetFaceWithMostAdj( nSurfId, nPartId, bCompare3Fc, dCosSideAng)
    end
    -- recupero le normali delle facce
    local vvtN = {}
+   local vPtC = {}
    for i = 1, nFacCnt do
-      local vtN = EgtSurfTmFacetNormVersor( nSurfId, i - 1, GDB_ID.ROOT)
+      local ptC, vtN = EgtSurfTmFacetCenter( nSurfId, i - 1, GDB_ID.ROOT)
       vvtN[i] = vtN ;
+      vPtC[i] = ptC
    end
    -- adiacenze e sottosquadra delle facce
    local vAdj = {}
    local vUcut = {}
    local vOrtho = {}
+   local nFacesWithUnderCut = 0
    for i = 1, nFacCnt do
       -- recupero le adiacenze del loop esterno
       local vFacAdj = EgtSurfTmFacetAdjacencies( nSurfId, i - 1)[1]
@@ -540,7 +548,7 @@ function BeamLib.GetFaceWithMostAdj( nSurfId, nPartId, bCompare3Fc, dCosSideAng)
          end
       end
       vAdj[i] = nCount
-      -- ne determino eventuale sottosquadra ( dal valore passato o - 3deg) e ortogonalità
+      -- ne determino eventuale sottosquadro ( dal valore passato o - 3deg) e ortogonalità, per facce adiacenti
       local bUcut = false
       local bOrtho = true
       for j = 1, #vFacAdj do
@@ -556,39 +564,58 @@ function BeamLib.GetFaceWithMostAdj( nSurfId, nPartId, bCompare3Fc, dCosSideAng)
             end
          end
       end
+      -- determino evenutale sottosquadro per facce non adiacenti
+      for j = 1, nFacCnt do
+         local bIsFaceAdjacent = false
+         for k = 1, #vFacAdj do
+            if j - 1 == vFacAdj[k] then
+               bIsFaceAdjacent = true
+               break
+            end
+         end
+         if not ( bIsFaceAdjacent or ( j == i)) then
+            -- recupero gruppo per geometria addizionale
+            local nAddGrpId = BeamLib.GetAddGroup( nPartId)
+            if not nAddGrpId then
+               EgtOutLog( 'Error : missing AddGroup')
+               return 0, 0, 0
+            end
+            -- verifico eventuale intersezione tra la faccia i-1 esima e la proiezione dell'altra sulla stessa
+            local nShadowFacetId = EgtCopySurfTmFacet( nSurfId, j - 1, nAddGrpId)
+            local nMasterContourId = EgtExtractSurfTmFacetLoops( nSurfId, i - 1, nAddGrpId)
+            local nShadowContourId = EgtExtractSurfTmLoops( nShadowFacetId, nAddGrpId)
+            EgtCutSurfTmPlane( nShadowFacetId, vPtC[i], -vvtN[i], false, GDB_RT.GLOB)
+            if nMasterContourId and nShadowContourId then
+               local nMasterSurfFlatId = EgtSurfFlatRegion( nAddGrpId, nMasterContourId)
+               local frOCSMaster = Frame3d( vPtC[i], vvtN[i])
+               EgtScale( nShadowContourId, frOCSMaster, 1, 1, 0, GDB_RT.GLOB)
+               EgtModifyCurveExtrusion( nShadowContourId, vvtN[i], GDB_RT.GLOB)
+               local dResV = vvtN[i] * vvtN[j]
+               EgtOffsetCurve( nShadowContourId, EgtIf( dResV > GEO.EPS_SMALL, - 50, 50) * GEO.EPS_SMALL)
+               local nShadowSurfFlatId = EgtSurfFlatRegion( nAddGrpId, nShadowContourId)
+               local bShadowSurfIsExternal = EgtSurfFrTestExternal( nMasterSurfFlatId, nShadowSurfFlatId)
+               if not bShadowSurfIsExternal then
+                  bUcut = true
+               end
+               EgtErase( { nMasterSurfFlatId, nShadowSurfFlatId})
+            end
+            EgtErase( { nMasterContourId, nShadowContourId, nShadowFacetId})
+         end
+      end
+      if bUcut then
+         nFacesWithUnderCut = nFacesWithUnderCut + 1
+      end
       vUcut[i] = bUcut
       vOrtho[i] = bOrtho
    end
-   -- se 4 facce tutte con adiacenza 2, allora è un tunnel
-   if nFacCnt == 4 then
-      if vAdj[1] == 2 and vAdj[2] == 2 and vAdj[3] == 2 and vAdj[4] == 2 then
-         -- se tutte le facce sono ortogonali tra loro esco con un flag che ne indica questa propietà
-         if vOrtho[1] == true and vOrtho[2] == true and vOrtho[3] == true and vOrtho[4] == true then
-            return -1, GEO.INFINITO, true
-         else         
-            return -1, GEO.INFINITO
-         end
-      end
-   end
-   -- se 3 facce con una che ha 2 adiacenze e le altre hanno 1 adiacenza, allora è una semi-fessura
-   if bCompare3Fc and nFacCnt == 3 then
-      local nCount2Adc = 0
-      local nCount1Adc = 0
-      -- ottengo il numero di facce con due adiacenze e il numero di facce con una adiacenza
-      for i = 1, #vAdj do
-         if vAdj[i] == 2 then
-            nCount2Adc = nCount2Adc + 1
-         elseif vAdj[i] == 1 then
-            nCount1Adc = nCount1Adc + 1
-         end
-      end
-      -- se il numero di adiacenze corrisponde
-      if nCount2Adc == 1 and nCount1Adc == 2 then
-         if vOrtho[1] == true and vOrtho[2] == true and vOrtho[3] == true then
-            return -1, GEO.INFINITO, true
-         else         
-            return -1, GEO.INFINITO
-         end
+   local bEveryFaceHasUndercut = ( nFacesWithUnderCut == nFacCnt)
+   -- tunnel o assimilabile
+   if Proc.Topology == 'Tunnel' or ( bCompare3Fc and Proc.Topology == 'Groove' and Proc.IsThrough) or ( Proc.Topology == 'Strip' and bEveryFaceHasUndercut) then
+      if Proc.IsParallel or Proc.AllRightAngles then
+         return -1, GEO.INFINITO, true
+      else
+         -- non gestito nella LapJoint
+         return -1, GEO.INFINITO
       end
    end
    -- recupero le facce non in sottosquadra e con il maggior numero di adiacenze
@@ -1058,5 +1085,67 @@ function BeamLib.CurveWithOnlyStraightLines( nPathInt)
    return true
 end
 
+---------------------------------------------------------------------
+-- Funzione per determinare se la faccia ha lati molto corti (trascurabili) ed è quindi approssimabile ad una 3 facce
+function BeamLib.Is3EdgesApprox( Proc, nFacet, nAddGrpId)
+   nAddGrpId = nAddGrpId or BeamLib.GetAddGroup( Proc.PartId)
+   local bResult = false
+   local nContourId, nContourCnt = EgtExtractSurfTmFacetLoops( Proc.Id, nFacet, nAddGrpId)
+   if not nContourId then return false end
+   EgtMergeCurvesInCurveCompo( nContourId)
+   -- recupero il numero effettivo di lati
+   local _, nEntityCount = EgtCurveDomain( nContourId)
+   local nEdges = nEntityCount
+   if nEntityCount and nEntityCount  == 3 then
+      bResult = true
+   -- rimuovo i lati molto corti dal conteggio totale
+   elseif nEntityCount then
+      for i = 1, nEntityCount do
+         local dLength = EgtCurveCompoLength( nContourId, i - 1)
+         if dLength < 15 then nEdges = nEdges - 1 end
+      end
+   end
+   if nEdges == 3 then bResult = true end
+   -- cancello tutti i contorni appena creati
+   EgtErase( EgtTableFill( nContourId, nContourCnt))
+   return bResult
+end
+
+-------------------------------------------------------------------------------------------------------------
+-- restituisce le facce della parte interessate dalla feature Proc
+function BeamLib.GetProcessAffectedFaces( Proc)
+   local nBoxSolidId = EgtGetFirstNameInGroup( Proc.PartId or GDB_ID.NULL, 'Box')
+   local b3Part = EgtGetBBoxGlob( nBoxSolidId, GDB_BB.STANDARD)
+   local vtFacesAffected = { Top = false, Bottom = false, Front = false, Back = false, Left = false, Right = false}
+   if Proc.Box and not Proc.Box:isEmpty() then
+      if Proc.Box:getMax():getZ() > b3Part:getMax():getZ() - 500 * GEO.EPS_SMALL then
+         vtFacesAffected.Top = true
+      end
+      if Proc.Box:getMin():getZ() < b3Part:getMin():getZ() + 500 * GEO.EPS_SMALL then
+         vtFacesAffected.Bottom = true
+      end
+      if Proc.Box:getMin():getY() < b3Part:getMin():getY() + 500 * GEO.EPS_SMALL then
+         vtFacesAffected.Front = true
+      end
+      if Proc.Box:getMax():getY() > b3Part:getMax():getY() - 500 * GEO.EPS_SMALL then
+         vtFacesAffected.Back = true
+      end
+      if Proc.Box:getMin():getX() < b3Part:getMin():getX() + 500 * GEO.EPS_SMALL then
+         vtFacesAffected.Left = true
+      end
+      if Proc.Box:getMax():getX() > b3Part:getMax():getX() - 500 * GEO.EPS_SMALL then
+         vtFacesAffected.Right = true
+      end
+   end
+   
+   return vtFacesAffected
+end
+
+-------------------------------------------------------------------------------------------------------------
+-- restituisce vero se la feature con box b3Proc taglia l'intera sezione della barra, rappresentata dalle sue dimensioni W e H
+function BeamLib.IsFeatureCuttingEntireSection( b3Proc, dRawW, dRawH)
+   return ((abs(b3Proc:getDimY() - dRawW) < 10 * GEO.EPS_SMALL or b3Proc:getDimY() > dRawW) and (abs(b3Proc:getDimZ() - dRawH) < 10 * GEO.EPS_SMALL or b3Proc:getDimZ() > dRawH))
+end
+
 -------------------------------------------------------------------------------------------------------------
 return BeamLib

LuaLibs/FacesBySaw.lua

@@ -1,4 +1,4 @@
--- FacesBySaw.lua by Egaltech s.r.l. 2023/06/14
+-- FacesBySaw.lua by Egaltech s.r.l. 2023/09/26
 -- Gestione taglio con lama di feature con una o due facce
 -- 2021/01/06 Cambiato limite per attacco Tg con lama e CalcLeadInOutGeom rinominata in CalcLeadInOutPerpGeom.
 -- 2021/02/03 In taglio lama si accettano anche due lati con deviazione minore di 20deg.
@@ -17,6 +17,7 @@
 -- 2023/04/20 Alcune modifiche per gestire tagli con faceuse parallelo.
 -- 2023/05/18 Imposto in ogni caso dVzLimDwnUp prima anche di bDownHead.
 -- 2023/06/14 Aggiunta MakeParallelOne e tolti parametri nForceWorkSide, nForceSCC da MakeOne.
+-- 2023/09/26 Piccola modifica per Turn su bInvert di MakeOne.
 
 -- Tabella per definizione modulo
 local FacesBySaw = {}
@@ -183,7 +184,7 @@ function FacesBySaw.MakeOne( nSurfId, nFacet, sCutting, dSawDiam, Par5, dVzLimDw
    local bInvert = bForceInvert
    -- se la lama ruota in senso antiorario inverto la direzione di lavorazione, per avere rotazione lama opposta a avanzamento
    if bInvert == nil then
-      if not BD.TURN then
+      if not BD.TURN or abs( ptP2:getY() - ptP1:getY()) < 250 then
          if bIsSawCCW ~= bDownUp then
             bInvert = ( ptP2:getZ() < ptP1:getZ() + 100 * GEO.EPS_SMALL)
          else

LuaLibs/FeatureTopology.lua

@@ -0,0 +1,244 @@
+-- FeatureTopology.lua by Egaltech s.r.l. 2023/09/12
+-- Libreria per classificazione topologica feature travi
+-- 2023/09/26 Aggiunte topologie Strip e Cut.
+-- 2023/09/27 Modificata GetFacesParallelToPart per tunnel, pocket e groove 3 lati
+
+-- Tabella per definizione modulo
+local FeatureTopology = {}
+
+-- Include
+require( 'EgtBase')
+
+-- Carico le librerie
+local BL = require( 'BeamLib')
+
+EgtOutLog( ' FeatureTopology started', 1)
+
+---------------------------------------------------------------------
+-- restituisce la matrice delle adiacenze di Proc dove i e j sono le facce e a(ij) è l'angolo tra di esse; 0 se nessuna adiacenza
+local function GetAdjacencyMatrix( Proc)
+   local vAdj = {}
+   for i = 1, Proc.Fct do
+      vAdj[i] = {}
+      for j = 1, Proc.Fct do
+         if i == j then
+            vAdj[i][j] = 0
+         else
+            _, _, _, vAdj[i][j] = EgtSurfTmFacetsContact( Proc.Id, i - 1, j - 1, GDB_ID.ROOT)
+            if not vAdj[i][j] then vAdj[i][j] = 0 end
+         end
+         j = j + 1
+      end
+      i = i + 1
+   end
+   return vAdj
+end
+
+---------------------------------------------------------------------
+-- restituisce gli id delle facce di Proc che hanno il numero di adiacenze nAdj
+function FeatureTopology.GetFacesWithGivenAdjacencyNumber( Proc, nAdj)
+   local vAdj = GetAdjacencyMatrix( Proc)
+   local vFacesWithGivenAdj = {}
+   for i = 1, Proc.Fct do
+      local nAdjCount = 0
+      for j = 1, Proc.Fct do
+         if vAdj[i][j] and vAdj[i][j] ~= 0 then
+            nAdjCount = nAdjCount + 1
+         end
+      end
+      if nAdjCount == nAdj then
+         table.insert( vFacesWithGivenAdj, i - 1)
+      end
+   end
+   return vFacesWithGivenAdj
+end
+
+---------------------------------------------------------------------
+-- restituisce true se Proc ha tutti gli angoli concavi (bAllConcave) e, nei casi in cui ha senso, se questi sono esattamente 90 deg (bAllRight)
+local function AreAllAnglesConcaveOrRight( Proc)
+   local vAdj = GetAdjacencyMatrix( Proc)
+   local bAllConcave, bAllRight = true, true
+   for i = 1, Proc.Fct do
+      for j = 1, Proc.Fct do
+         -- se trovo un angolo convesso restituisco falso e esco subito
+         if vAdj[i][j] and vAdj[i][j] > 0 then
+            bAllConcave = false
+            bAllRight = false
+            break
+         elseif vAdj[i][j] and vAdj[i][j] ~= 0 and vAdj[i][j] + 90 > 500 * GEO.EPS_ANG_SMALL then
+            bAllRight = false
+         end
+      end
+   end
+   -- se 1 faccia oppure 2 facce con angolo convesso non ha senso ritornare valori per bAllRight
+   if Proc.Fct < 2 or ( Proc.Fct == 2 and vAdj[1][2] > 0) then
+      return bAllConcave
+   else
+      return bAllConcave, bAllRight
+   end
+end
+
+---------------------------------------------------------------------
+-- restituisce true se almeno una delle dimensioni della feature è maggiore o uguale ad una delle dimensioni principali del pezzo (tolleranza 1 mm)
+local function IsAnyDimensionLongAsPart( Proc)
+   local bResult = false
+   local nBoxSolidId = EgtGetFirstNameInGroup( Proc.PartId or GDB_ID.NULL, 'Box')
+   local b3Solid = EgtGetBBoxGlob( nBoxSolidId, GDB_BB.STANDARD)
+   if Proc.Box:getDimX() > b3Solid:getDimX() - 1000 * GEO.EPS_SMALL or
+      Proc.Box:getDimY() > b3Solid:getDimY() - 1000 * GEO.EPS_SMALL or
+      Proc.Box:getDimZ() > b3Solid:getDimZ() - 1000 * GEO.EPS_SMALL then
+      bResult = true
+   end
+   return bResult
+end
+
+---------------------------------------------------------------------
+-- retituisce un vettore con gli indici (0 based) delle facce triangolari (o quasi) di Proc
+local function GetTriangularFaces( Proc)
+   local vTriangularFaces = {}
+   for i = 1, Proc.Fct do
+      if BL.Is3EdgesApprox( Proc, i - 1) then
+         table.insert( vTriangularFaces, i - 1)
+      end
+   end
+   return vTriangularFaces
+end
+
+---------------------------------------------------------------------
+-- restituisce un vettore contenente gli indici delle facce di Proc parallele ad una delle direzioni principali; il check varia in base alla famiglia topologica
+local function GetFacesParallelToPart( Proc, sFamily)
+   local vFacesParallelToPart = {}
+   for i = 0, Proc.Fct - 1 do
+      local vtN = EgtSurfTmFacetNormVersor( Proc.Id, i, GDB_ID.ROOT)
+      if sFamily == 'Bevel' or sFamily == 'DoubleBevel' or sFamily == 'Strip' or sFamily == 'Tunnel' or sFamily == 'Pocket' or ( sFamily == 'Groove' and Proc.IsThrough == true) then
+         local vTriangularFaces = GetTriangularFaces( Proc)
+         local bIsTriangularFace = false
+         -- verifico se la faccia è triangolare
+         for j = 1, #vTriangularFaces do
+            if i == vTriangularFaces[j] then
+               bIsTriangularFace = true
+            end
+         end
+         -- se faccia triangolare deve avere la normale parallela ad una direzione principale
+         if bIsTriangularFace then
+            if AreSameOrOppositeVectorApprox( vtN, X_AX()) or AreSameOrOppositeVectorApprox( vtN, Y_AX()) or AreSameOrOppositeVectorApprox( vtN, Z_AX()) then
+               table.insert( vFacesParallelToPart, i)
+            end
+         -- altrimenti deve avere una componente della normale nulla
+         else
+            if abs( vtN:getX()) < 10 * GEO.EPS_SMALL or abs( vtN:getY()) < 10 * GEO.EPS_SMALL or abs( vtN:getZ()) < 10 * GEO.EPS_SMALL then
+               table.insert( vFacesParallelToPart, i)
+            end
+         end
+      else
+         -- la normale deve essere parallela ad una direzione principale
+         if AreSameOrOppositeVectorApprox( vtN, X_AX()) or AreSameOrOppositeVectorApprox( vtN, Y_AX()) or AreSameOrOppositeVectorApprox( vtN, Z_AX()) then
+            table.insert( vFacesParallelToPart, i)
+         end
+      end
+   end
+   return vFacesParallelToPart
+end
+
+---------------------------------------------------------------------
+-- restituisce una stringa con il nome esteso della topologia della feature
+-- *famiglia-passante-angoli tutti concavi a 90deg-facce tutte parallele alle dimensioni principali-numero di facce*
+local function GetTopologyLongName( sFamily, bIsThrough, bAllRightAngles, bIsParallel, nNumberOfFaces)
+   -- feature passante o cieca
+   local sThrough = '_'
+   if bIsThrough ~= nil then sThrough = EgtIf( bIsThrough, 'Through', 'Blind') end
+   -- tutti gli angoli della feature sono retti oppure no
+   local sAllRightAngles = '_'
+   if bAllRightAngles ~= nil then sAllRightAngles = EgtIf( bAllRightAngles, 'RightAngles', 'NotRightAngles') end
+   -- tutte le dimensioni della feature sono parallele agli assi principali del pezzo oppure no
+   local sParallel = '_'
+   if bIsParallel ~= nil then sParallel = EgtIf( bIsParallel, 'Parallel', 'NotParallel') end
+
+   local sLongName = sFamily .. '-' .. sThrough .. '-' .. sAllRightAngles .. '-' .. sParallel .. '-' .. nNumberOfFaces
+   return sLongName
+end
+
+---------------------------------------------------------------------
+-- riconosce se Proc è una delle topologie standard e, in caso positivo, ne scrive le caratteristiche in campi specifici della Proc stessa restituendo true
+function FeatureTopology.Classify( Proc, b3Raw)
+   local bRecognized = false
+   local sFamily
+   local bIsThrough
+   local bAllRightAngles
+   local bIsParallel
+   local sLongName = ''
+
+   -- SE NON HA TUTTE LE FACCE PIANE RITORNARE NIL!!
+
+   local bAllAnglesConcave
+   bAllAnglesConcave, bAllRightAngles = AreAllAnglesConcaveOrRight( Proc)
+   local vTriangularFaces = GetTriangularFaces( Proc)
+   local bIsAnyDimensionLongAsPart = IsAnyDimensionLongAsPart( Proc)
+   local vFacesWithOneAdj = FeatureTopology.GetFacesWithGivenAdjacencyNumber( Proc, 1)
+   local vFacesWithTwoAdj = FeatureTopology.GetFacesWithGivenAdjacencyNumber( Proc, 2)
+   local vFacesWithThreeAdj = FeatureTopology.GetFacesWithGivenAdjacencyNumber( Proc, 3)
+   local vFacesWithFourAdj = FeatureTopology.GetFacesWithGivenAdjacencyNumber( Proc, 4)
+   local dRawW, dRawH = b3Raw:getDimY(), b3Raw:getDimZ()
+   local bIsFeatureCuttingEntireSection = BL.IsFeatureCuttingEntireSection( Proc.Box, dRawW, dRawH)
+
+   if Proc.IsOutline then
+      sFamily = 'OUTLINE'
+   elseif Proc.Prc == 40 then
+      sFamily = 'DRILLING'
+   elseif Proc.Fct == 1 and bIsAnyDimensionLongAsPart and bIsFeatureCuttingEntireSection then
+      sFamily = 'Cut'
+   elseif Proc.Fct == 1 and bIsAnyDimensionLongAsPart then
+      sFamily = 'Bevel'
+      bIsThrough = true
+   elseif Proc.Fct == 2 and bAllAnglesConcave and #vTriangularFaces == 1 then
+      sFamily = 'Bevel'
+      bIsThrough = false
+   elseif Proc.Fct == 2 and bAllAnglesConcave then
+      sFamily = 'Rabbet'
+      bIsThrough = true
+   elseif Proc.Fct == 2 and not bAllAnglesConcave and bIsAnyDimensionLongAsPart then
+      sFamily = 'DoubleBevel'
+      bIsThrough = true
+   elseif Proc.Fct == 3 and bAllAnglesConcave and #vFacesWithTwoAdj == 1 and #vTriangularFaces == 2 then
+      sFamily = 'Bevel'
+      bIsThrough = false
+   elseif Proc.Fct == 3 and bAllAnglesConcave and #vFacesWithTwoAdj == 1 and bIsAnyDimensionLongAsPart then
+      sFamily = 'Groove'
+      bIsThrough = true
+   elseif Proc.Fct == 3 and bAllAnglesConcave and #vFacesWithTwoAdj == 3 then
+      sFamily = 'Groove'
+      bIsThrough = false
+   elseif Proc.Fct == 4 and bAllAnglesConcave and #vFacesWithThreeAdj == 2 then
+      sFamily = 'Groove'
+      bIsThrough = false
+   elseif Proc.Fct == 4 and bAllAnglesConcave and #vFacesWithTwoAdj == 4 and bIsAnyDimensionLongAsPart then
+      sFamily = 'Tunnel'
+      bIsThrough = true
+   elseif Proc.Fct >= 4 and #vFacesWithOneAdj == 2 and bIsAnyDimensionLongAsPart then
+      sFamily = 'Strip'
+      bIsThrough = true
+   elseif Proc.Fct == 5 and bAllAnglesConcave and #vFacesWithFourAdj == 1 then
+      sFamily = 'Pocket'
+      bIsThrough = false
+   end
+   local vFacesParallelToPart = GetFacesParallelToPart( Proc, sFamily)
+   bIsParallel = ( #vFacesParallelToPart == Proc.Fct)
+
+   if sFamily == 'OUTLINE' or sFamily == 'DRILLING' then
+      Proc.Topology = sFamily
+      Proc.TopologyLongName = sFamily
+      bRecognized = true
+   elseif sFamily then
+      sLongName = GetTopologyLongName( sFamily, bIsThrough, bAllRightAngles, bIsParallel, Proc.Fct)
+      Proc.Topology, Proc.IsThrough, Proc.AllRightAngles, Proc.IsParallel, Proc.TopologyLongName = sFamily, bIsThrough, bAllRightAngles, bIsParallel, sLongName
+      bRecognized = true
+   else
+      Proc.Topology = 'OTHER'
+      Proc.TopologyLongName = 'OTHER'
+   end
+
+   return bRecognized
+end
+
+-------------------------------------------------------------------------------------------------------------
+return FeatureTopology

LuaLibs/ProcessCut.lua

@@ -21,6 +21,7 @@
 -- 2023/06/19 Aggiunti tagli speciali per evitare il rischio che il cubetto rimanga appoggiato al motore.
 -- 2023/08/01 Correzione su offset per taglio doppio di lato.
 -- 2023/08/01 In caso di pezzi alti su macchina tipo PF i tagli doppi di lato non vengono usati, si usa il metodo standard.
+-- 2023/09/12 In MakeFromTop abbassato a 590 mm il limite per convertire in LongCut.
 
 -- Tabella per definizione modulo
 local ProcessCut = {}
@@ -249,14 +250,15 @@ local function MakeFromTop( sCutting, Proc, nPhase, nRawId, nPartId, dOvmHead, b
    if bFromBottom == nil then
       bFromBottom = ( vtN:getZ() > 0.25 and b3Solid:getDimX() <  BD.LEN_SHORT_PART and not Proc.AdvTail and vtN:getX() < 0 and abs( vtN:getY()) < 0.259)
    end
+   local dMinFeatureLengthForLongCut = 590.000
    -- verifico se da considerare taglio lungo ( non da sotto, inclinato non più di 30deg, largo come la trave e abbastanza lungo)
    local bLongCut = ( not bDownCut and vtN:getZ() > 0.865 and
                       Proc.Box:getDimY() > b3Solid:getDimY() - 10 * GEO.EPS_SMALL and
-                      ( Proc.Box:getDimX() > 0.75 * b3Solid:getDimX() or Proc.Box:getDimX() > 600.000))
+                      ( Proc.Box:getDimX() > 0.75 * b3Solid:getDimX() or Proc.Box:getDimX() > dMinFeatureLengthForLongCut))
    -- verifico se da considerare taglio lungo dal lato, solo per macchine tipo PF1250, inclinato non più di 30deg
    local bLongCutFromSide = ( not bDownCut and ( BD.C_SIMM and BD.DOWN_HEAD and ( abs(vtN:getY()) > 0.865) and
                               Proc.Box:getDimZ() > b3Solid:getDimZ() - 10 * GEO.EPS_SMALL) and
-                              ( Proc.Box:getDimX() > 0.75 * b3Solid:getDimX() or Proc.Box:getDimX() > 600.000))
+                              ( Proc.Box:getDimX() > 0.75 * b3Solid:getDimX() or Proc.Box:getDimX() > dMinFeatureLengthForLongCut))
    -- se taglio lungo e Q04 = 1 allora lancio il processo dell'L10
    local bNoDicing = EgtGetInfo( Proc.Id, 'Q04', 'i') == 1
    if bNoDicing then

LuaLibs/ProcessDovetail.lua

@@ -1,6 +1,7 @@
 -- ProcessDovetail.lua by Egaltech s.r.l. 2021/02/04
 -- Gestione calcolo giunzione coda di rondine
 -- 2022/06/10 Aggiunto il parametro dOvmTail per gestire sovramateriali in coda diversi da OVM_MID (sezioni alte e larghe)
+-- 2023/09/26 Modificata chiamata a GetFaceWithMostAdj.
 
 -- Tabella per definizione modulo
 local ProcessDovetail = {}
@@ -280,7 +281,7 @@ function ProcessDovetail.IsHeadFeature( Proc, b3Raw, dCurrOvmH)
    end
    -- deve avere la normale principale diretta verso la testa
    local nPartId = EgtGetParent( EgtGetParent( Proc.Id) or GDB_ID.NULL)
-   local nFacInd, dElev, nFacInd2, dElev2 = BL.GetFaceWithMostAdj( Proc.Id, nPartId)
+   local nFacInd, dElev, nFacInd2, dElev2 = BL.GetFaceWithMostAdj( Proc, nPartId)
    local ptC, vtN = EgtSurfTmFacetCenter( Proc.Id, nFacInd, GDB_ID.ROOT)
    if vtN:getZ() < BD.NZ_MINA and nFacInd2 then
       ptC, vtN = EgtSurfTmFacetCenter( Proc.Id, nFacInd2, GDB_ID.ROOT)
@@ -341,7 +342,7 @@ function ProcessDovetail.IsTailFeature( Proc, b3Raw)
       end
    end
    -- deve avere la normale principale diretta verso la coda
-   local nFacInd, dElev, nFacInd2, dElev2 = BL.GetFaceWithMostAdj( Proc.Id, nPartId)
+   local nFacInd, dElev, nFacInd2, dElev2 = BL.GetFaceWithMostAdj( Proc, nPartId)
    local ptC, vtN = EgtSurfTmFacetCenter( Proc.Id, nFacInd, GDB_ID.ROOT)
    if vtN:getZ() < BD.NZ_MINA and nFacInd2 then
       ptC, vtN = EgtSurfTmFacetCenter( Proc.Id, nFacInd2, GDB_ID.ROOT)
@@ -933,10 +934,10 @@ local function MachSideFaces( Proc, nPartId, b3Raw, nFacetCnt, bForceUseRough,
          end
       else
          -- recupero la faccia con il maggior numero di adiacenze e l'elevazione relativa
-         nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc.Id, nPartId, false, sin(dSideAngle))
+         nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc, nPartId, false, sin(dSideAngle))
          if not nFacInd or nFacInd < 0 then
             -- provo eliminando i sottosquadra
-            nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc.Id, nPartId, false, -2)
+            nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc, nPartId, false, -2)
             if not nFacInd or nFacInd < 0 then
                local sErr = 'Error on process ' .. tostring( Proc.Id) .. ' MachSideFaces could not find reference face'
                EgtOutLog( sErr)

LuaLibs/ProcessDrill.lua

@@ -1,4 +1,4 @@
--- ProcessDrill.lua by Egaltech s.r.l. 2023/05/11
+-- ProcessDrill.lua by Egaltech s.r.l. 2023/09/26
 -- Gestione calcolo forature per Travi
 -- 2021/04/14 DS Corretta gestione Invert con testa da sotto. Migliorato calcolo ingombro portautensile.
 -- 2021/05/03 DS Su macchina con testa da sotto aggiunta possibilità di fare fori verticali metà sopra e l'altra sotto.
@@ -18,6 +18,7 @@
 -- 2022/12/28 Implementata gestione forature in doppio
 -- 2023/05/11 Se lavorazione in doppio e precedente no oppure di tipo diverso, forzo risalita a Zmax.
 -- 2023/07/28 Aggiunta gestione rinvio 90deg solo per forature da sotto esattamente verticali.
+-- 2023/09/26 Se errore in applicazione lavorazione si inverte e riprova solo se foratura singola su foro aperto.
 
 -- Tabella per definizione modulo
 local ProcessDrill = {}
@@ -166,7 +167,7 @@ function ProcessDrill.Classify( Proc, b3Raw)
       return false, false, false
    end
    local bOpen = ( Proc.Fce ~= 0)
-   local bFaceDown = ( ptCen:getZ() < b3Raw:getMin():getZ() + 2  and  (( not Proc.Head  and  not Proc.Tail) or (( BD.ROT90 or BD.DOWN_HEAD or BD.TURN) and ( Proc.Flg == 2 or Proc.Flg == -2))))
+   local bFaceDown = ( ptCen:getZ() < b3Raw:getMin():getZ() + 2  and  (( not Proc.Head  and  not Proc.Tail) or (( BD.ROT90 or BD.DOWN_HEAD or BD.TURN) and abs( Proc.Flg) == 2)))
    -- verifico se il foro è sotto e quindi va spostato o sopra o sul fianco
    if (( vtExtr:getZ() < BD.DRILL_VZ_MIN or bFaceDown) and ( not bOpen or Proc.Flg ~= 1)) then
       if not bOpen then
@@ -249,7 +250,7 @@ function ProcessDrill.Make( Proc, nPhase, nRawId, nPartId)
    -- verifico possibilità di foratura da sotto
    local bDownDrill = ( BD.TURN or ( BD.DOWN_HEAD and ML.FindDrilling( dDiam, nil, false, true)))
    -- se presente rinvio angolare
-   local bDrillAngTrasm = ( BD.ANG_TRASM and AreOppositeVectorApprox( vtExtr, Z_AX()))
+   local bDrillAngTrasm = ( BD.ANG_TRASM and AreOppositeVectorApprox( vtExtr, Z_AX()) and not bOpen)
    -- verifico che il foro non sia fattibile solo da sotto
    local bToInvert = ( vtExtr:getZ() < BD.DRILL_VZ_MIN and not bDownDrill) and not bDrillAngTrasm
    if bToInvert  and  ( not bOpen or Proc.Flg ~= 1) then
@@ -516,7 +517,7 @@ function ProcessDrill.Make( Proc, nPhase, nRawId, nPartId)
    EgtSetMachiningParam( MCH_MP.USERNOTES, sUserNotes)
    -- eseguo
    local bOk = ML.ApplyMachining( true, false)
-   if not bOk and bDownDrill then
+   if not bOk and bDownDrill and bOpen and abs( Proc.Flg) == 1 then
       if sType == 'Drill' or sType == 'Drill_H2' or sType == 'AngleDrill' then
          EgtSetMachiningParam( MCH_MP.INVERT, true)
       else

LuaLibs/ProcessLapJoint.lua

@@ -1,4 +1,4 @@
--- ProcessLapJoint.lua by Egaltech s.r.l. 2023/07/19
+-- ProcessLapJoint.lua by Egaltech s.r.l. 2023/09/26
 -- Gestione calcolo mezzo-legno per Travi
 -- 2019/10/08 Agg. gestione OpenPocket.
 -- 2021/01/24 Con sega a catena ora sempre impostato asse A.
@@ -72,6 +72,13 @@
 -- 2023/08/01 Aggiunta gestione rinvio angolare solo per tasche 5 lati, da sotto e perfettamente verticali.
 -- 2023/08/02 Piccola correzione al calcolo della distanza di sicurezza in base al gambo dell'utensile.
 -- 2023/08/02 Corretto incremento della profondità pari al raggio utensile per tunnel con sega a catena.
+-- 2023/08/10 Modificata scelta SCC per tasche in Y+/- in coda o quasi (anche per relativi antischeggia con lama).
+-- 2023/09/12 In MakeByPockets gestito correttamente il ritorno nil di GetUShapeWidth.
+-- 2023/09/26 In MakeLongMoreFaces migliorata gestione spezzatura di Tunnel.
+-- 2023/09/26 Modificata chiamata a GetFaceWithMostAdj.
+-- 2023/09/27 In MakeMoreFaces e Classify refactoring della gestione scanalatura chiusa.
+-- 2023/09/27 In MakeLongMoreFaces aggiunto Topology.Classify alla Proc creata.
+-- 2023/10/03 In MakePocket aggiunto messaggio in caso si rovini il pezzo successivo.
 
 -- Tabella per definizione modulo
 local ProcessLapJoint = {}
@@ -84,6 +91,7 @@ local Cut = require( 'ProcessCut')
 local DoubleCut = require( 'ProcessDoubleCut')
 local LongCut = require( 'ProcessLongCut')
 local Long2Cut = require( 'ProcessLongDoubleCut')
+local Topology = require( 'FeatureTopology')
 
 EgtOutLog( ' ProcessLapJoint started', 1)
 
@@ -632,7 +640,8 @@ function ProcessLapJoint.IsHeadFeature( Proc, b3Raw, dCurrOvmH)
       end
    end
    -- deve avere la normale principale diretta verso la testa
-   local nFacInd, dElev = BL.GetFaceWithMostAdj( Proc.Id, nPartId)
+   Topology.Classify( Proc, b3Raw)
+   local nFacInd, dElev = BL.GetFaceWithMostAdj( Proc, nPartId)
    local ptC, vtN = EgtSurfTmFacetCenter( Proc.Id, nFacInd, GDB_ID.ROOT)
    if vtN and vtN:getX() < 0.499 then
       return false
@@ -686,7 +695,8 @@ function ProcessLapJoint.IsTailFeature( Proc, b3Raw)
       end
    end
    -- deve avere la normale principale diretta verso la coda (oppure tunnel)
-   local nFacInd, dElev, nFacInd2, dElev2 = BL.GetFaceWithMostAdj( Proc.Id, nPartId)
+   Topology.Classify( Proc, b3Raw)
+   local nFacInd, dElev, nFacInd2, dElev2 = BL.GetFaceWithMostAdj( Proc, nPartId)
    local ptC, vtN = EgtSurfTmFacetCenter( Proc.Id, nFacInd, GDB_ID.ROOT)
    if vtN and vtN:getZ() < BD.NZ_MINA and nFacInd2 then
       ptC, vtN = EgtSurfTmFacetCenter( Proc.Id, nFacInd2, GDB_ID.ROOT)
@@ -735,7 +745,7 @@ function ProcessLapJoint.Classify( Proc, b3Raw)
       local nDeletedFace
       -- recupero la faccia con il maggior numero di adiacenze e minor elevazione
       local nPartId = EgtGetParent( EgtGetParent( Proc.Id) or GDB_ID.NULL)
-      local nFacInd, dElev, nFacInd2, dElev2 = BL.GetFaceWithMostAdj( Proc.Id, nPartId)
+      local nFacInd, dElev, nFacInd2, dElev2 = BL.GetFaceWithMostAdj( Proc, nPartId)
       if not nFacInd or nFacInd < 0 then
          if nFacInd == -1 then
             bClosedOrthoFaces = nFacInd2
@@ -750,45 +760,11 @@ function ProcessLapJoint.Classify( Proc, b3Raw)
             return true, false
          end
       end
-      -- se è una feature scanalatura (con 5 facce) e non è stata riconosciuta come fessura, eseguo altre verifiche
-      if Proc.Prc == 16 and Proc.Fct == 5 and not bClosedOrthoFaces then
-         -- recupero gruppo per geometria addizionale
-         local nAddGrpId = BL.GetAddGroup( nPartId)
-         if not nAddGrpId then
-            return false
-         end
-         -- dalla copia della superfice, ciclo eliminando una faccia per volta per verificare se trova fessura
-         for i = 1, Proc.Fct do
-            local nNewProc = EgtCopyGlob( Proc.Id, nAddGrpId) or GDB_ID.NULL
-            -- elimino una faccia
-            nDeletedFace = i - 1
-            if EgtSurfTmRemoveFacet( nNewProc, nDeletedFace) then
-               -- recupero la faccia con il maggior numero di adiacenze e l'elevazione relativa
-               nFacInd, dElev, nFacInd2, dElev2 = BL.GetFaceWithMostAdj( nNewProc, nPartId)
-               if not nFacInd or nFacInd < 0 then
-                  if nFacInd == -1 then
-                     bClosedOrthoFaces = nFacInd2
-                     EgtErase( nNewProc)
-                     break
-                  else
-                     EgtErase( nNewProc)
-                     return false
-                  end
-               end
-            -- altrimenti esco
-            else
-               EgtErase( nNewProc)
-               break
-            end
-         end
-         -- se riconosciuta fessura ricalcolo l'elevazione dalla faccia di fondo
-         if bClosedOrthoFaces then
-            nFacInd = nDeletedFace
-            -- rendo nulla la faccia opzionale perchè si tratta di una fessura
-            nFacInd2 = nil
-            dElev = BL.GetFaceElevation( Proc.Id, nFacInd)
-            bClosedOrthoFaces = false -- non lo setto come tunnel
-         end
+      -- se scanalatura chiusa lavoro la faccia di fondo
+      if Proc.Topology == 'Pocket' and ( Proc.IsParallel or Proc.AllRightAngles) and not bClosedOrthoFaces then
+         nFacInd = Topology.GetFacesWithGivenAdjacencyNumber( Proc, 4)[1]
+         nFacInd2 = nil
+         dElev = Proc.Face[ nFacInd + 1].Elevation
       end
       -- se facce formano un tunnel e sono ortogonali
       if bClosedOrthoFaces then
@@ -2679,7 +2655,7 @@ local function MakeByChainOrSaw( Proc, nPhase, nRawId, nPartId, nFacInd,
    local bIs3Faces = ( Proc.Fct == 3)
    if bIs3Faces then
       -- recupero la faccia con il maggior numero di adiacenze e l'elevazione relativa
-      local nFacInd1, dFacElev1, nFacInd2, dFacElev2 = BL.GetFaceWithMostAdj( Proc.Id, nPartId, bIs3Faces)
+      local nFacInd1, dFacElev1, nFacInd2, dFacElev2 = BL.GetFaceWithMostAdj( Proc, nPartId, bIs3Faces)
       if not nFacInd1 or nFacInd1 < 0 then
          if nFacInd1 == -1 then
             bOrthoFaces = nFacInd2
@@ -3102,8 +3078,9 @@ local function MakeAntiSplintBySaw( Proc, nFacet, vtN, b3Raw, nFacInd, bReduceDe
 end
 
 ---------------------------------------------------------------------
-local function MakePocket( Proc, nPartId, b3Solid, ptPs, tvtN, nFaceRef, sMchFind, nUseRoughTool, sMasterPocket, dPrevFaceElev, tDimAndRef, dAng, bOpenOutRaw, bLapJointAngTrasm)
-
+local function MakePocket( Proc, nPartId, b3Solid, ptPs, tvtN, nFaceRef, sMchFind, nUseRoughTool, sMasterPocket, dPrevFaceElev, tDimAndRef, dAng, bOpenOutRaw, bLapJointAngTrasm, nPhase, nRawId)
+   -- distanza dal pezzo successivo
+   local dDistToNextPiece = BL.GetDistanceToNextPart( nRawId, nPhase)
    -- calcolo l'elevazione dal punto medio
    local dElev
    local dLenIn, dLedOut = BL.GetPointDirDepth( nPartId, ptPs, tvtN[2])
@@ -3173,7 +3150,10 @@ local function MakePocket( Proc, nPartId, b3Solid, ptPs, tvtN, nFaceRef, sMchFin
       if AreSameVectorApprox( tvtN[2], Z_AX()) then
          nSCC = MCH_SCC.ADIR_YM
       elseif abs( tvtN[2]:getX()) < 0.1 then
-         nSCC = EgtIf( BL.IsPartFinalPhase( EgtGetCurrPhase()), MCH_SCC.ADIR_XM, MCH_SCC.ADIR_XP)
+         local bNearTail = ( Proc.Box:getMax():getX() < b3Solid:getCenter():getX() + 50 and
+                             Proc.Box:getMax():getX() - b3Solid:getMin():getX() < 1000)
+         local bVeryShortPart = ( BD.LEN_VERY_SHORT_PART and b3Solid:getDimX() < BD.LEN_VERY_SHORT_PART)
+         nSCC = EgtIf( BL.IsPartFinalPhase( EgtGetCurrPhase()) or ( bNearTail and not bVeryShortPart), MCH_SCC.ADIR_XM, MCH_SCC.ADIR_XP)
       elseif tvtN[2]:getY() > 0.1 then
          nSCC = MCH_SCC.ADIR_YP
       end
@@ -3223,6 +3203,11 @@ local function MakePocket( Proc, nPartId, b3Solid, ptPs, tvtN, nFaceRef, sMchFin
          return false, sErr
       end
    end
+   -- se la lavorazione rovina il pezzo successivo emetto warning
+   if ( not Proc.Tail or Proc.AdvTail) and Proc.AffectedFaces.Left and dDistToNextPiece < dDiamTool / 2 then
+      sWarn = 'Warning on pocketing ' .. sName .. ': machining can damage next piece. Minimum distance needed : ' .. ( 1 + dDiamTool / 2) .. ' mm'
+      EgtOutLog( sWarn)
+   end
    if BD.TURN then
       -- centro del pezzo
       local ptCen = ORIG()
@@ -3266,7 +3251,7 @@ local function MachineByMill( Proc, nPhase, nRawId, nPartId, b3Solid, tvtN, nBas
       -- se l'angolo dalla verticale si discosta di più dell'angolo limite impostato, utilizzo la svuotatura
       if cos( dDiffFromSqAng) < cos( dAngLimit) then
          -- applico la svuotatura
-         local bOk, sWarn, sTuuidPk, dDiamTool, dElev = MakePocket( Proc, nPartId, b3Solid, ptPs, tvtN, nSideFace, sMchFind, nUseRoughTool, sPocketing, dPrevFaceElev, tDimAndRef, dAng)
+         local bOk, sWarn, sTuuidPk, dDiamTool, dElev = MakePocket( Proc, nPartId, b3Solid, ptPs, tvtN, nSideFace, sMchFind, nUseRoughTool, sPocketing, dPrevFaceElev, tDimAndRef, dAng, nil, nil, nPhase, nRawId)
          if not bOk then
             -- se ho id utensile e diametro è perchè non ha fatto svuotatura perchè la faccia è più stretta del diametro utensile
             -- e provo ad inserire singola passata di testa
@@ -3475,7 +3460,7 @@ local function MakeByPockets( Proc, nPhase, nRawId, nPartId, nChamfer, dDepthCha
    end
    if bIs3Faces then
       -- recupero la faccia con il maggior numero di adiacenze e l'elevazione relativa
-      local nFacInd, dFacElev, nFacInd2, dFacElev2 = BL.GetFaceWithMostAdj( Proc.Id, nPartId, bIs3Faces)
+      local nFacInd, dFacElev, nFacInd2, dFacElev2 = BL.GetFaceWithMostAdj( Proc, nPartId, bIs3Faces)
       if not nFacInd or nFacInd < 0 then
          if nFacInd == -1 then
             bOrthoFaces = nFacInd2
@@ -3498,8 +3483,8 @@ local function MakeByPockets( Proc, nPhase, nRawId, nPartId, nChamfer, dDepthCha
          dDimMin = 2 * dDimMin
       -- se tre facce aumento il diametro di quanto possibile
       elseif Proc.Fct == 3 then
-         local nFacInd = BL.GetFaceWithMostAdj( Proc.Id, nPartId, false)
-         local dWidth = GetUShapeWidth( Proc, nFacInd or -1)
+         local nFacInd = BL.GetFaceWithMostAdj( Proc, nPartId, false)
+         local dWidth = GetUShapeWidth( Proc, nFacInd or -1) or dDimMin
          if dWidth > dDimMin then
             dDimMin = min( 2 * dDimMin, dWidth)
          end
@@ -3525,7 +3510,7 @@ local function MakeByPockets( Proc, nPhase, nRawId, nPartId, nChamfer, dDepthCha
                EgtInvertSurf( nSurfInt)
             end
          end
-         local nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc.Id, nPartId)
+         local nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc, nPartId)
          local vtN = EgtSurfTmFacetNormVersor( Proc.Id, nFacInd, GDB_ID.ROOT)
          local bMakeContour = false
          local sMilling
@@ -3872,7 +3857,7 @@ local function MakeByPockets( Proc, nPhase, nRawId, nPartId, nChamfer, dDepthCha
             -- se non completo e U, cerco di lavorare anche la faccia di fondo (con il massimo affondamento possibile)
             if not bComplete and bIs3Faces then
                -- recupero la faccia da lavorare
-               local nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc.Id, nPartId)
+               local nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc, nPartId)
                local vtN
                local bPock3rd = false
                if nFacInd then
@@ -3949,7 +3934,7 @@ local function MakeByPockets( Proc, nPhase, nRawId, nPartId, nChamfer, dDepthCha
 end
 
 ---------------------------------------------------------------------
-local function ManageAntiSplintBySaw( Proc, b3Raw, bIsU, vtN, nFacInd, sWarn, bMillDown, bReduceDepth)
+local function ManageAntiSplintBySaw( Proc, b3Raw, b3Solid, bIsU, vtN, nFacInd, sWarn, bMillDown, bReduceDepth)
 
    local bMadeASbyBld = false
    local nNumFac = EgtIf( bIsU, 2, 1)
@@ -3989,9 +3974,18 @@ local function ManageAntiSplintBySaw( Proc, b3Raw, bIsU, vtN, nFacInd, sWarn, bM
                if abs( vtN:getZ()) > 0.63 or abs( vtN:getY()) > 0.63 then
                   -- se superficie principale parallela al piano XZ
                   if nPrefSide == 0 then
-                     -- se facce praticamente verticali || mandrino a destra
+                     -- se facce praticamente verticali ||
                      if abs( vtNFc:getZ()) < 0.1 then
-                        nPrefSide = 2
+                        -- se FAST, mandrino a destra o sinistra come svuotatura
+                        if not BD.C_SIMM and not BD.TURN then
+                           local bNearTail = ( Proc.Box:getMax():getX() < b3Solid:getCenter():getX() + 50 and
+                                             Proc.Box:getMax():getX() - b3Solid:getMin():getX() < 1000)
+                           local bVeryShortPart = ( BD.LEN_VERY_SHORT_PART and b3Solid:getDimX() < BD.LEN_VERY_SHORT_PART)
+                           nPrefSide = EgtIf( BL.IsPartFinalPhase( EgtGetCurrPhase()) or ( bNearTail and not bVeryShortPart), 1, 2)
+                        -- altrimenti sempre a destra
+                        else
+                           nPrefSide = 2
+                        end
                      -- se facce inclinate \\ allora mandrino a destra (per essere verso l'alto)
                      elseif vtNFc:getX() * vtNFc:getZ() > 0 then
                         nPrefSide = 2
@@ -4488,7 +4482,7 @@ local function MakeMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead, bSinglePa
    local nBottomFace
    local sMchFindBackUp
    -- recupero la faccia con il maggior numero di adiacenze e l'elevazione relativa
-   nFacInd, dFacElev, nFacInd2, dFacElev2 = BL.GetFaceWithMostAdj( Proc.Id, nPartId)
+   nFacInd, dFacElev, nFacInd2, dFacElev2 = BL.GetFaceWithMostAdj( Proc, nPartId)
    if not nFacInd or nFacInd < 0 then
       if nFacInd == -1 then
          bClosedOrthoFaces = nFacInd2
@@ -4498,42 +4492,11 @@ local function MakeMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead, bSinglePa
          return false, sErr
       end
    end
-   -- se è una feature scanalatura (con 5 facce) e non è stata riconosciuta come fessura, eseguo altre verifiche
-   if Proc.Prc == 16 and Proc.Fct == 5 and not bClosedOrthoFaces then
-      -- dalla copia della superfice, ciclo eliminando una faccia per volta per verificare se trova fessura
-      for i = 1, Proc.Fct do
-         local nNewProc = EgtCopyGlob( Proc.Id, nAddGrpId) or GDB_ID.NULL
-         -- elimino una faccia
-         nBottomFace = i - 1
-         if EgtSurfTmRemoveFacet( nNewProc, nBottomFace) then
-            -- recupero la faccia con il maggior numero di adiacenze e l'elevazione relativa
-            nFacInd, dFacElev, nFacInd2, dFacElev2 = BL.GetFaceWithMostAdj( nNewProc, nPartId)
-            if not nFacInd or nFacInd < 0 then
-               if nFacInd == -1 then
-                  bClosedOrthoFaces = nFacInd2
-                  EgtErase( nNewProc)
-                  break
-               else
-                  EgtErase( nNewProc)
-                  local sErr = 'Error : MakeMoreFaces could not find reference face'
-                  EgtOutLog( sErr)
-                  return false, sErr
-               end
-            end
-         -- altrimenti esco
-         else
-            EgtErase( nNewProc)
-            break
-         end
-      end
-      -- se riconosciuta fessura ricalcolo l'elevazione dalla faccia di fondo
-      if bClosedOrthoFaces then
-         nFacInd = nBottomFace
-         -- rendo nulla la faccia opzionale perchè si tratta di una fessura
-         nFacInd2 = nil
-         dFacElev = BL.GetFaceElevation( Proc.Id, nFacInd)
-         bClosedOrthoFaces = false -- non setto come tunnel
-      end
+   -- se scanalatura chiusa lavoro la faccia di fondo
+   if Proc.Topology == 'Pocket' and ( Proc.IsParallel or Proc.AllRightAngles) and not bClosedOrthoFaces then
+      nFacInd = Topology.GetFacesWithGivenAdjacencyNumber( Proc, 4)[1]
+      nFacInd2 = nil
+      dFacElev = Proc.Face[ nFacInd + 1].Elevation
    end
    -- se proviene da divisione in parti lungo X ed è un tunnel non è fattibile
    if not bSinglePart and bClosedOrthoFaces then
@@ -5147,12 +5110,12 @@ local function MakeMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead, bSinglePa
                      local bMadeASbyBld = false
                      local bPassThrou = ( Proc.Box:getDimY() > b3Raw:getDimY() - 1 or Proc.Box:getDimZ() > b3Raw:getDimZ() - 1)
                      local bPassEdge = ((( bIsU and Proc.Fct == 3) or ( bIsL and Proc.Fct == 2)) and bSinglePart and Proc.Box:getDimX() < 0.9 * b3Raw:getDimX())
-                     local nFacIndOri = BL.GetFaceWithMostAdj( Proc.Id, nPartId)
+                     local nFacIndOri = BL.GetFaceWithMostAdj( Proc, nPartId)
                      local vtNOri = EgtSurfTmFacetNormVersor( Proc.Id, nFacIndOri, GDB_ID.ROOT)
                      if nChamfer < 2 and nQAntisplintResult == 1 and (( bIsU or bIsL) and ( bPassThrou or bPassEdge)) and ( vtNOri:getZ() > -0.087 or BD.DOWN_HEAD or BD.TURN) then
                         local bOk
                         local bSawDown = ( bMillDown and not bMillUp)
-                        bMadeASbyBld, bOk, sWarn = ManageAntiSplintBySaw( Proc, b3Raw, bIsU, vtNOri, nFacIndOri, sWarn, bSawDown, true)
+                        bMadeASbyBld, bOk, sWarn = ManageAntiSplintBySaw( Proc, b3Raw, b3Solid, bIsU, vtNOri, nFacIndOri, sWarn, bSawDown, true)
                         if not bOk then return false, sWarn end
                      end
                      -- in base al flag interno e al numero di facce e se ha forma ad U
@@ -5175,7 +5138,7 @@ local function MakeMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead, bSinglePa
                         -- 3 facce non a L ma con una faccia favorevole a Y, oppure
                         -- 4 facce
                         -- lancio la MakePocket ( lavorazione solo da un lato, fondo della tasca la faccia più favorevole a Y)
-                        if bForceSideMill and (( Proc.Fct == 3 and bIsL) or ( ( Proc.Fct == 3 or Proc.Fct == 2) and abs( vtNFacApproxY:getY()) >= 0.707 ) or Proc.Fct == 4) then
+                        if bForceSideMill and (( Proc.Fct == 3 and bIsL) or ( ( Proc.Fct == 3 or Proc.Fct == 2) and abs( vtNFacApproxY:getY()) >= 0.707 ) or ( Proc.Fct == 4 and Proc.Topology == 'Groove')) then
                            nFacInd, vtN, ptC = nFacApproxY, vtNFacApproxY, ptCFacApproxY
                            local tvtNx = {}
                            tvtNx[2] = vtN
@@ -5222,7 +5185,7 @@ local function MakeMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead, bSinglePa
                               dToolMaxDiam = min ( dFacDim1, dFacDim2, dToolTargetDiam)
                            end
                            local _, sPocketing = VerifyPocket( Proc, dToolMaxDiam, dFacElev, nil, sMyMchFind)
-                           bOk, sWarn2, sTuuidPk, dDiamTool, dDepth = MakePocket( Proc, nPartId, b3Solid, ptPs, tvtNx, nFacInd, sMyMchFind, nUseRoughTool, sPocketing, dMachDepth, nil, nil, bAllWithEndCap)
+                           bOk, sWarn2, sTuuidPk, dDiamTool, dDepth = MakePocket( Proc, nPartId, b3Solid, ptPs, tvtNx, nFacInd, sMyMchFind, nUseRoughTool, sPocketing, dMachDepth, nil, nil, bAllWithEndCap, nil, nPhase, nRawId)
                            if not bOk then return false, sWarn2 end
                            if sWarn2 then
                               if not sWarn then sWarn = '' end
@@ -5285,7 +5248,7 @@ local function MakeMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead, bSinglePa
                               sWarn = sErr
                               -- se ho antischeggia con fresa le inserisco
                            -- if nChamfer < 2 and nQAntisplintResult == 2 and ( bIsU or bIsL) then
-                              local nFacIndOri, dFacElevOri = BL.GetFaceWithMostAdj( Proc.Id, nPartId)
+                              local nFacIndOri, dFacElevOri = BL.GetFaceWithMostAdj( Proc, nPartId)
                               if nChamfer < 2 and nQAntisplintResult == 2 then
                                  local bOk, sWarn2 = ManageAntiSplintByMill( Proc, nPhase, nRawId, nPartId, b3Raw, 
                                                                            nFacIndOri, nAddGrpId, bMillDown, dDiamTool, nil, 
@@ -5386,12 +5349,12 @@ local function MakeMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead, bSinglePa
             local bMadeASbyBld = false
             local bPassThrou = ( Proc.Box:getDimY() > b3Raw:getDimY() - 1 or Proc.Box:getDimZ() > b3Raw:getDimZ() - 1)
             local bPassEdge = ((( bIsU and Proc.Fct == 3) or ( bIsL and Proc.Fct == 2)) and bSinglePart and Proc.Box:getDimX() < 0.9 * b3Raw:getDimX())
-            local nFacIndOri = BL.GetFaceWithMostAdj( Proc.Id, nPartId)
+            local nFacIndOri = BL.GetFaceWithMostAdj( Proc, nPartId)
             local vtNOri = EgtSurfTmFacetNormVersor( Proc.Id, nFacIndOri, GDB_ID.ROOT)
             if nChamfer < 2 and nQAntisplintResult == 1 and (( bIsU or bIsL) and ( bPassThrou or bPassEdge)) and ( vtNOri:getZ() > -0.087 or BD.DOWN_HEAD or BD.TURN) then
                local bOk
                local bSawDown = ( bMillDown and not bMillUp)
-               bMadeASbyBld, bOk, sWarn = ManageAntiSplintBySaw( Proc, b3Raw, bIsU, vtNOri, nFacIndOri, sWarn, bSawDown)
+               bMadeASbyBld, bOk, sWarn = ManageAntiSplintBySaw( Proc, b3Raw, b3Solid, bIsU, vtNOri, nFacIndOri, sWarn, bSawDown)
                if not bOk then return false, sWarn end
             end
             if nChamfer < 2 and nQAntisplintResult == 2 then
@@ -5432,7 +5395,7 @@ local function MakeMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead, bSinglePa
                   dMachiningDepth = dFacElev + dCollSic
                end
                local bOk, sWarn2
-               bOk, sWarn2, sTuuidPk, dDiamTool = MakePocket( Proc, nPartId, b3Solid, ptC, tvtNx, nFacInd, sMchFind, nUseRoughTool, sPocketing, dMachiningDepth, nil, nil, bAllWithEndCap, bLapJointAngTrasm)
+               bOk, sWarn2, sTuuidPk, dDiamTool = MakePocket( Proc, nPartId, b3Solid, ptC, tvtNx, nFacInd, sMchFind, nUseRoughTool, sPocketing, dMachiningDepth, nil, nil, bAllWithEndCap, bLapJointAngTrasm, nPhase, nRawId)
                if not bOk then return false, sWarn2 end
                if sWarn2 then
                   if not sWarn then sWarn = '' end
@@ -5530,62 +5493,104 @@ local function MakeLongMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead)
    -- ottengo la distanza tra la fine del pezzo e il pezzo successivo
    local dDistToNextPiece = BL.GetDistanceToNextPart( nRawId, nPhase)
    -- verifico se applicare gestione speciale delle giunzioni (U diretta come asse X)
-   local bAddEndCap = false
+   local bAddEndCapLeftSide, bAddEndCapRightSide = false, false
    local dAddLen = 0
    local bIsOpenU = ( Proc.Fct == 3 and not TestElleShape3( Proc))
    if bIsOpenU then
       local ptC0, vtN0 = EgtSurfTmFacetCenter( Proc.Id, 0, GDB_ID.ROOT)
       local ptC1, vtN1 = EgtSurfTmFacetCenter( Proc.Id, 1, GDB_ID.ROOT)
       local ptC2, vtN2 = EgtSurfTmFacetCenter( Proc.Id, 2, GDB_ID.ROOT)
-      if vtN0:getX() < 0.0175 and vtN1:getX() < 0.0175 and vtN2:getX() < 0.0175 then
+      if abs( vtN0:getX()) < 0.0175 and abs( vtN1:getX()) < 0.0175 and abs( vtN2:getX()) < 0.0175 then
          local dWidth = 0
          if vtN0 * vtN1 < -0.9998 then
-            bAddEndCap = true
+            bAddEndCapLeftSide = true
             dWidth = ( ptC1 - ptC0) * vtN0
          elseif vtN0 * vtN2 < -0.9998 then
-            bAddEndCap = true
+            bAddEndCapLeftSide = true
             dWidth = ( ptC2 - ptC0) * vtN0
          elseif vtN1 * vtN2 < -0.9998 then
-            bAddEndCap = true
+            bAddEndCapLeftSide = true
             dWidth = ( ptC2 - ptC1) * vtN1
          end
          dAddLen = min( dWidth, 100) / 2
       end
    end
+   -- verifico se applicare gestione speciale per tunnel
+   local nSurfBottomId
+   if Proc.Topology == 'Tunnel' then
+      bAddEndCapLeftSide = true
+      bAddEndCapRightSide = true
+      -- recupero centro e normale delle facce
+      local vtN = {}
+      for i = 1, Proc.Fct do
+         vtN[i] = EgtSurfTmFacetNormVersor( Proc.Id, i - 1, GDB_ID.ROOT)
+      end
+      -- calcolo l'orientamento del tunnel
+      local vtOrtho
+      local bAdj = EgtSurfTmFacetsContact( Proc.Id, 0, 1)
+      if bAdj then
+         vtOrtho = vtN[1] ^ vtN[2]
+      else
+         if Proc.Fct >= 3 then
+            vtOrtho = vtN[1] ^ vtN[3]
+         else
+            local sErr = 'Error : disconnected tunnel'
+            EgtOutLog( sErr)
+            return false, sErr
+         end
+      end
+      vtOrtho = EgtIf( abs( vtOrtho:getY()) > abs( vtOrtho:getZ()), -Y_AX(), -Z_AX())
+      local ptMaxBox = Proc.Box:getMax()
+      -- per creare le superfici di cap nei tunnel serve aggiungere una faccia di fondo
+      -- recupero il box della trave e lo ingrandisco di poco per essere sicuro di avere intersezione con il piano di fondo
+      local b3SolidExtended = EgtGetBBoxGlob( EgtGetFirstNameInGroup( nPartId, 'Box') or GDB_ID.NULL, GDB_BB.STANDARD)
+      b3SolidExtended:expand( 100 * GEO.EPS_SMALL)
+      -- superficie di fondo
+      nSurfBottomId = EgtSurfTmPlaneInBBox( nAddGrpId, ptMaxBox, vtOrtho, b3SolidExtended, GDB_ID.ROOT)
+      -- calcolo di quanto allargare le superfici
+      local dWidth = EgtIf( Proc.AffectedFaces.Front, Proc.Box:getDimZ(), Proc.Box:getDimY())
+      dAddLen = min( dWidth, 100) / 2
+   end
    -- la divido in parti lungo X
    local vAddId = {}
-   local bAllWithEndCap = bAddEndCap
+   local bAllWithEndCap = bAddEndCapLeftSide or bAddEndCapRightSide
    local nPart = max( ceil( Proc.Box:getDimX() / BD.LONGCUT_MAXLEN), 2)
    local dPartLen = Proc.Box:getDimX() / nPart
    local Xmin = Proc.Box:getMin():getX()
+   -- trimesh per la creazione degli EndCap
+   local AddIdCopy = EgtCopyGlob( Proc.Id, nAddGrpId) or GDB_ID.NULL
+   -- aggiunta di eventuale fondo alla trimesh per creare correttamente le superfici EndCap
+   if nSurfBottomId then EgtSurfTmAdd( AddIdCopy, nSurfBottomId) end
    for i = 1, nPart do
+      local nCapIdLeftSide, nCapIdRightSide
       -- eseguo divisione
       local AddId = EgtCopyGlob( Proc.Id, nAddGrpId) or GDB_ID.NULL
-      if i > 1 or bAddEndCap then
-         -- definizione del piano
+      if i > 1 or bAddEndCapLeftSide then
+         -- definizione del piano sinsitro
          local dAdd = EgtIf( i > 1, dAddLen, 0)
          local ptOn = Point3d( Xmin - dAdd + ( i - 1) * dPartLen, 0, 0)
          local vtN = -X_AX()
-         -- se richiesto, creazione tappo
-         local CapId
-         if bAddEndCap then
-            local nFirstId, nPnt, nCrv, nSrf = EgtPlaneSurfTmInters( ptOn, vtN, AddId, nAddGrpId, GDB_RT.GLOB)
+         -- se richiesto, creazione tappo sinistro
+         if bAddEndCapLeftSide then
+            local nFirstId, nPnt, nCrv, nSrf = EgtPlaneSurfTmInters( ptOn, vtN, AddIdCopy, nAddGrpId, GDB_RT.GLOB)
             if nPnt == 0 and nCrv == 1 and nSrf == 0 then
                EgtCloseCurveCompo( nFirstId)
-               CapId = EgtSurfTmByFlatContour( nAddGrpId, nFirstId)
-               if not CapId then bAllWithEndCap = false end
-            end
-            if nFirstId then
-               for nId = nFirstId, nFirstId + nPnt + nCrv + nSrf - 1 do
-                  EgtErase( nId)
-               end
+               nCapIdLeftSide = EgtSurfTmByFlatContour( nAddGrpId, nFirstId)
+               if not nCapIdLeftSide then bAllWithEndCap = false end
             end
+          if nFirstId then
+             for nId = nFirstId, nFirstId + nPnt + nCrv + nSrf - 1 do
+                EgtErase( nId)
+             end
+          end
          end
-         -- taglio della superficie
-         EgtCutSurfTmPlane( AddId, ptOn, vtN, true, GDB_RT.GLOB)
-         -- se esiste, aggiunta del tappo
-         if CapId then
-            AddId = EgtSurfTmBySewing( nAddGrpId, { AddId, CapId})
+         if i > 1 then
+            -- taglio della superficie lato sinistro
+            EgtCutSurfTmPlane( AddId, ptOn, vtN, true, GDB_RT.GLOB)
+         end
+         -- se esiste, aggiunta del tappo sinistro
+         if nCapIdLeftSide then
+            AddId = EgtSurfTmBySewing( nAddGrpId, { AddId, nCapIdLeftSide})
             -- se prima spezzatura, allungamento per non lasciare archi
             if i == 1 then
                local b3Box = EgtGetBBoxGlob( AddId, GDB_BB.STANDARD)
@@ -5595,21 +5600,45 @@ local function MakeLongMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead)
          end
       end
       if i < nPart then
-         local ptOn = Point3d( Xmin + i * dPartLen, 0, 0)
+         -- definizione del piano destro
+         local dAdd = EgtIf( i < nPart, dAddLen, 0)
+         local ptOn = Point3d( Xmin + dAdd + i * dPartLen, 0, 0)
          local vtN = X_AX()
+         -- taglio della superficie lato destro
          EgtCutSurfTmPlane( AddId, ptOn, vtN, true, GDB_RT.GLOB)
+         -- se richiesto, creazione tappo destro
+         if bAddEndCapRightSide then
+            local nFirstId, nPnt, nCrv, nSrf = EgtPlaneSurfTmInters( ptOn, vtN, AddIdCopy, nAddGrpId, GDB_RT.GLOB)
+            if nPnt == 0 and nCrv == 1 and nSrf == 0 then
+               EgtCloseCurveCompo( nFirstId)
+               nCapIdRightSide = EgtSurfTmByFlatContour( nAddGrpId, nFirstId)
+               if not nCapIdRightSide then bAllWithEndCap = false end
+            end
+           if nFirstId then
+              for nId = nFirstId, nFirstId + nPnt + nCrv + nSrf - 1 do
+                 EgtErase( nId)
+              end
+           end
+         end
+      end
+      -- se esiste, aggiunta del tappo destro
+      if nCapIdRightSide then
+         AddId = EgtSurfTmBySewing( nAddGrpId, { AddId, nCapIdRightSide})
       end
       EgtSetName( AddId, 'AddPart_' .. tostring( Proc.Id) .. '_' .. tostring( i))
       -- eseguo inserimento in modo da ordinare da X+ a X-
       table.insert( vAddId, 1, AddId)
    end
+   EgtErase( AddIdCopy)
    -- applico le lavorazioni sulle diverse parti
    local sWarn
    local bPrevBhSideMill
    for i = 1, #vAddId do
       local b3Box = EgtGetBBoxGlob( vAddId[i], GDB_BB.STANDARD)
       local nFct = EgtSurfTmFacetCount( vAddId[i])
-      local AddProc = { Id = vAddId[i], Grp = Proc.Grp, Prc = Proc.Prc, Box = b3Box, TotBox = Proc.Box, Fct = nFct, Flg = Proc.Flg}
+      local AddProc = { Id = vAddId[i], Grp = Proc.Grp, Prc = Proc.Prc, Box = b3Box, TotBox = Proc.Box, Fct = nFct, Flg = Proc.Flg, PartId = Proc.PartId}
+      Topology.Classify( AddProc, b3Raw)
+      AddProc.AffectedFaces = BL.GetProcessAffectedFaces( AddProc)
       -- lasciare il false nel sesto parametro (perchè internamente viene verificato se diverso da nil)
       local bOk, sMyWarn 
       bOk, sMyWarn, bPrevBhSideMill = MakeMoreFaces( AddProc, nPhase, nRawId, nPartId, dOvmHead, bAllWithEndCap, bPrevBhSideMill, bAllWithEndCap)

LuaLibs/ProcessLongCut.lua

@@ -29,6 +29,7 @@
 -- 2023/05/03 Corretto SCC in caso di asse utensile allineato con Z.
 -- 2023/05/19 Migliorato calcolo e verifica affondamento per lavorazione con lama con codolo in mezzo.
 -- 2023/08/01 Ammesso uso lama da sotto fino a N +3deg in verticale.
+-- 2023/09/26 Modificata chiamata a GetFaceWithMostAdj.
 
 -- Tabella per definizione modulo
 local ProcessLongCut = {}
@@ -437,7 +438,7 @@ end
 -- Applicazione della lavorazione
 local function MakeByPocketing( Proc, nPhase, nRawId, nPartId)
    -- recupero la faccia con il maggior numero di adiacenze e l'elevazione relativa
-   local nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc.Id, nPartId)
+   local nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc, nPartId)
    -- cerco la svuotatura opportuna
    local sPocketing = ML.FindPocketing( 'OpenPocket', Proc.Box:getDimX())
    if not sPocketing then

LuaLibs/ProcessMortise.lua

@@ -1,4 +1,4 @@
--- ProcessMortise.lua by Egaltech s.r.l. 2023/07/21
+-- ProcessMortise.lua by Egaltech s.r.l. 2023/08/10
 -- Gestione calcolo mortase per Travi
 -- 2021/07/20 Aggiunta gestione rinvio angolare su FAST.
 -- 2021/12/01 Se frontale aggiungo taglio con Grp e Proc di vero taglio (per aggiornare ingombro di testa /coda).
@@ -12,6 +12,7 @@
 -- 2023/05/11 Se lavorazione in doppio e precedente no oppure di tipo diverso, forzo risalita a Zmax.
 -- 2023/06/27 Esclusa la ricerca lati aperti per le mortase passanti, per le quali il contorno finisce sulla faccia e il riconoscimento sbaglia.
 -- 2023/07/21 Correzioni per mortise passanti con curva Aux che guarda in basso.
+-- 2023/08/10 Modificata scelta SCC per tasche in Y+/- in coda o quasi.
 
 -- Tabella per definizione modulo
 local ProcessMortise = {}
@@ -469,11 +470,14 @@ function ProcessMortise.Make( Proc, nPhase, nRawId, nPartId, dCurrOvmH)
    local nSCC = MCH_SCC.NONE
    if bPockAngTrasm then
       nSCC = MCH_SCC.ADIR_NEAR
-   elseif not BD.C_SIMM then
-      if abs( vtN:getX()) < GEO.EPS_SMALL and abs( vtN:getY()) < GEO.EPS_SMALL then
+   elseif not BD.C_SIMM and not BD.TURN then
+      if AreSameOrOppositeVectorApprox( vtN, Z_AX()) then
          nSCC = MCH_SCC.ADIR_YM
-      elseif abs( vtN:getX()) < GEO.EPS_SMALL then
-         nSCC = MCH_SCC.ADIR_XP
+      elseif abs( vtN:getX()) < 0.1 then
+         local bNearTail = ( Proc.Box:getMax():getX() < b3Solid:getCenter():getX() and
+                             Proc.Box:getMax():getX() - b3Solid:getMin():getX() < 1000)
+         local bVeryShortPart = ( BD.LEN_VERY_SHORT_PART and b3Solid:getDimX() < BD.LEN_VERY_SHORT_PART)
+         nSCC = EgtIf( BL.IsPartFinalPhase( EgtGetCurrPhase()) or ( bNearTail and not bVeryShortPart), MCH_SCC.ADIR_XM, MCH_SCC.ADIR_XP)
       elseif vtN:getY() < GEO.EPS_SMALL then
          nSCC = EgtIf( bRevertSide, MCH_SCC.ADIR_YP, MCH_SCC.ADIR_YM)
       else
@@ -551,11 +555,14 @@ function ProcessMortise.Make( Proc, nPhase, nRawId, nPartId, dCurrOvmH)
       EgtSetMachiningParam( MCH_MP.INVERT, EgtIf( bCW, true, false))
       -- imposto posizione braccio porta testa
       local nSCC = MCH_SCC.NONE
-      if not BD.C_SIMM then
-         if abs( vtN:getX()) < GEO.EPS_SMALL and abs( vtN:getY()) < GEO.EPS_SMALL then
+      if not BD.C_SIMM and not BD.TURN then
+         if AreSameOrOppositeVectorApprox( vtN, Z_AX()) then
             nSCC = MCH_SCC.ADIR_YM
-         elseif abs( vtN:getX()) < GEO.EPS_SMALL then
-            nSCC = MCH_SCC.ADIR_XP
+         elseif abs( vtN:getX()) < 0.1 then
+            local bNearTail = ( Proc.Box:getMax():getX() < b3Solid:getCenter():getX() and
+                                Proc.Box:getMax():getX() - b3Solid:getMin():getX() < 1000)
+            local bVeryShortPart = ( BD.LEN_VERY_SHORT_PART and b3Solid:getDimX() < BD.LEN_VERY_SHORT_PART)
+            nSCC = EgtIf( BL.IsPartFinalPhase( EgtGetCurrPhase()) or ( bNearTail and not bVeryShortPart), MCH_SCC.ADIR_XM, MCH_SCC.ADIR_XP)
          elseif vtN:getY() < GEO.EPS_SMALL then
             nSCC = MCH_SCC.ADIR_YP
          else

LuaLibs/ProcessRidgeLap.lua

@@ -1,4 +1,4 @@
--- ProcessRidgeLap.lua by Egaltech s.r.l. 2023/05/18
+-- ProcessRidgeLap.lua by Egaltech s.r.l. 2023/08/10
 -- Gestione calcolo mezzolegno di testa per Travi
 
 -- Tabella per definizione modulo
@@ -195,7 +195,7 @@ function ProcessRidgeLap.Make( Proc, nPhase, nRawId, nPartId, dOvmHead)
    -- aggiornamento ingombro di testa o coda
    if Proc.Head then
       local dHCI = 0
-      if abs( vtN[vFaceOrd[3]]:getZ()) > 0.1 then
+      if vtN[vFaceOrd[3]]:getZ() > 0.1 then
          local b3Fac1 = EgtSurfTmGetFacetBBoxGlob( Proc.Id, vFaceOrd[1] - 1, GDB_BB.STANDARD)
          if b3Fac1 then dHCI = b3Raw:getMax():getX() - dOvmHead - b3Fac1:getMin():getX() end
       else
@@ -204,7 +204,7 @@ function ProcessRidgeLap.Make( Proc, nPhase, nRawId, nPartId, dOvmHead)
       BL.UpdateHCING( nRawId, dHCI)
    elseif Proc.Tail then
       local dTCI = 0
-      if abs( vtN[vFaceOrd[3]]:getZ()) > 0.1 then
+      if vtN[vFaceOrd[3]]:getZ() > 0.1 then
          local b3Fac1 = EgtSurfTmGetFacetBBoxGlob( Proc.Id, vFaceOrd[1] - 1, GDB_BB.STANDARD)
          if b3Fac1 then dTCI = b3Fac1:getMax():getX() - b3Solid:getMin():getX() end
       else

LuaLibs/ProcessTenon.lua

@@ -1,4 +1,4 @@
--- ProcessTenon.lua by Egaltech s.r.l. 2023/03/27
+-- ProcessTenon.lua by Egaltech s.r.l. 2023/09/26
 -- Gestione calcolo tenone per Travi
 -- 2021/10/04 Corretto calcolo HCING per pezzi piccoli.
 -- 2022/02/15 Aggiornata VerifyOrientation per macchine con testa da sotto.
@@ -10,6 +10,7 @@
 -- 2022/12/20 Aggiunta gestione smusso. Migliorata scelta tra sopra e sotto per il punto di inizio. In VerifyOrientation aggiunta trave medio alta.
 -- 2023/01/20 Modificata scelta lato di attacco per pezzo piccolo su macchine con pinza speciale (pinza 5).
 -- 2023/03/27 Migliorate condizioni scelta pretaglio con lama o fresa; gestione unificata con ProcessDtTenon. Migliorato calcolo altezza tenone.
+-- 2023/09/26 Per macchina TURN si assegna SCC per privilegiare accesso dal lato corto della trave.
 
 -- Tabella per definizione modulo
 local ProcessTenon = {}
@@ -325,7 +326,13 @@ function ProcessTenon.Make( Proc, nPhase, nRawId, nPartId, dOvmHead)
       end
       -- imposto posizione braccio porta testa
       local nSCC = MCH_SCC.NONE
-      if not BD.C_SIMM and not BD.TURN then
+      if BD.TURN then
+         if b3Solid:getDimY() >= b3Solid:getDimZ() then
+            nSCC = MCH_SCC.ADIR_ZP
+         else
+            nSCC = MCH_SCC.ADIR_YM
+         end
+      elseif not BD.C_SIMM then
          nSCC = MCH_SCC.ADIR_YM
          if abs( vtExtr:getY()) > 0.088 then
             nSCC = EgtIf( vtExtr:getX() < GEO.EPS_SMALL, MCH_SCC.ADIR_XM, MCH_SCC.ADIR_XP)

LuaLibs/ProcessTyroleanDovetail.lua

@@ -1,6 +1,7 @@
 -- ProcessTyroleanDovetail.lua by Egaltech s.r.l. 2022/03/21
 -- Gestione calcolo giunzione tirolese
 -- 2022/06/10 Aggiunto il parametro dOvmTail per gestire sovramateriali in coda diversi da OVM_MID (sezioni alte e larghe)
+-- 2023/09/26 Modificata chiamata a GetFaceWithMostAdj.
 
 -- Tabella per definizione modulo
 local ProcessTyroleanDovetail = {}
@@ -510,10 +511,10 @@ local function MakeMachByMill( Proc, nPhase, nRawId, nPartId, dOvmHead, b3Raw, b
       dStepmach = dMaxMat * 0.5
    end
    -- recupero la faccia con il maggior numero di adiacenze e l'elevazione relativa
-   local nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc.Id, nPartId, false, sin(dSideAngle))
+   local nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc, nPartId, false, sin(dSideAngle))
    if not nFacInd or nFacInd < 0 then
       -- provo eliminando i sottosquadra
-      nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc.Id, nPartId, false, -2)
+      nFacInd, dFacElev = BL.GetFaceWithMostAdj( Proc, nPartId, false, -2)
       if not nFacInd or nFacInd < 0 then
          local sErr = 'Error on process ' .. tostring( Proc.Id) .. ' MakeMachByMill could not find reference face'
          EgtOutLog( sErr)

UpdateLog.txt

@@ -1,6 +1,21 @@
            ==== Beam Update Log ====
 
-Versione 2.5g4 (07/08/2023)
+Versione 2.5i2 (27/09/2023)
+- Added : aggiunto riconoscimento sottosquadro da facce non adiacenti
+- Fixed : in LapJoint corretto caso tunnel lungo non lavorato correttamente
+- Fixed : in foratura quando errore in applicazione lavorazione si inverte e riprova solo se singola su foro aperto
+- Modif : in lavorazione tenoni per Turn si assegna SCC per privilegiare accesso dal lato corto della trave.
+
+Versione 2.5i1 (12/09/2023)
+- Fixed : in LapJoint gestito correttamente il ritorno nil di GetUShapeWidth [Ticket #1354]
+- Modif : in Cut abbassato a 590 mm il limite per convertire in LongCut [Ticket #1448].
+
+Versione 2.5h2 (11/08/2023)
+- Fixed : corrette forature con aggregato quando lato aperto
+- Fixed : in RidgeLap corretto ingombro lavorazione in testa e coda quando rivolta verso il basso
+- Modif : in LapJoint e Mortise modificata posizione braccio per FAST quando vicino alla coda della trave.
+
+Versione 2.5h1 (07/08/2023)
 - Fixed : tagli doppi di lato non effettuati se macchina tipo PF e pezzo alto [Ticket #1400]
 - Fixed : in split e headcut corretto calcolo allungamenti/accorciamenti pezzi alti per contemplare anche taglio singolo [Ticket #1227]
 - Fixed : in lapjoint piccola correzione al calcolo della distanza di sicurezza in base al gambo dell'utensile [Ticket #1227]

Version.lua

@@ -1,6 +1,6 @@
--- Version.lua by Egaltech s.r.l. 2023/08/01
+-- Version.lua by Egaltech s.r.l. 2023/09/26
 -- Gestione della versione di Beam
 
 NAME = 'Beam'
-VERSION = '2.5g4'
+VERSION = '2.5i2'
 MIN_EXE = '2.5c1'