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 edgesFilling = geompy.ExtractShapes(fillingFaceExterne, geompy.ShapeType["EDGE"], False)
37 for i_aux, edge in enumerate(edgesFilling):
38 edgepeau = geompy.GetInPlace(facePeau, edge)
39 name = "edgepeau{}".format(i_aux)
40 geomPublishInFather(initLog.debug, facePeau,edgepeau, name)
41 logging.debug("edgepeau %s", geompy.ShapeInfo(edgepeau))
42 if geompy.ShapeInfo(edgepeau)['EDGE'] > 1:
43 logging.debug(" EDGES multiples")
44 l_edges = geompy.ExtractShapes(edgepeau, geompy.ShapeType["EDGE"], False)
46 edgesListees += l_edges
48 logging.debug(" EDGE")
49 edgesBords.append(edgepeau)
50 edgesListees.append(edgepeau)
52 groupEdgesBordPeau = geompy.CreateGroup(facePeau, geompy.ShapeType["EDGE"])
53 geompy.UnionList(groupEdgesBordPeau, edgesBords)
55 if aretesVivesC is not None:
56 logging.debug("identification des bords vifs par GetInPlace")
57 bordsVifs = geompy.GetInPlace(facePeau, aretesVivesC)
59 logging.debug("pas d'identification des bords vifs par GetInPlace: test par distance")
61 arvives = geompy.ExtractShapes(aretesVivesC, geompy.ShapeType["EDGE"], False)
62 l_edges = geompy.ExtractShapes(facePeau, geompy.ShapeType["EDGE"], False)
64 vxs = geompy.ExtractShapes(edge, geompy.ShapeType["VERTEX"], False)
66 dist = geompy.MinDistance(vxs[0], arete)
67 dist += geompy.MinDistance(vxs[1], arete)
68 logging.debug("test distance bord face peau - arete vive: %s",dist)
73 bordsVifs = geompy.CreateGroup(facePeau,geompy.ShapeType["EDGE"])
75 geompy.AddObject(bordsVifs, geompy.GetSubShapeID(facePeau, edge))
77 if bordsVifs is not None:
78 geomPublishInFather(initLog.always, facePeau, bordsVifs, "bordsVifs")
79 groupEdgesBordPeau = geompy.CutGroups(groupEdgesBordPeau, bordsVifs)
81 if len(aretesVivesCoupees) > 0:
82 grpc = geompy.MakeCompound(aretesVivesCoupees)
83 grptmp = geompy.GetInPlace(facePeau, grpc)
84 if grptmp is not None:
85 grpnew = geompy.CutGroups(bordsVifs, grptmp) # ce qui est nouveau dans bordsVifs
88 if grpnew is not None:
89 edv = geompy.ExtractShapes(grpnew, geompy.ShapeType["EDGE"], False)
90 aretesVivesCoupees += edv
92 logging.debug("aretesVivesCoupees %s",aretesVivesCoupees)
93 geomPublishInFather(initLog.always, facePeau, groupEdgesBordPeau , "EdgesBords")
95 return groupEdgesBordPeau, bordsVifs