src/Tools/blocFissure/gmu/identifieEdgesPeau_c.py

   1 # -*- coding: utf-8 -*-
   2 # Copyright (C) 2014-2021  EDF R&D
   3 #
   4 # This library is free software; you can redistribute it and/or
   5 # modify it under the terms of the GNU Lesser General Public
   6 # License as published by the Free Software Foundation; either
   7 # version 2.1 of the License, or (at your option) any later version.
   8 #
   9 # This library is distributed in the hope that it will be useful,
  10 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  11 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  12 # Lesser General Public License for more details.
  13 #
  14 # You should have received a copy of the GNU Lesser General Public
  15 # License along with this library; if not, write to the Free Software
  16 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  17 #
  18 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  19 #
  20 """edges de la face de peau partagées avec la face de fissure"""
  21
  22 import logging
  23
  24 from . import initLog
  25
  26 from .geomsmesh import geompy
  27 from .geomsmesh import geomPublishInFather
  28
  29 from .substractSubShapes import substractSubShapes
  30
  31 def identifieEdgesPeau_c(verticesPipePeau, facePeau, edgesListees, verticesCircPeau):
  32   """edges de la face de peau partagées avec la face de fissure"""
  33   logging.info('start')
  34   logging.info("Traitement des arêtes de '%s'", facePeau.GetName())
  35   logging.info('Nombre de sommets : len(verticesPipePeau) = %d', len(verticesPipePeau))
  36   #logging.info('verticesPipePeau = %s', verticesPipePeau)
  37   logging.info('Nombre de sommets : len(verticesCircPeau) = %d', len(verticesCircPeau))
  38   #logging.info('verticesCircPeau = %s', verticesCircPeau)
  39
  40   edgesPeau = geompy.ExtractShapes(facePeau, geompy.ShapeType["EDGE"], False)
  41   logging.info('Nombre total d arêtes de la peau : len(edgesPeau) = %d', len(edgesPeau))
  42
  43   edges = substractSubShapes(facePeau, edgesPeau, edgesListees)
  44
  45   edgesFissurePeau = list()
  46 # --- au moins une extrémité du pipe sur cette face de peau
  47   if verticesPipePeau:
  48     # En premier, les edges associés aux extrémités du pipe
  49     edgesFissurePeau = [None for _ in range(len(verticesCircPeau))]
  50     for edge in edges:
  51       for i_aux, vertex in enumerate(verticesCircPeau):
  52         logging.info(".. distance %s", geompy.MinDistance(vertex, edge))
  53         if ( ( geompy.MinDistance(vertex, edge) < 1.e-3 ) and ( edge not in edgesFissurePeau ) ):
  54           edgesFissurePeau[i_aux] = edge
  55           name = "edgeFissurePeau_{}".format(i_aux)
  56           logging.info("... entrée de %s à la place %d", edge, i_aux)
  57           geomPublishInFather(initLog.debug, facePeau, edge, name)
  58     # Ensuite, on ajoute les edges manquantes
  59     for edge in edges:
  60       if edge not in edgesFissurePeau:
  61         logging.info("... ajout")
  62         edgesFissurePeau.append(edge)
  63
  64   else:
  65     for i_aux, edge in enumerate(edges):
  66       edgesFissurePeau.append(edge)
  67       name = "edgeFissurePeau{}".format(i_aux)
  68       geomPublishInFather(initLog.debug, facePeau, edge, name)
  69
  70   logging.info('==> Nombre d arêtes : len(edgesFissurePeau) = %d', len(edgesFissurePeau))
  71
  72   return edgesFissurePeau