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 bord de la face de peau"""
26 from .geomsmesh import geompy
27 from .geomsmesh import geomPublishInFather
29 def identifieEdgesPeau_b(facePeau, edgesListees, \
30 fillingFaceExterne, aretesVivesC, aretesVivesCoupees, \
32 """edges de bord de la face de peau"""
35 # Liste des arêtes de bord
36 edgesFilling = geompy.ExtractShapes(fillingFaceExterne, geompy.ShapeType["EDGE"], False)
37 logging.info('Détermination des arêtes de bord à partir des %d arêtes de fillingFaceExterne', len(edgesFilling))
39 for i_aux, edge in enumerate(edgesFilling):
40 edgepeau = geompy.GetInPlace(facePeau, edge)
41 name = "edgepeau_{}".format(i_aux)
42 geomPublishInFather(initLog.debug, facePeau,edgepeau, name)
43 logging.debug("edgepeau %s", geompy.ShapeInfo(edgepeau))
44 if geompy.ShapeInfo(edgepeau)['EDGE'] > 1:
45 logging.debug(" EDGES multiples")
46 l_edges = geompy.ExtractShapes(edgepeau, geompy.ShapeType["EDGE"], False)
47 edgesBords.extend(l_edges)
48 edgesListees.extend(l_edges)
50 logging.debug(" EDGE")
51 edgesBords.append(edgepeau)
52 edgesListees.append(edgepeau)
53 logging.info('==> Nombre d arêtes de bord : len(edgesBords) = %d', len(edgesBords))
55 groupEdgesBordPeau = geompy.CreateGroup(facePeau, geompy.ShapeType["EDGE"])
56 geompy.UnionList(groupEdgesBordPeau, edgesBords)
58 if aretesVivesC is not None:
59 logging.info("identification des bords vifs par GetInPlace")
60 bordsVifs = geompy.GetInPlace(facePeau, aretesVivesC)
62 logging.debug("pas d'identification des bords vifs par GetInPlace: test par distance")
64 arvives = geompy.ExtractShapes(aretesVivesC, geompy.ShapeType["EDGE"], False)
65 l_edges = geompy.ExtractShapes(facePeau, geompy.ShapeType["EDGE"], False)
67 vxs = geompy.ExtractShapes(edge, geompy.ShapeType["VERTEX"], False)
69 dist = geompy.MinDistance(vxs[0], arete)
70 dist += geompy.MinDistance(vxs[1], arete)
71 logging.debug("test distance bord face peau - arete vive: %s",dist)
76 logging.info('==> Nombre d arêtes de bord vif : len(edvifs) = %d', len(edvifs))
77 bordsVifs = geompy.CreateGroup(facePeau,geompy.ShapeType["EDGE"])
79 geompy.AddObject(bordsVifs, geompy.GetSubShapeID(facePeau, edge))
81 if bordsVifs is not None:
82 geomPublishInFather(initLog.always, facePeau, bordsVifs, "bordsVifs")
83 groupEdgesBordPeau = geompy.CutGroups(groupEdgesBordPeau, bordsVifs)
85 if len(aretesVivesCoupees) > 0:
86 grpc = geompy.MakeCompound(aretesVivesCoupees)
87 grptmp = geompy.GetInPlace(facePeau, grpc)
88 if grptmp is not None:
89 grpnew = geompy.CutGroups(bordsVifs, grptmp) # ce qui est nouveau dans bordsVifs
92 if grpnew is not None:
93 edv = geompy.ExtractShapes(grpnew, geompy.ShapeType["EDGE"], False)
94 aretesVivesCoupees += edv
96 logging.debug("aretesVivesCoupees %s",aretesVivesCoupees)
97 geomPublishInFather(initLog.always, facePeau, groupEdgesBordPeau , "EdgesBords")
99 return groupEdgesBordPeau, bordsVifs