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))
37 edgesPeau = geompy.ExtractShapes(facePeau, geompy.ShapeType["EDGE"], False)
38 logging.info('Nombre total d arêtes de la peau : len(edgesPeau) = %d', len(edgesPeau))
40 edges = substractSubShapes(facePeau, edgesPeau, edgesListees)
42 l_edge_cercle = list()
43 # --- au moins une extrémité du pipe sur cette face de peau : arêtes sui les prenent toutes
47 for groupe in verticesCircPeau:
49 for id_vertex in geompy.GetObjectIDs(groupe):
50 vertex = geompy.GetSubShape(facePeau, [id_vertex])
51 distance = geompy.MinDistance(vertex, edge)
52 logging.info(".. distance %s", distance)
53 if ( distance > 1.e-3 ):
57 logging.info("... entrée de %s dans l_edge_cercle", edge)
58 l_edge_cercle.append(edge)
60 name = "edgeFissurePeauCercle_{}".format(i_aux)
61 geomPublishInFather(initLog.info, facePeau, edge, name)
63 # --- Les arêtes n'appartenant pas aux cercles
64 edgesFissurePeau = list()
65 for i_aux, edge in enumerate(edges):
66 if edge not in l_edge_cercle:
67 edgesFissurePeau.append(edge)
68 name = "edgeFissurePeau_{}".format(i_aux)
69 geomPublishInFather(initLog.debug, facePeau, edge, name)
71 logging.info('==> Nombre d arêtes : len(edgesFissurePeau) = %d', len(edgesFissurePeau))
73 return edgesFissurePeau