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)