src/TEST_PY/test_v6/makeCylinderTop.py

   1
   2 import hexablock
   3
   4 # ================================================ rename
   5 def rename (elt, pref, indices) :
   6     if elt == None : return
   7     elt.setName (pref + indices)
   8 # ================================================ begin
   9 doc = hexablock.addDocument ("model")
  10 nr = 3
  11 na = 8
  12 nh = 3
  13
  14 grid = doc.makeCylinderTop (nr, na, nh)
  15 for i in range (nr+1) :
  16    for j in range (na) :
  17       for k in range (nh+1) :
  18           name  = "(%d,%d,%d)" % (i,j,k)
  19           vertex =  grid.getVertexIJK (i,j,k)
  20           ei     =  grid.getEdgeI (i,j,k)
  21           ej     =  grid.getEdgeJ (i,j,k)
  22           ek     =  grid.getEdgeK (i,j,k)
  23           rename (vertex, 'v', name)
  24           rename (ei,    'ei', name)
  25           rename (ej,    'ej', name)
  26           rename (ek,    'ek', name)
  27
  28 doc.save    ("cylinder")
  29 doc.saveVtk ("cylinder.vtk")
  30 doc.addLaws (1, True)
  31 mesh_hexas = hexablock.mesh(doc)