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                           nro_cas=None):
  32   """edges de bord de la face de peau"""
  33   logging.info('start')
  34
  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))
  38   edgesBords = list()
  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)
  49     else:
  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))
  54
  55   groupEdgesBordPeau = geompy.CreateGroup(facePeau, geompy.ShapeType["EDGE"])
  56   geompy.UnionList(groupEdgesBordPeau, edgesBords)
  57   bordsVifs = None
  58   if aretesVivesC is not None:
  59     logging.info("identification des bords vifs par GetInPlace")
  60     bordsVifs = geompy.GetInPlace(facePeau, aretesVivesC)
  61     if bordsVifs is None:
  62       logging.debug("pas d'identification des bords vifs par GetInPlace: test par distance")
  63       edvifs = list()
  64       arvives = geompy.ExtractShapes(aretesVivesC, geompy.ShapeType["EDGE"], False)
  65       l_edges = geompy.ExtractShapes(facePeau, geompy.ShapeType["EDGE"], False)
  66       for edge in l_edges:
  67         vxs = geompy.ExtractShapes(edge, geompy.ShapeType["VERTEX"], False)
  68         for arete in arvives:
  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)
  72           if ( dist < 0.001 ):
  73             edvifs.append(edge)
  74             break
  75       if edvifs:
  76         logging.info('==> Nombre d arêtes de bord vif : len(edvifs) = %d', len(edvifs))
  77         bordsVifs = geompy.CreateGroup(facePeau,geompy.ShapeType["EDGE"])
  78         for edge in edvifs:
  79           geompy.AddObject(bordsVifs, geompy.GetSubShapeID(facePeau, edge))
  80
  81   if bordsVifs is not None:
  82     geomPublishInFather(initLog.always, facePeau, bordsVifs, "bordsVifs")
  83     groupEdgesBordPeau = geompy.CutGroups(groupEdgesBordPeau, bordsVifs)
  84     grptmp = None
  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
  90     else:
  91       grpnew = bordsVifs
  92     if grpnew is not None:
  93       edv = geompy.ExtractShapes(grpnew, geompy.ShapeType["EDGE"], False)
  94       aretesVivesCoupees += edv
  95
  96   logging.debug("aretesVivesCoupees %s",aretesVivesCoupees)
  97   geomPublishInFather(initLog.always, facePeau, groupEdgesBordPeau , "EdgesBords")
  98
  99   return groupEdgesBordPeau, bordsVifs