]> SALOME platform Git repositories - modules/smesh.git/blob - src/Tools/blocFissure/gmu/facesCirculaires.py
Salome HOME
Allow saving groups with non-ascii names.
[modules/smesh.git] / src / Tools / blocFissure / gmu / facesCirculaires.py
1 # -*- coding: utf-8 -*-
2
3 import logging
4 from .geomsmesh import geompy
5 from .geomsmesh import geomPublish
6 from .geomsmesh import geomPublishInFather
7 from . import initLog
8
9 # -----------------------------------------------------------------------------
10 # --- TORE
11 ## --- faces 1/2 circulaires et leur centre, edges de ces faces dans le plan de fissure
12
13 def facesCirculaires(bloc, tore):
14   """
15   Extraction des faces demi circulaires à l'intersection du tore partitionné et de la paroi,
16   de leur centre, les edges de ces faces situees dans le plan de fissure et un booleen par edge,
17   indiquant son sens (normal / reversed).
18   @param bloc : bloc defaut
19   @param tore : le tore partitionné et coupé
20   @return (faces, centres, edges, reverses)
21   """
22   logging.info("start")
23
24   faces = geompy.GetShapesOnShape(bloc, tore, geompy.ShapeType["FACE"], GEOM.ST_ON)
25
26   geomPublishInFather(initLog.debug, tore, faces[0], 'face0' )
27   geomPublishInFather(initLog.debug, tore, faces[1], 'face1' )
28   geomPublishInFather(initLog.debug, tore, faces[2], 'face2' )
29   geomPublishInFather(initLog.debug, tore, faces[3], 'face3' )
30
31   centres = [None, None, None, None]
32   [v1,centres[0],v3] = geompy.ExtractShapes(faces[0], geompy.ShapeType["VERTEX"], True)
33   [v1,centres[1],v3] = geompy.ExtractShapes(faces[1], geompy.ShapeType["VERTEX"], True)
34   [v1,centres[2],v3] = geompy.ExtractShapes(faces[2], geompy.ShapeType["VERTEX"], True)
35   [v1,centres[3],v3] = geompy.ExtractShapes(faces[3], geompy.ShapeType["VERTEX"], True)
36
37   geomPublishInFather(initLog.debug, faces[0], centres[0], 'centre0' )
38   geomPublishInFather(initLog.debug, faces[1], centres[1], 'centre1' )
39   geomPublishInFather(initLog.debug, faces[2], centres[2], 'centre2' )
40   geomPublishInFather(initLog.debug, faces[3], centres[3], 'centre3' )
41
42   alledges = [None, None, None, None]
43   alledges[0] = geompy.ExtractShapes(faces[0], geompy.ShapeType["EDGE"], True)
44   alledges[1] = geompy.ExtractShapes(faces[1], geompy.ShapeType["EDGE"], True)
45   alledges[2] = geompy.ExtractShapes(faces[2], geompy.ShapeType["EDGE"], True)
46   alledges[3] = geompy.ExtractShapes(faces[3], geompy.ShapeType["EDGE"], True)
47
48   dicoedge = {}
49   edges = []
50   reverses = []
51   for i in range(len(alledges)):
52     edgesface = alledges[i]
53     lenef = []
54     for j in range(len(edgesface)):
55       props = geompy.BasicProperties(edgesface[j])
56       lenef.append(props[0])
57       pass
58     maxlen = max(lenef)
59     for j in range(len(edgesface)):
60       if lenef[j] < maxlen:
61         edgid = geompy.GetSubShapeID(tore, edgesface[j])
62         if not (edgid in dicoedge):
63           dicoedge[edgid] = edgesface[j]
64           edges.append(edgesface[j])
65           named = 'edge_' + str(i) + '_' +str(j)
66           geomPublishInFather(initLog.debug, faces[i], edgesface[j], named)
67           vertices = geompy.ExtractShapes(edgesface[j], geompy.ShapeType["VERTEX"], False)
68           #firstVertex = geompy.GetFirstVertex(edgesface[j])
69           if geompy.GetSubShapeID(tore, vertices[0]) != geompy.GetSubShapeID(tore, centres[i]):
70             reverses.append(1)
71             #print 'reversed ' + str(edgid)
72           else:
73             reverses.append(0)
74             #print 'normal' + str(edgid)
75           pass
76         pass
77       pass
78     pass
79
80   return faces, centres, edges, reverses
81