-- ProcessLongDoubleCut.lua by Egaltech s.r.l. 2020/02/02
-- Gestione calcolo doppio taglio longitudinale per Travi

-- Tabella per definizione modulo
local ProcessLong2Cut = {}

-- Include
require( 'EgtBase')
local BL = require( 'BeamLib')

EgtOutLog( ' ProcessLongDoubleCut started', 1)

-- Dati
local BD = require( 'BeamData')
local ML = require( 'MachiningLib')

---------------------------------------------------------------------
-- Riconoscimento della feature
function ProcessLong2Cut.Identify( Proc)
   return ( Proc.Grp == 0  and  Proc.Prc == 12)
end

---------------------------------------------------------------------
-- Classificazione della feature
function ProcessLong2Cut.Classify( Proc)
   return true, false
end

---------------------------------------------------------------------
-- Classificazione della feature, calcolo il numero di facce longitudinali e limite.
-- Ritorna due valori :
--  il primo indica il numero di facce longitudinali
--  il secondo indica quante e quali facce limite sono presenti con questa codifica
--    0: nessuna faccia limite; 1: faccia/e limite a destra (X-); 2: faccia/e limite a sinistra (X+); 3: due facce limite.
function ProcessLong2Cut.IdentifyFaces( Proc, bDefineFaces)
   -- verifico le normali delle facce
   local nFaceLong = 0
   local nFaceLimit = 0
   local tFaceLong = {}
   local tFaceLimit = {}
   local tptC = {}
   local tvtN = {}
   local nFacetCnt = EgtSurfTmFacetCount( Proc.Id)
   local nLong = 0
   local nLimit = 0
   local bLimitPos
   local bLimitNeg

   for i = 1, nFacetCnt do
      local vtN = EgtSurfTmFacetNormVersor( Proc.Id, i-1, GDB_ID.ROOT)
      -- se è una faccia limite a sinistra
      if vtN:getX() > GEO.EPS_SMALL then
         if not bLimitPos then 
            nFaceLimit = nFaceLimit + 2
            bLimitPos = true
         end
         -- se devo restituire una tabella con le facce limite
         if bDefineFaces then
            nLimit = nLimit + 1
            tFaceLimit[nLimit] = i-1
         end
      -- se è una faccia limite a destra
      elseif vtN:getX() < -GEO.EPS_SMALL then
         if not bLimitNeg then 
            nFaceLimit = nFaceLimit + 1
            bLimitNeg = true
         end
         -- se devo restituire una tabella con le facce limite
         if bDefineFaces then
            nLimit = nLimit + 1
            tFaceLimit[nLimit] = i-1
         end
      -- altrimenti è una faccia longitudinale
      else
         nFaceLong = nFaceLong + 1
         -- se devo restituire una tabella con le facce longitudinali
         if bDefineFaces then
            nLong = nLong + 1
            tFaceLong[nLong] = i-1
            tptC[nLong], tvtN[nLong] = EgtSurfTmFacetCenter( Proc.Id, (i-1), GDB_ID.ROOT)
         end
      end
   end
   return nFaceLong, nFaceLimit, tFaceLong, tFaceLimit, tptC, tvtN
end

---------------------------------------------------------------------
local function MakeSideFace( nId, nFac, nSide, sMilling, dToolDiam)
   -- inserisco la lavorazione
   local sNameF = 'L2C_' .. ( EgtGetName( nId) or tostring( nId)) .. '_' .. tostring( nFac)
   local nMchFId = EgtAddMachining( sNameF, sMilling)
   if not nMchFId then
      local sErr = 'Error adding machining ' .. sNameF .. '-' .. sMilling
      EgtOutLog( sErr)
      return false, sErr
   end
   -- aggiungo geometria
   EgtSetMachiningGeometry( {{ nId, nFac}})
   -- lato di lavoro e inversione
   EgtSetMachiningParam( MCH_MP.INVERT, true)
   EgtSetMachiningParam( MCH_MP.WORKSIDE, MCH_MILL_WS.LEFT)
   -- uso della faccia
   local nFaceUse = MCH_MILL_FU.PARAL_DOWN
   if nSide == -2 then
      nFaceUse = MCH_MILL_FU.PARAL_BACK
   elseif nSide == 2 then
      nFaceUse = MCH_MILL_FU.PARAL_FRONT
   end
   EgtSetMachiningParam( MCH_MP.FACEUSE, nFaceUse)
   -- annullo offset radiale
   EgtSetMachiningParam( MCH_MP.OFFSR, 0)
   -- attacco e uscita
   EgtSetMachiningParam( MCH_MP.LIPERP, 0)
   EgtSetMachiningParam( MCH_MP.LITANG, dToolDiam / 2 + 30)
   EgtSetMachiningParam( MCH_MP.LOPERP, 0)
   EgtSetMachiningParam( MCH_MP.LOTANG, dToolDiam / 2 + 30)
   -- eseguo
   if not EgtApplyMachining( true, false) then
      local _, sErr = EgtGetLastMachMgrError()
      EgtSetOperationMode( nMchFId, false)
      return false, sErr
   end
end

---------------------------------------------------------------------
-- Applicazione della lavorazione
local function MakeByPocketing( Proc, nPhase, nRawId, nPartId, nFacInd, dFacElev)
   -- cerco la svuotatura opportuna
   local sPocketing = ML.FindPocketing( 'OpenPocket', Proc.Box:getDimX())
   if not sPocketing then
      local sErr = 'Error on process ' .. tostring( Proc.Id) .. ' pocketing not found in library'
      EgtOutLog( sErr)
      return false, sErr
   end
   -- recupero i dati dell'utensile
   local dMaxDepth = 0
   if EgtMdbSetCurrMachining( sPocketing) then
      local sTuuid = EgtMdbGetCurrMachiningParam( MCH_MP.TUUID)
      if EgtTdbSetCurrTool( EgtTdbGetToolFromUUID( sTuuid) or '') then
         dMaxDepth = ( EgtTdbGetCurrToolMaxDepth() or dMaxDepth)
      end
   end
   -- inserisco la lavorazione
   local sName = 'Pock_' .. ( EgtGetName( Proc.Id) or tostring( Proc.Id))
   local nMchFId = EgtAddMachining( sName, sPocketing)
   if not nMchFId then
      local sErr = 'Error adding machining ' .. sName .. '-' .. sPocketing
      EgtOutLog( sErr)
      return false, sErr
   end
   -- aggiungo geometria
   EgtSetMachiningGeometry( {{ Proc.Id, nFacInd}})
   -- imposto uso faccia
   EgtSetMachiningParam( MCH_MP.FACEUSE, MCH_MILL_FU.ORTHO_CONT)
   -- imposto attacco per tasca aperta
   EgtSetMachiningParam( MCH_MP.SUBTYPE, MCH_POCK_SUB.SPIRALIN)
   -- se elevazione superiore a massimo affondamento della fresa, riduco opportunamente
   local sWarn
   if dFacElev > dMaxDepth + 10 * GEO.EPS_SMALL then
      EgtSetMachiningParam( MCH_MP.DEPTH, dMaxDepth - dFacElev)
      dFacElev = dMaxDepth
      sWarn = 'Warning in process ' .. tostring( Proc.Id) .. ' : elevation bigger than max tool depth'
      EgtOutLog( sWarn)
   end
   -- imposto elevazione e dichiaro non si generano sfridi per VMill
   local sNotes = 'MaxElev=' .. EgtNumToString( dFacElev, 1) .. ';'
   sNotes = sNotes .. 'VMRS=0;'
   EgtSetMachiningParam( MCH_MP.USERNOTES, sNotes)
   -- eseguo
   if not EgtApplyMachining( true, false) then
      -- provo ad allargare leggermente la tasca
      EgtSetMachiningParam( MCH_MP.OFFSR, -0.1)
      if not EgtApplyMachining( true, false) then
         local _, sErr = EgtGetLastMachMgrError()
         EgtSetOperationMode( nMchFId, false)
         return false, sErr
      end
   end
   return true, sWarn
end

---------------------------------------------------------------------
-- Applicazione della lavorazione
function ProcessLong2Cut.Make( Proc, nPhase, nRawId, nPartId)
   -- dati delle facce
   local nFaceLimit = 0
   local tFaceLong = {}
   local tFaceLimit = {}
   local ptC = {}
   local vtN = {}
   -- carico i dati delle face già inserite nelle opportune tabelle
   _, nFaceLimit, tFaceLong, tFaceLimit, ptC, vtN = ProcessLong2Cut.IdentifyFaces( Proc, true)

   local dLen = EgtGetBBoxGlob( Proc.Id, GDB_BB.STANDARD):getDimX()
   -- verifico posizione (+1=sopra, -1=sotto, 0=sui lati)
   local nSide = 0
   if vtN[1]:getZ() > -10 * GEO.EPS_SMALL and vtN[2]:getZ() > -10 * GEO.EPS_SMALL then
      nSide = 1
   elseif vtN[1]:getZ() < -10 * GEO.EPS_SMALL and vtN[2]:getZ() < -10 * GEO.EPS_SMALL then
      nSide = -1
   end
   -- angolo diedro per stabilire se taglio convesso
   local bInt, ptP1, ptP2, dAng = EgtSurfTmFacetsContact( Proc.Id, tFaceLong[1], tFaceLong[2], GDB_ID.ROOT)
   local bConvex
   local bOrtho
   local ptM
   if bInt then
      bConvex = ( dAng >= 0)
      bOrtho = ( abs( dAng + 90) < 1)
      ptM = ( ptP1 + ptP2) / 2
   else
      bConvex = true
      bOrtho = false
      ptM = ( ptC[1] + ptC[2]) / 2
   end
   local ptRef = ( ptC[1] + ptC[2]) / 2
   -- analisi del taglio
   local vOrd = {}
   local vFaceUse = {}
   if nSide == 1 then
      vOrd = EgtIf( ptC[1]:getY() < ptRef:getY(), { 1, 2}, { 2, 1})
      vFaceUse = { BL.GetNearestOrthoOpposite( ptC[1] - ptM), BL.GetNearestOrthoOpposite( ptC[2] - ptM)}
   elseif nSide == -1 then
      vOrd = EgtIf( ptC[1]:getY() < ptM:getY(), { 1, 2}, { 2, 1})
      vFaceUse = { BL.GetNearestParalOpposite( ptC[1] - ptM), BL.GetNearestParalOpposite( ptC[2] - ptM)}
   else
      local bFront = ( vtN[1]:getY() < 0)
      if bFront then
         vOrd = EgtIf( ptC[1]:getZ() < ptM:getZ(), { 1, 2}, { 2, 1})
         vFaceUse = { BL.GetNearestOrthoOpposite( ptC[1] - ptM), BL.GetNearestOrthoOpposite( ptC[2] - ptM)}
      else
         vOrd = EgtIf( ptC[1]:getZ() < ptM:getZ(), { 2, 1}, { 1, 2})
         vFaceUse = { BL.GetNearestOrthoOpposite( ptC[1] - ptM), BL.GetNearestOrthoOpposite( ptC[2] - ptM)}
      end
   end
   local vWidth = {}
   _, _, vWidth[1] = BL.GetFaceHvRefDim( Proc.Id, tFaceLong[1])
   _, _, vWidth[2] = BL.GetFaceHvRefDim( Proc.Id, tFaceLong[2])

   local sWarn
   -- Se ho 2 face longitudinali e non limitato da altre facce e da sopra e richiesto con doppio taglio di lama
   if #tFaceLong == 2 and Proc.Fct == 2 and nSide == 1 and EgtGetInfo( Proc.Id, 'Q01', 'i') == 1 then
      -- recupero la lavorazione
      local sCutting = ML.FindCutting( 'HeadSide')
      if not sCutting then
         local sErr = 'Error on process ' .. tostring( Proc.Id) .. ' sawing not found in library'
         EgtOutLog( sErr)
         return false, sErr
      end
      -- recupero i dati dell'utensile
      local dToolDiam = 0
      local dMaxDepth = 0
      local dToolThick = 0
      if EgtMdbSetCurrMachining( sCutting) then
         local sTuuid = EgtMdbGetCurrMachiningParam( MCH_MP.TUUID)
         if EgtTdbSetCurrTool( EgtTdbGetToolFromUUID( sTuuid) or '') then
            dToolDiam = EgtTdbGetCurrToolParam( MCH_TP.DIAM) or dToolDiam
            dMaxDepth = EgtTdbGetCurrToolMaxDepth() or dMaxDepth
            dToolThick = EgtTdbGetCurrToolParam(MCH_TP.THICK) or dToolThick
         end
      end

      -- determino numero di parti
      local dStartAccDist = BD.LONGCUT_ENDLEN
      local dEndAccDist = BD.LONGCUT_ENDLEN
      local nC = ceil( ( dLen - dStartAccDist - dEndAccDist) / BD.LONGCUT_MAXLEN)
      local dC = 0
      if nC > 0 then
         dC = dLen / ( nC + 2)
         dStartAccDist = dC
         dEndAccDist = dC
         nC = nC + 2
      else
         if dLen > min( dStartAccDist, dEndAccDist) then
            nC = 2
            dStartAccDist = dLen/2
            dEndAccDist = dStartAccDist
         else
            nC = 1
            dStartAccDist = 0
            dEndAccDist = 0
         end
      end

      local nM = 0
      -- se convesso lavoro ogni faccia in due metà lasciando attaccata la parte centrale
      if bConvex then
         -- si percorrono i lati alto e basso della faccia
         for i = 1, nC do
            -- Posizione braccio portatesta
            local nSCC = EgtIf( ( i == 1 or i == nC - 1), MCH_SCC.ADIR_XP, MCH_SCC.ADIR_XM)
            -- ciclo sulle due facce
            for j = 1, #vOrd do
               -- determino se lavorazione da davanti o da dietro
               local bFront = ( vtN[vOrd[j]]:getY() < 0)
               -- determino l'utilizzo della faccia
               local nFaceUse  = EgtIf( abs( vtN[vOrd[j]]:getY()) > 0.01, MCH_MILL_FU.ORTHO_TOP, EgtIf( bFront, MCH_MILL_FU.ORTHO_BACK, MCH_MILL_FU.ORTHO_FRONT))
               local nFaceUse2 = EgtIf( abs( vtN[vOrd[j]]:getY()) > 0.01, MCH_MILL_FU.ORTHO_DOWN, EgtIf( bFront, MCH_MILL_FU.ORTHO_FRONT, MCH_MILL_FU.ORTHO_BACK))
               -- ciclo sulle passate
               local dOffset = ( vWidth[vOrd[j]] + BD.DIM_STRIP_SMALL) / 2 ;
               local dLioTang = 0
               local dLioPerp = ( vWidth[vOrd[j]] - BD.DIM_STRIP_SMALL) / 2 + BD.CUT_SIC ;
               for k = 1, 2 do
                  -- inserisco le parti di lavorazione
                  nM = nM + 1
                  local sNameF = 'L2C_' .. ( EgtGetName( Proc.Id) or tostring( Proc.Id)) .. '_' .. tostring( nM) .. '_' .. tostring(j)
                  local nMchFId = EgtAddMachining( sNameF, sCutting)
                  if not nMchFId then
                     local sErr = 'Error adding machining ' .. sNameF .. '-' .. sCutting
                     EgtOutLog( sErr)
                     return false, sErr
                  end
                  -- aggiungo geometria
                  EgtSetMachiningGeometry( {{ Proc.Id, vOrd[j]-1}})
                  -- limito opportunamente la lavorazione
                  local dSal = EgtIf( i == nC, 0, - dEndAccDist - ( nC - i - 1) * dC)
                  local dEal = EgtIf( i == 1, 0, - dStartAccDist - ( i - 2) * dC)
                  if ( not bFront and k == 1) or ( bFront and k == 2) then
                     dSal, dEal = dEal, dSal
                  end
                  EgtSetMachiningParam( MCH_MP.STARTADDLEN, dSal)
                  EgtSetMachiningParam( MCH_MP.ENDADDLEN, dEal)
                  -- imposto offset radiale
                  EgtSetMachiningParam( MCH_MP.OFFSR, dOffset)
                  -- imposto attacco/uscita
                  EgtSetMachiningParam( MCH_MP.LITANG, dLioTang)
                  EgtSetMachiningParam( MCH_MP.LIPERP, dLioPerp)
                  EgtSetMachiningParam( MCH_MP.LOTANG, dLioTang)
                  EgtSetMachiningParam( MCH_MP.LOPERP, dLioPerp)
                  -- imposto posizione braccio porta testa per non ingombrare agli estremi
                  EgtSetMachiningParam( MCH_MP.SCC, nSCC)
                  -- imposto uso della faccia
                  EgtSetMachiningParam( MCH_MP.FACEUSE, EgtIf( k == 1, nFaceUse, nFaceUse2))
                  -- eseguo
                  if not EgtApplyMachining( true, false) then
                     local _, sErr = EgtGetLastMachMgrError()
                     EgtSetOperationMode( nMchFId, false)
                     return false, sErr
                  end
               end
            end
         end
      -- altrimenti concavo
      else
         -- si percorrono i lati alto e basso della faccia
         for i = 1, nC do
            -- Posizione braccio portatesta
            local nSCC = EgtIf( ( i == 1 or i == nC - 1), MCH_SCC.ADIR_XP, MCH_SCC.ADIR_XM)
            -- ciclo sulle due facce
            for j = 1, #vOrd do
               -- determino se lavorazione da davanti o da dietro
               local bFront = ( vtN[vOrd[j]]:getY() < 0)
               -- determino l'utilizzo della faccia
               local nFaceUse  = EgtIf( abs( vtN[vOrd[j]]:getY()) > 0.01, MCH_MILL_FU.ORTHO_TOP, EgtIf( bFront, MCH_MILL_FU.ORTHO_BACK, MCH_MILL_FU.ORTHO_FRONT))
               local nFaceUse2 = EgtIf( abs( vtN[vOrd[j]]:getY()) > 0.01, MCH_MILL_FU.ORTHO_DOWN, EgtIf( bFront, MCH_MILL_FU.ORTHO_FRONT, MCH_MILL_FU.ORTHO_BACK))
               -- ciclo sulle passate
               local dOffset
               local dLioPerp
               if j == 1 then -- il primo taglio lo faccio completo
                  dOffset = 0
                  dLioPerp = vWidth[vOrd[j]] + BD.CUT_SIC
               else -- il secondo ridotto della distanza minima e della componente spessore della lama
                  dOffset = BD.DIM_STRIP_SMALL + ((dToolThick* vtN[vOrd[1]]) - (dToolThick* vtN[vOrd[1]]) * vtN[vOrd[2]] * vtN[vOrd[2]]):len()
                  dLioPerp = vWidth[vOrd[j]] - dOffset + BD.CUT_SIC
               end
               local dLioTang = 0
               
               local sNameF = 'L2C_' .. ( EgtGetName( Proc.Id) or tostring( Proc.Id)) .. '_' .. tostring(j)
               local nMchFId = EgtAddMachining( sNameF, sCutting)
               if not nMchFId then
                  local sErr = 'Error adding machining ' .. sNameF .. '-' .. sCutting
                  EgtOutLog( sErr)
                  return false, sErr
               end
               -- aggiungo geometria
               EgtSetMachiningGeometry( {{ Proc.Id, vOrd[j]-1}})
               -- limito opportunamente la lavorazione
               local dSal = EgtIf( i == nC, 0, - dEndAccDist - ( nC - i - 1) * dC)
               local dEal = EgtIf( i == 1, 0, - dStartAccDist - ( i - 2) * dC)
               if ( bFront) then
                  dSal, dEal = dEal, dSal
               end
               EgtSetMachiningParam( MCH_MP.STARTADDLEN, dSal)
               EgtSetMachiningParam( MCH_MP.ENDADDLEN, dEal)
               -- imposto offset radiale
               EgtSetMachiningParam( MCH_MP.OFFSR, dOffset)
               -- imposto attacco/uscita
               EgtSetMachiningParam( MCH_MP.LITANG, dLioTang)
               EgtSetMachiningParam( MCH_MP.LIPERP, dLioPerp)
               EgtSetMachiningParam( MCH_MP.LOTANG, dLioTang)
               EgtSetMachiningParam( MCH_MP.LOPERP, dLioPerp)
               -- imposto posizione braccio porta testa per non ingombrare agli estremi
               EgtSetMachiningParam( MCH_MP.SCC, nSCC)
               -- imposto uso della faccia
               EgtSetMachiningParam( MCH_MP.FACEUSE, nFaceUse2)
               -- eseguo
               if not EgtApplyMachining( true, false) then
                  local _, sErr = EgtGetLastMachMgrError()
                  EgtSetOperationMode( nMchFId, false)
                  return false, sErr
               end
            end
         end
      end

   -- Se non è sotto : lavorazione Long2Cut
   elseif nSide ~= - 1 then

      -- recupero la lavorazione
      local sMilling = ML.FindMilling( 'Long2Cut')
      if not sMilling then
         local sErr = 'Error on process ' .. tostring( Proc.Id) .. ' milling not found in library'
         EgtOutLog( sErr)
         return false, sErr
      end
      -- recupero i dati dell'utensile
      local dToolDiam = 0
      local dMaxDepth = 0
      if EgtMdbSetCurrMachining( sMilling) then
         local sTuuid = EgtMdbGetCurrMachiningParam( MCH_MP.TUUID)
         if EgtTdbSetCurrTool( EgtTdbGetToolFromUUID( sTuuid) or '') then
            dToolDiam = EgtTdbGetCurrToolParam( MCH_TP.DIAM) or dToolDiam
            dMaxDepth = EgtTdbGetCurrToolMaxDepth() or dMaxDepth
         end
      end
      -- se chiuso e corto, applico svuotatura con fresa opportuna
      if nFaceLimit == 3 and Proc.Box:getDimX() < 2 * dToolDiam then
         -- ottengo le facce con il maggior  numero di adiacenze e l'elevazione relativa
         local nFacInd, dFacElev, nFacInd2, dFacElev2 = BL.GetFaceWithMostAdj( Proc.Id, nPartId)
         if nFacInd > 0 and dFacElev then
            local bOk, sErr = MakeByPocketing( Proc, nPhase, nRawId, nPartId, nFacInd, dFacElev)
            if not bOk then return bOk, sErr end
         end
         if nFacInd2 and nFacInd2 > 0 and dFacElev2 then
            return MakeByPocketing( Proc, nPhase, nRawId, nPartId, nFacInd2, dFacElev2)
         end
      end
      -- determino gli estremi
      local dStartDist = 0
      local dStartAccDist = BD.LONGCUT_ENDLEN
      local bStartFixed = true
      -- se ho facce limite a destra
      if nFaceLimit % 2 == 1 then
         dStartDist = dToolDiam / 2
         dStartAccDist = BD.LONGCUT_MAXLEN
         bStartFixed = false
      end
      local dEndDist = 0
      local dEndAccDist = BD.LONGCUT_ENDLEN
      local bEndFixed = true
      -- se ho facce limite a sinistra
      if nFaceLimit >= 2  then
         dEndDist = dToolDiam / 2
         dEndAccDist = BD.LONGCUT_MAXLEN
         bEndFixed = false
      end 
      -- determino il numero di parti in cui dividere la lavorazione
      local nC = ceil( ( dLen - dStartAccDist - dEndAccDist) / BD.LONGCUT_MAXLEN)
      local dC = 0
      if nC > 0 then
         if bStartFixed and bEndFixed then
            dC = ( dLen - dStartAccDist - dEndAccDist) / nC
         elseif bStartFixed then
            dC = ( dLen - dStartAccDist) / ( nC + 1)
            dEndAccDist = dC
         elseif bEndFixed then
            dC = ( dLen - dEndAccDist) / ( nC + 1)
            dStartAccDist = dC
         else
            dC = dLen / ( nC + 2)
            dStartAccDist = dC
            dEndAccDist = dC
         end
         nC = nC + 2
      else
         if dLen > min( dStartAccDist, dEndAccDist) then
            nC = 2
            if bStartFixed and not bEndFixed then
               dStartAccDist = min( dStartAccDist, dLen/2)
               dEndAccDist = dLen - dStartAccDist
            elseif not bStartFixed and bEndFixed then
               dEndAccDist = min( dEndAccDist, dLen/2)
               dStartAccDist = dLen - dEndAccDist
            else
               dStartAccDist = dLen/2
               dEndAccDist = dStartAccDist
            end
         else
            nC = 1
            dStartAccDist = 0
            dEndAccDist = 0
         end
      end
      -- ciclo sulle parti
      local nM = 0
      for j = 1, nC do
         -- se facce ortogonali (concave), lavoro solo quella con versore maggiormente verso l'alto
         local nIni, nFin = 1, 2
         if bOrtho then
            if vtN[vOrd[1]]:getZ() >= vtN[vOrd[2]]:getZ() then
               nFin = 1
            else
               nIni = 2
            end
         end
         -- su entrambe le facce
         for i = nIni, nFin do
            -- Limitazioni della lavorazione
            local nPos = EgtIf( i == 1, j, nC - j + 1)
            local dSal = EgtIf( nPos == 1, - EgtIf( i == nIni, dStartDist, dEndDist), - EgtIf( i == nIni, dStartAccDist, dEndAccDist) - ( nPos - 2) * dC)
            local dEal = EgtIf( nPos == nC, - EgtIf( i == nIni, dEndDist, dStartDist), - EgtIf( i == nIni, dEndAccDist, dStartAccDist) - ( nC - nPos - 1) * dC)
            -- Posizione braccio portatesta
            local nSCC = EgtIf( ( j == 1 or j == nC - 1), MCH_SCC.ADIR_XP, MCH_SCC.ADIR_XM)
            -- Verifico se da fare di fianco perchè normale troppo verso il basso (minore di -30deg)
            local bSide = ( vtN[vOrd[i]]:getZ() < -0.5)
            -- ciclo sulle passate
            local nO = 1
            local dStep = 0
            local dAgg = EgtIf( bConvex, 2 * BD.CUT_EXTRA, BD.CUT_EXTRA)
            local dLargh = vWidth[vOrd[i]]
            if not bSide and dLargh > 0.8 * dToolDiam then 
               nO = ceil( dLargh / ( 0.6 * dToolDiam))
               if nO > 1 then
                  dStep = dLargh / nO
               end
            end
            for k = 1, nO do
               -- inserisco le parti di lavorazione
               nM = nM + 1
               local sNameF = 'L2C_' .. ( EgtGetName( Proc.Id) or tostring( Proc.Id)) .. '_' .. tostring( nM)
               local nMchFId = EgtAddMachining( sNameF, sMilling)
               if not nMchFId then
                  local sErr = 'Error adding machining ' .. sNameF .. '-' .. sMilling
                  EgtOutLog( sErr)
                  return false, sErr
               end
               -- aggiungo geometria
               EgtSetMachiningGeometry( {{ Proc.Id, tFaceLong[vOrd[i]]}})
               -- limito opportunamente la lavorazione
               EgtSetMachiningParam( MCH_MP.STARTADDLEN, dSal)
               EgtSetMachiningParam( MCH_MP.ENDADDLEN, dEal)
               -- imposto posizione braccio porta testa per non ingombrare agli estremi
               EgtSetMachiningParam( MCH_MP.SCC, nSCC)
               -- imposto uso faccia
               EgtSetMachiningParam( MCH_MP.FACEUSE, EgtIf( bSide, MCH_MILL_FU.PARAL_DOWN, vFaceUse[vOrd[i]]))
               -- imposto lato di lavoro e inversione
               EgtSetMachiningParam( MCH_MP.WORKSIDE, MCH_MILL_WS.LEFT)
               EgtSetMachiningParam( MCH_MP.INVERT, EgtIf( bSide, true, false))
               -- imposto offset radiale (nullo se concavo o di lato)
               local dOffsR = EgtIf( k < nO, ( nO - k) * dStep, EgtIf( bConvex, - BD.CUT_EXTRA, 0))
               if bSide then dOffsR = 0 end
               EgtSetMachiningParam( MCH_MP.OFFSR, dOffsR)
               -- se di lato, imposto offset longitudinale
               if bSide then
                  EgtSetMachiningParam( MCH_MP.OFFSL, - BD.CUT_EXTRA)
               end
               -- attacco e uscita
               local dLioPerp = EgtIf( bSide, dToolDiam, dLargh - dOffsR + BD.CUT_SIC)
               EgtSetMachiningParam( MCH_MP.LIPERP, dLioPerp)
               EgtSetMachiningParam( MCH_MP.LOPERP, dLioPerp)
               -- eseguo
               if not EgtApplyMachining( true, false) then
                  local _, sErr = EgtGetLastMachMgrError()
                  EgtSetOperationMode( nMchFId, false)
                  return false, sErr
               end
            end
         end
         -- se facce principali convesse, eventuale lavorazione della faccia limitante l'inizio (a destra)
         if bConvex and j == 1 and not bStartFixed then
            local vtIni = -X_AX() 
            for i = 1, #tFaceLimit do
               local _, vtN = EgtSurfTmFacetCenter( Proc.Id, tFaceLimit[i], GDB_ID.ROOT)
               if vtIni * vtN > 0 then
                  MakeSideFace( Proc.Id, tFaceLimit[i], nSide, sMilling, dToolDiam)
               end
            end
         end
      end
      -- se facce principali convesse, eventuale lavorazione della faccia limitante la fine (a sinistra)
      if bConvex and not bEndFixed then
         local vtFin = X_AX() 
         for i = 1, #tFaceLimit do
            local _, vtN = EgtSurfTmFacetCenter( Proc.Id, tFaceLimit[i], GDB_ID.ROOT)
            if vtFin * vtN > 0 then
               MakeSideFace( Proc.Id, tFaceLimit[i], nSide, sMilling, dToolDiam)
            end
         end
      end
   -- altrimenti è sotto : lavorazione Long2CutDown
   else
      -- recupero la lavorazione
      local sMilling = ML.FindMilling( 'Long2CutDown')
      if not sMilling then
         local sErr = 'Error on process ' .. tostring( Proc.Id) .. ' milling not found in library'
         EgtOutLog( sErr)
         return false, sErr
      end
      -- recupero i dati dell'utensile
      local dToolDiam = 0
      local dMaxDepth = 0
      local dThDiam = 0
      if EgtMdbSetCurrMachining( sMilling) then
         local sTuuid = EgtMdbGetCurrMachiningParam( MCH_MP.TUUID)
         if EgtTdbSetCurrTool( EgtTdbGetToolFromUUID( sTuuid) or '') then
            dToolDiam = EgtTdbGetCurrToolParam( MCH_TP.DIAM) or dToolDiam
            dMaxDepth = EgtTdbGetCurrToolMaxDepth() or dMaxDepth
            dThDiam = EgtTdbGetCurrToolThDiam() or dThDiam 
         end
      end
      -- aggiuntivo sull'affondamento
      local dAgg = EgtIf( bConvex, BD.CUT_EXTRA, 0)
      -- determino gli estremi
      local dStartDist = 0
      local dStartAccDist = BD.LONGCUT_ENDLEN
      local bStartFixed = true
      -- se ho facce limite a destra
      if nFaceLimit % 2 == 1 then
         dStartDist = dToolDiam / 2
         dStartAccDist = BD.LONGCUT_MAXLEN
         bStartFixed = false
      end
      local dEndDist = 0
      local dEndAccDist = BD.LONGCUT_ENDLEN
      local bEndFixed = true
      -- se ho facce limite a sinistra
      if nFaceLimit >= 2  then
         dEndDist = dToolDiam / 2
         dEndAccDist = BD.LONGCUT_MAXLEN
         bEndFixed = false
      end 
      -- determino il numero di parti in cui dividere la lavorazione
      local nC = ceil( ( dLen - dStartAccDist - dEndAccDist) / BD.LONGCUT_MAXLEN)
      local dC = 0
      if nC > 0 then
         if bStartFixed and bEndFixed then
            dC = ( dLen - dStartAccDist - dEndAccDist) / nC
         elseif bStartFixed then
            dC = ( dLen - dStartAccDist) / ( nC + 1)
            dEndAccDist = dC
         elseif bEndFixed then
            dC = ( dLen - dEndAccDist) / ( nC + 1)
            dStartAccDist = dC
         else
            dC = dLen / ( nC + 2)
            dStartAccDist = dC
            dEndAccDist = dC
         end
         nC = nC + 2
      else
         if dLen > min( dStartAccDist, dEndAccDist) then
            nC = 2
            if bStartFixed and not bEndFixed then
               dStartAccDist = min( dStartAccDist, dLen/2)
               dEndAccDist = dLen - dStartAccDist
            elseif not bStartFixed and bEndFixed then
               dEndAccDist = min( dEndAccDist, dLen/2)
               dStartAccDist = dLen - dEndAccDist
            else
               dStartAccDist = dLen/2
               dEndAccDist = dStartAccDist
            end
         else
            nC = 1
            dStartAccDist = 0
            dEndAccDist = 0
         end
      end
      -- ciclo sulle parti
      local nM = 0
      for j = 1, nC do
         -- su entrambe le facce
         for i = 1, 2 do
            -- Limitazioni della lavorazione
            local nPos = EgtIf( i == 1, j, nC - j + 1)
            local dSal = EgtIf( nPos == 1, - EgtIf( i == 1, dStartDist, dEndDist), - EgtIf( i == 1, dStartAccDist, dEndAccDist) - ( nPos - 2) * dC)
            local dEal = EgtIf( nPos == nC, - EgtIf( i == 1, dEndDist, dStartDist), - EgtIf( i == 1, dEndAccDist, dStartAccDist) - ( nC - nPos - 1) * dC)
            -- Posizione braccio portatesta
            local nSCC = MCH_SCC.ADIR_XP
            -- inserisco le parti di lavorazione
            nM = nM + 1
            local sNameF = 'L2C_' .. ( EgtGetName( Proc.Id) or tostring( Proc.Id)) .. '_' .. tostring( nM)
            local nMchFId = EgtAddMachining( sNameF, sMilling)
            if not nMchFId then
               local sErr = 'Error adding machining ' .. sNameF .. '-' .. sMilling
               EgtOutLog( sErr)
               return false, sErr
            end
            -- aggiungo geometria
            EgtSetMachiningGeometry( {{ Proc.Id, tFaceLong[vOrd[i]]}})
            -- limito opportunamente la lavorazione
            EgtSetMachiningParam( MCH_MP.STARTADDLEN, dSal)
            EgtSetMachiningParam( MCH_MP.ENDADDLEN, dEal)
            -- imposto posizione braccio porta testa per non ingombrare agli estremi
            EgtSetMachiningParam( MCH_MP.SCC, nSCC)
            -- imposto uso faccia
            EgtSetMachiningParam( MCH_MP.FACEUSE, vFaceUse[vOrd[i]])
            -- imposto lato di lavoro e inversione
            EgtSetMachiningParam( MCH_MP.WORKSIDE, MCH_MILL_WS.LEFT)
            EgtSetMachiningParam( MCH_MP.INVERT, true)
             -- verifico massimo affondamento (tengo conto dell'inclinazione utensile e della pinza con estremità conica)
            local dCollSic = max( BD.COLL_SIC, ( dThDiam - dToolDiam) / 2 * abs( vtN[vOrd[i]]:getY() / vtN[vOrd[i]]:getZ()) - 3)
            if vWidth[vOrd[i]] + dAgg > dMaxDepth - dCollSic then
               sWarn = 'Warning in LongDoubleCut : depth (' .. EgtNumToString( vWidth[vOrd[i]] + dAgg, 1) .. ') bigger than max tool depth (' .. EgtNumToString( dMaxDepth - dCollSic, 1) .. ')'
            end
            local dDepth = min( dMaxDepth - dCollSic, vWidth[vOrd[i]] + dAgg)
            EgtSetMachiningParam( MCH_MP.DEPTH, dDepth)
            -- eseguo
            if not EgtApplyMachining( true, false) then
               local _, sErr = EgtGetLastMachMgrError()
               EgtSetOperationMode( nMchFId, false)
               return false, sErr
            end
         end
      end
   end

   return true, sWarn
end   

---------------------------------------------------------------------
return ProcessLong2Cut