t_before = time.time()
print "threads=",bdc1.getThreads()
bdc1.run()
t_after = time.time()
calctime = t_after-t_before
print " done in ", calctime," s"
clpoints = bdc1.getCLPoints()
print len(clpoints), " cl points evaluated"
print "rendering...",
# draw the CL-points
camvtk.drawCLPointCloud(myscreen, clpoints)
print "done"
myscreen.camera.SetPosition(3, 23, 15)
myscreen.camera.SetFocalPoint(4, 5, 0)
w2if = vtk.vtkWindowToImageFilter()
w2if.SetInput(myscreen.renWin)
lwr = vtk.vtkPNGWriter()
lwr.SetInput( w2if.GetOutput() )
w2if.Modified()
lwr.SetFileName("tux1.png")
#lwr.Write()
path = ocl.Path() # create an empty path object
# add Line objects to the path in this loop
for n in range(0, Ny):
y = ymin + n * dy
p1 = ocl.Point(0, y, 0) # start-point of line
p2 = ocl.Point(9, y, 0) # end-point of line
l = ocl.Line(p1, p2) # line-object
path.append(l) # add the line to the path
print(" set the path for pdf ")
pdc.setPath(path)
print(" run the calculation ")
t_before = time.time()
pdc.run() # run drop-cutter on the path
t_after = time.time()
print("run took ", t_after - t_before, " s")
print("get the results ")
clp = pdc.getCLPoints() # get the cl-points from pdf
print(" render the CL-points")
camvtk.drawCLPointCloud(myscreen, clp)
#myscreen.addActor( camvtk.PointCloud(pointlist=clp, collist=ccp) )
myscreen.camera.SetPosition(3, 23, 15)
myscreen.camera.SetFocalPoint(5, 5, 0)
myscreen.render()
print(" All done.")
myscreen.iren.Start()
t_before = time.time()
bdc.run() # run the actual drop-cutter
dc_calls = bdc.getCalls()
t_after = time.time()
calctime = t_after - t_before
print " BDC done in ", calctime, " s", dc_calls, " dc-calls"
dropcutter_time = calctime
clpoints = bdc.getCLPoints()
#print len(clpoints), " cl points evaluated"
print "rendering raw CL-points."
# draw the CL-points
camvtk.drawCLPointCloud(myscreen, clpoints)
# filter
print "filtering. before filter we have", len(clpoints), "cl-points"
t_before = time.time()
f = ocl.LineCLFilter()
f.setTolerance(0.001)
for p in clpoints:
f.addCLPoint(p)
f.run()
t_after = time.time()
calctime = t_after - t_before
print " done in ", calctime, " s"
clp2 = f.getCLPoints()
print "after filtering we have", len(clp2), "cl-points"
print "filtering. before filter we have", len(aclp_lifted), "cl-points"
t_before = time.time()
f = ocl.LineCLFilter()
f.setTolerance(0.001)
for p in aclp_lifted:
p2 = ocl.CLPoint(p.x, p.y, p.z)
f.addCLPoint(p2)
f.run()
t_after = time.time()
calctime = t_after - t_before
print " done in ", calctime, " s"
cl_filtered = f.getCLPoints()
aclp_lifted2 = []
for p in cl_filtered:
p2 = ocl.Point(p.x, p.y, p.z) + ocl.Point(0, 0, 1)
aclp_lifted2.append(p2)
print " render the CL-points"
camvtk.drawCLPointCloud(myscreen, clp)
camvtk.drawCLPointCloud(myscreen, aclp_lifted)
camvtk.drawCLPointCloud(myscreen, aclp_lifted2)
#myscreen.addActor( camvtk.PointCloud(pointlist=clp, collist=ccp) )
myscreen.camera.SetPosition(3, 23, 15)
myscreen.camera.SetFocalPoint(5, 5, 0)
myscreen.render()
print " All done."
myscreen.iren.Start()
#cutter.edgeDrop(cl,t)
#cutter.facetDrop(cl,t)
#c2.vertexDrop(cl2,t)
cutter.dropCutter(cl, t)
c2.dropCutter(cl2, t)
n = n + 1
if (n % int(len(clpoints) / 10)) == 0:
print n / int(len(clpoints) / 10), " ",
print "done."
print "rendering..."
print " len(clpoints)=", len(clpoints)
print " len(ccl2pts)=", len(cl2pts)
print "rendering clpoints...",
camvtk.drawCLPointCloud(myscreen, clpoints)
print "done."
cl2ptsofs = []
for p in cl2pts:
p.z = p.z + offset
cl2ptsofs.append(p)
print "rendering offset clpoints...",
camvtk.drawCLPointCloud(myscreen, cl2ptsofs)
print "done."
origo = camvtk.Sphere(center=(0, 0, 0), radius=0.1, color=camvtk.blue)
origo.SetOpacity(0.2)
myscreen.addActor(origo)
myscreen.camera.SetPosition(0.5, 3, 2)
myscreen.camera.SetFocalPoint(0.5, 0.5, 0)
myscreen.render()
print "filtering. before filter we have", len(aclp_lifted),"cl-points"
t_before = time.time()
f = ocl.LineCLFilter()
f.setTolerance(0.001)
for p in aclp_lifted:
p2 = ocl.CLPoint(p.x,p.y,p.z)
f.addCLPoint(p2)
f.run()
t_after = time.time()
calctime = t_after-t_before
print " done in ", calctime," s"
cl_filtered = f.getCLPoints()
aclp_lifted2=[]
for p in cl_filtered:
p2 = ocl.Point(p.x,p.y,p.z) + ocl.Point(0,0,1)
aclp_lifted2.append(p2)
print " render the CL-points"
camvtk.drawCLPointCloud(myscreen, clp)
camvtk.drawCLPointCloud(myscreen, aclp_lifted)
camvtk.drawCLPointCloud(myscreen, aclp_lifted2)
#myscreen.addActor( camvtk.PointCloud(pointlist=clp, collist=ccp) )
myscreen.camera.SetPosition(3, 23, 15)
myscreen.camera.SetFocalPoint(5, 5, 0)
myscreen.render()
print " All done."
myscreen.iren.Start()
