1 # -*- coding: utf-8 -*-
2 # Copyright (C) 2014-2023 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 """Recherche edges communes entre une face inside et (faces onside, edges pipe et fond débouchante)"""
26 from .geomsmesh import geompy
27 from .geomsmesh import geomPublishInFather
29 def identifieElementsDebouchants_d(iedf,\
30 partitionPeauFissFond, edgesFondIn, edgesPipeIn, \
31 facesInside, facesOnside, \
32 ptPeau, facesFissExt, edgesFissExtPeau, edgesFissExtPipe):
33 """Recherche edges communes entre une face inside et (faces onside, edges pipe et fond débouchante)"""
35 for ifa, face in enumerate(facesInside):
36 logging.debug("recherche edges communes entre une face inside et (faces onside, edges pipe et fond débouchante)")
41 edgesPeauFis = geompy.GetSharedShapesMulti([geompy.MakeCompound(facesOnside), face], geompy.ShapeType["EDGE"])
42 logging.debug(" faces onside %s",edgesPeauFis)
43 edgesPipeFis = geompy.GetSharedShapesMulti([geompy.MakeCompound(edgesPipeIn), face], geompy.ShapeType["EDGE"])
44 logging.debug(" edgesPipeIn %s", edgesPipeFis)
45 edgesPipeFnd = geompy.GetSharedShapesMulti([geompy.MakeCompound(edgesFondIn), face], geompy.ShapeType["EDGE"])
46 logging.debug(" edgesFondIn %s ", edgesPipeFnd)
48 logging.debug(" pb edges communes %s %s %s",edgesPeauFis, edgesPipeFis, edgesPipeFnd)
49 if (len(edgesPeauFis) > 0) and (len(edgesPipeFis) > 0) and (len(edgesPipeFnd) == 0):
50 dist = geompy.MinDistance(geompy.MakeCompound(edgesPeauFis), ptPeau)
51 logging.debug(" test distance extrémité reference %s", dist)
52 if dist < 1.e-3: # c'est la face de fissure externe associée
53 logging.debug(" face %s inside ajoutée", ifa)
54 facesFissExt.append(face)
55 name="faceFissExt%d"%iedf
56 geomPublishInFather(initLog.debug, partitionPeauFissFond, face, name)
58 for edpe in edgesPeauFis:
59 for edpi in edgesPipeFis:
60 dist = geompy.MinDistance(edpe, edpi)
62 edgesFissExtPeau.append(edpe)
63 name="edgesFissExtPeau%d"%iedf
64 geomPublishInFather(initLog.debug, partitionPeauFissFond, edpe, name)
65 edgesFissExtPipe.append(edpi)
66 name="edgesFissExtPipe%d"%iedf
67 geomPublishInFather(initLog.debug, partitionPeauFissFond, edpi, name)