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