Salome HOME
debug
authorPaul RASCLE <paul.rascle@edf.fr>
Sun, 21 Sep 2014 19:49:53 +0000 (21:49 +0200)
committerPaul RASCLE <paul.rascle@edf.fr>
Sun, 21 Sep 2014 19:49:53 +0000 (21:49 +0200)
src/Tools/blocFissure/CasTests/execution_Cas.py
src/Tools/blocFissure/gmu/construitEdgesRadialesDebouchantes.py
src/Tools/blocFissure/gmu/construitFissureGenerale.py
src/Tools/blocFissure/gmu/construitMaillagePipe.py
src/Tools/blocFissure/gmu/identifieElementsGeometriquesPeau.py
src/Tools/blocFissure/gmu/identifieFacesEdgesFissureExterne.py
src/Tools/blocFissure/gmu/mailleAretesEtJonction.py
src/Tools/blocFissure/gmu/mailleFacesFissure.py
src/Tools/blocFissure/gmu/mailleFacesPeau.py

index 9936248f31cbd7c94552e2693b82fe5b69149341..580861961fbe9542882c30c2cca449ac96c22b42 100644 (file)
@@ -131,7 +131,7 @@ runall = True
 if runall:
   torun = [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,]
 else: #prob 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27   
-  torun = [ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,]
+  torun = [ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0,]
   
 for i in range(len(problemes)):
   if torun[i]:
index 29430917716386f1cbd6b808f65cee6d6df2dbd6..91b1b725c87af6781ad83e17c5f25eedc4567723 100644 (file)
@@ -3,6 +3,7 @@
 import logging
 
 from geomsmesh import geompy
+import GEOM
 from sortEdges import sortEdges
 
 def construitEdgesRadialesDebouchantes(idisklim, idiskout, gptsdisks, raydisks,
@@ -48,6 +49,9 @@ def construitEdgesRadialesDebouchantes(idisklim, idiskout, gptsdisks, raydisks,
   for i, nappes in enumerate(listNappes):
     if facesDebouchantes[i]:
       for k, face in enumerate(facesPipePeau):
+        #logging.debug('i, k, face, nappes[0] %s %s %s %s', i, k, face, nappes[0])
+        #geompy.addToStudy(nappes[0], 'lanappe')
+        #geompy.addToStudy(face, 'laface')
         edge = geompy.MakeSection(face, nappes[0])
         if geompy.NbShapes(edge, geompy.ShapeType["EDGE"]) > 0:
           idFacesDebouchantes[i] = k
index f6858132f469a8aaec7bb8b1a8e771f69bf1b61e..d4936d0ae5f303e9353a062ca547cc39bdf31616 100644 (file)
@@ -154,8 +154,8 @@ def construitFissureGenerale(maillagesSains,
   nbFacesFilling = len(partitionsPeauFissFond)
   
   ptEdgeFond = [ []  for i in range(nbFacesFilling)] # pour chaque face [points edge fond de fissure aux débouchés du pipe]
-  facesPipePeau = [ []  for i in range(nbFacesFilling)] # pour chaque face [faces du pipe débouchantes]
-  edgeRadFacePipePeau = [ []  for i in range(nbFacesFilling)] # pour chaque face [edge radiale des faces du pipe débouchantes ]
+  fsPipePeau = [ []  for i in range(nbFacesFilling)] # pour chaque face [faces du pipe débouchantes]
+  edRadFPiPo = [ []  for i in range(nbFacesFilling)] # pour chaque face [edge radiale des faces du pipe débouchantes ]
   fsFissuExt = [ []  for i in range(nbFacesFilling)] # pour chaque face [faces de fissure externes au pipe]
   edFisExtPe = [ []  for i in range(nbFacesFilling)] # pour chaque face [edge en peau des faces de fissure externes (pas subshape facePeau)]
   edFisExtPi = [ []  for i in range(nbFacesFilling)] # pour chaque face [edge commun au pipe des faces de fissure externes]
@@ -174,8 +174,8 @@ def construitFissureGenerale(maillagesSains,
                                                                       facesDefaut, centreFondFiss, rayonPipe,
                                                                       aretesVivesCoupees)      
       ptEdgeFond[ifil] = dataPPFF['endsEdgeFond']
-      facesPipePeau[ifil] = dataPPFF['facesPipePeau']
-      edgeRadFacePipePeau[ifil] = dataPPFF['edgeRadFacePipePeau']
+      fsPipePeau[ifil] = dataPPFF['facesPipePeau']
+      edRadFPiPo[ifil] = dataPPFF['edgeRadFacePipePeau']
       fsFissuExt[ifil] = dataPPFF['facesFissExt']
       edFisExtPe[ifil] = dataPPFF['edgesFissExtPeau']
       edFisExtPi[ifil] = dataPPFF['edgesFissExtPipe']
@@ -187,6 +187,11 @@ def construitFissureGenerale(maillagesSains,
       edFissPeau[ifil] = dataPPFF['edgesFissurePeau']
       ptFisExtPi[ifil] = dataPPFF['verticesPipePeau']
 
+  facesPipePeau = []
+  edgeRadFacePipePeau = []
+  for ifil in range(nbFacesFilling):
+    facesPipePeau += fsPipePeau[ifil]
+    edgeRadFacePipePeau += edRadFPiPo[ifil]
   
   for i, avc in enumerate(aretesVivesCoupees):
     name = "areteViveCoupee%d"%i
index d5f395db33790a71ac68f72f37d71f793499dd34..e08ba47c343c2046005a8347966984560b5e2fd8 100644 (file)
@@ -4,6 +4,7 @@ import logging
 
 from geomsmesh import geompy
 from geomsmesh import smesh
+import SMESH
 
 def construitMaillagePipe(gptsdisks, idisklim, nbsegCercle, nbsegRad):
   """
index a2a19001c3935c3f1f19d882c544aca56878a9e5..c60682c513b17832f08b1f6fc5340602f08937bf 100644 (file)
@@ -54,7 +54,7 @@ def identifieElementsGeometriquesPeau(ifil, partitionPeauFissFond, edgesPipeFiss
                                                                          edgesPeauFondIn, fillingFaceExterne, aretesVivesC, aretesVivesCoupees)
   
   dataPPFF = dict(endsEdgeFond        = endsEdgeFond,        # pour chaque face [points edge fond de fissure aux débouchés du pipe]
-                  facesPipePeau       = facesPipePeau,       # pour chaque face [face
+                  facesPipePeau       = facesPipePeau,       # pour chaque face [faces du pipe débouchantes]
                   edgeRadFacePipePeau = edgeRadFacePipePeau, # pour chaque face [edge radiale des faces du pipe débouchantes ]
                   facesFissExt        = facesFissExt,        # pour chaque face [faces de fissure externes au pipe]
                   edgesFissExtPeau    = edgesFissExtPeau,    # pour chaque face [edge en peau des faces de fissure externes (pas subshape facePeau)]
index 0b4c780b183a8e2f9426c80d6d549cd3d8f4507f..856cfae483e8a0354ea6ae17a85afc4caa77538a 100644 (file)
@@ -13,7 +13,7 @@ def identifieFacesEdgesFissureExterne(fsFissuExt, edFisExtPe, edFisExtPi, edgesP
   facesFissExt = []
   edgesFissExtPeau = []
   edgesFissExtPipe = []
-  for ifil in len(fsFissuExt): # TODO: éliminer les doublons (comparer tous les vertices triés, avec mesure de distance ?)
+  for ifil in range(len(fsFissuExt)): # TODO: éliminer les doublons (comparer tous les vertices triés, avec mesure de distance ?)
     facesFissExt += fsFissuExt[ifil]
     edgesFissExtPeau += edFisExtPe[ifil]
     edgesFissExtPipe += edFisExtPi[ifil]
index bd070af8e5902a723daa0f5dfd6bc6db3462442f..179124308406918a1b1cf30e57c645d7a8ffc223 100644 (file)
@@ -4,6 +4,7 @@ import logging
 
 from geomsmesh import geompy
 from geomsmesh import smesh
+import SMESH
 
 from putName import putName
   
@@ -30,6 +31,7 @@ def mailleAretesEtJonction(internalBoundary, aretesVivesCoupees, dmoyen):
 
   # --- maillage des éventuelles arêtes vives entre faces reconstruites
   
+  grpAretesVives = None
   if len(aretesVivesCoupees) > 0:
     aretesVivesC = geompy.MakeCompound(aretesVivesCoupees)
     meshAretesVives = smesh.Mesh(aretesVivesC)
index 94bf8c5a5ead539286df5b9e2e1bd6d23b6a7cf1..90bd8701271531570522c4ecf221a1032cb5fca9 100644 (file)
@@ -5,6 +5,7 @@ import logging
 from geomsmesh import geompy
 from geomsmesh import smesh
 from salome.smesh import smeshBuilder
+import SMESH
 
 from putName import putName
   
index 7b25c43acb4ad836b52f6bad1b5286e3f9742c36..5ff5170bedd6655b61e4a2f813fa5013e96981d3 100644 (file)
@@ -5,6 +5,7 @@ import logging
 from geomsmesh import geompy
 from geomsmesh import smesh
 from salome.smesh import smeshBuilder
+import SMESH
 
 from putName import putName