1 # -*- coding: utf-8 -*-
2 # Copyright (C) 2014-2021 EDF R&D
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.
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.
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
18 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
20 """edges de la face de peau partagées avec la face de fissure"""
26 from .geomsmesh import geompy
27 from .geomsmesh import geomPublishInFather
29 from .substractSubShapes import substractSubShapes
31 def identifieEdgesPeau_c(verticesPipePeau, facePeau, edgesListees, verticesCircPeau):
32 """edges de la face de peau partagées avec la face de fissure"""
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)
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))
43 edges = substractSubShapes(facePeau, edgesPeau, edgesListees)
45 edgesFissurePeau = list()
46 # --- au moins une extrémité du pipe sur cette face de peau
48 # En premier, les edges associés aux extrémités du pipe
49 edgesFissurePeau = [None for _ in range(len(verticesCircPeau))]
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
60 if edge not in edgesFissurePeau:
61 logging.info("... ajout")
62 edgesFissurePeau.append(edge)
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)
70 logging.info('==> Nombre d arêtes : len(edgesFissurePeau) = %d', len(edgesFissurePeau))
72 return edgesFissurePeau