doc/examples/modifying_meshes_ex17.py

   1 # Uniting two Triangles
   2
   3 import salome
   4 salome.salome_init_without_session()
   5
   6 from salome.smesh import smeshBuilder
   7
   8 smesh_builder = smeshBuilder.New()
   9
  10 # create an empty mesh structure
  11 mesh = smesh_builder.Mesh()
  12
  13 # create the following mesh:
  14 # .----.----.----.
  15 # |   /|   /|   /|
  16 # |  / |  / |  / |
  17 # | /  | /  | /  |
  18 # |/   |/   |/   |
  19 # .----.----.----.
  20
  21 bb = [0, 0, 0, 0]
  22 tt = [0, 0, 0, 0]
  23 ff = [0, 0, 0, 0, 0, 0]
  24
  25 bb[0] = mesh.AddNode( 0., 0., 0.)
  26 bb[1] = mesh.AddNode(10., 0., 0.)
  27 bb[2] = mesh.AddNode(20., 0., 0.)
  28 bb[3] = mesh.AddNode(30., 0., 0.)
  29
  30 tt[0] = mesh.AddNode( 0., 15., 0.)
  31 tt[1] = mesh.AddNode(10., 15., 0.)
  32 tt[2] = mesh.AddNode(20., 15., 0.)
  33 tt[3] = mesh.AddNode(30., 15., 0.)
  34
  35 ff[0] = mesh.AddFace([bb[0], bb[1], tt[1]])
  36 ff[1] = mesh.AddFace([bb[0], tt[1], tt[0]])
  37 ff[2] = mesh.AddFace([bb[1], bb[2], tt[2]])
  38 ff[3] = mesh.AddFace([bb[1], tt[2], tt[1]])
  39 ff[4] = mesh.AddFace([bb[2], bb[3], tt[3]])
  40 ff[5] = mesh.AddFace([bb[2], tt[3], tt[2]])
  41
  42 # delete the diagonal bb[1] - tt[2]
  43 print("\nUnite two triangles ... ", end=' ')
  44 res = mesh.DeleteDiag(bb[1], tt[2])
  45 if not res: print("failed!")
  46 else:       print("done.")