X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FTools%2FblocFissure%2Fgmu%2FelimineExtremitesPipe.py;h=36927198062b3e7fd3cdc9db8e32faa47578a6a2;hb=f44efbc02056a2734440b736789b8cfa380430bb;hp=b584c0c91c87b23c2daf04d1152c205b9668abda;hpb=71de383cb614e31ded6aa6343c052fbd8335b472;p=modules%2Fsmesh.git diff --git a/src/Tools/blocFissure/gmu/elimineExtremitesPipe.py b/src/Tools/blocFissure/gmu/elimineExtremitesPipe.py index b584c0c91..369271980 100644 --- a/src/Tools/blocFissure/gmu/elimineExtremitesPipe.py +++ b/src/Tools/blocFissure/gmu/elimineExtremitesPipe.py @@ -17,6 +17,7 @@ # # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com # +"""Recherche des points en trop (externes au volume à remailler)""" import logging @@ -24,33 +25,34 @@ from .geomsmesh import geompy from .whichSideVertex import whichSideVertex def elimineExtremitesPipe(ptEdgeFond, facesDefaut, centres, gptsdisks, nbsegCercle): - """ - recherche des points en trop (externes au volume à remailler) - - on associe chaque extrémité du pipe à une face filling + """Recherche des points en trop (externes au volume à remailler) + + - on associe chaque extrémité du pipe à une face filling - on part des disques aux extrémités du pipe - pour chaque disque, on prend les vertices de géométrie, on marque leur position relative à la face. - on s'arrete quand tous les noeuds sont dedans """ - + logging.info('start') - pt0 = centres[0] - pt1 = centres[-1] + point_0 = centres[0] + point_1 = centres[-1] idFillingFromBout = [None, None] # contiendra l'index du filling pour les extrémités 0 et 1 nbFacesFilling = len(ptEdgeFond) for ifil in range(nbFacesFilling): - for ipt, pt in enumerate(ptEdgeFond[ifil]): # il y a un ou deux points débouchant sur cette face - if geompy.MinDistance(pt,pt0) < geompy.MinDistance(pt,pt1): # TODO: trouver plus fiable pour les cas tordus... + for point in ptEdgeFond[ifil]: # il y a un ou deux points débouchant sur cette face + if geompy.MinDistance(point,point_0) < geompy.MinDistance(point,point_1): # TODO: trouver plus fiable pour les cas tordus... idFillingFromBout[0] = ifil - else: + else: idFillingFromBout[1] = ifil logging.debug("association bouts du pipe - faces de filling: %s", idFillingFromBout) - + logging.debug("recherche des disques de noeuds complètement internes") - idisklim = [] # indices des premier et dernier disques internes - idiskout = [] # indices des premier et dernier disques externes + idisklim = list() # indices des premier et dernier disques internes + idiskout = list() # indices des premier et dernier disques externes for bout in range(2): + if bout == 0: idisk = -1 inc = 1 @@ -59,26 +61,28 @@ def elimineExtremitesPipe(ptEdgeFond, facesDefaut, centres, gptsdisks, nbsegCerc idisk = len(gptsdisks) inc = -1 numout = len(gptsdisks) + inside = False outside = True while not inside: - idisk = idisk + inc + idisk += inc logging.debug("examen disque %s", idisk) gptdsk = gptsdisks[idisk] inside = True - for k in range(nbsegCercle): - points = gptdsk[k] - for j, pt in enumerate(points): - side = whichSideVertex(facesDefaut[idFillingFromBout[bout]], pt) + for n_seg in range(nbsegCercle): + points = gptdsk[n_seg] + for point in points: + side = whichSideVertex(facesDefaut[idFillingFromBout[bout]], point) if side < 0: if outside: # premier point detecté dedans outside = False numout = idisk -inc # le disque précédent était dehors else: - inside = False # ce point est dehors + inside = False # ce point est dehors if not inside and not outside: break + idisklim.append(idisk) # premier et dernier disques internes idiskout.append(numout) # premier et dernier disques externes - - return (idFillingFromBout, idisklim, idiskout) \ No newline at end of file + + return (idFillingFromBout, idisklim, idiskout)