+
+# ======================================================= test_cone
+def test_cone () :
+
+ doc = hexablock.addDocument ("cyl:hexas")
+
+ orig = doc.addVertex (0, 0, 0)
+ ori2 = doc.addVertex (0, 0, 10)
+ vz = doc.addVector (0, 0, 1)
+ vx = doc.addVector (1 ,0, 0)
+
+ nr = 2
+ na = 8
+ nl = 1
+ c1 = doc.makeCylindrical (orig, vx,vz, 2.0, 360.0, 5.0, nr,na,nl, False)
+ c2 = doc.makeCylindrical (ori2, vx,vz, 1.0, 360.0, 5.0, nr,na,nl, False)
+
+ vh0 = c2.getVertexIJK (0, 0, 0)
+ vh1 = c2.getVertexIJK (1, 0, 0)
+
+ vb0 = c1.getVertexIJK (0, 0, nl)
+ vb1 = c1.getVertexIJK (1, 0, nl)
+
+ qcible = c2.getQuadIJ (0, 0, 0)
+ qliste = []
+ for ni in range (nr) :
+ for nj in range (na) :
+ quad = c1.getQuadIJ (ni, nj, nl)
+ qliste.append (quad)
+ quad .setColor (3);
+
+ qcible.setColor (3);
+ vb0.setColor (5);
+ vh0.setColor (5);
+ vb1.setColor (5);
+ vh1.setColor (5);
+
+ doc.joinQuads (qliste, qcible, vb0, vh0, vb1, vh1, 1)
+
+ doc.saveVtk ("cone.vtk")
+ return doc
+