def xfiber(cutter,xvals,t,zh,myscreen):
for x in xvals:
f1 = ocl.Point(x,-5.5,zh) # start point of fiber
f2 = ocl.Point(x,5.5,zh) # end point of fiber
f = ocl.Fiber( f1, f2)
i = ocl.Interval()
cutter.pushCutter(f,i,t)
f.addInterval(i)
drawFiber_clpts(myscreen, f, camvtk.lblue)
def yfiber(cutter,yvals,t,zh,myscreen):
for y in yvals:
f1 = ocl.Point(-5.5,y,zh) # start point of fiber
f2 = ocl.Point(5.5,y,zh) # end point of fiber
f = ocl.Fiber( f1, f2)
i = ocl.Interval()
cutter.pushCutter(f,i,t)
f.addInterval(i)
drawFiber_clpts(myscreen, f, camvtk.red)
def xfiber(xvals, t, zh, myscreen):
for x in xvals:
f1 = ocl.Point(x, -0.5, zh) # start point of fiber
f2 = ocl.Point(x, 1.5, zh) # end point of fiber
f = ocl.Fiber(f1, f2)
i = ocl.Interval()
#cutter.vertexPush(f,i,t)
#cutter.facetPush(f,i,t)
#cutter.edgePush(f,i,t)
cutter.pushCutter(f, i, t)
f.addInterval(i)
drawFiber_clpts(myscreen, f, camvtk.lblue)
def yfiber(yvals, t, zh, myscreen):
for y in yvals:
f1 = ocl.Point(-0.5, y, zh) # start point of fiber
f2 = ocl.Point(1.5, y, zh) # end point of fiber
f = ocl.Fiber(f1, f2)
i = ocl.Interval()
#cutter.vertexPush(f,i,t)
#cutter.facetPush(f,i,t)
#cutter.edgePush(f,i,t)
cutter.pushCutter(f, i, t)
f.addInterval(i)
drawFiber_clpts(myscreen, f, camvtk.red)
def xfiber(xvals,s,zh,myscreen):
for x in xvals:
f1 = ocl.Point(x,-20,zh) # start point of fiber
f2 = ocl.Point(x,+20,zh) # end point of fiber
f = ocl.Fiber( f1, f2)
for t in s.getTriangles():
i = ocl.Interval()
#cutter.vertexPush(f,i,t)
#cutter.facetPush(f,i,t)
#cutter.edgePush(f,i,t)
cutter.pushCutter(f,i,t)
f.addInterval(i)
drawFiber_clpts(myscreen, f, camvtk.lblue)
#cutter = ocl.BallCutter(diameter, length)
#cutter = ocl.BullCutter(diameter, diameter/5, length)
#cutter = ocl.ConeCutter(diameter, angle, length)
print cutter
range=6
Nmax = 500
yvals = [float(n-float(Nmax)/2)/Nmax*range for n in xrange(0,Nmax+1)]
xvals = [float(n-float(Nmax)/2)/Nmax*range for n in xrange(0,Nmax+1)]
zh = 0.1071567
for y in yvals:
f1 = ocl.Point(-0.5,y,zh) # start point of fiber
f2 = ocl.Point(1.5,y,zh) # end point of fiber
f = ocl.Fiber( f1, f2)
i = ocl.Interval()
cutter.pushCutter(f,i,t)
#cutter.vertexPush(f,i,t)
#cutter.facetPush(f,i,t)
#cutter.edgePush(f,i,t)
f.addInterval(i)
drawFiber(myscreen, f, camvtk.red)
for x in xvals:
f1 = ocl.Point(x,-0.5,zh) # start point of fiber
f2 = ocl.Point(x,1.5,zh) # end point of fiber
f = ocl.Fiber( f1, f2)
i = ocl.Interval()
cutter.pushCutter(f,i,t)
#cutter.vertexPush(f,i,t)
#cutter.facetPush(f,i,t)
import ocl
import camvtk
import time
import vtk
import datetime
import math
if __name__ == "__main__":
print(ocl.version())
i = ocl.Interval()
print(i)
i = ocl.Interval(0.123, 0.456)
print(i)
cc = ocl.CCPoint(1, 2, 3)
cc.type = ocl.CCType.VERTEX
i.updateLower(-0.222, cc)
print(i)
i.updateLower(0, cc)
print(i)
i.updateUpper(0, cc)
print(i)
i.updateUpper(2, cc)
print(i)
