Merge branch 'release/2.5l3'

- modifiche per centrare i pezzi in Y sulla tavola (da flag BeamData.CENTER_BEAM)
- modifiche per rendere più veloci i calcoli di FeatureTopology
- scorporata CalcHeadTailMachBeforeIntersDrillings da CollectFeatures
- migliorato calcolo attacco/uscite lame in FacesBySaw.

LuaLibs/BeamExec.lua

@@ -1,4 +1,4 @@
--- BeamExec.lua by Egaltech s.r.l. 2023/11/08
+-- BeamExec.lua by Egaltech s.r.l. 2023/12/26
 -- Libreria esecuzione lavorazioni per Travi
 -- 2019/07/11 Aggiunta gestione stato rotazione di feature per TS3.
 -- 2019/09/04 Corretto controllo feature di testa e coda con sovramateriale di testa elevato.
@@ -51,6 +51,8 @@
 -- 2022/09/26 In ClassifyTopology aggiunto passaggio del parametro nRawId
 -- 2023/10/24 Aggiunta scrittura parametro BARLEN nelle info del mach group
 -- 2023/11/08 Aggiunta gestione processi Variant.
+-- 2023/11/30 Migliorato il calcolo elevazione con l'utilizzo della nuova funzione EgtSurfTmFacetElevationInBBox.
+-- 2023/12/26 Modifiche per centrare i pezzi in Y sulla tavola.
 
 -- Tabella per definizione modulo
 local BeamExec = {}
@@ -255,6 +257,40 @@ local function IsTailFeature( Proc, b3Raw, dCurrOvmH, dCurrOvmT)
    return false
 end
 
+-------------------------------------------------------------------------------------------------------------
+local function NeedTopologyFeature( Proc)
+   -- richiedono calcolo topologia
+   if Cut.Identify( Proc) then
+      return true
+   end
+   if DoubleCut.Identify( Proc) then
+      return true
+   end
+   if LongCut.Identify( Proc) then
+      return true
+   end
+   if Long2Cut.Identify( Proc) then
+      return true
+   end
+   if SawCut.Identify( Proc) then
+      return true
+   end
+   if RidgeLap.Identify( Proc) then
+      return true
+   end
+   if LapJoint.Identify( Proc) then
+      return true
+   end
+   if FrenchRidgeLap.Identify( Proc) then
+      return true
+   end
+   if Chamfer.Identify( Proc) then
+      return true
+   end
+   -- tutte le altre non richiedono calcolo topologia
+   return false
+end
+
 -------------------------------------------------------------------------------------------------------------
 local function CollectFeatures( PartId, b3Raw, dCurrOvmH, dCurrOvmT)
    local dRawW = b3Raw:getDimY()
@@ -264,8 +300,6 @@ local function CollectFeatures( PartId, b3Raw, dCurrOvmH, dCurrOvmT)
    local LayerId = {}
    LayerId[1] = BL.GetAddGroup( PartId)
    LayerId[2] = EgtGetFirstNameInGroup( PartId or GDB_ID.NULL, 'Processings')
-   local nMachineBeforeIntersectingDrillingsIdHead, dMachineBeforeIntersectingDrillingsXHead, nMachineBeforeIntersectingDrillingsIdTail, dMachineBeforeIntersectingDrillingsXTail = nil, GEO.INFINITO, nil, -GEO.INFINITO
-   local b3MachineBeforeIntersectingDrillingsBoxHead, b3MachineBeforeIntersectingDrillingsBoxTail = nil, nil
    for nInd = 1, 2 do
       local ProcId = EgtGetFirstInGroup( LayerId[nInd] or GDB_ID.NULL)
       while ProcId do
@@ -301,28 +335,24 @@ local function CollectFeatures( PartId, b3Raw, dCurrOvmH, dCurrOvmT)
                end
                Proc.Box = EgtGetBBoxGlob( ProcId, GDB_BB.STANDARD)
                if b3Raw then
+                  local b3Solid = EgtGetBBoxGlob( EgtGetFirstNameInGroup( PartId, 'Box') or GDB_ID.NULL, GDB_BB.STANDARD)
                   -- 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)}
+                  if not Proc.Face then
+                     Proc.Face = {}
+                     for i = 1, Proc.Fct do
+                        Proc.Face[i] = { Id = i - 1, VtN = EgtSurfTmFacetNormVersor( Proc.Id, i - 1, GDB_ID.ROOT )}
+                        if Proc.Fct < 10 then
+                           Proc.Face[i].Elevation = EgtSurfTmFacetElevationInBBox( Proc.Id, i - 1, b3Solid, true, GDB_ID.ROOT)
+                        end
+                     end
                   end
                end
                if Proc.Box and not Proc.Box:isEmpty() then 
+                  Proc.NeedTopology = NeedTopologyFeature( Proc)
                   Proc.Head = IsHeadFeature( Proc, b3Raw, dCurrOvmH)
                   Proc.Tail, Proc.AdvTail = IsTailFeature( Proc, b3Raw, dCurrOvmH, dCurrOvmT)
-                  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
-                        b3MachineBeforeIntersectingDrillingsBoxHead = Proc.Box
-                     elseif Proc.Tail and Proc.Box:getCenter():getX() > dMachineBeforeIntersectingDrillingsXTail then
-                        dMachineBeforeIntersectingDrillingsXTail = Proc.Box:getCenter():getX()
-                        nMachineBeforeIntersectingDrillingsIdTail = Proc.Id
-                        b3MachineBeforeIntersectingDrillingsBoxTail = Proc.Box
-                     end
-                  end
                   table.insert( vProc, Proc)
                   -- se foro
                   if Drill.Identify( Proc) then
@@ -376,11 +406,38 @@ local function CollectFeatures( PartId, b3Raw, dCurrOvmH, dCurrOvmT)
          ProcId = EgtGetNext( ProcId)
       end
    end
-   local vMachineBeforeIntersectingDrillings = {
-      Head = { Id = nMachineBeforeIntersectingDrillingsIdHead, Box = b3MachineBeforeIntersectingDrillingsBoxHead},
-      Tail = { Id = nMachineBeforeIntersectingDrillingsIdTail, Box = b3MachineBeforeIntersectingDrillingsBoxTail}
+   return vProc
+end
+
+-------------------------------------------------------------------------------------------------------------
+local function CalcHeadTailMachBeforeIntersDrillings( vProc, b3Raw)
+   local nHeadId
+   local dHeadX = GEO.INFINITO
+   local nTailId
+   local dTailX = -GEO.INFINITO
+   local b3HeadBox
+   local b3TailBox
+   for i = 1, #vProc do
+      local Proc = vProc[i]
+      if Proc.Box and not Proc.Box:isEmpty() then 
+         if Proc.Fct == 1 and BL.IsFeatureCuttingEntireSection( Proc.Box, b3Raw:getDimY(), b3Raw:getDimZ()) and ( Proc.Head or Proc.Tail) and Proc.Prc ~= 340 and Proc.Prc ~= 350 then
+            if Proc.Head and Proc.Box:getCenter():getX() < dHeadX then
+               dHeadX = Proc.Box:getCenter():getX()
+               nHeadId = Proc.Id
+               b3HeadBox = Proc.Box
+            elseif Proc.Tail and Proc.Box:getCenter():getX() > dTailX then
+               dTailX = Proc.Box:getCenter():getX()
+               nTailId = Proc.Id
+               b3TailBox = Proc.Box
+            end
+         end
+      end
+   end
+   local vMachBeforeIntersDrillings = {
+      Head = { Id = nHeadId, Box = b3HeadBox},
+      Tail = { Id = nTailId, Box = b3TailBox}
    }
-   return vProc, vMachineBeforeIntersectingDrillings
+   return vMachBeforeIntersDrillings
 end
 
 -------------------------------------------------------------------------------------------------------------
@@ -494,7 +551,8 @@ function BeamExec.ProcessBeams( dRawW, dRawH, dRawL, dOvmHead, vBeam, bMachGroup
    -- Area tavola
    local b3Tab = EgtGetTableArea()
    -- Calcolo posizione estremo TR/BR della tavola rispetto a sua origine in BL
-   BD.OriXR = Point3d( b3Tab:getDimX(), EgtIf( BD.RIGHT_LOAD, 0, b3Tab:getDimY()), 0)
+   local dPosY = EgtIf( BD.CENTER_BEAM, ( b3Tab:getDimY() + dRawW * EgtIf( BD.RIGHT_LOAD, -1, 1)) / 2, EgtIf( BD.RIGHT_LOAD, 0, b3Tab:getDimY()))
+   BD.OriXR = Point3d( b3Tab:getDimX(), dPosY, 0)
    BD.PosXR = EgtIf( BD.RIGHT_LOAD, MCH_CR.BR, MCH_CR.TR)
    
    -- Impostazione dell'attrezzaggio di default
@@ -1172,7 +1230,7 @@ local function ClassifyTopology( vProc, nRawId)
          nRecognized = nRecognized + 1
       end
    end
-   
+
    return nRecognized
 end
 
@@ -1850,10 +1908,11 @@ function BeamExec.ProcessFeatures()
       local dCurrOvmH = EgtGetInfo( nRawId, 'HOVM', 'd') or 0
       local dCurrOvmT = EgtGetInfo( nRawId, 'TOVM', 'd') or 0
       -- recupero le feature di lavorazione della trave
-      local vProc, vMachineBeforeIntersectingDrillings = CollectFeatures( nPartId, b3Raw, dCurrOvmH, dCurrOvmT)
+      local vProc = CollectFeatures( nPartId, b3Raw, dCurrOvmH, dCurrOvmT)
       -- verifica presenza forature influenzate da lavorazioni di testa o coda
       if BD.IMPROVE_HEAD_TAIL_DRILLINGS then
-         SetDrillingsToMachineAfterHeadOrTailCut( vProc, vMachineBeforeIntersectingDrillings)
+         vMachBeforeIntersDrillings = CalcHeadTailMachBeforeIntersDrillings( vProc, b3Raw)
+         SetDrillingsToMachineAfterHeadOrTailCut( vProc, vMachBeforeIntersDrillings)
       end
       -- verifica presenza di feature specchiate per eventuali lavorazioni simultanee
       if BD.TWO_EQUAL_HEADS or BD.DOWN_HEAD then

LuaLibs/BeamLib.lua

@@ -1,4 +1,4 @@
--- BeamLib.lua by Egaltech s.r.l. 2023/03/22
+-- BeamLib.lua by Egaltech s.r.l. 2023/12/15
 -- Libreria globale per Travi
 -- 2020/07/28 Corretto calcolo attacchi e uscite di lame per non uscire dalla faccia sotto.
 -- 2020/08/18 Aggiunto a GetNearestParalOpposite e GetNearestOrthoOpposite parametro opzionale vtNorm.
@@ -26,6 +26,9 @@
 -- 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.
+-- 2023/11/30 Calcolo elevazione velocizzato e centralizzato tramite la funzione GetFaceElevation. Se l'elevazione è già calcolata la recupera da Proc, altrimenti la calcola al momento.
+-- 2023/12/15 Correzione a GetFaceWithMostAdj per calcoli con facce che fanno ombra.
+-- 2023/08/01 Correzione a Is3EdgesApprox per casi in cui nAddGroup non esista ancora.
 
 -- Tabella per definizione modulo
 local BeamLib = {}
@@ -440,54 +443,22 @@ function BeamLib.GetBoxFaceNorm( b3Box, ptP, vtV)
 end
 
 ---------------------------------------------------------------------
-function BeamLib.GetFaceElevation( nSurfId, nFac, nPartId)
-   -- centro e normale della faccia
-   local ptC, vtN = EgtSurfTmFacetCenter( nSurfId, nFac, GDB_ID.ROOT)
-   if not ptC or not vtN then return 0 end
-   -- riferimento OCS della faccia per ricavare elevazione rispetto alla faccia della superficie
-   local frOCS = Frame3d( ptC, vtN)
-   local b3Box = EgtGetBBoxRef( nSurfId, GDB_BB.STANDARD, frOCS)
-   local dElev = b3Box:getMax():getZ()
-   -- se definito identificativo di pezzo
-   if nPartId then
-      -- se superficie con più facce
-      if EgtSurfTmFacetCount( nSurfId) > 1 then
-         -- determino elevazione del centro faccia rispetto al box del pezzo
-         local _, dCenElev = BeamLib.GetPointDirDepth( nPartId, ptC, vtN)
-         if dCenElev and dCenElev > dElev then
-            dElev = dCenElev
-         end
-      -- altrimenti superficie ad una sola faccia
-      else
-         -- determino elevazione box del pezzo rispetto alla faccia
-         local b3Solid = EgtGetBBoxRef( EgtGetFirstNameInGroup( nPartId, 'Box') or GDB_ID.NULL, GDB_BB.STANDARD, frOCS)
-         local dSolidElev = b3Solid:getMax():getZ()
-         if b3Solid and dSolidElev > dElev then
-            dElev = dSolidElev
-         end
-      end
+function BeamLib.GetFaceElevation( procOrProcId, nFacet, b3Solid)
+   local Proc, nProcId
+   if type( procOrProcId) == "table" then
+      Proc = procOrProcId
+      nProcId = Proc.Id
+   else
+      nProcId = procOrProcId
+   end
+   local dElevation
+   if not Proc or not Proc.Face or not Proc.Face[nFacet + 1].Elevation then
+      dElevation = EgtSurfTmFacetElevationInBBox( nProcId, nFacet, b3Solid, true, GDB_ID.ROOT)
+   else
+      dElevation = Proc.Face[nFacet + 1].Elevation
    end
-   return dElev
-end
 
----------------------------------------------------------------------
-function BeamLib.GetOtherFaceElevation( nSurfId, nOtherSurfId, nOtherFac, nPartId)
-   -- centro e normale della faccia
-   local ptC, vtN = EgtSurfTmFacetCenter( nOtherSurfId, nOtherFac, GDB_ID.ROOT)
-   if not ptC or not vtN then return 0 end
-   -- riferimento OCS della faccia per ricavare elevazione rispetto alla faccia della superficie
-   local frOCS = Frame3d( ptC, vtN)
-   local b3Box = EgtGetBBoxRef( nSurfId, GDB_BB.STANDARD, frOCS)
-   local dElev = b3Box:getMax():getZ()
-   -- se definito identificativo di pezzo
-   if nPartId then
-      -- determino elevazione del centro faccia rispetto al box del pezzo
-      local _, dCenElev = BeamLib.GetPointDirDepth( nPartId, ptC, vtN)
-      if dCenElev and dCenElev > dElev then
-         dElev = dCenElev
-      end
-   end
-   return dElev
+   return dElevation
 end
 
 ---------------------------------------------------------------------
@@ -515,6 +486,7 @@ end
 ---------------------------------------------------------------------
 function BeamLib.GetFaceWithMostAdj( Proc, nPartId, bCompare3Fc, dCosSideAng)
    local nSurfId = Proc.Id
+   local b3Solid = EgtGetBBoxGlob( EgtGetFirstNameInGroup( nPartId, 'Box') or GDB_ID.NULL, GDB_BB.STANDARD)
    -- verifica che la superficie non sia chiusa e quindi non lavorabile
    if EgtSurfIsClosed( nSurfId) then
       return
@@ -590,14 +562,19 @@ function BeamLib.GetFaceWithMostAdj( Proc, nPartId, bCompare3Fc, dCosSideAng)
                local frOCSMaster = Frame3d( vPtC[i], vvtN[i])
                EgtScale( nShadowContourId, frOCSMaster, 1, 1, 0, GDB_RT.GLOB)
                EgtModifyCurveExtrusion( nShadowContourId, vvtN[i], GDB_RT.GLOB)
+               local nShadowSurfFlatId
+               -- se non troppo ortogonale (almeno 0.5deg di delta)
                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
+               if abs( dResV) > 0.087 then
+                  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( nShadowSurfFlatId)
                end
-               EgtErase( { nMasterSurfFlatId, nShadowSurfFlatId})
+               EgtErase( nMasterSurfFlatId)
             end
             EgtErase( { nMasterContourId, nShadowContourId, nShadowFacetId})
          end
@@ -645,7 +622,7 @@ function BeamLib.GetFaceWithMostAdj( Proc, nPartId, bCompare3Fc, dCosSideAng)
    local dMinElev, dMinElev2 = GEO.INFINITO, GEO.INFINITO
    local dtElev = {}
    for i = 1, #nFacInd do
-      local dElev =  BeamLib.GetFaceElevation( nSurfId, nFacInd[i], nPartId)
+      local dElev = BeamLib.GetFaceElevation( Proc, nFacInd[i], b3Solid)
       table.insert( dtElev, dElev)
       if dElev < dMinElev then
          if dMinElev < dMinElev2 then
@@ -1096,6 +1073,10 @@ 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)
+   if not nAddGrpId then
+      local nEdges = #(EgtSurfTmFacetAdjacencies( Proc.Id, nFacet)[1])
+      return ( nEdges == 3)
+   end
    local bResult = false
    local nContourId, nContourCnt = EgtExtractSurfTmFacetLoops( Proc.Id, nFacet, nAddGrpId)
    if not nContourId then return false end

LuaLibs/FaceByPocket.lua

@@ -1,5 +1,6 @@
 -- FaceByPocket.lua by Egaltech s.r.l. 2023/04/04
 -- Gestione svuotatura di feature con una faccia
+-- 2023/11/30 Calcolo elevazione velocizzato e centralizzato tramite la funzione GetFaceElevation.
 
 -- Tabella per definizione modulo
 local FaceByPocket = {}
@@ -78,7 +79,7 @@ function FaceByPocket.Make( Proc, nSurfId, nFacet, sPocketing, nPartId, b3Solid,
    end
    -- dati della faccia
    local ptC, vtN = EgtSurfTmFacetCenter( nSurfId, nFacet, GDB_ID.ROOT)
-   local dElev = BL.GetFaceElevation( nSurfId, nFacet, nPartId)
+   local dElev = BL.GetFaceElevation( nSurfId, nFacet, b3Solid)
    -- determino numero e valore degli step di lavorazione
    local nSurfStep = ceil( dElev / dMaxDepth)
    local dSurfStep = dElev / nSurfStep

LuaLibs/FacesBySaw.lua

@@ -1,4 +1,4 @@
--- FacesBySaw.lua by Egaltech s.r.l. 2023/09/26
+-- FacesBySaw.lua by Egaltech s.r.l. 2023/12/30
 -- 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.
@@ -20,6 +20,9 @@
 -- 2023/09/26 Piccola modifica per Turn su bInvert di MakeOne.
 -- 2023/10/24 In MakeOne migliorata gestione taglio con percorso bilinea. Aggiunta funzione GetNameSolidFaceIncludingLine.
 -- 2023/11/14 In MakeOne migliorato calcolo scelta soluzione per macchina TURN
+-- 2023/11/28 In MakeTwo raffinamento calcolo vtRef per casi dubbi.
+-- 2023/12/06 In CalcLeadInOutPerpGeom gestito caso in cui la geometria della feature esce dal grezzo.
+-- 2023/12/30 Modifiche in CalcLeadInOutPerpGeom e CalcLeadInOutTangGeom con uso di EgtCAvToolPosBox.
 
 -- Tabella per definizione modulo
 local FacesBySaw = {}
@@ -37,6 +40,7 @@ local ML = require( 'MachiningLib')
 
 ---------------------------------------------------------------------
 function MakeParallelOne( nSurfId, nFacet, sCutting, dSawDiam, nFaceUse, dVzLimDwnUp, dCutExtra, dCutSic, dCutOffset, dAccStart, dAccEnd, sNotes, b3Raw, bForceInvert)
+   EgtOutLog( 'FacesBySaw.MakeParallelOne', 3)
    -- dati della faccia
    local ptC, vtN = EgtSurfTmFacetCenter( nSurfId, nFacet, GDB_ID.ROOT)
    -- accetto solo facce perpendicolari all'asse X della trave
@@ -179,6 +183,7 @@ function FacesBySaw.MakeOne( nSurfId, nFacet, sCutting, dSawDiam, Par5, dVzLimDw
       return MakeParallelOne( nSurfId, nFacet, sCutting, dSawDiam, Par5, dVzLimDwnUp, dCutExtra, dCutSic, dCutOffset, dAccStart, dAccEnd, sNotes, b3Raw, bForceInvert)
    end
    -- la lama ha asse perpendicolare alla faccia
+   EgtOutLog( 'FacesBySaw.MakeOne', 3)
    -- dati della faccia
    local ptC, vtN = EgtSurfTmFacetCenter( nSurfId, nFacet, GDB_ID.ROOT)
    -- risolvo parametro ambiguo
@@ -191,6 +196,7 @@ function FacesBySaw.MakeOne( nSurfId, nFacet, sCutting, dSawDiam, Par5, dVzLimDw
       nOrthoOpposite = Par5
       vtOrthO = BL.GetVersRef( Par5)
    end
+   EgtOutLog( 'VtOrthO='..tostring( vtOrthO)..' FaceUse='..tostring( nOrthoOpposite), 3)
    -- verifico se testa da sotto oppure se lavorazione sotto con testa da sopra
    if not dVzLimDwnUp then dVzLimDwnUp = BL.GetNzLimDownUp( b3Raw, vtN, vtOrthO) end
    local bDownHead = ( dVzLimDwnUp and dVzLimDwnUp < - 1.5)
@@ -230,7 +236,7 @@ function FacesBySaw.MakeOne( nSurfId, nFacet, sCutting, dSawDiam, Par5, dVzLimDw
    local dAllStart = 0
    local dAllEnd = 0
    local bIsBiLinea = false
-   local bCosAngleL1L2 = 0
+   local dCosAngleL1L2 = 0
    local dDist1 = dist( ptP1, ptPm)
    local dDist2 = dist( ptP2, ptPm)
    -- verifico se la bilinea si trova sul bordo del solido, quindi è una geometria aperta
@@ -249,10 +255,8 @@ function FacesBySaw.MakeOne( nSurfId, nFacet, sCutting, dSawDiam, Par5, dVzLimDw
       local dCosMax = 0.951       -- cos( 18°)
       local dLenMin = 30
       local dLenMax = max( 0.5 * dSawDiam * 0.17365 + 1, 2 * dLenMin)
-      --if vtTg1 * vtTg2 < dCosMax or ( dDist1 < dLenMax and dDist1 > dLenMin) or ( dDist2 < dLenMax and dDist2 > dLenMin) then
-
-      bCosAngleL1L2 = vtTg1 * vtTg2
-      if bCosAngleL1L2 < dCosMax then
+      dCosAngleL1L2 = vtTg1 * vtTg2
+      if dCosAngleL1L2 < dCosMax then
          local dOrtho1 = abs( vtTg1 * vtOrthO)
          local dOrtho2 = abs( vtTg2 * vtOrthO)
          if dOrtho1 < dOrtho2 or ( abs( dOrtho1 - dOrtho2) < 0.1 and dDist1 > 4 * dDist2) then
@@ -299,43 +303,46 @@ function FacesBySaw.MakeOne( nSurfId, nFacet, sCutting, dSawDiam, Par5, dVzLimDw
    -- parametri di attacco/uscita
    local b3Box = BBox3d( b3Raw)
    b3Box:expand( dCutSic)
+   local ptPa1 = ptP1 + dAccStart * vtTg
+   local ptPa2 = ptP2 + dAccEnd * vtTg
    -- attacco perpendicolare
-   local dLiTang, dLiPerp, dLoTang, dLoPerp, vtLio = FacesBySaw.CalcLeadInOutPerpGeom( ptP1, ptP2, vtV1, vtV2, vtN, dSawDiam/2, vtRef, dCutExtra, b3Box)
+   local dLiTang, dLiPerp, dLoTang, dLoPerp, vtLio = FacesBySaw.CalcLeadInOutPerpGeom( ptPa1, ptPa2, vtV1, vtV2, vtN, dSawDiam / 2, vtRef, dCutExtra, b3Box)
    local dLenLi = sqrt( dLiTang * dLiTang + dLiPerp * dLiPerp)
    local dLenLo = sqrt( dLoTang * dLoTang + dLoPerp * dLoPerp)
    -- attacco tangente
-   local dLi2Tang, dLi2Perp, dLo2Tang, dLo2Perp = FacesBySaw.CalcLeadInOutTangGeom( ptP1, ptP2, vtN, dSawDiam/2, vtRef, dCutExtra, b3Box)
+   local dLi2Tang, dLi2Perp, dLo2Tang, dLo2Perp = FacesBySaw.CalcLeadInOutTangGeom( ptPa1, ptPa2, vtN, dSawDiam / 2, vtRef, dCutExtra, b3Box)
    local dLenLi2 = abs( dLi2Tang)
    local dLenLo2 = abs( dLo2Tang)
 
    -- se il lato non lavorato della bilinea è aperto, setto entrata/uscita con la stessa direzione
    if bIsBiLinea then
       -- angolo tra le due linee
-      local dAlpha = acos( abs( bCosAngleL1L2))
+      local dCosAlpha = dCosAngleL1L2
+      local dSinAlpha = sqrt( max( 1 - dCosAlpha * dCosAlpha, 0))
       -- se ho accorciato ingresso, setto componente tangente e perpendicolare sul percorso di entrata
       if abs( dAllStart) > 100 * GEO.EPS_SMALL and bIsL1OnFace then
          -- controllo prima che il secondo lato non sia già incluso nella lavorazione del primo
-         local dDistPtTang = cos( dAlpha) * dDist1
-         local dDistPtPerp = abs( sin( dAlpha) * dDist1)
+         local dDistPtTang = dCosAlpha * dDist1
+         local dDistPtPerp = dSinAlpha * dDist1
          local dDistToCenter = 0.5 * dSawDiam - dDistPtPerp
          local dDistPointToCenter = sqrt( dDistPtTang * dDistPtTang + dDistToCenter * dDistToCenter)
          -- se distanza al punto è maggiore del raggio lama, significa che non ho già lavorato, quindi calcolo entrata opportunamente
          if dDistPointToCenter > 0.5 * dSawDiam then
-            dLiTang = -dAllStart * cos( dAlpha)
-            dLiPerp = sin( dAlpha) * dDist1
+            dLiTang = -dAllStart * dCosAlpha
+            dLiPerp = dDist1 * dSinAlpha 
          end
       end
       -- se ho accorciato uscita, setto componente tangente e perpendicolare sul percorso di uscita
       if abs( dAllEnd) > 100 * GEO.EPS_SMALL and bIsL2OnFace then
          -- controllo prima che il secondo lato non sia già incluso nella lavorazione del primo
-         local dDistPtTang = cos( dAlpha) * dDist2
-         local dDistPtPerp = abs( sin( dAlpha) * dDist2)
+         local dDistPtTang = dCosAlpha * dDist2
+         local dDistPtPerp = dSinAlpha * dDist2
          local dDistToCenter = 0.5 * dSawDiam - dDistPtPerp
          local dDistPointToCenter = sqrt( dDistPtTang * dDistPtTang + dDistToCenter * dDistToCenter)
          -- se distanza al punto è maggiore del raggio lama, significa che non ho già lavorato, quindi calcolo uscita opportunamente
          if dDistPointToCenter > 0.5 * dSawDiam then
-            dLoTang = -dAllEnd * cos( dAlpha)
-            dLoPerp = sin( dAlpha) * dDist2
+            dLoTang = -dAllEnd * dCosAlpha
+            dLoPerp = dDist2 * dSinAlpha 
          end
       end
    end
@@ -354,7 +361,7 @@ function FacesBySaw.MakeOne( nSurfId, nFacet, sCutting, dSawDiam, Par5, dVzLimDw
       dLiTang, dLiPerp, dLoTang, dLoPerp = dLi2Tang, dLi2Perp, dLo2Tang, dLo2Perp
       if BD.TURN then
          local dMove = dist( ptP1, ptP2)
-         dLoTang = -( dLiTang + dAllStart - dAccStart + dAllEnd - dAccEnd + dMove)
+         dLoTang = -( dLiTang - dAccStart - dAccEnd + dMove)
          dLoPerp = BD.COLL_SIC
       end
    end
@@ -438,6 +445,7 @@ end
 
 ---------------------------------------------------------------------
 function FacesBySaw.MakeTwo( Proc, nPhase, nRawId, nPartId, dOvmHead, sCutType, bUpdateIng, bDownHead)
+   EgtOutLog( 'FacesBySaw.MakeTwo', 3)
    -- bUpdateIng : parametro opzionale con default true
    if bUpdateIng == nil then bUpdateIng = true end 
    -- recupero l'ingombro del grezzo di appartenenza
@@ -485,14 +493,14 @@ function FacesBySaw.MakeTwo( Proc, nPhase, nRawId, nPartId, dOvmHead, sCutType,
    vtRef[1] = vtN[1] ^ vtTg
    if vtRef[1] * vtN[2] < 0 then vtRef[1] = - vtRef[1] end
    vtRef[1]:normalize()
-   if abs( vtTg:getZ()) < 0.708 and abs( vtRef[1]:getZ()) < 0.577 and abs( abs( vtRef[1]:getX()) - abs( vtRef[1]:getY())) < 0.1 then
+   if abs( vtTg:getZ()) < 0.708 and abs( vtRef[1]:getZ()) < 0.577 and abs( abs( vtRef[1]:getX()) - abs( vtRef[1]:getY())) < 0.05 then
       vtRef[1] = ptC[1] - ptM
       vtRef[1]:normalize()
    end
    vtRef[2] = vtN[2] ^ vtTg
    if vtRef[2] * vtN[1] < 0 then vtRef[2] = - vtRef[2] end
    vtRef[2]:normalize()
-   if abs( vtTg:getZ()) < 0.708 and abs( vtRef[2]:getZ()) < 0.577 and abs( abs( vtRef[2]:getX()) - abs( vtRef[2]:getY())) < 0.1 then
+   if abs( vtTg:getZ()) < 0.708 and abs( vtRef[2]:getZ()) < 0.577 and abs( abs( vtRef[2]:getX()) - abs( vtRef[2]:getY())) < 0.05 then
       vtRef[2] = ptC[2] - ptM
       vtRef[2]:normalize()
    end
@@ -690,97 +698,49 @@ function FacesBySaw.CalcLeadInOutPerpGeom( ptP1, ptP2, vtV1, vtV2, vtN, dRad, vt
    -- Sistema di riferimento intrinseco al taglio
    local vtX = vtTg ^ vtN
    local frFace = Frame3d( ptP1, vtX, vtTg, vtN)
-   EgtOutLog( 'Vref=' .. tostring( vtRef) .. ' V1=' .. tostring( vtV1) .. ' V2=' .. tostring( vtV2), 3)
+   local bRight = ( vtX * vtRef > 0)
+   EgtOutLog( 'LioPerp --> Vref=' .. tostring( vtRef) .. ' V1=' .. tostring( vtV1) .. ' V2=' .. tostring( vtV2), 3)
+
    -- Versore di attacco e uscita
    local dCos1 = vtV1 * vtRef
    local dCos2 = vtV2 * vtRef
    local vtLio
    if abs( dCos1 - dCos2) < 0.001 then
-      if abs( vtV1:getZ()) < abs( vtV2:getZ()) then
-         vtLio = vtV1
-      else
-         vtLio = vtV2
-      end
+      vtLio = EgtIf( abs( vtV1:getZ()) < abs( vtV2:getZ()), vtV1, vtV2)
    elseif dCos1 > dCos2 then
       vtLio = vtV1
    else
       vtLio = vtV2
    end
-   local bRight = ( vtX * vtLio > 0)
    -- Versore di attacco e uscita nel riferimento intrinseco al taglio
    local vtLioL = Vector3d( vtLio) ;   vtLioL:toLoc( frFace)
-   -- Spostamento punti per effetto dell'extra o della deficienza di taglio
-   ptP1 = ptP1 + vtX * ( EgtIf( bRight, - dCutExtra, dCutExtra))
-   ptP2 = ptP2 + vtX * ( EgtIf( bRight, - dCutExtra, dCutExtra))
-   -- Non va considerata l'uscita dalla faccia sotto, pertanto va abbassata
-   -- 2021/02/26 Abilito anche uscita sotto
-   local b3MyBox = BBox3d( b3Box) ; --b3MyBox:Add( b3MyBox:getMin() - 1000 * Z_AX())
-   -- Attacco
-   local dLiTang = 10000
-   local dLiPerp = 10000
-   local bLiOk, _, vLiPar = EgtLineBoxInters( ptP1, vtLio, b3MyBox)
-   if bLiOk and #vLiPar > 0 then
-      -- con la prima faccia di uscita
-      local dLen = vLiPar[#vLiPar]
-      local ptInt = ptP1 + vtLio * dLen
-      local vtFN = BL.GetBoxFaceNorm( b3MyBox, ptInt, vtLio)
-      EgtOutLog( 'LiFaceNorm=' .. tostring( vtFN), 3)
-      local dAddLen = dRad * ( sqrt( 1 - ( vtN * vtFN) * ( vtN * vtFN)) - abs( vtX * vtFN)) / ( vtLio * vtFN)
-      local dLiLen = dLen + dAddLen
-      EgtOutLog( 'LeadIn Dist=' .. EgtNumToString( dLiLen), 3)
-      dLiTang = - dLiLen * vtLioL:getY()
-      dLiPerp = EgtIf( bRight, dLiLen, - dLiLen) * vtLioL:getX()
-      -- verifico se miglioro calcolando con faccia successiva
-      local b3Mod = BBox3d( b3MyBox) ;  b3Mod:Add( ptInt + 1000 * vtFN)
-      local bLiOk2, _, vLiPar2 = EgtLineBoxInters( ptP1, vtLio, b3Mod)
-      if bLiOk2 and #vLiPar2 > 0 then
-         local dLen2 = vLiPar2[#vLiPar2]
-         local vtFN2 = BL.GetBoxFaceNorm( b3Mod, ptP1 + vtLio * dLen2, vtLio)
-         EgtOutLog( 'LiFaceNorm2=' .. tostring( vtFN2), 3)
-         local dAddLen2 = dRad * ( sqrt( 1 - ( vtN * vtFN2) * ( vtN * vtFN2)) - abs( vtX * vtFN2)) / ( vtLio * vtFN2)
-         local dLiLen2 = dLen2 + dAddLen2
-         EgtOutLog( 'LeadIn Dist2=' .. EgtNumToString( dLiLen2), 3)
-         local dLiTang2 = - dLiLen2 * vtLioL:getY()
-         local dLiPerp2 = EgtIf( bRight, dLiLen2, - dLiLen2) * vtLioL:getX()
-         if dLiLen2 < dLiLen then
-            dLiTang = dLiTang2
-            dLiPerp = dLiPerp2
-         end
-      end
-   end
-   -- Lunghezza di uscita
-   local dLoTang = 10000
-   local dLoPerp = 10000
-   local bLoOk, _, vLoPar = EgtLineBoxInters( ptP2, vtLio, b3MyBox)
-   if bLoOk and #vLoPar > 0 then
-      -- con la prima faccia di uscita
-      local dLen = vLoPar[#vLoPar]
-      local ptInt = ptP2 + vtLio * dLen
-      local vtFN = BL.GetBoxFaceNorm( b3MyBox, ptInt, vtLio)
-      EgtOutLog( 'LoFaceNorm=' .. tostring( vtFN), 3)
-      local dAddLen = dRad * ( sqrt( 1 - ( vtN * vtFN) * ( vtN * vtFN)) - abs( vtX * vtFN)) / ( vtLio * vtFN)
-      local dLoLen = dLen + dAddLen
-      EgtOutLog( 'LeadOut Dist=' .. EgtNumToString( dLoLen), 3)
-      dLoTang = dLoLen * vtLioL:getY()
-      dLoPerp = EgtIf( bRight, dLoLen, - dLoLen) * vtLioL:getX()
-      -- verifico se miglioro calcolando con faccia successiva
-      local b3Mod = BBox3d( b3MyBox) ;  b3Mod:Add( ptInt + 1000 * vtFN)
-      local bLoOk2, _, vLoPar2 = EgtLineBoxInters( ptP2, vtLio, b3Mod)
-      if bLoOk2 and #vLoPar2 > 0 then
-         local dLen2 = vLoPar2[#vLoPar2]
-         local vtFN2 = BL.GetBoxFaceNorm( b3Mod, ptP2 + vtLio * dLen2, vtLio)
-         EgtOutLog( 'LoFaceNorm2=' .. tostring( vtFN2), 3)
-         local dAddLen2 = dRad * ( sqrt( 1 - ( vtN * vtFN2) * ( vtN * vtFN2)) - abs( vtX * vtFN2)) / ( vtLio * vtFN2)
-         local dLoLen2 = dLen2 + dAddLen2
-         EgtOutLog( 'LeadOut Dist2=' .. EgtNumToString( dLoLen2), 3)
-         local dLoTang2 = dLoLen2 * vtLioL:getY()
-         local dLoPerp2 = EgtIf( bRight, dLoLen2, - dLoLen2) * vtLioL:getX()
-         if dLoLen2 < dLoLen then
-            dLoTang = dLoTang2
-            dLoPerp = dLoPerp2
-         end
-      end
-   end
+
+   -- Centro lama a inizio fine percorso
+   local dOffs = ( -dCutExtra + dRad) * EgtIf( bRight, 1, -1)
+   local ptCen1 = ptP1 + dOffs * vtX
+   local ptCen2 = ptP2 + dOffs * vtX
+
+   -- Box del pezzo espansi ortogonalmente alle tre direzioni canoniche X, Y e Z
+   local b3BoxX = BBox3d( b3Box) ; b3BoxX:expand( 0, 1000, 1000)
+   local b3BoxY = BBox3d( b3Box) ; b3BoxY:expand( 1000, 0, 1000)
+   local b3BoxZ = BBox3d( b3Box) ; b3BoxZ:expand( 1000, 1000, 0)
+
+   -- Calcolo elevazione dei due centri
+   local dSawThick = 6
+   EgtCAvSetSawTool( dSawThick, 2 * dRad, dSawThick, 0, 0)
+   local dElev1 = min( EgtCAvToolPosBox( ptCen1 + dSawThick * vtN, vtN, b3BoxX, vtLio),
+                       EgtCAvToolPosBox( ptCen1 + dSawThick * vtN, vtN, b3BoxY, vtLio),
+                       EgtCAvToolPosBox( ptCen1 + dSawThick * vtN, vtN, b3BoxZ, vtLio))
+   EgtOutLog( 'Elev1=' .. EgtNumToString( dElev1), 3)
+   local dLiTang = -dElev1 * vtLioL:getY()
+   local dLiPerp = EgtIf( bRight, dElev1, -dElev1) * vtLioL:getX()
+   local dElev2 = min( EgtCAvToolPosBox( ptCen2 + dSawThick * vtN, vtN, b3BoxX, vtLio),
+                       EgtCAvToolPosBox( ptCen2 + dSawThick * vtN, vtN, b3BoxY, vtLio),
+                       EgtCAvToolPosBox( ptCen2 + dSawThick * vtN, vtN, b3BoxZ, vtLio))
+   EgtOutLog( 'Elev2=' .. EgtNumToString( dElev2), 3)
+   local dLoTang = dElev2 * vtLioL:getY()
+   local dLoPerp = EgtIf( bRight, dElev2, -dElev2) * vtLioL:getX()
+
    return dLiTang, dLiPerp, dLoTang, dLoPerp, vtLio
 end
 
@@ -792,72 +752,35 @@ function FacesBySaw.CalcLeadInOutTangGeom( ptP1, ptP2, vtN, dRad, vtRef, dCutExt
    -- Sistema di riferimento intrinseco al taglio
    local vtX = vtTg ^ vtN
    local frFace = Frame3d( ptP1, vtX, vtTg, vtN)
-   if ( vtX * vtRef < 0) then
-      vtX = - vtX
-   end
-   EgtOutLog( 'Vref=' .. tostring( vtRef) .. ' V1=' .. tostring( vtV1) .. ' V2=' .. tostring( vtV2), 3)
-   -- Spostamento punti per effetto dell'extra o della deficienza di taglio
-   ptP1 = ptP1 - vtX * dCutExtra
-   ptP2 = ptP2 - vtX * dCutExtra
-   -- Non va considerata l'uscita dalla faccia sotto, pertanto va abbassata 
-   -- 2021/02/26 Abilito anche uscita sotto
-   local b3MyBox = BBox3d( b3Box) ; --b3MyBox:Add( b3MyBox:getMin() - 1000 * Z_AX())
-   -- Attacco
-   local dLiTang = 10000
+   local bRight = ( vtX * vtRef > 0)
+   EgtOutLog( 'LioTang --> Vref=' .. tostring( vtRef), 3)
+
+   -- Centro lama a inizio fine percorso
+   local dOffs = ( -dCutExtra + dRad) * EgtIf( bRight, 1, -1)
+   local ptCen1 = ptP1 + dOffs * vtX
+   local ptCen2 = ptP2 + dOffs * vtX
+
+   -- Box del pezzo espansi ortogonalmente alle tre direzioni canoniche X, Y e Z
+   local b3BoxX = BBox3d( b3Box) ; b3BoxX:expand( 0, 1000, 1000)
+   local b3BoxY = BBox3d( b3Box) ; b3BoxY:expand( 1000, 0, 1000)
+   local b3BoxZ = BBox3d( b3Box) ; b3BoxZ:expand( 1000, 1000, 0)
+
+   -- Calcolo elevazione dei due centri
+   local dSawThick = 6
+   EgtCAvSetSawTool( dSawThick, 2 * dRad, dSawThick, 0, 0)
+   local dElev1 = min( EgtCAvToolPosBox( ptCen1 + dSawThick * vtN, vtN, b3BoxX, -vtTg),
+                       EgtCAvToolPosBox( ptCen1 + dSawThick * vtN, vtN, b3BoxY, -vtTg),
+                       EgtCAvToolPosBox( ptCen1 + dSawThick * vtN, vtN, b3BoxZ, -vtTg))
+   EgtOutLog( 'Elev1=' .. EgtNumToString( dElev1), 3)
+   local dLiTang = dElev1
    local dLiPerp = 0
-   local bLiOk, _, vLiPar = EgtLineBoxInters( ptP1, vtTg, b3MyBox)
-   if bLiOk and #vLiPar > 0 then
-      local dLen = vLiPar[1]
-      local ptInt = ptP1 + vtTg * dLen
-      local vtFN = BL.GetBoxFaceNorm( b3MyBox, ptInt, -vtTg)
-      EgtOutLog( 'LiFaceNorm=' .. tostring( vtFN), 3)
-      local dAddLen = dRad * ( sqrt( 1 - ( vtN * vtFN) * ( vtN * vtFN)) - ( vtX * vtFN)) / ( vtTg * vtFN)
-      local dLiLen = dLen + dAddLen
-      EgtOutLog( 'LeadIn Dist=' .. EgtNumToString( dLiLen), 3)
-      dLiTang = - dLiLen
-      -- verifico se miglioro calcolando con faccia successiva
-      local b3Mod = BBox3d( b3MyBox) ;  b3Mod:Add( ptInt + 1000 * vtFN)
-      local bLiOk2, _, vLiPar2 = EgtLineBoxInters( ptP1, vtTg, b3Mod)
-      if bLiOk2 and #vLiPar2 > 0 then
-         local dLen2 = vLiPar2[1]
-         local vtFN2 = BL.GetBoxFaceNorm( b3Mod, ptP1 + vtTg * dLen2, -vtTg)
-         EgtOutLog( 'LiFaceNorm2=' .. tostring( vtFN2), 3)
-         local dAddLen2 = dRad * ( sqrt( 1 - ( vtN * vtFN2) * ( vtN * vtFN2)) - ( vtX * vtFN2)) / ( vtTg * vtFN2)
-         local dLiLen2 = dLen2 + dAddLen2
-         EgtOutLog( 'LeadIn Dist2=' .. EgtNumToString( dLiLen2), 3)
-         if -dLiLen2 < -dLiLen then
-            dLiTang = - dLiLen2
-         end
-      end
-   end
-   -- Lunghezza di uscita
-   local dLoTang = 10000
+   local dElev2 = min( EgtCAvToolPosBox( ptCen2 + dSawThick * vtN, vtN, b3BoxX, vtTg),
+                       EgtCAvToolPosBox( ptCen2 + dSawThick * vtN, vtN, b3BoxY, vtTg),
+                       EgtCAvToolPosBox( ptCen2 + dSawThick * vtN, vtN, b3BoxZ, vtTg))
+   EgtOutLog( 'Elev2=' .. EgtNumToString( dElev2), 3)
+   local dLoTang = dElev2
    local dLoPerp = 0
-   local bLoOk, _, vLoPar = EgtLineBoxInters( ptP2, vtTg, b3MyBox)
-   if bLoOk and #vLoPar > 0 then
-      local dLen = vLoPar[#vLoPar]
-      local ptInt = ptP2 + vtTg * dLen
-      local vtFN = BL.GetBoxFaceNorm( b3MyBox, ptInt, vtTg)
-      EgtOutLog( 'LoFaceNorm=' .. tostring( vtFN), 3)
-      local dAddLen = dRad * ( sqrt( 1 - ( vtN * vtFN) * ( vtN * vtFN)) - ( vtX * vtFN)) / ( vtTg * vtFN)
-      local dLoLen = dLen + dAddLen
-      EgtOutLog( 'LeadOut Dist=' .. EgtNumToString( dLoLen), 3)
-      dLoTang = dLoLen
-      -- verifico se miglioro calcolando con faccia successiva
-      local b3Mod = BBox3d( b3MyBox) ;  b3Mod:Add( ptInt + 1000 * vtFN)
-      local bLoOk2, _, vLoPar2 = EgtLineBoxInters( ptP2, vtTg, b3Mod)
-      if bLoOk2 and #vLoPar2 > 0 then
-         local dLen2 = vLoPar2[#vLoPar2]
-         local vtFN2 = BL.GetBoxFaceNorm( b3Mod, ptP2 + vtTg * dLen2, vtTg)
-         EgtOutLog( 'LoFaceNorm2=' .. tostring( vtFN2), 3)
-         local dAddLen2 = dRad * ( sqrt( 1 - ( vtN * vtFN2) * ( vtN * vtFN2)) - ( vtX * vtFN2)) / ( vtTg * vtFN2)
-         local dLoLen2 = dLen2 + dAddLen2
-         EgtOutLog( 'LeadOut Dist2=' .. EgtNumToString( dLoLen2), 3)
-         if dLoLen2 < dLoLen then
-            dLoTang = dLoLen2
-         end
-      end
-   end
+
    return dLiTang, dLiPerp, dLoTang, dLoPerp
 end
 

LuaLibs/FeatureTopology.lua

@@ -21,30 +21,33 @@ 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 = {}
+   -- essendo la matrice simmetrica a diagonale nulla, ne calcolo solo la metà superiore
    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
+      for j = i + 1, Proc.Fct do
+         _, _, _, 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
+   end
+   -- riempio di conseguenza il resto della matrice
+   for i = 1, Proc.Fct do
+      vAdj[i][i] = 0
+      for j = i + 1, Proc.Fct do
+         vAdj[j][i] = vAdj[i][j]
       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)
+function FeatureTopology.GetFacesWithGivenAdjacencyNumber( Proc, vAdj, nAdj)
+   if not vAdj then vAdj = GetAdjacencyMatrix( Proc) end
+   local nFct = #( vAdj or {}) 
    local vFacesWithGivenAdj = {}
-   for i = 1, Proc.Fct do
+   for i = 1, nFct do
       local nAdjCount = 0
-      for j = 1, Proc.Fct do
+      for j = 1, nFct do
          if vAdj[i][j] and vAdj[i][j] ~= 0 then
             nAdjCount = nAdjCount + 1
          end
@@ -58,11 +61,11 @@ 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 function AreAllAnglesConcaveOrRight( vAdj)
    local bAllConcave, bAllRight = true, true
-   for i = 1, Proc.Fct do
-      for j = 1, Proc.Fct do
+   local nFct = #( vAdj or {})
+   for i = 1, nFct do
+      for j = 1, nFct do
          -- se trovo un angolo convesso restituisco falso e esco subito
          if vAdj[i][j] and vAdj[i][j] > 0 then
             bAllConcave = false
@@ -74,7 +77,7 @@ local function AreAllAnglesConcaveOrRight( Proc)
       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
+   if nFct < 2 or ( nFct == 2 and vAdj[1][2] > 0) then
       return bAllConcave
    else
       return bAllConcave, bAllRight
@@ -165,36 +168,42 @@ 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)
 
-   if not Proc.AffectedFaces then Proc.AffectedFaces = BL.GetProcessAffectedFaces( Proc) end
+   -- la feature deve avere geometria
    if not Proc.Box or Proc.Box:isEmpty() then
       return false
    end
 
-   local bRecognized = false
-   local sFamily
-   local bIsThrough
-   local bAllRightAngles
-   local bIsParallel
-   local sLongName = ''
+   -- se non richiesto (default true, quindi nil vale true), esco
+   if Proc.NeedTopology == false then
+      Proc.Topology = 'SPECIAL'
+      Proc.TopologyLongName = Proc.Topology
+      return true
+   end
 
-   -- SE NON HA TUTTE LE FACCE PIANE RITORNARE NIL!!
+   -- se già calcolato, esco
+   if Proc.Topology then
+      return true
+   end
 
-   local bAllAnglesConcave
-   bAllAnglesConcave, bAllRightAngles = AreAllAnglesConcaveOrRight( Proc)
+   -- calcoli
+   if not Proc.AffectedFaces then Proc.AffectedFaces = BL.GetProcessAffectedFaces( Proc) end
+
+   local vAdj = GetAdjacencyMatrix( Proc)
+   local bAllAnglesConcave, bAllRightAngles = AreAllAnglesConcaveOrRight( vAdj)
    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 vFacesWithOneAdj = FeatureTopology.GetFacesWithGivenAdjacencyNumber( Proc, vAdj, 1)
+   local vFacesWithTwoAdj = FeatureTopology.GetFacesWithGivenAdjacencyNumber( Proc, vAdj, 2)
+   local vFacesWithThreeAdj = FeatureTopology.GetFacesWithGivenAdjacencyNumber( Proc, vAdj, 3)
+   local vFacesWithFourAdj = FeatureTopology.GetFacesWithGivenAdjacencyNumber( Proc, vAdj, 4)
+   local dRawW = b3Raw:getDimY()
+   local dRawH = 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
+   -- assegnazione tipologia
+   local sFamily
+   local bIsThrough
+   if Proc.Fct == 1 and bIsAnyDimensionLongAsPart and bIsFeatureCuttingEntireSection then
       sFamily = 'Cut'
    elseif Proc.Fct == 1 and bIsAnyDimensionLongAsPart then
       sFamily = 'Bevel'
@@ -233,23 +242,24 @@ function FeatureTopology.Classify( Proc, b3Raw)
       sFamily = 'Pocket'
       bIsThrough = false
    end
-   local vFacesParallelToPart = GetFacesParallelToPart( Proc, sFamily, bIsThrough)
-   bIsParallel = ( #vFacesParallelToPart == Proc.Fct)
 
-   if sFamily == 'OUTLINE' or sFamily == 'DRILLING' then
+   -- verifico se facce parallele a quelle della trave
+   local vFacesParallelToPart = GetFacesParallelToPart( Proc, sFamily, bIsThrough)
+   local bIsParallel = ( #vFacesParallelToPart == Proc.Fct)
+
+   -- assegnazioni
+   if sFamily 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
+      Proc.TopologyLongName = GetTopologyLongName( sFamily, bIsThrough, bAllRightAngles, bIsParallel, Proc.Fct)
+      Proc.IsThrough = bIsThrough
+      Proc.AllRightAngles = bAllRightAngles
+      Proc.IsParallel = bIsParallel
+      return true
    else
       Proc.Topology = 'OTHER'
-      Proc.TopologyLongName = 'OTHER'
+      Proc.TopologyLongName = Proc.Topology
+      return false
    end
-
-   return bRecognized
 end
 
 -------------------------------------------------------------------------------------------------------------

LuaLibs/ProcessCut.lua

@@ -23,6 +23,7 @@
 -- 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.
 -- 2023/10/27 In MakeFromTop corretto massimo materiale in caso di lavorazione da sotto.
+-- 2023/12/07 Correzione in Classify in scelta ribaltamento trave quando si è in condizioni downUp.
 
 -- Tabella per definizione modulo
 local ProcessCut = {}
@@ -67,10 +68,22 @@ function ProcessCut.Classify( Proc, b3Raw)
    local dNzLimDwnUp = BL.GetNzLimDownUp( b3Raw)
    local bDownCut = ( vtN:getZ() <= dNzLimDwnUp)
    if bDownCut then
+      -- recupero i dati della lama
+      local sCutType = EgtIf( Proc.Head, 'HeadSide', 'TailSide')
+      local sCutting = ML.FindCutting( sCutType, true, false)
+      local dMaxDepth = 0
+      local dMaxMat = 0
+      if EgtMdbSetCurrMachining( sCutting) then
+         local sTuuid = EgtMdbGetCurrMachiningParam( MCH_MP.TUUID)
+         if EgtTdbSetCurrTool( EgtTdbGetToolFromUUID( sTuuid) or '') then
+            dMaxDepth = EgtTdbGetCurrToolMaxDepth() or dMaxDepth
+            dMaxMat = EgtTdbGetCurrToolParam( MCH_TP.MAXMAT) or dMaxMat
+         end
+      end
       -- calcolo l'ingombro orizzontale della faccia
       local _, DimH, DimV = BL.GetFaceHvRefDim( Proc.Id, 0)
-      -- confronto questo ingombro con il doppio della massima dimensione del DiceCut (impossibile lavorare sotto da sopra se più di 2 tagli)
-      if DimH > 2 * BD.MAX_DIM_DICE then
+      -- confronto questo ingombro con il doppio della massima dimensione del DiceCut (impossibile lavorare sotto da sopra se più di 2 tagli oppure se tipo PF, taglio inclinato in Y e non taglio singolo orizzontale)
+      if DimH > 2 * BD.MAX_DIM_DICE  or ( BD.C_SIMM and abs( vtN:getY()) > 0.1 and dMaxMat < DimH + BD.CUT_EXTRA) then
          return true, true
       end
    end

LuaLibs/ProcessDoubleCut.lua

@@ -2,6 +2,7 @@
 -- Gestione calcolo doppi tagli di lama per Travi
 -- 2022/08/29 Implementata la fresatura dal lato per tagli problematici per PF1250.
 -- 2023/03/22 Eliminata SetOpenSide locale, si usa quella di libreria.
+-- 2023/12/07 Correzione in Classify in scelta ribaltamento trave quando si è in condizioni downUp
 
 -- Tabella per definizione modulo
 local ProcessDoubleCut = {}
@@ -54,10 +55,22 @@ function ProcessDoubleCut.Classify( Proc, b3Raw)
          local dNzLimDwnUp = BL.GetNzLimDownUp( b3Raw, vtN)
          local bDownCut = ( vtN:getZ() <= dNzLimDwnUp)
          if bDownCut then
+            -- recupero i dati della lama
+            local sCutType = EgtIf( Proc.Head, 'HeadSide', 'TailSide')
+            local sCutting = ML.FindCutting( sCutType, true, false)
+            local dMaxDepth = 0
+            local dMaxMat = 0
+            if EgtMdbSetCurrMachining( sCutting) then
+               local sTuuid = EgtMdbGetCurrMachiningParam( MCH_MP.TUUID)
+               if EgtTdbSetCurrTool( EgtTdbGetToolFromUUID( sTuuid) or '') then
+                  dMaxDepth = EgtTdbGetCurrToolMaxDepth() or dMaxDepth
+                  dMaxMat = EgtTdbGetCurrToolParam( MCH_TP.MAXMAT) or dMaxMat
+               end
+            end
             -- calcolo l'ingombro orizzontale della faccia
-            local _, DimH, DimV = BL.GetFaceHvRefDim( Proc.Id, nFac)
-            -- confronto questo ingombro con il doppio della massima dimensione del DiceCut (impossibile lavorare sotto da sopra se più di 2 tagli)
-            if DimH > 2 * BD.MAX_DIM_DICE then
+            local _, DimH = BL.GetFaceHvRefDim( Proc.Id, nFac)
+            -- confronto questo ingombro con il doppio della massima dimensione del DiceCut (impossibile lavorare sotto da sopra se più di 2 tagli oppure se tipo PF, taglio inclinato in Y e non taglio singolo orizzontale)
+            if DimH > 2 * BD.MAX_DIM_DICE or ( BD.C_SIMM and abs( vtN:getY()) > 0.1 and dMaxMat < DimH + BD.CUT_EXTRA) then
                return true, true
             end
          end

LuaLibs/ProcessFreeContour.lua

@@ -35,6 +35,7 @@ local Q_RADIAL_OFFSET = 'Q06' -- d, valido solo per pocket
 function ProcessFreeContour.Identify( Proc)
    return ( ( Proc.Grp == 0 or Proc.Grp == 3 or Proc.Grp == 4)  and  Proc.Prc == 250)
 end
+
 ---------------------------------------------------------------------
 -- Verifica se feature di testa
 function ProcessFreeContour.IsHeadFeature( Proc, b3Raw, dCurrOvmH)

LuaLibs/ProcessLapJoint.lua

@@ -84,6 +84,9 @@
 -- 2023/11/03 Aggiunta MakeStaircaseStep per lavorazione con fresa + lama di gradini scala. Smusso opzionale.
 -- 2023/11/06 In MakeStaircaseStep gestito ritorno.
 -- 2023/11/06 Migliorie per lavorazione con AngularTransmission con FAST.
+-- 2023/11/30 Calcolo elevazione velocizzato e centralizzato tramite la funzione GetFaceElevation.
+-- 2023/12/01 In MakeMoreFaces, nelle OpenPocket, ammesso utensile con diametro fino a 3 volte la dimensione della tasca (era 2).
+-- 2023/12/06 In VerifySideMillAsSaw, se SIDEDEPTH non definita, viene calcolata.
 
 -- Tabella per definizione modulo
 local ProcessLapJoint = {}
@@ -547,12 +550,14 @@ local function VerifySideMillAsSaw( Proc, nAddGrpId, vtN, dDiam, dFacElev)
    local dToolDiamOnSide = 0
    local dThickMillOnSide = 0
    local dMaxDepthOnSide = 0
+   local dTDiamStem = 50
    if EgtMdbSetCurrMachining( sMillingOnSide) then
       local sTuuid = EgtMdbGetCurrMachiningParam( MCH_MP.TUUID)
       if EgtTdbSetCurrTool( EgtTdbGetToolFromUUID( sTuuid) or '') then
          dToolDiamOnSide = EgtTdbGetCurrToolParam( MCH_TP.DIAM) or dToolDiamOnSide
          dThickMillOnSide = EgtTdbGetCurrToolParam( MCH_TP.MAXMAT) or dThickMillOnSide
-         dMaxDepthOnSide = EgtTdbGetCurrToolValInNotes( MCH_TP.USERNOTES, 'SIDEDEPTH', 'd') or dMaxDepthOnSide
+         dTDiamStem = EgtTdbGetCurrToolParam( MCH_TP.STEMDIAM) or dTDiamStem
+         dMaxDepthOnSide = min( EgtTdbGetCurrToolValInNotes( MCH_TP.USERNOTES, 'SIDEDEPTH', 'd') or 999, 0.5 * ( dToolDiamOnSide - dTDiamStem)) or dMaxDepthOnSide
       end
    end
    local bSpecialMillOnSide = ( dDiam > dThickMillOnSide - 10 * GEO.EPS_SMALL and dFacElev < dMaxDepthOnSide + 10 * GEO.EPS_SMALL)
@@ -770,7 +775,7 @@ function ProcessLapJoint.Classify( Proc, b3Raw)
       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]
+         nFacInd = Topology.GetFacesWithGivenAdjacencyNumber( Proc, nil, 4)[1]
          nFacInd2 = nil
          dElev = Proc.Face[ nFacInd + 1].Elevation
       end
@@ -1589,7 +1594,7 @@ local function MakeRoundCleanContour( Proc, nPhase, nRawId, nPartId, b3Raw,
       -- distanza di sicurezza per evitare collisioni
       dCollSic = BL.CalcCollisionSafety( vtN1)
       -- elevazione massima della faccia
-      dMaxElev = BL.GetFaceElevation( Proc.Id, nFacInd, nPartId)
+      dMaxElev = BL.GetFaceElevation( Proc, nFacInd, b3Solid)
       -- ciclo tutta la tabella
       local tPaths = {}
       local ptIniPath
@@ -1634,7 +1639,7 @@ local function MakeRoundCleanContour( Proc, nPhase, nRawId, nPartId, b3Raw,
       -- distanza di sicurezza per evitare collisioni
       dCollSic = BL.CalcCollisionSafety( vtOrtho)
       -- calcolo elevazione dalla faccia trasversale aggiunta
-      local dSurfIntElev = BL.GetOtherFaceElevation( Proc.Id, nSurfInt, 0)
+      local dSurfIntElev = BL.GetFaceElevation( nSurfInt, 0, b3Solid)
       if bDoubleSide then
          dMaxElev = dSurfIntElev
       else
@@ -1645,13 +1650,13 @@ local function MakeRoundCleanContour( Proc, nPhase, nRawId, nPartId, b3Raw,
             -- se non è possibile svuotare completamente da una sola parte
             if dTMaxDepth <= ( dMaxElev + BD.CUT_EXTRA + dCollSic) then
                bDoubleSide = true
-               dMaxElev = BL.GetOtherFaceElevation( Proc.Id, nSurfInt, 0)
+               dMaxElev = BL.GetFaceElevation( nSurfInt, 0, b3Solid)
             else
-               dExtraDepth = dMaxElev - BL.GetOtherFaceElevation( Proc.Id, nSurfInt, 0)
+               dExtraDepth = dMaxElev - BL.GetFaceElevation( nSurfInt, 0, b3Solid)
             end
          -- altrimenti non è stata fatta completamente calcolo la distanza tra faccia aggiunta e profondità superficie
          else
-            dExtraDepth = dMaxElev - BL.GetOtherFaceElevation( Proc.Id, nSurfInt, 0)
+            dExtraDepth = dMaxElev - BL.GetFaceElevation( nSurfInt, 0, b3Solid)
          end
       end
       -- normale alla faccia aggiunta
@@ -4240,7 +4245,7 @@ local function ManageAntiSplintByMill( Proc, nPhase, nRawId, nPartId, b3Raw,
       dCheckDepth = 0.5 * dDepth
    else
       if not nPathInt then
-         dCheckDepth = BL.GetFaceElevation( Proc.Id, nFacInd)
+         dCheckDepth = BL.GetFaceElevation( Proc, nFacInd, b3Solid)
       else
          dCheckDepth = dDepth
       end
@@ -4289,7 +4294,7 @@ local function ManageAntiSplintByMill( Proc, nPhase, nRawId, nPartId, b3Raw,
          end
          local ptIniPath
          local nMaxLen = 0
-         dMaxElevMaster = BL.GetFaceElevation( Proc.Id, nFacInd)
+         dMaxElevMaster = BL.GetFaceElevation( Proc, nFacInd, b3Solid)
          -- ciclo tutta la tabella
          for i = 1, #tFacAdjMain do
             -- le 2 facce di contatto devono essere perpendicolari o non sottosquadra rispetto alla faccia di fondo
@@ -4348,13 +4353,13 @@ local function ManageAntiSplintByMill( Proc, nPhase, nRawId, nPartId, b3Raw,
                -- se non è possibile svuotare completamente da una sola parte
                if dMaxDepth <= ( dMaxElevMaster + BD.CUT_EXTRA + dCollSic) then
                   bDoubleSide = true
-                  dMaxElevMaster = BL.GetOtherFaceElevation( Proc.Id, nSurfInt, 0)
+                  dMaxElevMaster = BL.GetFaceElevation( nSurfInt, 0, b3Solid)
                else
-                  dExtraDepth = dMaxElevMaster - BL.GetOtherFaceElevation( Proc.Id, nSurfInt, 0)
+                  dExtraDepth = dMaxElevMaster - BL.GetFaceElevation( nSurfInt, 0, b3Solid)
                end
             -- altrimenti non è stata fatta completamente calcolo la distanza tra faccia aggiunta e profondità superficie
-            else
-               dExtraDepth = dMaxElevMaster - BL.GetOtherFaceElevation( Proc.Id, nSurfInt, 0)
+            else 
+               dExtraDepth = dMaxElevMaster - BL.GetFaceElevation( nSurfInt, 0, b3Solid)
             end
          end
          vtN1 = Vector3d(vtOrtho)
@@ -4502,7 +4507,7 @@ local function MakeMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead, bSinglePa
    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]
+      nFacInd = Topology.GetFacesWithGivenAdjacencyNumber( Proc, nil, 4)[1]
       nFacInd2 = nil
       dFacElev = Proc.Face[ nFacInd + 1].Elevation
    end
@@ -4897,7 +4902,7 @@ local function MakeMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead, bSinglePa
                         return false, sErr
                      end
                      nFacInd = nFacAdj
-                     dFacElev = BL.GetFaceElevation( Proc.Id, nFacInd)
+                     dFacElev = BL.GetFaceElevation( Proc, nFacInd, b3Solid)
                      ptC, vtN = EgtSurfTmFacetCenter( Proc.Id, nFacInd, GDB_ID.ROOT)
                      rfFac, dH, dV = EgtSurfTmFacetMinAreaRectangle( Proc.Id, nFacInd, GDB_ID.ROOT)
                   end
@@ -4935,7 +4940,7 @@ local function MakeMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead, bSinglePa
                   -- per rifinire gli angoli premio utensile diam 25 o da BD
                   dDiam = min( 2 * dDiam + 10 * GEO.EPS_SMALL, BD.MAXDIAM_POCK_CORNER)
                else
-                  dDiam = 2 * dDiam
+                  dDiam = 3 * dDiam
                end
             end
             sMchFindBackUp = sMchFind
@@ -5159,7 +5164,7 @@ local function MakeMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead, bSinglePa
                            local tvtNx = {}
                            tvtNx[2] = vtN
                            local ptPs = ptC
-                           dFacElev = BL.GetFaceElevation( Proc.Id, nFacInd)
+                           dFacElev = BL.GetFaceElevation( Proc, nFacInd, b3Solid)
                            dCollSic = BL.CalcCollisionSafety( tvtNx[2])
                            local dMachDepth = dFacElev + dCollSic
                            local frFacRec, dFacDim1, dFacDim2 = EgtSurfTmFacetMinAreaRectangle( Proc.Id, nFacInd, GDB_ID.ROOT)
@@ -5181,7 +5186,7 @@ local function MakeMoreFaces( Proc, nPhase, nRawId, nPartId, dOvmHead, bSinglePa
                                  end
                                  vAdj[i] = nCount
                                  if vAdj[i] == Proc.Fct - 1 and i ~= nFacInd + 1 then
-                                    dV = BL.GetFaceElevation( Proc.Id, i - 1)
+                                    dV = BL.GetFaceElevation( Proc, i - 1, b3Solid)
                                     vtN2 = EgtSurfTmFacetNormVersor( Proc.Id, i - 1, GDB_ID.ROOT)
                                  end
                               end

LuaLibs/ProcessLongCut.lua

@@ -31,6 +31,7 @@
 -- 2023/08/01 Ammesso uso lama da sotto fino a N +3deg in verticale.
 -- 2023/09/26 Modificata chiamata a GetFaceWithMostAdj.
 -- 2023/10/24 Migliorata spezzatura taglio passante con due spezzoni
+-- 2023/11/30 Calcolo elevazione velocizzato e centralizzato tramite la funzione GetFaceElevation.
 
 -- Tabella per definizione modulo
 local ProcessLongCut = {}
@@ -959,7 +960,7 @@ function ProcessLongCut.Make( Proc, nPhase, nRawId, nPartId, bCustUseBlade, nCus
       if bForceUseBladeOnNotThruFace and ( bLimXmin or bLimXmax) then
          local nCountMilHead = 0
          -- determino la massima elevazione
-         local dElev = BL.GetFaceElevation( Proc.Id, 0, nPartId)
+         local dElev = BL.GetFaceElevation( Proc, 0, b3Solid)
          -- recupero la lavorazione
          local sMilling = ML.FindMilling( 'Long2Cut_H2', dElev, nil, nil, nil, not bCanUseUnderBlade, bCanUseUnderBlade)
          if not sMilling then
@@ -1209,7 +1210,7 @@ function ProcessLongCut.Make( Proc, nPhase, nRawId, nPartId, bCustUseBlade, nCus
    -- se non è sotto e non uso fresa di fianco: lavorazione Long2Cut
    elseif ( nSide ~= - 1 or BD.DOWN_HEAD) and nUseMillOnSide == 0 then
       -- determino la massima elevazione
-      local dElev = BL.GetFaceElevation( Proc.Id, 0, nPartId)
+      local dElev = BL.GetFaceElevation( Proc, 0, b3Solid)
       -- recupero la lavorazione
       local bDownHead = ( nSide == - 1)
       sMchType = EgtIf( bDownHead, 'Long2Cut_H2', 'Long2Cut')
@@ -1553,7 +1554,7 @@ function ProcessLongCut.Make( Proc, nPhase, nRawId, nPartId, bCustUseBlade, nCus
       local dDistToEnd = dToolDiam / 2
       local dDistToEndDn = dToolDiamDn / 2
       -- calcolo l'elevazione della faccia principale
-      local dFacElev = BL.GetFaceElevation( Proc.Id, 0, nPartId)
+      local dFacElev = BL.GetFaceElevation( Proc, 0, b3Solid)
       -- se fresa di fianco o da sotto calcolo quanto l'utensile può andare vicino all'estremo con l'elevazione della faccia minore del raggio utensile
       if ( nUseMillOnSide <= 1 or nSide == -1) and dFacElev < dToolDiam / 2 then
          dDistToEnd = sqrt( dFacElev * ( dToolDiam - dFacElev))

LuaLibs/ProcessLongDoubleCut.lua

@@ -21,6 +21,7 @@
 -- 2023/10/24 Migliorata spezzatura taglio passante con due spezzoni
 -- 2023/10/25 Se effettivamente un taglio longitudinale e lama non taglia completamente, limito la lavorazione. Altrimenti esco.
 -- 2023/11/24 Aggiunta Q05 per utilizzo lama anche in feature cieche conme per LongCut.
+-- 2023/11/30 Calcolo elevazione velocizzato e centralizzato tramite la funzione GetFaceElevation.
 
 -- Tabella per definizione modulo
 local ProcessLong2Cut = {}
@@ -903,8 +904,8 @@ function ProcessLong2Cut.Make( Proc, nPhase, nRawId, nPartId, bForcedBladeMaster
             if bForceUseBladeOnNotThruFace then
                -- determino la massima elevazione
                local dElev = 0
-               local dFacElev1 = EgtSurfTmFacetElevationInBBox( Proc.Id, tFaceLong[1], b3Solid, GDB_ID.ROOT)
-               local dFacElev2 = EgtSurfTmFacetElevationInBBox( Proc.Id, tFaceLong[2], b3Solid, GDB_ID.ROOT)
+               local dFacElev1 = BL.GetFaceElevation( Proc, tFaceLong[1], b3Solid)
+               local dFacElev2 = BL.GetFaceElevation( Proc, tFaceLong[2], b3Solid)
                dElev = max( dFacElev1, dFacElev2)
                -- recupero la lavorazione
                local sMilling
@@ -1079,8 +1080,8 @@ function ProcessLong2Cut.Make( Proc, nPhase, nRawId, nPartId, bForcedBladeMaster
    elseif ( nSide ~= -1 or BD.DOWN_HEAD or BD.TURN) and nUseMillOnSide == 0 then
       -- determino la massima elevazione
       local dElev = 0
-      local dFacElev1 = EgtSurfTmFacetElevationInBBox( Proc.Id, tFaceLong[1], b3Solid, GDB_ID.ROOT)
-      local dFacElev2 = EgtSurfTmFacetElevationInBBox( Proc.Id, tFaceLong[2], b3Solid, GDB_ID.ROOT)
+      local dFacElev1 = BL.GetFaceElevation( Proc, tFaceLong[1], b3Solid)
+      local dFacElev2 = BL.GetFaceElevation( Proc, tFaceLong[2], b3Solid)
       -- se facce concave e a 90 gradi, prendo l'elevazione minima
       if bOrtho then
          dElev = min( dFacElev1, dFacElev2)
@@ -1450,8 +1451,8 @@ function ProcessLong2Cut.Make( Proc, nPhase, nRawId, nPartId, bForcedBladeMaster
       local dDistToEnd = dToolDiam / 2
       local dDistToEndDn = dToolDiamDn / 2
          -- calcolo l'elevazione della faccia principale
-      local dFacElev1 = BL.GetFaceElevation( Proc.Id, tFaceLong[vOrd[1]], nPartId)
-      local dFacElev2 = BL.GetFaceElevation( Proc.Id, tFaceLong[vOrd[2]], nPartId)
+      local dFacElev1 = BL.GetFaceElevation( Proc, tFaceLong[vOrd[1]], b3Solid)
+      local dFacElev2 = BL.GetFaceElevation( Proc, tFaceLong[vOrd[2]], b3Solid)
       -- se fresa di fianco o da sotto calcolo quanto l'utensile può andare vicino al limite se l'elevazione della faccia è minore del raggio utensile
       if nUseMillOnSide <= 1 or nSide == -1 then
          local dFacElev = max( dFacElev1, dFacElev2)

LuaLibs/ProcessProfCamb.lua

@@ -9,6 +9,7 @@
 -- 2022/05/24 Aggiunta fresatura da sotto su macchine con testa da sotto.
 -- 2022/05/28 Aggiunto calcolo svuotatura da modulo di libreria.
 -- 2022/11/03 Correzione per riconoscimento testa da sotto su fresatura.
+-- 2023/11/30 Calcolo elevazione velocizzato e centralizzato tramite la funzione GetFaceElevation.
 
 -- Tabella per definizione modulo
 local ProcessProfCamb = {}
@@ -285,7 +286,7 @@ function ProcessProfCamb.Make( Proc, nPhase, nRawId, nPartId, dCurrOvmH)
       -- aggiungo piano di sgrossatura e lo lavoro
       local ptStart, vtNP = GetSawCutData( AuxId, vtNF)
       local AddId = EgtSurfTmPlaneInBBox( nAddGrpId, ptStart, vtNP, b3Solid, GDB_RT.GLOB)
-      if AddId and BL.GetFaceElevation( AddId, 0, nPartId) > 20.0 then
+      if AddId and BL.GetFaceElevation( AddId, 0, b3Solid) > 20.0 then
          EgtSetName( AddId, 'AddCut_' .. tostring( Proc.Id))
          EgtSetInfo( AddId, 'TASKID', Proc.TaskId)
          -- se macchina PF o ONE e profilo sopra davanti e pezzo alto applico svuotatura

LuaLibs/ProcessProfConcave.lua

@@ -7,6 +7,7 @@
 -- 2022/05/28 Aggiunto calcolo svuotatura da modulo di libreria.
 -- 2023/05/25 Sistemazione SCC per TURN.
 -- 2023/07/31 Correzione e semplificazione di ModifySideInvertLead per invert.
+-- 2023/11/30 Calcolo elevazione velocizzato e centralizzato tramite la funzione GetFaceElevation.
 
 -- Tabella per definizione modulo
 local ProcessProfConcave = {}
@@ -289,7 +290,7 @@ function ProcessProfConcave.Make( Proc, nPhase, nRawId, nPartId, dCurrOvmH)
       -- aggiungo piano di sgrossatura e lo lavoro
       local ptStart, vtNP = GetSawCutData( AuxId, vtN)
       local AddId = EgtSurfTmPlaneInBBox( nAddGrpId, ptStart, vtNP, b3Solid, GDB_RT.GLOB)
-      if AddId and BL.GetFaceElevation( AddId, 0, nPartId) > 20.0 then
+      if AddId and BL.GetFaceElevation( AddId, 0, b3Solid) > 20.0 then
          EgtSetName( AddId, 'AddCut_' .. tostring( Proc.Id))
          EgtSetInfo( AddId, 'TASKID', Proc.TaskId)
          -- se macchina PF o ONE e profilo sopra davanti e pezzo alto applico svuotatura

LuaLibs/ProcessProfConvex.lua

@@ -7,6 +7,7 @@
 -- 2022/05/28 Aggiunto calcolo svuotatura da modulo di libreria.
 -- 2023/05/25 Sistemazione SCC per TURN.
 -- 2023/07/31 Correzione e semplificazione di ModifySideInvertLead per invert.
+-- 2023/11/30 Calcolo elevazione velocizzato e centralizzato tramite la funzione GetFaceElevation.
 
 -- Tabella per definizione modulo
 local ProcessProfConvex = {}
@@ -289,7 +290,7 @@ function ProcessProfConvex.Make( Proc, nPhase, nRawId, nPartId, dCurrOvmH)
       -- aggiungo piano di sgrossatura e lo lavoro
       local ptStart, vtNP = GetSawCutData( AuxId, vtN)
       local AddId = EgtSurfTmPlaneInBBox( nAddGrpId, ptStart, vtNP, b3Solid, GDB_RT.GLOB)
-      if AddId and BL.GetFaceElevation( AddId, 0, nPartId) > 20.0 then
+      if AddId and BL.GetFaceElevation( AddId, 0, b3Solid) > 20.0 then
          EgtSetName( AddId, 'AddCut_' .. tostring( Proc.Id))
          EgtSetInfo( AddId, 'TASKID', Proc.TaskId)
          -- se macchina PF o ONE e profilo sopra davanti e pezzo alto applico svuotatura

LuaLibs/ProcessProfFront.lua

@@ -5,6 +5,7 @@
 -- 2022/02/02 Aggiunta funzione OnlyChamfer.
 -- 2022/05/24 Aggiunta fresatura da sotto su macchine con testa da sotto.
 -- 2022/05/28 Aggiunto calcolo svuotatura da modulo di libreria.
+-- 2023/11/30 Calcolo elevazione velocizzato e centralizzato tramite la funzione GetFaceElevation.
 
 -- Tabella per definizione modulo
 local ProcessProfFront = {}
@@ -236,7 +237,7 @@ function ProcessProfFront.Make( Proc, nPhase, nRawId, nPartId, dCurrOvmH)
       -- aggiungo piano di sgrossatura e lo lavoro
       local ptStart, vtNP = GetSawCutData( AuxId, vtN)
       local AddId = EgtSurfTmPlaneInBBox( nAddGrpId, ptStart, vtNP, b3Solid, GDB_RT.GLOB)
-      if AddId and BL.GetFaceElevation( AddId, 0, nPartId) > 20.0 then
+      if AddId and BL.GetFaceElevation( AddId, 0, b3Solid) > 20.0 then
          EgtSetName( AddId, 'AddCut_' .. tostring( Proc.Id))
          EgtSetInfo( AddId, 'TASKID', Proc.TaskId)
          -- se macchina PF o ONE e profilo sopra davanti e pezzo alto applico svuotatura

LuaLibs/ProcessProfHead.lua

@@ -7,6 +7,7 @@
 -- 2022/05/24 Aggiunta fresatura da sotto su macchine con testa da sotto.
 -- 2022/05/28 Aggiunto calcolo svuotatura da modulo di libreria.
 -- 2023/07/31 Correzione e semplificazione di ModifySideInvertLead per invert.
+-- 2023/11/30 Calcolo elevazione velocizzato e centralizzato tramite la funzione GetFaceElevation.
 
 -- Tabella per definizione modulo
 local ProcessProfHead = {}
@@ -271,7 +272,7 @@ function ProcessProfHead.Make( Proc, nPhase, nRawId, nPartId, dCurrOvmH)
       -- aggiungo piano di sgrossatura e lo lavoro
       local ptStart, vtNP = GetSawCutData( AuxId, vtN)
       local AddId = EgtSurfTmPlaneInBBox( nAddGrpId, ptStart, vtNP, b3Solid, GDB_RT.GLOB)
-      if AddId and BL.GetFaceElevation( AddId, 0, nPartId) > 20.0 then
+      if AddId and BL.GetFaceElevation( AddId, 0, b3Solid) > 20.0 then
          EgtSetName( AddId, 'AddCut_' .. tostring( Proc.Id))
          EgtSetInfo( AddId, 'TASKID', Proc.TaskId)
          -- se macchina PF o ONE e profilo sopra davanti e pezzo alto applico svuotatura

LuaLibs/ProcessRoundArch.lua

@@ -1,4 +1,4 @@
--- ProcessRoundArch.lua by Egaltech s.r.l. 2023/05/18
+-- ProcessRoundArch.lua by Egaltech s.r.l. 2023/12/07
 -- Gestione calcolo archi per Travi
 
 -- Tabella per definizione modulo
@@ -118,10 +118,13 @@ function ProcessRoundArch.Make( Proc, nPhase, nRawId, nPartId, dOvmHead)
    end
    -- recupero i dati dell'utensile
    local dMaxDepth = 0
+   local bCW = true
    if EgtMdbSetCurrMachining( sMilling) then
       local sTuuid = EgtMdbGetCurrMachiningParam( MCH_MP.TUUID)
       if EgtTdbSetCurrTool( EgtTdbGetToolFromUUID( sTuuid) or '') then
          dMaxDepth = EgtTdbGetCurrToolMaxDepth() or dMaxDepth
+         local dSpeed = EgtMdbGetCurrMachiningParam( MCH_MP.SPEED) or 0
+         bCW = ( dSpeed >= 0)
       end
    end
    -- ne verifico la lunghezza per eventuale spezzatura e lavorazione in doppio
@@ -180,6 +183,10 @@ function ProcessRoundArch.Make( Proc, nPhase, nRawId, nPartId, dOvmHead)
          end
       end
    end
+   -- se lavorazione in doppio senza codolo, aggiungo sovrapposizione
+   if nDouble == 2 and dDimStrip < 10 * GEO.EPS_SMALL then
+      dDepth = dDepth + BD.MILL_OVERLAP / 2
+   end
    -- se utensile orizzontale verso Y+, non in doppio e codolo da lasciare, devo invertire per lavorare sempre da Y-
    if vtExtr:getY() > 0.707 and  nDouble == 1 and bStripOnSide then
       bToolInv = true
@@ -188,7 +195,7 @@ function ProcessRoundArch.Make( Proc, nPhase, nRawId, nPartId, dOvmHead)
    for i = 1, nStep do
       for j = 1, nDouble do
          -- inserisco la lavorazione
-         local sName = 'Free_' .. ( EgtGetName( Proc.Id) or tostring( Proc.Id))
+         local sName = 'Arch_' .. ( EgtGetName( Proc.Id) or tostring( Proc.Id))
          local nMchId = EgtAddMachining( sName, sMilling)
          if not nMchId then
             local sErr = 'Error adding machining ' .. sName .. '-' .. sMilling
@@ -212,20 +219,19 @@ function ProcessRoundArch.Make( Proc, nPhase, nRawId, nPartId, dOvmHead)
             EgtSetMachiningParam( MCH_MP.TOOLINVERT, true)
          end
          -- se seconda passata, inverto direzione di lavoro
-         if j == 2 then
-            EgtSetMachiningParam( MCH_MP.INVERT, true)
+         if ( j == 1 and not bToolInv) or ( j == 2 and bToolInv) then
+            EgtSetMachiningParam( MCH_MP.INVERT, EgtIf( bCW, true, false))
+         else
+            EgtSetMachiningParam( MCH_MP.INVERT, EgtIf( bCW, false, true))
          end
          -- assegno affondamento
          EgtSetMachiningParam( MCH_MP.DEPTH, dDepth)
-         -- assegno lato di lavoro
-         if bToolInv then
-            EgtSetMachiningParam( MCH_MP.WORKSIDE, MCH_MILL_WS.LEFT)
-         else
-            EgtSetMachiningParam( MCH_MP.WORKSIDE, MCH_MILL_WS.RIGHT)
-         end
+         -- assegno lato di lavoro (per avere lavorazione climb)
+         EgtSetMachiningParam( MCH_MP.WORKSIDE, EgtIf( bCW, MCH_MILL_WS.LEFT, MCH_MILL_WS.RIGHT))
          -- posizione braccio porta testa
          local nSCC = MCH_SCC.NONE
          if not BD.C_SIMM and not BD.TURN then
+            nSCC = MCH_SCC.ADIR_XP
             if Proc.Head then
                nSCC = MCH_SCC.ADIR_XP
             elseif Proc.Tail then

LuaLibs/ProcessVariant.lua

@@ -1,6 +1,7 @@
 -- ProcessVariant.lua by Egaltech s.r.l. 2023/11/08
 -- Gestione calcolo Feature Custom (Variant) per Travi
 -- 2023/11/08 Creazione modulo.
+-- 2023/11/30 Calcolo elevazione velocizzato e centralizzato tramite la funzione GetFaceElevation.
 
 -- Tabella per definizione modulo
 local ProcessVariant = {}
@@ -72,6 +73,7 @@ end
 ---------------------------------------------------------------------
 -- Applicazione della lavorazione
 local function MakePocket( Proc, nRawId, b3Raw, nPartId)
+   local b3Solid = EgtGetBBoxGlob( EgtGetFirstNameInGroup( nPartId, 'Box') or GDB_ID.NULL, GDB_BB.STANDARD)
    -- recupero gli indici delle facce da lavorare
    local vFace = EgtGetInfo( Proc.Id, 'Faces', 'vi')
    if not vFace or #vFace == 0 then
@@ -83,7 +85,7 @@ local function MakePocket( Proc, nRawId, b3Raw, nPartId)
    local dTotDepth = 0 
    local vDepth = {}
    for i = 1, #vFace do
-      vDepth[i] = BL.GetFaceElevation( Proc.Id, vFace[1], nPartId)
+      vDepth[i] = BL.GetFaceElevation( Proc, vFace[1], b3Solid)
       dTotDepth = max( dTotDepth, vDepth[i])
    end
    local vtN = EgtSurfTmFacetNormVersor( Proc.Id, vFace[1], GDB_ID.ROOT)

UpdateLog.txt

@@ -1,5 +1,22 @@
            ==== Beam Update Log ====
 
+Versione 2.5l3 (08/01/2023)
+- Modif : modifiche per centrare i pezzi in Y sulla tavola (da flag BeamData.CENTER_BEAM)
+- Modif : resi più veloci i calcoli di FeatureTopology
+- Modif : migliorato calcolo attacco/uscite lame in FacesBySaw.
+
+Versione 2.5l2 (20/12/2023)
+- Modif : in Cut e DoubleCut modifica in scelta ribaltamento trave quando si è in condizioni downUp
+- Modif : migliorato e esteso il funzionamento di SideMillAsSaw
+- Modif : in LapJoint con più facce, nelle OpenPocket, ammesso utensile con diametro fino a 3 volte la dimensione della tasca (era 2)
+- Modif : in RoundArch aggiunta gestione senso di rotazione fresa per lavorare sempre in concordanza
+- Fixed : correzione a GetFaceWithMostAdj per calcoli con facce che fanno ombra
+
+Versione 2.5l1 (01/12/2023)
+- Added : in LongDoubleCut aggiunta Q05 per utilizzo lama anche in feature cieche conme per LongCut
+- Modif : Migliorato e velocizzato il calcolo elevazione (EgtSurfTmFacetElevationInBBox)
+- Fixed : in MakeTwo di FacesBySaw raffinamento calcolo vtRef per casi dubbi.
+
 Versione 2.5k3 (28/11/2023)
 - Added : in Long2Cut gestite lavorazioni lama con Q05 come in LongCut.
 
@@ -15,26 +32,26 @@ Versione 2.5k1 (02/11/2023)
 
 Versione 2.5j1 (23/10/2023)
 - Added : in Process (lancio da EgtCAM5 quindi uso praticamente solo interno) aggiunta possibilità di tenere ordine pezzi in barra come da selezione
+- Modif : in Topology aggiunti casi al check solo direzione principale
+- Modif : nella ricerca lavorazione da usare controllo di utensile attivo sostituito con controllo utensile presente nel setup corrente
 - Fixed : piccola correzione a Classify di FeatureTopology per casi senza geometria
 - Fixed : in HeadCut e Split corretto calcolo allungamenti/accorciamenti per evitare lunghezze del percorso negative
-- Fixed : in Lapjoint -> MakePocket aggiunto messaggio in caso si rovini il pezzo successivo
-- Modif : in Topology aggiunti casi al check solo direzione principale
-- Modif : nella ricerca lavorazione da usare controllo di utensile attivo sostituito con controllo utensile presente nel setup corrente.
+- Fixed : in Lapjoint -> MakePocket aggiunto messaggio in caso si rovini il pezzo successivo.
 
 Versione 2.5i2 (27/09/2023)
 - Added : aggiunto riconoscimento sottosquadro da facce non adiacenti
+- Modif : in lavorazione tenoni per Turn si assegna SCC per privilegiare accesso dal lato corto della trave.
 - 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].
+- Modif : in Cut abbassato a 590 mm il limite per convertire in LongCut [Ticket #1448]
+- Fixed : in LapJoint gestito correttamente il ritorno nil di GetUShapeWidth [Ticket #1354].
 
 Versione 2.5h2 (11/08/2023)
+- Modif : in LapJoint e Mortise modificata posizione braccio per FAST quando vicino alla coda della trave
 - 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.
+- Fixed : in RidgeLap corretto ingombro lavorazione in testa e coda quando rivolta verso il basso.
 
 Versione 2.5h1 (07/08/2023)
 - Fixed : tagli doppi di lato non effettuati se macchina tipo PF e pezzo alto [Ticket #1400]

Version.lua

@@ -1,6 +1,6 @@
--- Version.lua by Egaltech s.r.l. 2023/11/16
+-- Version.lua by Egaltech s.r.l. 2024/01/08
 -- Gestione della versione di Beam
 
 NAME = 'Beam'
-VERSION = '2.5k5'
-MIN_EXE = '2.5c1'
+VERSION = '2.5l3'
+MIN_EXE = '2.5l3'