[bos #40653][CEA] New mesh import export formats with meshio.

[modules/smesh.git] / src / Tools / blocFissure / gmu / peauInterne.py

diff --git a/src/Tools/blocFissure/gmu/peauInterne.py b/src/Tools/blocFissure/gmu/peauInterne.py
index 6c241a75ce883946ce0fd04bd8e04e82c1281828..defb410012e48989a84638b6b37ef971f0429339 100644 (file)
--- a/src/Tools/blocFissure/gmu/peauInterne.py
+++ b/src/Tools/blocFissure/gmu/peauInterne.py
@@ -1,5 +1,5 @@
 # -*- coding: utf-8 -*-
-# Copyright (C) 2014-2021  EDF R&D
+# Copyright (C) 2014-2024  EDF
 #
 # This library is free software; you can redistribute it and/or
 # modify it under the terms of the GNU Lesser General Public
@@ -29,8 +29,10 @@ from .geomsmesh import smesh
 from .fissError import fissError
 from .listOfExtraFunctions import lookForCorner
 from .fusionMaillageAttributionDefaut import fusionMaillageDefaut
+from .putName import putName
 
-def peauInterne(fichierMaillage, shapeDefaut, nomZones):
+def peauInterne(fichierMaillage, shapeDefaut, nomZones, \
+                nro_cas=None):
   """Retrouve les groupes de défaut dans le maillage sain modifié par CreateHoleSkin (CreeZoneDefautMaillage)
 
   On récupère le volume et la peau de la zone de défaut, les éventuelles faces et arêtes internes de cette zone.
@@ -44,31 +46,20 @@ def peauInterne(fichierMaillage, shapeDefaut, nomZones):
   #for grp in groups:
     #print ("\t{}".format(grp.GetName()))
   zoneDefaut = None
-  for grp in groups:
-    logging.debug("groupe %s",grp.GetName())
-    if grp.GetName() == nomZones + "_vol":
-      zoneDefaut = grp
-      break
-  #print ("zoneDefaut = {}".format(zoneDefaut))
   zoneDefaut_skin = None
-  for grp in groups:
-    if grp.GetName() == nomZones + "_skin":
-      zoneDefaut_skin = grp
-      break
-  #print ("zoneDefaut_skin = {}".format(zoneDefaut_skin))
   zoneDefaut_internalFaces = None
-  for grp in groups:
-    if grp.GetName() == nomZones + "_internalFaces":
-      zoneDefaut_internalFaces = grp
-      break
-  #print ("zoneDefaut_internalFaces = {}".format(zoneDefaut_internalFaces))
   zoneDefaut_internalEdges = None
   for grp in groups:
-    if grp.GetName() == nomZones + "_internalEdges":
+    nom = grp.GetName()
+    logging.debug("groupe %s",nom)
+    if ( nom == nomZones + "_vol" ):
+      zoneDefaut = grp
+    elif ( nom == nomZones + "_skin" ):
+      zoneDefaut_skin = grp
+    elif ( nom == nomZones + "_internalFaces" ):
+      zoneDefaut_internalFaces = grp
+    elif ( nom == nomZones + "_internalEdges" ):
       zoneDefaut_internalEdges = grp
-      break
-  #print ("zoneDefaut_internalEdges = {}".format(zoneDefaut_internalEdges))
-
   # --- Le groupe ZoneDefaut ne doit contenir que des Hexaèdres"
 
   info = maillageSain.GetMeshInfo(zoneDefaut)
@@ -95,13 +86,17 @@ def peauInterne(fichierMaillage, shapeDefaut, nomZones):
   _, maillageSain, DefautBoundary = maillageSain.MakeBoundaryElements( SMESH.BND_2DFROM3D, 'DefBound', '', 0, [ zoneDefaut ])
   internal = maillageSain.GetMesh().CutListOfGroups( [ DefautBoundary ], [ zoneDefaut_skin ], 'internal' )
   internalBoundary = smesh.CopyMesh( internal, 'internalBoundary', 0, 0)
+  putName(internalBoundary, "internalBoundary", i_pref=nro_cas)
 
   maillageDefautCible = smesh.CopyMesh(zoneDefaut_skin, 'maillageCible', 0, 0)
+  putName(maillageDefautCible, "maillageCible", i_pref=nro_cas)
   listOfCorner = lookForCorner(maillageDefautCible)
   texte = "listOfCorner = {}".format(listOfCorner)
   logging.debug(texte)
   if listOfCorner:
     logging.info("présence de coins à la surface externe de la zone à reconstruire")
-    zoneDefaut_skin, internalBoundary = fusionMaillageDefaut(maillageSain, maillageDefautCible, internalBoundary, zoneDefaut_skin, shapeDefaut, listOfCorner)
+    zoneDefaut_skin, internalBoundary = \
+          fusionMaillageDefaut(maillageSain, maillageDefautCible, internalBoundary, zoneDefaut_skin, shapeDefaut, listOfCorner, \
+          nro_cas)
 
   return maillageSain, internalBoundary, zoneDefaut, zoneDefaut_skin, zoneDefaut_internalFaces, zoneDefaut_internalEdges