]> SALOME platform Git repositories - modules/smesh.git/blob - doc/salome/examples/modifying_meshes_ex15.py
Salome HOME
Fix for issue '[EDF]: Broken links in the Developer's user guide in GEOM and SMESH'
[modules/smesh.git] / doc / salome / examples / modifying_meshes_ex15.py
1 # Moving Nodes
2
3
4 import salome
5 salome.salome_init()
6
7 from salome.geom import geomBuilder
8 geompy = geomBuilder.New()
9
10 import SMESH, SALOMEDS
11 from salome.smesh import smeshBuilder
12 smesh =  smeshBuilder.New()
13
14
15 box = geompy.MakeBoxDXDYDZ(200, 200, 200)
16
17 mesh = smesh.Mesh( box )
18 mesh.Segment().AutomaticLength(0.1)
19 mesh.Quadrangle()
20 mesh.Compute()
21
22 # find node at (0,0,0) which is located on a geom vertex
23 node000 = None
24 for vId in geompy.SubShapeAllIDs( box, geompy.ShapeType["VERTEX"]):
25     if node000: break
26     nodeIds = mesh.GetSubMeshNodesId( vId, True )
27     for node in nodeIds:
28         xyz = mesh.GetNodeXYZ( node )
29         if xyz[0] == 0 and xyz[1] == 0 and xyz[2] == 0 :
30             node000 = node
31             pass
32         pass
33     pass
34
35 if not node000:
36     raise Exception("node000 not found")
37
38 # find node000 using a dedicated function 
39 n = mesh.FindNodeClosestTo( -1,-1,-1 )
40 if not n == node000:
41     raise Exception("FindNodeClosestTo() returns " + str( n ) + " != " + str( node000 ))
42
43 # move node000 to a new location
44 x,y,z = -10, -10, -10
45 n = mesh.MoveNode( n,x,y,z )
46 if not n:
47     raise Exception("MoveNode() returns " + n)
48
49 # check the coordinates of the node000
50 xyz = mesh.GetNodeXYZ( node000 )
51 if not ( xyz[0] == x and xyz[1] == y and xyz[2] == z) :
52     raise Exception("Wrong coordinates: " + str( xyz ) + " != " + str( [x,y,z] ))