doc/test_doc/make_symmetry/make_sym_plane.py

   1 ######## Make Symmetry Plane ########
   2
   3 import hexablock
   4
   5
   6 # ======================================================= make_grid
   7 def make_grid (doc) :
   8
   9     ori  = doc.addVertex ( 0, 0, 0)
  10     vz   = doc.addVector ( 0, 0, 1)
  11     vx   = doc.addVector ( 1 ,0, 0)
  12
  13     dr = 1
  14     da = 360
  15     dl = 1
  16
  17     nr = 1
  18     na = 6
  19     nl = 1
  20     grid = doc.makeCylindrical (ori, vx,vz, dr,da,dl, nr,na,nl, False)
  21     #####  doc .saveVtk ("cyl_grid.vtk")
  22     return grid
  23
  24
  25 # ======================================================= test_sym_plane
  26 def test_sym_plane () :
  27
  28     doc  = hexablock.addDocument ("default")
  29     grid = make_grid (doc)
  30
  31     orig = grid.getVertex(3)
  32     dir   = doc.addVector (0, 1, 0);
  33     grid2 = doc.makeSymmetryPlane (grid, orig, dir)
  34
  35     #####  doc .saveVtk ("sym_plane.vtk")
  36     return doc
  37
  38 # ================================================================= Begin
  39 doc = test_sym_plane  ()