print "dir=", dir
(cl,ccp,cp,d) = cutter.intersect_sphere_line(dir, edge)
print "cp=", cp
print "ccp=", ccp
print "d=", d
print "cl=", cl
exit()
if False:
cutter = SphericalCutter(1, Point(-10,0,0))
edge = Line(Point(0,-5,1),Point(3,+5,1))
dir = Point(1,-0.2,0)
print "cutter=", cutter
print "edge=", edge
print "dir=", dir
(cl,ccp,cp,d) = cutter.intersect_cylinder_line(dir, edge)
print "cp=", cp
print "ccp=", ccp
print "d=", d
print "cl=", cl
dir = Point(0,0,-1)
c = SphericalCutter(1, Point(0,0,6))
#c = SphericalCutter(1, Point(-2.2,0.2,6))
#c = SphericalCutter(1, Point(-1.7,0.5,6))
print "c=", c
t = Triangle(Point(-2,0,2), Point(2,1,3), Point(2,-1,4))
#t = Triangle(Point(-2,0,2), Point(2,-1,4), Point(2,1,3))
#t = Triangle(Point(2,0,4), Point(2,-1,2), Point(2,1,2))
#t = Triangle(Point(2,0,2), Point(-2,1,2), Point(-2,-1,2))
print "t=", t
print "dir=", dir
(cl, ccp, cp, d) = cutter.intersect_sphere_line(dir, edge)
print "cp=", cp
print "ccp=", ccp
print "d=", d
print "cl=", cl
exit()
if False:
cutter = SphericalCutter(1, Point(-10, 0, 0))
edge = Line(Point(0, -5, 1), Point(3, +5, 1))
dir = Point(1, -0.2, 0)
print "cutter=", cutter
print "edge=", edge
print "dir=", dir
(cl, ccp, cp, d) = cutter.intersect_cylinder_line(dir, edge)
print "cp=", cp
print "ccp=", ccp
print "d=", d
print "cl=", cl
dir = Point(0, 0, -1)
c = SphericalCutter(1, Point(0, 0, 6))
#c = SphericalCutter(1, Point(-2.2,0.2,6))
#c = SphericalCutter(1, Point(-1.7,0.5,6))
print "c=", c
t = Triangle(Point(-2, 0, 2), Point(2, 1, 3), Point(2, -1, 4))
#t = Triangle(Point(-2,0,2), Point(2,-1,4), Point(2,1,3))
#t = Triangle(Point(2,0,4), Point(2,-1,2), Point(2,1,2))
#t = Triangle(Point(2,0,2), Point(-2,1,2), Point(-2,-1,2))
print "t=", t
