doc/salome/examples/modifying_meshes_ex16.py

   1 # Diagonal Inversion
   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 # create an empty mesh structure
  17 mesh = smesh.Mesh()
  18
  19 # create the following mesh:
  20 # .----.----.----.
  21 # |   /|   /|   /|
  22 # |  / |  / |  / |
  23 # | /  | /  | /  |
  24 # |/   |/   |/   |
  25 # .----.----.----.
  26
  27 bb = [0, 0, 0, 0]
  28 tt = [0, 0, 0, 0]
  29 ff = [0, 0, 0, 0, 0, 0]
  30
  31 bb[0] = mesh.AddNode( 0., 0., 0.)
  32 bb[1] = mesh.AddNode(10., 0., 0.)
  33 bb[2] = mesh.AddNode(20., 0., 0.)
  34 bb[3] = mesh.AddNode(30., 0., 0.)
  35
  36 tt[0] = mesh.AddNode( 0., 15., 0.)
  37 tt[1] = mesh.AddNode(10., 15., 0.)
  38 tt[2] = mesh.AddNode(20., 15., 0.)
  39 tt[3] = mesh.AddNode(30., 15., 0.)
  40
  41 ff[0] = mesh.AddFace([bb[0], bb[1], tt[1]])
  42 ff[1] = mesh.AddFace([bb[0], tt[1], tt[0]])
  43 ff[2] = mesh.AddFace([bb[1], bb[2], tt[2]])
  44 ff[3] = mesh.AddFace([bb[1], tt[2], tt[1]])
  45 ff[4] = mesh.AddFace([bb[2], bb[3], tt[3]])
  46 ff[5] = mesh.AddFace([bb[2], tt[3], tt[2]])
  47
  48 # inverse the diagonal bb[1] - tt[2]
  49 print "\nDiagonal inversion ... ",
  50 res = mesh.InverseDiag(bb[1], tt[2])
  51 if not res: print "failed!"
  52 else:       print "done."
  53
  54 salome.sg.updateObjBrowser(1)