doc/test_doc/quadRevolution/revolution_quad.py

   1 # -*- coding: utf-8 -*-
   2 #### Quad Revolution Test #####
   3
   4 import hexablock
   5
   6
   7 doc = hexablock.addDocument("Quad Revolution Test")
   8
   9 center = doc.addVertex(0, 0, 0)
  10 vz     = doc.addVector(0, 0, 1)
  11
  12 v1     = doc.addVertex (10, 0, 0)
  13 v2     = doc.addVertex (11, 0, 0)
  14 v3     = doc.addVertex (11, 0, 2)
  15 v4     = doc.addVertex (10, 0, 2)
  16 quad   = doc.addQuadVertices (v1, v2, v3, v4)
  17 doc.saveVtk("revolution1.vtk")
  18
  19 angle = 180
  20 nbLayers = 8
  21 grid0 = doc.revolutionQuadUni  (quad, center, vz, angle, nbLayers)
  22 doc.saveVtk ("revolution2.vtk")
  23
  24 nr = 1
  25 na = 6
  26 nl = 1
  27
  28 grid = doc.makeCylinderTop (nr,na,nl)
  29
  30 liste = []
  31 for nx in range(nr):
  32     for ny in range(na):
  33         cell = grid.getQuadIJ (nx, ny, nl)
  34         liste.append(cell)
  35
  36 center = doc.addVertex(0, -10, 0)
  37 axis = doc.addVector (1, 0, 0)
  38 angle = 180
  39 nbLayers = 9
  40 grid1 = doc.revolutionQuadsUni  (liste, center, axis, angle, nbLayers)
  41 if grid1.isValid():
  42     doc.saveVtk ("revolution3.vtk")