src/Tools/blocFissure/gmu/identifieEdgesPeau_b.py

   1 # -*- coding: utf-8 -*-
   2 # Copyright (C) 2014-2021  EDF R&D
   3 #
   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.
   8 #
   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.
  13 #
  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
  17 #
  18 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  19 #
  20 """edges de bord de la face de peau"""
  21
  22 import logging
  23
  24 from . import initLog
  25
  26 from .geomsmesh import geompy
  27 from .geomsmesh import geomPublishInFather
  28
  29 def identifieEdgesPeau_b(facePeau, edgesListees, \
  30                          fillingFaceExterne, aretesVivesC, aretesVivesCoupees):
  31   """edges de bord de la face de peau"""
  32   logging.info('start')
  33
  34   edgesFilling = geompy.ExtractShapes(fillingFaceExterne, geompy.ShapeType["EDGE"], False)
  35   edgesBords = list()
  36   for i, edge in enumerate(edgesFilling):
  37     edgepeau = geompy.GetInPlace(facePeau, edge)
  38     name = "edgepeau%d"%i
  39     geomPublishInFather(initLog.debug, facePeau,edgepeau, name)
  40     logging.debug("edgepeau %s", geompy.ShapeInfo(edgepeau))
  41     if geompy.ShapeInfo(edgepeau)['EDGE'] > 1:
  42       logging.debug("  EDGES multiples")
  43       l_edges = geompy.ExtractShapes(edgepeau, geompy.ShapeType["EDGE"], False)
  44       edgesBords += l_edges
  45       edgesListees += l_edges
  46     else:
  47       logging.debug("  EDGE")
  48       edgesBords.append(edgepeau)
  49       edgesListees.append(edgepeau)
  50
  51   groupEdgesBordPeau = geompy.CreateGroup(facePeau, geompy.ShapeType["EDGE"])
  52   geompy.UnionList(groupEdgesBordPeau, edgesBords)
  53   bordsVifs = None
  54   if aretesVivesC is not None:
  55     logging.debug("identification des bords vifs par GetInPlace")
  56     bordsVifs = geompy.GetInPlace(facePeau, aretesVivesC)
  57     if bordsVifs is None:
  58       logging.debug("pas d'identification des bords vifs par GetInPlace: test par distance")
  59       edvifs = list()
  60       arvives = geompy.ExtractShapes(aretesVivesC, geompy.ShapeType["EDGE"], False)
  61       l_edges = geompy.ExtractShapes(facePeau, geompy.ShapeType["EDGE"], False)
  62       for edge in l_edges:
  63         vxs = geompy.ExtractShapes(edge, geompy.ShapeType["VERTEX"], False)
  64         for arete in arvives:
  65           dist = geompy.MinDistance(vxs[0], arete)
  66           dist += geompy.MinDistance(vxs[1], arete)
  67           logging.debug("test distance bord face peau - arete vive: %s",dist)
  68           if ( dist < 0.001 ):
  69             edvifs.append(edge)
  70             break
  71       if len(edvifs) >0:
  72         bordsVifs = geompy.CreateGroup(facePeau,geompy.ShapeType["EDGE"])
  73         for edge in edvifs:
  74           geompy.AddObject(bordsVifs, geompy.GetSubShapeID(facePeau, edge))
  75
  76   if bordsVifs is not None:
  77     geomPublishInFather(initLog.debug, facePeau, bordsVifs, "bordsVifs")
  78     groupEdgesBordPeau = geompy.CutGroups(groupEdgesBordPeau, bordsVifs)
  79     grptmp = None
  80     if len(aretesVivesCoupees) > 0:
  81       grpc = geompy.MakeCompound(aretesVivesCoupees)
  82       grptmp = geompy.GetInPlace(facePeau, grpc)
  83     if grptmp is not None:
  84       grpnew = geompy.CutGroups(bordsVifs, grptmp) # ce qui est nouveau dans bordsVifs
  85     else:
  86       grpnew = bordsVifs
  87     if grpnew is not None:
  88       edv = geompy.ExtractShapes(grpnew, geompy.ShapeType["EDGE"], False)
  89       aretesVivesCoupees += edv
  90
  91   logging.debug("aretesVivesCoupees %s",aretesVivesCoupees)
  92   geomPublishInFather(initLog.debug, facePeau, groupEdgesBordPeau , "EdgesBords")
  93
  94   return groupEdgesBordPeau, bordsVifs