doc/test_doc/make_symmetry/make_sym_point.py

   1 ######## Make Symmetry Point ########
   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_point
  26 def test_sym_point () :
  27
  28     doc  = hexablock.addDocument ("default")
  29     grid = make_grid (doc)
  30
  31     orig = grid.getVertex(3)
  32     grid2 = doc.makeSymmetryPoint (grid, orig)
  33
  34     #####  doc .saveVtk ("sym_point.vtk")
  35     return doc
  36
  37
  38 # ================================================================= Begin
  39 doc = test_sym_point  ()
  40