-- WFeatureTopology.lua by Egaltech s.r.l. 2023/06/23
-- Libreria per classificazione topologica feature pareti

-- Tabella per definizione modulo
local WFeatureTopology = {}

-- Include
require( 'EgtBase')

-- Carico le librerie
local WL = require( 'WallLib')

EgtOutLog( ' WFeatureTopology started', 1)

---------------------------------------------------------------------
local function GetAdjacencyMatrix(Proc )
   local vAdj = {}
   for i = 1, Proc.Fct do
      vAdj[i] = {}
      for j = 1, Proc.Fct do
         if i == j then
            vAdj[i][j] = 0
         else
            _, _, _, vAdj[i][j] = EgtSurfTmFacetsContact( Proc.Id, i - 1, j - 1, GDB_ID.ROOT)
         end
         j = j + 1
      end
      i = i + 1
   end
   return vAdj
end

---------------------------------------------------------------------
local function AreAllAnglesConcaveOrRight( Proc)
   if Proc.Fct < 2 then return nil end
   local vAdj = GetAdjacencyMatrix( Proc)
   local bAllConcave, bAllRight = true, true
   for i = 1, Proc.Fct do
      for j = 1, Proc.Fct do
         if vAdj[i][j] and vAdj[i][j] > 0 then
            bAllConcave = false
            bAllRight = false
            break
         elseif vAdj[i][j] and vAdj[i][j] ~= 0 and vAdj[i][j] + 90 > GEO.EPS_SMALL then
            bAllRight = false
         end
      end
   end
   return bAllConcave, bAllRight
end

---------------------------------------------------------------------
local function IsAnyDimensionLongAsPart( Proc )
   local bResult = false
   local b3Solid = EgtGetBBoxGlob( Proc.PartId, GDB_BB.STANDARD)
   if Proc.Box:getDimX() > b3Solid:getDimX() - 1000 * GEO.EPS_SMALL or
      Proc.Box:getDimY() > b3Solid:getDimY() - 1000 * GEO.EPS_SMALL or
      Proc.Box:getDimZ() > b3Solid:getDimZ() - 1000 * GEO.EPS_SMALL then
      bResult = true
   end
   return bResult
end

---------------------------------------------------------------------
local function GetTriangularFaces( Proc)
   local vTriangularFaces = {}
   for i = 1, Proc.Fct do
      if WL.Is3EdgesApprox( Proc, i - 1) then
         table.insert( vTriangularFaces, i - 1)
      end
   end
   return vTriangularFaces
end

---------------------------------------------------------------------
local function GetFacesParallelToPart( Proc, sFamily)
   local vFacesParallelToPart = {}
   for i = 0, Proc.Fct - 1 do
      local vtN = EgtSurfTmFacetNormVersor( Proc.Id, i, GDB_ID.ROOT)
      if sFamily == 'Bevel' then
         local vTriangularFaces = GetTriangularFaces( Proc)
         local bIsTriangularFace = false
         -- verifico se la faccia è triangolare
         for j = 1, #vTriangularFaces do
            if i == vTriangularFaces[j] then
               bIsTriangularFace = true
            end
         end
         -- se faccia triangolare deve avere la normale parallela ad una direzione principale
         if bIsTriangularFace then
            if AreSameOrOppositeVectorApprox( vtN, X_AX()) or AreSameOrOppositeVectorApprox( vtN, Y_AX()) or AreSameOrOppositeVectorApprox( vtN, Z_AX()) then
               table.insert( vFacesParallelToPart, i)
            end
         -- altrimenti deve avere una componente della normale nulla
         else
            if vtN:getX() < 10 * GEO.EPS_SMALL or vtN:getY() < 10 * GEO.EPS_SMALL or vtN:getZ() < 10 * GEO.EPS_SMALL then
               table.insert( vFacesParallelToPart, i)
            end
         end
      else
         -- la normale deve essere parallela ad una direzione principale
         if AreSameOrOppositeVectorApprox( vtN, X_AX()) or AreSameOrOppositeVectorApprox( vtN, Y_AX()) or AreSameOrOppositeVectorApprox( vtN, Z_AX()) then
            table.insert( vFacesParallelToPart, i)
         end
      end
   end
   return vFacesParallelToPart
end

---------------------------------------------------------------------
local function GetTopologyLongName( sFamily, bIsThrough, bAllRightAngles, bIsParallel, nNumberOfFaces)
   -- feature passante o cieca
   local sThrough = '_'
   if bIsThrough ~= nil then sThrough = EgtIf( bIsThrough, 'Through', 'Blind') end
   -- tutti gli angoli della feature sono retti oppure no
   local sAllRightAngles = '_'
   if bAllRightAngles ~= nil then sAllRightAngles = EgtIf( bAllRightAngles, 'RightAngles', 'NotRightAngles') end
   -- tutte le dimensioni della feature sono parallele agli assi principali del pezzo oppure no
   local sParallel = '_'
   if bIsParallel ~= nil then sParallel = EgtIf( bIsParallel, 'Parallel', 'NotParallel') end

   local sLongName = sFamily .. '-' .. sThrough .. '-' .. sAllRightAngles .. '-' .. sParallel .. '-' .. nNumberOfFaces
   return sLongName
end

---------------------------------------------------------------------
function WFeatureTopology.Classify( Proc)
   local bRecognized = false
   local sFamily
   local bIsThrough
   local bAllRightAngles
   local bIsParallel
   local sLongName = ''

   -- SE NON HA TUTTE LE FACCE PIANE RITORNARE NIL!!
   local bAllAnglesConcave
   bAllAnglesConcave, bAllRightAngles = AreAllAnglesConcaveOrRight( Proc)
   local vTriangularFaces = GetTriangularFaces( Proc)
   local bIsAnyFaceTriangular = #vTriangularFaces > 0
   local bIsAnyDimensionLongAsPart = IsAnyDimensionLongAsPart( Proc)
   if Proc.Fct == 1 and bIsAnyDimensionLongAsPart then
      sFamily = 'Bevel'
      bIsThrough = true
   elseif Proc.Fct == 2 and bAllAnglesConcave and bIsAnyFaceTriangular then
      sFamily = 'Bevel'
      bIsThrough = false
   elseif Proc.Fct == 2 and bAllAnglesConcave and bIsAnyDimensionLongAsPart then
      sFamily = 'Rabbet'
      bIsThrough = true
   end
   local vFacesParallelToPart = GetFacesParallelToPart( Proc, sFamily)
   bIsParallel = ( #vFacesParallelToPart == Proc.Fct)

   if 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
   end

   return bRecognized
end

-------------------------------------------------------------------------------------------------------------
return WFeatureTopology